Oxidation-reduction process of Arabidopsis thaliana roots induced by bisphenol compounds based on RNA-seq analysis.

Bisphenol compounds (BPs) have various industrial uses and can enter the environment through various sources. To evaluate the ecotoxicity of BPs and identify potential gene candidates involved in the plant toxicity, Arabidopsis thaliana was exposed to bisphenol A (BPA), BPB, BPE, BPF, and BPS at 1, 3, 10 mg/L for a duration of 14 days, and their growth status were monitored. At day 14, roots and leaves were collected for internal BPs exposure concentration detection, RNA-seq (only roots), and morphological observations. As shown in the results, exposure to BPs significantly disturbed root elongation, exhibiting a trend of stimulation at low concentration and inhibition at high concentration. Additionally, BPs exhibited pronounced generation of reactive oxygen species, while none of the pollutants caused significant changes in root morphology. Internal exposure concentration analysis indicated that BPs tended to accumulate in the roots, with BPS exhibiting the highest level of accumulation. The results of RNA-seq indicated that the shared 211 differently expressed genes (DEGs) of these 5 exposure groups were enriched in defense response, generation of precursor metabolites, response to organic substance, response to oxygen-containing, response to hormone, oxidation-reduction process and so on. Regarding unique DEGs in each group, BPS was mainly associated with the redox pathway, BPB primarily influenced seed germination, and BPA, BPE and BPF were primarily involved in metabolic signaling pathways. Our results provide new insights for BPs induced adverse effects on Arabidopsis thaliana and suggest that the ecological risks associated with BPA alternatives cannot be ignored.

Identification of deregulated lncRNAs in Alzheimer's disease: an integrated gene co-expression network analysis of hippocampus and fusiform gyrus RNA-seq datasets.

Introduction: The deregulation of lncRNAs expression has been associated with neuronal damage in Alzheimer's disease (AD), but how or whether they can influence its onset is still unknown. We investigated 2 RNA-seq datasets consisting, respectively, of the hippocampal and fusiform gyrus transcriptomic profile of AD patients, matched with non-demented controls. Methods: We performed a differential expression analysis, a gene correlation network analysis (WGCNA) and a pathway enrichment analysis of two RNA-seq datasets. Results: We found deregulated lncRNAs in common between hippocampus and fusiform gyrus and deregulated gene groups associated to functional pathways related to neurotransmission and memory consolidation. lncRNAs, co-expressed with known AD-related coding genes, were identified from the prioritized modules of both brain regions. Discussion: We found common deregulated lncRNAs in the AD hippocampus and fusiform gyrus, that could be considered common signatures of AD pathogenesis, providing an important source of information for understanding the molecular changes of AD.

Isolation, expansion, and adoptive transfer of human ILC2s for the treatment of mice bearing liquid and solid tumors.

Type 2 innate lymphoid cells (ILC2s) are crucial in regulating immune responses and various physiological processes, including tissue repair, metabolic homeostasis, inflammation, and cancer surveillance. Here, we present a protocol that outlines the isolation, expansion, and adoptive transfer of human ILC2s from peripheral blood mononuclear cells for an in vivo lineage tracking experiment in a mouse model. Additionally, we detail the steps involved in the adoptive transfer of human ILC2s to recipient mice bearing human liquid or solid tumors. For complete details on the use and execution of this protocol, please refer to Li et al.1.

The PAH1-encoded phosphatidate phosphatase of Yarrowia lipolytica differentially affects gene expression and lipid biosynthesis.

Yarrowia lipolytica is a model oleaginous yeast with a strong capacity for lipid accumulation, yet its lipid metabolic pathways and regulatory mechanisms remain largely unexplored. The PAH1-encoded phosphatidate (PA) phosphatase governs lipid biosynthesis by its enzymatic activity and regulating the transcription of genes involved in phospholipid biosynthesis. In this work, we examined the effect of the loss of Pah1 (i.e., pah1Δ) on cell metabolism in cells growing in low- and high-glucose media. Multi-omics analyses revealed the global effect of the pah1Δ mutation on lipid and central carbon metabolism. Lipidomics analyses showed that the pah1Δ mutation caused a massive decrease in the masses of triacylglycerol (TAG) and diacylglycerol (DAG), and these effects were independent of glucose concentration in the media. Conversely, phospholipid levels declined in low-glucose media but increased in high-glucose media. The loss of Pah1 affected the expression of genes involved in key pathways of glucose metabolism, such as glycolysis, citric acid cycle, oxidative phosphorylation, and the pentose phosphate pathway, and these effects were more pronounced in high-glucose media. In lipid biosynthesis, the genes catalyzing phosphatidylcholine (PC) synthesis from phosphatidylethanolamine (PE) were upregulated within the CDP-DAG pathway. In contrast, PC synthesis through the Kennedy pathway was downregulated. The ethanolamine branch of the Kennedy pathway that synthesizes PE was also upregulated in pah1Δ. Interestingly, we noted a massive increase in the levels of lysophospholipids, consistent with the upregulation of genes involved in lipid turnover. Overall, this work identified novel regulatory roles of Pah1 in lipid biosynthesis and gene expression.

Snow alga Sanguina aurantia as revealed through de novo genome assembly and annotation.

To thrive on melting alpine and polar snow, some Chlorophytes produce an abundance of astaxanthin, causing red blooms, often dominated by genus Sanguina. The red cells have not been cultured, but we recently grew a green biciliate conspecific with Sanguina aurantia from a sample of watermelon snow. This culture provided source material for Oxford Nanopore Technology and Illumina sequencing. Our assembly pipeline exemplifies the value of a hybrid long- and short-read approach for the complexities of working with a culture grown from a field sample. Using bioinformatic tools we separated assembled contigs into two genomic pools based on a difference in GC content (57.5% and 55.1%). We present the data as two assemblies of S. aurantia variants but explore other possibilities. High-throughput chromatin conformation capture analysis (Hi-C sequencing) was used to scaffold the assemblies into a 96 MB genome designated 'A' and a 102 MB genome designated 'B'. Both assemblies are highly contiguous: genome A consists of 38 scaffolds with an N50 of 5.4 Mb while genome B has 50 scaffolds with an N50 of 6.4 Mb. RNA-sequencing was used to improve gene annotation.

Regeneration of tunic cuticle is suppressed in edible ascidian Halocynthia roretzi contracting soft tunic syndrome.

Soft tunic syndrome is an infectious disease caused by the flagellate Azumiobodo hoyamushi, which severely damages the aquaculture of the edible ascidian Halocynthia roretzi. Tunic is a cellulosic extracellular matrix entirely covering the body in ascidians and other tunicates, and its dense cuticle layer covers the tunic surface as a physical barrier against microorganisms. When the tunic of intact H. roretzi individuals was cut into strips, electron-dense fibers (DFs) appeared on the cut surface of the tunic matrix and aggregated to regenerate a new cuticular layer in seawater within a few days. DF formation was partially or completely inhibited in individuals with soft tunic syndrome, and DF formation was also inhibited by the presence of some proteases, indicating the involvement of proteolysis in the process of tunic softening as well as cuticle regeneration. Using pure cultures of the causative flagellate A. hoyamushi, the expression of protease genes and secretion of some proteases were confirmed by RNA-seq analysis and a 4-methylcoumaryl-7-amide substrate assay. Some of these proteases may degrade proteins in the tunic matrix. These findings suggest that the proteases of A. hoyamushi is the key to understanding the mechanisms of cuticular regeneration inhibition and tunic softening.

Transcriptomic changes and gene fusions during the progression from Barrett's esophagus to esophageal adenocarcinoma.

The incidence of esophageal adenocarcinoma (EAC) has surged by 600% in recent decades, with a dismal 5-year survival rate of just 15%. Barrett's esophagus (BE), affecting about 2% of the population, raises the risk of EAC by 40-fold. Despite this, the transcriptomic changes during the BE to EAC progression remain unclear. Our study addresses this gap through comprehensive transcriptomic profiling to identify key mRNA signatures and genomic alterations, such as gene fusions. We performed RNA-sequencing on BE and EAC tissues from 8 individuals, followed by differential gene expression, pathway and network analysis, and gene fusion prediction. We identified mRNA changes during the BE-to-EAC transition and validated our results with single-cell RNA-seq datasets. We observed upregulation of keratin family members in EAC and confirmed increased levels of keratin 14 (KRT14) using immunofluorescence. More differentiated BE marker genes are downregulated during progression to EAC, suggesting undifferentiated BE subpopulations contribute to EAC. We also identified several gene fusions absent in paired BE and normal esophagus but present in EAC. Our findings are critical for the BE-to-EAC transition and have the potential to promote early diagnosis, prevention, and improved treatment strategies for EAC.

Mapping of dynamic quantitative trait loci for plant height in a RIL population of foxtail millet (Setaria italica L.).

Plant height (PH) is a crucial trait for strengthening lodging resistance and boosting yield in foxtail millet. To identify quantitative trait loci (QTL) and candidate genes associated with PH, we first developed a genetic map using a recombinant inbred line (RIL) population derived from a cross between Aininghuang and Jingu 21. Then, PH phenotyping data and four variations of best linear unbiased prediction (BLUP) were collected from nine environments and three development stages. Next, QTL mapping was conducted using both unconditional and conditional QTL methods. Subsequently, candidate genes were predicted via transcriptome analysis of parental samples at three developmental stages. The results revealed that the genetic map, based on re-sequencing, consisted of 4,360 bin markers spanning 1,016.06 cM with an average genetic distance of 0.23 cM. A total of 19 unconditional QTL, accounting for 5.23%-35.36% of the phenotypic variation explained (PVE), which included 7 major and 4 stable QTL, were identified. Meanwhile, 13 conditional QTL, explaining 5.88%-40.35% of PVE, including 5 major and 3 stable QTL, were discovered. Furthermore, four consistent and stable QTL were identified. Finally, eight candidate genes were predicted through RNA-seq and weighted gene co-expression network analysis (WGCNA). Those findings provide a crucial foundation for understanding the genetic mechanisms underlying PH development and facilitate molecular marker-assisted breeding of ideal plant types in foxtail millet.

Transcriptomic analysis of primary nasal epithelial cells reveals altered interferon signalling in preterm birth survivors at one year of age.

Introduction: Many survivors of preterm birth (<37 weeks gestation) have lifelong respiratory deficits, the drivers of which remain unknown. Influencers of pathophysiological outcomes are often detectable at the gene level and pinpointing these differences can help guide targeted research and interventions. This study provides the first transcriptomic analysis of primary nasal airway epithelial cells in survivors of preterm birth at approximately 1 year of age. Methods: Nasal airway epithelial brushings were collected, and primary cell cultures established from term (>37 weeks gestation) and very preterm participants (≤32 weeks gestation). Ex vivo RNA was collected from brushings with sufficient cell numbers and in vitro RNA was extracted from cultured cells, with bulk RNA sequencing performed on both the sample types. Differential gene expression was assessed using the limma-trend pipeline and pathway enrichment identified using Reactome and GO analysis. To corroborate gene expression data, cytokine concentrations were measured in cell culture supernatant. Results: Transcriptomic analysis to compare term and preterm cells revealed 2,321 genes differentially expressed in ex vivo samples and 865 genes differentially expressed in cultured basal cell samples. Over one third of differentially expressed genes were related to host immunity, with interferon signalling pathways dominating the pathway enrichment analysis and IRF1 identified as a hub gene. Corroboration of disrupted interferon release showed that concentrations of IFN-α2 were below measurable limits in term samples but elevated in preterm samples [19.4 (76.7) pg/ml/µg protein, p = 0.03]. IFN-γ production was significantly higher in preterm samples [3.3 (1.5) vs. 9.4 (17.7) pg/ml/µg protein; p = 0.01] as was IFN-ß [7.8 (2.5) vs. 13.6 (19.5) pg/ml/µg protein, p = 0.01]. Conclusion: Host immunity may be compromised in the preterm nasal airway epithelium in early life. Altered immune responses may lead to cycles of repeated infections, causing persistent inflammation and tissue damage which can have significant impacts on long-term respiratory function.

First report of potato virus H infecting tomato (Solanum lycopersicum) in China.

Potato virus H (PVH), belonging to the genus Carlavirus in the family Betaflexiviridae, was initially discovered in potato plants in Inner Mongolia, China (Li et al., 2013). Subsequently, it was documented to infect pepino, a perennial shrub of the Solanaceae family like potatoes (Abouelnasr et al., 2014). Tomato (Solanum lycopersicum L.), a major global crop, faces threats from various plant viruses. In an open field survey in Yunnan, China during July 2023, tomatoes (cultivar: Liangsi) showed typical virus symptoms: leaf yellowing, curling, mottling, and fruit with abnormal shape and color. Eleven symptomatic tomato samples were collected for high-throughput sequencing to identify the potential pathogen. RNA sequencing libraries were prepared using the TruSeq RNA sample prep kit (Illumina, San Diego, CA, USA), followed by RNA-seq sequencing on an Illumina HiSeq4000 platform (LC Sciences, USA). Approximately 77,928,560 paired-end reads (150-bp each) were generated. After quality control, 75,808,296 reads were retained and subjected to de novo assembly using Trinity (version 2.8.5). The assembled contigs, ranging from 198 nt to 15865 nt, were used as queries to search against the NCBI non-redundant protein sequence database (NR) or nucleotide sequence database (NT) to detect the potential pathogens using BLASTx and BLASTn program with a cutoff e-value of 10-5. As a consequence, certain contigs were assigned to 3 plant viruses, including PVH (the highest RdRp blastx identity to UAD82396.1: 97.8%), Capsicum chlorosis virus (CaCV, the highest RdRp blastx identity to APQ31267.1: 98.4%), and southern tomato virus (STV, the highest CP-RdRp fusion protein blastx identity to QOW17541.1: 99.74%). The presence of the identified 3 viruses was subsequently screened in the 11 tomato samples originally collected from the corresponding field. Notably, the specific detection primers for the PVH genome was designed from the newly assembled PVH genome (Forward primer: 5'- ATAGTTGTGCACTGTGTGCCTG-3'; Reverse primer: 5'-GCTTAAGGTTCTTAGCGTATTC-3'), targeting ~1.1kb. Consequently, PVH was detected in 3 out of 11 samples: 2 leaf samples and 1 fruit sample, with one leaf sample showing a single infection. The complete genome sequence of PVH in tomatoes (PVH-tomato) was successfully obtained by assembling nine overlapping regions spanning the entire PVH-tomato genome, following the RT-PCR and the 5' RACE and 3' RACE approaches, and deposited in NCBI nucleotide database with accession number OR397130.1Phylogenetic analysis based on the full genome sequences of PVH-tomato and other publicly available PVH isolates revealed that PVH-tomato was closely related to a PVH isolate found in potatoes in Yunnan (blastn similarity: 97.76%) (Fig. S1A). To test PVH-tomato infectivity and pathogenicity, four healthy Nicotiana benthamiana and four healthy tomato plants were mechanically inoculated with PVH-infected leaf sap; controls used sap from healthy plants. Three weeks post-inoculation, all N. benthamiana (4/4) and three tomato plants (3/4) were PVH-positive by RT-PCR. Symptoms were milder in N. benthamiana, and only two tomato plants (2/4) showed leaf curling. No PVH was detected in control samples (Figure S1B, S1C). Sanger sequencing confirmed the amplicons' expected length of 1093 bp. Previously, PVH was documented only in potato and pepino. This is the first report of tomatoes as natural PVH hosts and PVH infecting N. benthamiana under lab conditions.

A molecular pathway for cancer cachexia-induced muscle atrophy revealed at single-nucleus resolution.

Cancer cachexia is a prevalent and often fatal wasting condition that cannot be fully reversed with nutritional interventions. Muscle atrophy is a central component of the syndrome, but the mechanisms whereby cancer leads to skeletal muscle atrophy are not well understood. We performed single-nucleus multi-omics on skeletal muscles from a mouse model of cancer cachexia and profiled the molecular changes in cachexic muscle. Our results revealed the activation of a denervation-dependent gene program that upregulates the transcription factor myogenin. Further studies showed that a myogenin-myostatin pathway promotes muscle atrophy in response to cancer cachexia. Short hairpin RNA inhibition of myogenin or inhibition of myostatin through overexpression of its endogenous inhibitor follistatin prevented cancer cachexia-induced muscle atrophy in mice. Our findings uncover a molecular basis of muscle atrophy associated with cancer cachexia and highlight potential therapeutic targets for this disorder.

Nitric oxide mitigates the phytotoxicity of imidazolium-based ionic liquids in Arabidopsis.

Ionic liquids (ILs) have many beneficial properties that are extensively used in various fields. Despite their utility, the phytotoxic aspects of ILs are poorly known. This is especially true at the transcriptomic level and the role of nitric oxide (NO) in this process. Herein, we studied the mechanism by which endogenous NO reduces the toxicity of ILs in Arabidopsis. We examined the effects of two imidazolium-based ILs (IILs) on three Arabidopsis lines, each characterized by distinct endogenous NO levels, using a combination of physiological and transcriptomics methods. IILs impaired seed germination, seedling development, chlorophyll content, and redox homeostasis in Arabidopsis. Notably, 1,3-dibutyl imidazole bromide had greater toxicity than 1-butyl-3-methylimidazolium chloride. Nox1, a mutant with an elevated NO level, had enhanced resistance, while nia1nia2, a mutant with a diminished NO level, had increased susceptibility compared to the wild type. RNA sequencing results suggested that NO mitigates IILs-induced phytotoxicity by modulating the metabolism of chlorophyll and secondary metabolites, and by bolstering the antioxidant defense system. These findings illustrate the complex molecular networks that respond to IIL stress and reveal the potential of endogenous NO as a mitigating factor in plant stress physiology.

Comparison of primordial germ cell differences at different developmental time points in chickens.

Objective: Recently, the application in the field of germplasm resource conservation has become an important application of primordial germ cells (PGCs). However, due to the lack of deep understanding of the biological characteristics of PGCs at different time points, there is no systematic scheme for the selection of PGCs at which time points in practical application, which affects the practical application effect of PGCs. This study aims to clarify the differences in PGCs during development. Methods: Here, migration experiment, EdU proliferation assay and cell apoptosis assay were conducted to compare the differences in the migration ability, the proliferation ability and the recovery efficiency among female and male PGCs at E3.5, E4.5 and E5.5, which were explained by the following transcriptome sequencing analysis. Results: We found that there were larger differences between female and male PGCs at different embryonic ages, while smaller differences between female and male PGCs at the same embryonic age. Further comparison showed that the cell migration ability of female and male PGCs decreased gradually during development, so female and male PGCs at E3.5 are more suitable for in vitro allotransplantation. At the same time, the proliferation ability of PGCs gradually decreased during development, and cell adhesion and extracellular matrix communication were weakened, indicating that female and male PGCs of E3.5 are more suitable for in vitro long-term culture cell line establishment. Interestingly, female and male PGCs at E5.5 showed strong DNA damage repair ability, thus more suitable for in vitro long-term cryopreservation. Conclusion: This study provides a theoretical basis for systematically selecting PGCs at suitable developmental time points as cell materials for efficient utilization by analyzing the characteristics of female and male PGCs at different developmental time points based on transcriptome.

Identification of a Novel ECM Remodeling Macrophage Subset in AKI to CKD Transition by Integrative Spatial and Single-Cell Analysis.

The transition from acute kidney injury (AKI) to chronic kidney disease (CKD) is a critical clinical issue. Although previous studies have suggested macrophages as a key player in promoting inflammation and fibrosis during this transition, the heterogeneity and dynamic characterization of macrophages are still poorly understood. Here, we used integrated single-cell RNA sequencing and spatial transcriptomic to characterize the spatiotemporal heterogeneity of macrophages in murine AKI-to-CKD model of unilateral ischemia-reperfusion injury. A marked increase in macrophage infiltration at day 1 was followed by a second peak at day 14 post AKI. Spatiotemporal profiling revealed that injured tubules and macrophages co-localized early after AKI, whereas in late chronic stages had spatial proximity to fibroblasts. Further pseudotime analysis revealed two distinct lineages of macrophages in this transition: renal resident macrophages differentiated into the pro-repair subsets, whereas infiltrating monocyte-derived macrophages contributed to chronic inflammation and fibrosis. A novel macrophage subset, extracellular matrix remodeling-associated macrophages (EAMs) originating from monocytes, linked to renal fibrogenesis and communicated with fibroblasts via insulin-like growth factors (IGF) signalling. In sum, our study identified the spatiotemporal dynamics of macrophage heterogeneity with a unique subset of EAMs in AKI-to-CKD transition, which could be a potential therapeutic target for preventing CKD development.

Analysis of transcriptional response in haploid and diploid Schizosaccharomyces pombe under genotoxic stress.

Whole genome duplications are implicated in genome instability and tumorigenesis. Human and yeast polyploids exhibit increased replication stress and chromosomal instability, both hallmarks of cancer. In this study, we investigate the transcriptional response of Schizosaccharomyces pombe to increased ploidy generally, and in response to treatment with the genotoxin methyl methanesulfonate (MMS). We find that treatment of MMS induces upregulation of genes involved in general response to genotoxins, in addition to cell cycle regulatory genes. Downregulated genes are enriched in transport and sexual reproductive pathways. We find that the diploid response to MMS is muted compared to the haploid response, although the enriched pathways remain largely the same. Overall, our data suggests that the global S. pombe transcriptome doubles in response to increased ploidy but undergoes modest transcriptional changes in both unperturbed and genotoxic stress conditions.

Transcriptional dynamic changes in energy metabolism, protein synthesis and cell cycle regulation reveal the biological adaptation mechanisms of juvenile Acrossocheilus wenchowensis under acute temperature changes.

In recent years, frequent acute temperature changes have posed a serious threat to the physiology and survival of fish. This study utilized RNA-Seq technology to analyze the transcriptional dynamics in the muscle tissues of Acrossocheilus wenchowensis under various acute temperature conditions (16◦C, 20◦C, 24◦C, 28◦C and 32◦C). Through comprehensive analysis, we identified 11509 differentially expressed genes (DEGs), a gene set (profiles 19) that was significantly up-regulated with increasing temperature, and two weighted gene co-expression network analysis (WGCNA) modules that were significantly correlated with acute temperature changes. Furthermore, we identified 28 transcription factors that are pivotal in oxidative stress and energy metabolism under acute temperature changes. Our results showed that, compared to the control group (24°C), KEGG functional enrichment analysis revealed significant enrichment of DEGs in the cell cycle, DNA replication, and p53 signaling pathway, with an overall trend of suppressed expression. This indicates that maintaining cell stability and reducing cell damage is an effective adaptive mechanism for A. wenchowensis to cope with acute temperature changes. Through STEM analysis and the black WGCNA module associated with high-temperature stress, we identified significant up-regulation of pathways and hub genes related to energy metabolism including oxidative phosphorylation, TCA cycle, purine metabolism, and glutathione metabolism, as well as the central roles of signal transduction pathways such as MAPK signaling pathway and AMPK signaling pathway, which synergistically regulate energy production. Under acute low-temperature stress, the turquoise WGCNA module highlighted significant up-regulation of hub genes associated with Ribosomal and Spliceosomal pathways related to protein synthesis and processing, as well as activation of calcium signaling pathways, which plays an important role in maintaining cellular function during low-temperature adaptation. These findings provide a critical theoretical and molecular basis for the adaptation of eurythermal fish to rapid temperature changes.

Hepatic BMAL1 and HIF1α regulate a time-dependent hypoxic response and prevent hepatopulmonary-like syndrome.

The transcriptional response to hypoxia is temporally regulated, yet the molecular underpinnings and physiological implications are unknown. We examined the roles of hepatic Bmal1 and Hif1α in the circadian response to hypoxia in mice. We found that the majority of the transcriptional response to hypoxia is dependent on either Bmal1 or Hif1α, through shared and distinct roles that are daytime determined. We further show that hypoxia-inducible factor (HIF)1α accumulation upon hypoxia is temporally regulated and Bmal1 dependent. Unexpectedly, mice lacking both hepatic Bmal1 and Hif1α are hypoxemic and exhibit increased mortality upon hypoxic exposure in a daytime-dependent manner. These mice display mild liver dysfunction with pulmonary vasodilation likely due to extracellular signaling regulated kinase (ERK) activation, endothelial nitric oxide synthase, and nitric oxide accumulation in lungs, suggestive of hepatopulmonary syndrome. Our findings indicate that hepatic BMAL1 and HIF1α are key time-dependent regulators of the hypoxic response and can provide molecular insights into the pathophysiology of hepatopulmonary syndrome.

The underlying mechanisms by which boron mitigated copper toxicity in Citrus sinensis leaves revealed by integrated analysis of transcriptome, metabolome, and physiology.

Both copper (Cu) excess and boron (B) deficiency are often observed in some citrus orchard soils. The molecular mechanisms by which B alleviates excessive Cu in citrus are poorly understood. Seedlings of sweet orange (Citrus sinensis (L.) Osbeck cv. Xuegan) were treated with 0.5 (Cu0.5) or 350 (Cu350 or Cu excess) µM CuCl2 and 2.5 (B2.5) or 25 (B25) µM HBO3 for 24 weeks. Thereafter, this study examined the effects of Cu and B treatments on gene expression levels revealed by RNA-Seq, metabolite profiles revealed by a widely targeted metabolome, and related physiological parameters in leaves. Cu350 upregulated 564 genes and 170 metabolites, and downregulated 598 genes and 58 metabolites in leaves of 2.5 µM B-treated seedlings (LB2.5), but it only upregulated 281 genes and 100 metabolites, and downregulated 136 genes and 40 metabolites in leaves of 25 µM B-treated seedlings (LB25). Cu350 decreased the concentrations of sucrose and total soluble sugars, and increased the concentrations of starch, glucose, fructose, and total nonstructural carbohydrates (TNC) in LB2.5, but it only increased the glucose concentration in LB25. Further analysis demonstrated that B addition reduced the oxidative damage and alterations in primary and secondary metabolisms caused by Cu350; and alleviated the impairment of Cu350 to photosynthesis and cell wall metabolism, thus improving leaf growth. LB2.5 exhibited some adaptive responses to Cu350 to meet the increasing need for the dissipation of excessive excitation energy (EEE) and the detoxification of reactive oxygen species (reactive aldehydes) and Cu. Cu350 increased photorespiration, xanthophyll cycle-dependent thermal dissipation, nonstructural carbohydrate accumulation, and secondary metabolite biosynthesis and abundances; and upregulated tryptophan metabolism and related metabolite abundances, and some antioxidant-related gene expression, and some antioxidant abundances. Additionally, this study identified some metabolic pathways, metabolites, and genes that might lead to Cu tolerance in leaves.

Screening of reliable reference genes for the normalization of RT-qPCR in chicken liver tissues and LMH cells.

The liver plays a vital role in lipid synthesis and metabolism in poultry. To study the functional genes more effectively, it is essential to screen of reliable reference genes in the chicken liver, including females, males, embryos, as well as the Leghorn Male Hepatoma (LMH) cell line. Traditional reference gene screening involves selecting commonly used housekeeping genes (HKGs) for RT-qPCR experiments and using different algorithms to identify the most stable ones. However, this approach is limited in selecting the best reference gene from a small pool of HKGs. High-throughput sequencing technology may offer a solution to this limitation. This study aimed to identify the most consistently expressed genes by utilizing multiple published RNA-seq data of chicken liver and LMH cells. Subsequently, the stability of the newly identified reference genes was assessed in comparison to previously validated stable poultry liver expressed reference genes and the commonly employed HKGs using RT-qPCR. The findings indicated that there is a higher degree of similarity in stable expression genes between female and male liver (such as LSM14A and CDC40). In embryonic liver, the optimal new reference genes were SUDS3, TRIM33, and ERAL1. For LMH cells, the optimal new reference genes were ALDH9A1, UGGT1, and C21H1orf174. However, it is noteworthy that most HKGs did not exhibit stable expression across multiple samples, indicating potential instability under diverse conditions. Furthermore, RT-qPCR experiments proved that the stable expression genes identified from RNA-seq data outperformed commonly used HKGs and certain validated reference genes specific to poultry liver. Over all, this study successfully identified new stable reference genes in chicken liver and LMH cells using RNA-seq data, offering researchers a wider range of reference gene options for RT-qPCR in diverse situations.

Combined transcriptome and metabolome analysis identifies triterpenoid-induced defense responses in Myzus persicae Sülzer-infested peach.

The defense response of peach (Prunus persica) to insect attack involves changes in gene expression and metabolites. Piercing/sucking insects such as green peach aphid cause direct damage by obtaining phloem nutrients and indirect damage by spreading plant viruses. To investigate the response of peach trees to aphids, the leaf transcriptome and metabolome of two genotypes with different sensitivities to green peach aphid (GPA, Myzus persicae) were studied. The transcriptome analysis of infected peach leaves showed two different response patterns. The gene expression of aphid-susceptible peach plants infected by aphids was more similar to that of the control plants, while the gene expression of aphid-resistant peach plants infected by aphids showed strongly induced changes in gene expression compared with the response in the control plants. Furthermore, gene transcripts in defense-related pathways, including plant-pathogen interaction, MAPK signaling, and several metabolic pathways, were more strongly enriched upon aphid infestation. Untargeted secondary metabolite profiling confirmed that aphid treatment induced larger changes in aphid-resistant peaches than in aphid-susceptible peaches. Consistent with transcriptomic alterations, nine triterpenoids showed extremely significant GPA-induced accumulation in aphid-resistant peaches, whereas triterpenoid abundance remained predominantly unchanged or undetected in aphid- susceptible peaches. Furthermore, some types of transcription factors (including WRKYs, ERFs, NACs, etc.) were more strongly induced upon GPA infestation in aphid-resistant peaches but not in aphid-susceptible peaches. Aphid feeding-dependent transcriptome and metabolite profiles provide the foundation for understanding the molecular mechanisms underlying the response of peach to aphid infestation. These results suggested that accumulation of specialized triterpenoids and the corresponding pathway transcripts may play a key role in peach GPA resistance.