SERE: single-parameter quality control and sample comparison for RNA-Seq

Comparative transcriptome analysis provides comprehensive insight into the molecular mechanisms of heat adaption in Plutella xylostella

Plutella xylostella is one of the most destructive pests for cruciferous vegetables, and is adaptability to different environmental stressors. However, we still know little about the molecular mechanisms of how P. xylostella adapt to thermal stress. Here, the comparative transcriptome analysis was conducted from the samples of control (27 °C, CK) and heat treatment (40 °C, 40 T) P. xylostella. The results showed 1253 genes were differentially expressed, with 624 and 629 genes up- and down-regulated respectively. The annotation analysis demonstrated that "Energy production and conversion", "Protein processing in endoplasmic reticulum", "Peroxisome" and "Tyrosine metabolism" pathways were significantly enriched. Additionally, we found the expression levels of heat shock protein genes (Hsps), cuticle related genes and mitochondrial genes were significantly up-regulated in 40 T insects, suggesting their vital roles in improving adaption to heat stress. Importantly, the SOD activity and MDA content of P. xylostella were both identified to be increased under high temperature stress, indicating the elevated antioxidant reactions might be involved in response to heat stress. In conclusion, the present study offered us an overview of gene expression changes after 40 °C treatments, and found some critical pathways and genes of P. xylostella might play the critical roles in resisting heat stress.

Transcriptomic and western blot characterisation of the human CLEFF4 clone, a new rapid cell line replacement for the Caco2 model

The CLEFF4 sub clone from stock late passage Caco2 cells has a unique property of being able to develop polarised cell monolayers with high P-gp expression and tight junctions much quicker than the original cell line. Instead of being useful for transport studies 21-24 days after initiating culture, the CLEFF4 cell line matures in 5-6 days with tight junctions surpassing that of 3 week old Caco2 cells in that time frame [1]. This has enabled the CLEFF4 cell line to provide measures of apparent permeability for potential drug candidates, so important for pre-clinical drug development, 4 times faster than the original cell line. RNA samples were collected and analysed at days 4 and 7 of culture over a 3 year period and had full RNA transcriptome analysed by the ranaseq.eu open bioinformatics platform. Protein was also collected from day 4 to day 22 of culture. Differential expression data from the FASTQ files have shown significant differences in expression in multiple genes involved with drug efflux, tight junctions, phase 2 metabolism and growth factors, which have been confirmed from protein determination that may hold the key to understanding accelerated human cell maturation. These gene expression results may be significant for other tissues beyond the gastrointestinal tract, and potentially for accelerated cell growth for the new field of laboratory grown tissues for organ replacement. The data also confirms the different genetic expression in CLEFF4 cells compared to Caco2 and the stable nature of the different expression over many years.

Computational analysis of peripheral blood RNA sequencing data unravels disrupted immune patterns in Alzheimer's disease

The central nervous system (CNS) and the immune system collectively coordinate cellular functionalities, sharing common developmental mechanisms. Immunity-related molecules exert an influence on brain development, challenging the conventional view of the brain as immune-privileged. Chronic inflammation emerges as a key player in the pathophysiology of Alzheimer's disease (AD), with increased stress contributing to the disease progression and potentially exacerbating existing symptoms. In this study, the most significant gene signatures from selected RNA-sequencing (RNA-seq) data from AD patients and healthy individuals were obtained and a functional analysis and biological interpretation was conducted, including network and pathway enrichment analysis. Important evidence was reported, such as enrichment in immune system responses and antigen processes, as well as positive regulation of T-cell mediated cytotoxicity and endogenous and exogenous peptide antigen, thus indicating neuroinflammation and immune response participation in disease progression. These findings suggest a disturbance in the immune infiltration of the peripheral immune environment, providing new challenges to explore key biological processes from a molecular perspective that strongly participate in AD development.

Changes in transcriptome regulations of a marine rotifer Brachionus plicatilis under methylmercury stress

Mercury (Hg), a heavy metal pollutant worldwide, can be transformed into methylmercury (MeHg) by various aquatic microorganisms in water, thus accumulating along the aquatic food chain and posing a particular challenge to human health. Zooplankton plays a crucial role in aquatic ecosystems and serves as a major component of the food chain. To evaluate the effects of MeHg on the rotifer Brachionus plicatilis and reveal the underlying mechanism of these effects, we exposed B. plicatilis to MeHg by either direct immersion or by feeding with MeHg-poisoned Chlorella pyrenoidesa, respectively, and conducted a transcriptomic analysis. The results showed that B. plicatilis directly exposed to MeHg by immersion showed significant enrichment of the glutathione metabolism pathway for detoxification of MeHg. In addition, the exposure to MeHg by feeding induced a significant enrichment of lysosome and notch signaling pathways of rotifers, supporting the hypothesis that MeHg can induce autophagy dysfunction in cells and disturb the nervous system of rotifers. In two different routes of MeHg exposure, the pathway of cytochrome P450 in rotifers showed significant enrichment for resisting MeHg toxicity. Our results suggest further studies on the potential mechanism and biological responses of MeHg toxicity in other links of the aquatic food chain.

Non-Coding RNAs Potentially Involved in Pyrethroid Resistance of Anopheles funestus Population in Western Kenya

Backgrounds

The resurgence of Anopheles funestus, a dominant vector of human malaria in western Kenya was partly attributed to insecticide resistance. However, evidence on the molecular basis of pyrethroid resistance in western Kenya is limited. Noncoding RNAs (ncRNAs) form a vast class of RNAs that do not code for proteins and are ubiquitous in the insect genome. Here, we demonstrated that multiple ncRNAs could play a potential role in An. funestusresistance to pyrethroid in western Kenya.

Materials and Methods

Anopheles funestus mosquitoes were sampled by aspiration methods in Bungoma, Teso, Siaya, Port Victoria and Kombewa in western Kenya. The F1 progenies were exposed to deltamethrin (0.05%), permethrin (0.75%), DDT (4%) and pirimiphos-methyl (0.25%) following WHO test guidelines. A synergist assay using piperonyl butoxide (PBO) (4%) was conducted to determine cytochrome P450s' role in pyrethroid resistance. RNA-seq was conducted on a combined pool of specimens that were resistant and unexposed, and the results were compared with those of the FANG susceptible strain. This approach aimed to uncover the molecular mechanisms underlying pyrethroid resistance.

Results

Pyrethroid resistance was observed in all the sites with an average mortality rate of 57.6%. Port Victoria had the highest level of resistance to permethrin (MR=53%) and deltamethrin (MR=11%) pyrethroids. Teso had the lowest level of resistance to permethrin (MR=70%) and deltamethrin (MR=87%). Resistance to DDT was observed only in Kombewa (MR=89%) and Port Victoria (MR=85%). A full susceptibility to P-methyl (0.25%) was observed in all the sites. PBO synergist assay revealed high susceptibility (>98%) to the pyrethroids in all the sites except for Port Victoria (MR=96%, n=100). Whole transcriptomic analysis showed that most of the gene families associated with pyrethroid resistance comprised non-coding RNAs (67%), followed by imipenemase (10%), cytochrome P450s (6%), cuticular proteins (5%), olfactory proteins (4%), glutathione S-transferases (3%), UDP-glycosyltransferases (2%), ATP-binding cassettes (2%) and carboxylesterases(1%).

Conclusions

This study unveils the molecular basis of insecticide resistance in An. funestus in western Kenya, highlighting for the first time the potential role of non-coding RNAs in pyrethroid resistance. Targeting non-coding RNAs for intervention development could help in insecticide resistance management.

The Dsup coordinates grain development and abiotic stress in rice

DNA damage is a serious threat to all living organisms and may be induced by environmental stressors. Previous studies have revealed that the tardigrade (Ramazzotius varieornatus) DNA damage suppressor protein Dsup has protective effects in human cells and tobacco. However, whether Dsup provides radiation damage protection more widely in crops is unclear. To explore the effects of Dsup in other crops, stable Dsup overexpression lines through Agrobacterium-mediated transformation were generated and their agronomic traits were deeply investigated. In this study, the overexpression of Dsup not only enhanced the DNA damage resistance at the seeds and seedlings stages, they also exhibited grain size enlargement and starch granule structure and cell size alteration by the scanning electron microscopy observation. Notably, the RNA-seq revealed that the Dsup plants increased radiation-related and abiotic stress-related gene expression in comparison to wild types, suggesting that Dsup is capable to coordinate normal growth and abiotic stress resistance in rice. Immunoprecipitation enrichment with liquid chromatography-tandem mass spectrometry (IP-LC-MS) assays uncovered 21 proteins preferably interacting with Dsup in plants, suggesting that Dsup binds to transcription and translation related proteins to regulate the homeostasis between DNA protection and plant development. In conclusion, our data provide a detailed agronomic analysis of Dsup plants and potential mechanisms of Dsup function in crops. Our findings provide novel insights for the breeding of crop radiation resistance.

Physiological, Cytological, and Transcriptomic Analysis of Magnesium Protoporphyrin IX Methyltransferase Mutant Reveal Complex Genetic Regulatory Network Linking Chlorophyll Synthesis and Chloroplast Development in Rice

Functional defects in key genes for chlorophyll synthesis usually cause abnormal chloroplast development, but the genetic regulatory network for these key genes in regulating chloroplast development is still unclear. Magnesium protoporphyrin IX methyltransferase (ChlM) is a key rate-limiting enzyme in the process of chlorophyll synthesis. Physiological analysis showed that the chlorophyll and carotenoid contents were significantly decreased in the chlm mutant. Transmission electron microscopy demonstrated that the chloroplasts of the chlm mutant were not well developed, with poor, loose, and indistinct thylakoid membranes. Hormone content analysis found that jasmonic acid, salicylic acid, and auxin accumulated in the mutant. A comparative transcriptome profiling identified 1534 differentially expressed genes (DEGs) between chlm and the wild type, including 876 up-regulated genes and 658 down-regulated genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that these DEGs were highly involved in chlorophyll metabolism, chloroplast development, and photosynthesis. Protein-protein interaction network analysis found that protein translation played an essential role in the ChlM gene-regulated process. Specifically, 62 and 6 DEGs were annotated to regulate chlorophyll and carotenoid metabolism, respectively; 278 DEGs were predicted to be involved in regulating chloroplast development; 59 DEGs were found to regulate hormone regulatory pathways; 192 DEGs were annotated to regulate signal pathways; and 49 DEGs were putatively identified as transcription factors. Dozens of these genes have been well studied and reported to play essential roles in chlorophyll accumulation or chloroplast development, providing direct evidence for the reliability of the role of the identified DEGs. These findings suggest that chlorophyll synthesis and chloroplast development are actively regulated by the ChlM gene. And it is suggested that hormones, signal pathways, and transcription regulation were all involved in these regulation processes. The accuracy of transcriptome data was validated by quantitative real-time PCR (qRT-PCR) analysis. This study reveals a complex genetic regulatory network of the ChlM gene regulating chlorophyll synthesis and chloroplast development. The ChlM gene's role in retrograde signaling was discussed. Jasmonic acid, salicylic acid, or their derivatives in a certain unknown state were proposed as retrograde signaling molecules in one of the signaling pathways from the chloroplast to nucleus.

Pyramiding Breeding of Low-Glutelin-Content Indica Rice with Good Quality and Resistance

Low-glutelin-content rice, a type of functional rice with glutelin levels below 4%, is an essential dietary supplement for chronic kidney disease (CKD) patients. Developing low-glutelin-content rice varieties is crucial to catering to the growing CKD population. In this study, we aimed to create a new low-glutelin indica rice variety with excellent agronomic traits. To achieve this, we employed a combination of molecular-marker-assisted selection and traditional breeding techniques. The cultivars W3660, Wushansimiao (WSSM), and Nantaixiangzhan (NTXZ) were crossbred, incorporating the Lgc-1, Pi-2, Xa23, and fgr alleles into a single line. The result of this breeding effort was "Yishenxiangsimiao", a new indica rice variety that inherits the desirable characteristics of its parent lines. Yishenxiangsimiao (YSXSM) possesses not only a low glutelin content but also dual resistance to blast and bacterial blight (BB). It exhibits high-quality grains with a fragrant aroma. This new low-glutelin indica cultivar not only ensures a stable food supply for CKD patients but also serves as a healthy dietary option for the general public. We also performed RNA-seq of these rice varieties to investigate their internal gene expression differences. The YSXSM exhibited a higher biotic-resistance gene expression in comparison to NTXZ. In summary, we successfully developed a novel low-glutelin indica rice variety, "Yishenxiangsimiao", with superior agronomic traits. This rice variety addresses the dietary needs of CKD patients and offers a nutritious choice for all consumers.

Genome-Wide Identification of the CDPK Gene Family and Their Involvement in Taproot Cracking in Radish

Taproot cracking, a severe and common physiological disorder, markedly reduces radish yield and commercial value. Calcium-dependent protein kinase (CDPK) plays a pivotal role in various plant developmental processes; however, its function in radish taproot cracking remains largely unknown. Here, 37 RsCDPK gene members were identified from the long-read radish genome "QZ-16". Phylogenetic analysis revealed that the CDPK members in radish, tomato, and Arabidopsis were clustered into four groups. Additionally, synteny analysis identified 13 segmental duplication events in the RsCDPK genes. Analysis of paraffin-embedded sections showed that the density and arrangement of fleshy taproot cortex cells are important factors that affect radish cracking. Transcriptome sequencing of the fleshy taproot cortex revealed 5755 differentially expressed genes (DEGs) (3252 upregulated and 2503 downregulated) between non-cracking radish "HongYun" and cracking radish "505". These DEGs were significantly enriched in plant hormone signal transduction, phenylpropanoid biosynthesis, and plant-pathogen interaction KEGG pathways. Furthermore, when comparing the 37 RsCDPK gene family members and RNA-seq DEGs, we identified six RsCDPK genes related to taproot cracking in radish. Soybean hairy root transformation experiments showed that RsCDPK21 significantly and positively regulates root length development. These findings provide valuable insights into the relationship between radish taproot cracking and RsCDPK gene function.

Hormonal Regulation and Stimulation Response of Jatropha curcas L. Homolog Overexpression on Tobacco Leaf Growth by Transcriptome Analysis

The Flowering locus T (FT) gene encodes the florigen protein, which primarily regulates the flowering time in plants. Recent studies have shown that FT genes also significantly affect plant growth and development. The FT gene overexpression in plants promotes flowering and suppresses leaf and stem development. This study aimed to conduct a transcriptome analysis to investigate the multiple effects of Jatropha curcas L. homolog (JcFT) overexpression on leaf growth in tobacco plants. The findings revealed that JcFT overexpression affected various biological processes during leaf development, including plant hormone levels and signal transduction, lipid oxidation metabolism, terpenoid metabolism, and the jasmonic-acid-mediated signaling pathway. These results suggested that the effects of FT overexpression in plants were complex and multifaceted, and the combination of these factors might contribute to a reduction in the leaf size. This study comprehensively analyzed the effects of JcFT on leaf development at the transcriptome level and provided new insights into the function of FT and its homologous genes.

Large-scale differentiation of iPSC-derived motor neurons from ALS and control subjects

Using induced pluripotent stem cells (iPSCs) to understand the mechanisms of neurological disease holds great promise; however, there is a lack of well-curated lines from a large array of participants. Answer ALS has generated over 1,000 iPSC lines from control and amyotrophic lateral sclerosis (ALS) patients along with clinical and whole-genome sequencing data. The current report summarizes cell marker and gene expression in motor neuron cultures derived from 92 healthy control and 341 ALS participants using a 32-day differentiation protocol. This is the largest set of iPSCs to be differentiated into motor neurons, and characterization suggests that cell composition and sex are significant sources of variability that need to be carefully controlled for in future studies. These data are reported as a resource for the scientific community that will utilize Answer ALS data for disease modeling using a wider array of omics being made available for these samples.

Sucrose promotes branch-thorn occurrence of Lycium ruthenicum through dual effects of energy and signal

Lycium ruthenicum is an important eco-economic thorny shrub. In this study, L. ruthenicum plants of a clone showed two types of 'less leaves without thorn' and 'more leaves with thorns' under the same condition after transplanting. Microscopic observation revealed that apical buds of the thornless (Thless) and thorny (Thorny) branches should be selected as materials for further study. RNA-Seq analysis showed that KEGG pathway of Starch and sucrose metabolism and DEGs of SUT13, SUS, TPP and TPS were significantly up-regulated in the Thorny. The results of qRT-PCR confirmed the accuracy and credibility of the RNA-Seq. The content of sucrose in the Thorny was significantly higher than that in the Thless, but the content of trehalose-6-phosphate was opposite. Leaf-clipping treatments reduced sucrose content and inhibited the occurrence/development of branch-thorns; exogenous sucrose of 16 g/L significantly promoted the occurrence and growth of branch-thorns and the promotion effects were significantly higher than those treated with non-metabolizable sucrose analogs (isomaltolose, melitose). These findings suggested that sucrose might play a dual role of energy and signal in the occurrence of branch-thorns. Higher sucrose supply in apical buds from more leaves promoted the occurrence of branch-thorns via lower content of trehalose-6-phosphate and higher expression levels of SUS, TPP and TPS, whereas less leaves inhibited the occurrence. Molecular hypothesis model of leaf number/sucrose supply regulating the occurrence of branch-thorns in L. ruthenicum was established in the study, which provides foundation for breeding both thornless L. ruthenicum and thornless types of other species.

Transcriptome changes of Apis mellifera female embryos with fem gene knockout by CRISPR/Cas9

The sex of honey bees is decided by a regulatory cascade comprising of csd, fem and Amdsx. In order to further identify other genes involved in sex determination and differentiation of honey bees in the early stages of embryo development, the CRISPR/Cas9 method was used to knock out fem gene in the embryonic stage of diploid western honey bees, and RNA-seq was used to analyze gene expression changes in the embryo after fem knockout. Finally, we found that the bees had undergone gender changes due to fem knockout. A total of 155 differentially expressed genes (DEGs) were obtained, with 48 up-regulated and 107 down-regulated DEGs in the mutant group compared to the control group. Of them, many genes are related to sex development or differentiation. In addition, 1502 differentially expressed alternative splicing events (DEASEs) related to 1011 genes, including the main honey bee sex-determining genes csd, tra2, fem, and Amdsx, were identified between the mutant group and control group, indicating that fem regulates alternative splicing of a large number of downstream genes. Our results provide valuable clues for further investigating the molecular mechanism of sex determination and differentiation in honey bees.

Cyclo-glycylproline attenuates hydrogen peroxide-induced cellular damage mediated by the MDM2-p53 pathway in human neural stem cells

Cyclo-glycylproline (cGP), a cyclic dipeptide containing a condensation bond between glycine and proline, is produced by the cyclization of the N-terminal tripeptide of insulin-like growth factor-1. Previous studies have shown that cGP administration exerts a neuroprotective effect and enhances the regenerative ability in rats with ischemic brain injury. The efficacy of cGP is medicated by regulating the bioavailability of insulin-like growth factor-1 (IGF-1), however, the molecular mechanisms underlying the neuroprotective effects of cGP on brain damage remains to be elucidated. In the current study, we investigated the cGP-mediated molecular mechanism in human fetal neural stem cells (hfNSCs) exposed to oxidative stress, which is a key factor affecting the development of several brain diseases, including traumatic brain injury and Parkinson's disease. We found that cGP treatment attenuated oxidative stress-induced cell death in cultured hfNSCs in a dose-dependent manner. Transcriptome analysis revealed that under oxidative stress conditions, p53-mediated signaling was activated, accompanied by upregulation of mouse double minute 2 homolog (MDM2), a p53-specific E3 ubiquitin ligase, in cGP-treated hfNSCs. By using a comprehensive protein phosphorylation array, we found that cGP induced the activation of Akt signaling pathway, which enhanced the expression of MDM2, in hfNSCs exposed to oxidative stress. Moreover, the MDM2 inhibitor nutlin-3 inhibited the protective effect of cGP on oxidative stress-induced cell death and apoptosis. Therefore, cGP attenuates oxidative stress-induced cell death mediated by the interplay between IGF-1 signaling and the MDM2-p53 pathway in human NSCs. We revealed the molecular mechanism underlying cGP-induced neuroprotective properties in a model of brain damage.

Basic helix-loop-helix transcription factor PxbHLH02 enhances drought tolerance in Populus (Populus simonii × P. nigra)

The basic helix-loop-helix (bHLH) transcription factors (TFs) are involved in plant morphogenesis and various abiotic and biotic stress responses. However, further exploration is required of drought-responsive bHLH family members and their detailed regulatory mechanisms in Populus. Two bHLH TF genes, PxbHLH01/02, were identified in Populus simonii × P. nigra and cloned. The aim of this study was to examine the role of bHLH TFs in drought tolerance in P. simonii × P. nigra. The results showed that the amino acid sequences of the two genes were homologous to Arabidopsis thaliana UPBEAT1 (AtUPB1) and overexpression of PxbHLH01/02 restored normal root length in the AtUPB1 insertional mutant (upb1-1). The PxbHLH01/02 gene promoter activity analysis suggested that they were involved in stress responses and hormone signaling. Furthermore, Arabidopsis transgenic lines overexpressing PxbHLH01/02 exhibited higher stress tolerance compared with the wild-type. Populus simonii × P. nigra overexpressing PxbHLH02 increased drought tolerance and exhibited higher superoxide dismutase and peroxidase activities, lower H2O2 and malondialdehyde content, and lower relative conductivity. The results of transcriptome sequencing (RNA-seq) and quantitative real-time PCR suggested that the response of PxbHLH02 to drought stress was related to abscisic acid (ABA) signal transduction. Overall, the findings of this study suggest that PxbHLH02 from P. simonii × P. nigra functions as a positive regulator of drought stress responses by regulating stomatal aperture and promoting ABA signal transduction.

Transcriptome analysis and exploration of genes involved in the biosynthesis of secoiridoids in Gentiana rhodantha

Gentiana rhodantha is a medicinally important perennial herb used as traditional Chinese and ethnic medicines. Secoiridoids are one of the major bioactive compounds in G. rhodantha. To better understand the secoiridoid biosynthesis pathway, we generated transcriptome sequences from four organs (root, leaf, stem and flower), followed by the de novo sequence assembly. We verified 8-HGO (8-hydroxygeraniol oxidoreductase), which may encode key enzymes of the secoiridoid biosynthesis by qRT-PCR. The mangiferin, swertiamarin and loganic acid contents in root, stem, leaf, and flower were determined by HPLC. The results showed that there were 47,871 unigenes with an average length of 1,107.38 bp. Among them, 1,422 unigenes were involved in 25 standard secondary metabolism-related pathways in the KEGG database. Furthermore, we found that 1,005 unigenes can be divided into 66 transcription factor (TF) families, with no family members exhibiting significant organ-specificity. There were 54 unigenes in G. rhodantha that encoded 17 key enzymes of the secoiridoid biosynthetic pathway. The qRT-PCR of the 8-HGO and HPLC results showed that the relative expression and the mangiferin, swertiamarin, and loganic acid contents of the aerial parts were higher than in the root. Six types of SSR were identified by SSR analysis of unigenes: mono-nucleoside repeat SSR, di-nucleoside repeat SSR, tri-nucleoside repeat SSR, tetra-nucleoside repeat SSR, penta-nucleoside repeat SSR, and hexa-nucleoside repeat SSR. This report not only enriches the Gentiana transcriptome database but helps further study the function and regulation of active component biosynthesis of G. rhodantha.

Integrated metabolomics and transcriptomics analysis of roots of Bupleurum chinense and B. scorzonerifolium, two sources of medicinal Chaihu

Radix Bupleuri (Chaihu in Chinese) is a traditional Chinese medicine commonly used to treat colds and fevers. The root metabolome and transcriptome of two cultivars of B. chinense (BCYC and BCZC) and one of B. scorzonerifolium (BSHC) were determined and analyzed. Compared with BSHC, 135 and 194 differential metabolites were identified in BCYC and BCZC, respectively, which were mainly fatty acyls, organooxygen metabolites. A total of 163 differential metabolites were obtained between BCYC and BCZC, including phenolic acids and lipids. Compared with BSHC, 6557 and 5621 differential expression genes (DEGs) were found in BCYC and BSHC, respectively, which were annotated into biosynthesis of unsaturated fatty acid and fatty acid metabolism. A total of 4,880 DEGs existed between the two cultivars of B. chinense. The abundance of flavonoids in B. scorzonerifolium was higher than that of B. chinense, with the latter having higher saikosaponin A and saikosaponin D than the former. Pinobanksin was the most major flavonoid which differ between the two cultivars of B. chinense. The expression of chalcone synthase gene was dramatically differential, which had a positive correlation with the biosynthesis of pinobanksin. The present study laid a foundation for further research on biosynthesis of flavonoids and terpenoids of Bupleurum L.

Comparative transcriptomic analysis-based identification of the regulation of foreign proteins with different stabilities expressed in Pichia pastoris

Introduction

The industrial yeast Pichia pastoris is widely used as a cell factory to produce proteins, chemicals and advanced biofuels. We have previously constructed P. pastoris strains that overexpress protein disulfide isomerase (PDI), which is a kind of molecular chaperone that can improve the expression of an exogenous protein when they are co-expressed. Chicken cystatin (cC) is a highly thermostable cysteine protease inhibitor and a homologous protein of human cystatin C (HCC). Wild-type cC and the two mutants, I66Q and ΔW (a truncated cC lacking the á-helix 2) represent proteins with different degrees of stability.

Methods

Wild-type cC, I66Q and ΔW were each overexpressed in P. pastoris without and with the coexpression of PDI and their extracellular levels were determined and compared. Transcriptomic profiling was performed to compare the changes in the main signaling pathways and cell components (other than endoplasmic reticulum quality control system represented by molecular chaperones) in P. pastoris in response to intracellular folding stress caused by the expression of exogenous proteins with different stabilities. Finally, hub genes hunting was also performed.

Results and discussion

The coexpression of PDI was able to increase the extracellular levels of both wild-type cC and the two mutants, indicating that overexpression of PDI could prevent the misfolding of unstable proteins or promote the degradation of the misfolded proteins to some extent. For P. pastoris cells that expressed the I66Q or ΔW mutant, GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses of the common DEGs in these cells revealed a significant upregulation of the genes involved in protein processing, but a significant downregulation of the genes enriched in the Ribosome, TCA and Glycolysis/Gluconeogenesis pathways. Hub genes hunting indicated that the most downregulated ribosome protein, C4QXU7 in this case, might be an important target protein that could be manipulated to increase the expression of foreign proteins, especially proteins with a certain degree of instability.

Conclusion

These findings should shed new light on our understanding of the regulatory mechanism in yeast cells that responds to intracellular folding stress, providing valuable information for the development of a convenient platform that could improve the efficiency of heterologous protein expression in P. pastoris.

Transcriptomics Analysis Reveals a More Refined Regulation Mechanism of Methylation in a Drought-Tolerant Variety of Potato

Whether DNA methylation modification affects the gene transcription and expression of potatoes under drought stress is still unknown. In this study, we used comparative transcriptomics to explore the expression pattern of related genes of the drought-tolerant variety Qingshu 9 (Q) and the drought-sensitive variety Atlantic (A) under drought stress and DNA methylation inhibitor treatment. The results showed that there was a significant difference in the number of DEGs between the two varieties' responses to mannitol and 5-azad C, especially when they were co-treated with two reagents, and the gene expression of Q was more sensitive to mannitol after two hours. Furthermore, we found that these differentially expressed genes (DEGs) were significantly enriched in DNA replication, transcription, translation, carbohydrate metabolism, photosynthesis, signal transduction, and glutathione metabolism. These results indicate that the difference in the background of methylation leads to the difference in drought resistance of the two varieties. The complexity of the DNA methylation of variety Q might be higher than that of variety A, and the method of methylation regulation is more refined. This study systematically expands the understanding of the molecular mechanism wherein DNA methylation regulates the response to drought stress.

Unpacking COVID-19 Systems Biology in Lung and Whole Blood with Transcriptomics and miRNA Regulators

COVID-19 is a systemic disease affecting tissues and organs, including and beyond the lung. Apart from the current pandemic context, we also have vastly inadequate knowledge of consequences of repeated exposures to SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2), the virus causing COVID-19, in multiple organ systems and the whole organism scales when the disease evolves from a pandemic to an endemic state. This calls for a systems biology and systems medicine approach and unpacking the effects of COVID-19 in lung as well as other tissues. We report here original findings from transcriptomics analyses and differentially expressed genes (DEGs) in lung samples from 60 patients and 27 healthy controls, and in whole blood samples from 255 patients and 103 healthy individuals. A total of 11 datasets with RNA-seq transcriptomic data were obtained from the Gene Expression Omnibus and the European Nucleotide Archive. The identified DEGs were used to construct protein interaction and functional networks and to identify related pathways and miRNAs. We found 35 DEGs common between lung and the whole blood, and importantly, 2 novel genes, namely CYP1B1 and TNFAIP6, which have not been previously implicated with COVID-19. We also identified four novel miRNA potential regulators, hsa-mir-192-5p, hsa-mir-221-3p, hsa-mir-4756-3p, and hsa-mir-10a-5p, implicated in lung or other diseases induced by coronaviruses. In summary, these findings offer new molecular leads and insights to unpack COVID-19 systems biology in a whole organism context and might inform future antiviral drug, diagnostics, and vaccine discovery efforts.

Transcriptome Analysis to Identify Responsive Genes under Sublethal Concentration of Bifenazate in the Diamondback Moth, Plutella xylostella (Linnaeus, 1758) (Lepidoptera: Plutellidae)

Bifenazate is a novel acaricide that has been widely used to control spider mites. Interestingly, we found bifenazate had a biological activity against the diamondback moth (Plutella xylostella), one of the most economically important pests on crucifer crops around the world. However, the molecular mechanisms underlying the response of P. xylostella to bifenazate treatment are not clear. In this study, we first estimated the LC₃₀ dose of bifenazate for third-instar P. xylostella larvae. Then, in order to identify genes that respond to the treatment of this insecticide, the comparative transcriptome profiles were used to analyze the gene expression changes in P. xylostella larvae after exposure to LC₃₀ of bifenazate. In total, 757 differentially expressed genes (DEGs) between bifenazate-treated and control P. xylostella larvae were identified, in which 526 and 231 genes were up-regulated and down-regulated, respectively. The further Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that the xenobiotics metabolisms pathway was significantly enriched, with ten detoxifying enzyme genes (four P450s, five glutathione S-transferases (GSTs), and one UDP-Glucuronosyltransferase (UGT)) were up-regulated, and their expression patterns were validated by qRT-PCR as well. Interestingly, the present results showed that 17 cuticular protein (CP) genes were also remarkably up-regulated, including 15 CPR family genes. Additionally, the oxidative phosphorylation pathway was found to be activated with eight mitochondrial genes up-regulated in bifenazate-treated larvae. In contrast, we found some genes that were involved in tyrosine metabolism and purine pathways were down-regulated, indicating these two pathways of bifenazate-exposed larvae were significantly inhibited. In conclusion, the present study would help us to better understand the molecular mechanisms of sublethal doses of bifenazate detoxification and action in P. xylostella.

Characterization of Abscisic Acid and Ethylene in Regulating the White Blush in Fresh-Cut Carrots

The surface of fresh-cut carrots is apt to white blush, however the physiological and molecular mechanism for this process is not yet fully understood. In this study, exogenous abscisic acid (ABA) and ethylene separately promoted and inhibited the white-blush formation after three days after treatment, respectively. Metabolome analysis found that white-blush components mainly consist of p-hydroxyphenyl lignin and guaiacyl lignin. Transcriptome analysis found an increase in the whiteness values was consistent with the higher expression of genes encoding O-methyltransferase, trans-anol O-methyltransferase, bergaptol O-methyltransferase, caffeic acid 3-O-methyltransferase, phenylalanine ammonia-lyase, and ferulate-5-hydroxylase, together with the lower expression of genes encoding cinnamic acid 4-hydroxylase caffeoyl-CoA O-methyltransferase and 5-O-(4-coumaroyl)-D-quinate 3′-monooxygenase. In conclusion, ABA plays an important role in lignin biosynthesis essential to the formation of white blush in fresh-cut carrots. This is the first report that uncovers the physiological and molecular causes of white blush in fresh-cut carrots, providing a basis for white-blush control in fresh-cut carrots.

Integrated transcriptional and metabolomics signature pattern of pigmented wheat to insight the seed pigmentation and other associated features

Anthocyanin biosynthesis in plants is complex, especially in a polyploid monocot wheat plant. Using whole-genome sequencing, transcriptomics, and LC-MS/MS, we investigated anthocyanin pigmentation patterns in (black, blue, and purple) colored wheat seeds. According to differential gene expression profiling, 2AS-MYC, 7DL-MYB, and WD40 regulatory genes control purple pericarp coloration, 4DL-MYC, 2AS-MYC, 7DL-MYB, WD40 control blue aleurone coloration, and 4DL-MYC, 7DL-MYB, WD40 controls black aleurone color. We hypothesized that at least one MYC and MYB isoform is sufficient to regulate the anthocyanin synthesis in pericarp or aleurone. Transcriptomics and metabolomics revealed that the purple pericarp trait is associated with acylated anthocyanins compared to blue aleurone. Based upon the reduced expressions of the genes belonging to the 4D, SSR molecular marker mapping, variant calling using genome sequencing, and IGV browser gene structure visualization, it was inferred that the advanced black and blue wheat lines were substitution lines (4E{4D}), with very small recombinations. Pericarp anthocyanin pigmentation is controlled by a mutation in chromosome 2AS of purple wheat, and environmental variations influence pigmented pericarp trait. The expression patterns of anthocyanin structural and other genes varied in different colored wheat, corroborating differences in agronomical metrics. Ovate seed shape trait in black and blue wheat dragged with 4E chromosome.

SALT AND ABA RESPONSE ERF1 improves seed germination and salt tolerance by repressing ABA signaling in rice

Rice (Oryza sativa) germination and seedling establishment, particularly in increasingly saline soils, are critical to ensure successful crop yields. Seed vigor, which determines germination and seedling growth, is a complex trait affected by exogenous (environmental) and endogenous (hormonal) factors. Here, we used genetic and biochemical analyses to uncover the role of an APETALA2-type transcription factor, SALT AND ABA RESPONSE ERF1 (OsSAE1), as a positive regulator of seed germination and salt tolerance in rice by repressing the expression of ABSCISIC ACID-INSENSITIVE5 (OsABI5). ossae1 knockout lines exhibited delayed seed germination, enhanced sensitivity to abscisic acid (ABA) during germination and in early seedling growth, and reduced seedling salt tolerance. OsSAE1 overexpression lines exhibited the converse phenotype, with increased seed germination and salt tolerance. In vivo and in vitro assays indicated that OsSAE1 binds directly to the promoter of OsABI5, a major downstream component of the ABA signaling pathway and acts as a major regulator of seed germination and stress response. Genetic analyses revealed that OsABI5-mediated ABA signaling functions downstream of OsSAE1. This study provides important insights into OsSAE1 regulation of seed vigor and salt tolerance and facilitates the practical use of OsSAE1 in breeding salt-tolerant varieties suitable for direct seeding cultivation.

OsQHB Improves Salt Tolerance by Scavenging Reactive Oxygen Species in Rice

Soil salinity is a major environmental stress that restricts the growth and yield of crops. Mining the key genes involved in the balance of rice salt tolerance and yield will be extremely important for us to cultivate salt-tolerance rice varieties. In this study, we report a WUSCHEL-related homeobox (WOX) gene, quiescent-center-specific homeobox (OsQHB), positively regulates yield-related traits and negatively regulates salt tolerance in rice. Mutation in OsQHB led to a decrease in plant height, tiller number, panicle length, grain length and grain width, and an increase in salt tolerance. Transcriptome and qPCR analysis showed that reactive oxygen species (ROS) scavenging-related genes were regulated by OsQHB. Moreover, the osqhb mutants have higher ROS-scavenging enzymes activities and lower accumulation of ROS and malondialdehyde (MDA) under salt stress. Thus, our findings provide new insights into the role of rice WOX gene family in rice development and salt tolerance, and suggest that OsQHB is a valuable target for improving rice production in environments characterized by salt stress.

Comparative transcriptomics provide new insights into the mechanisms by which foliar silicon alleviates the effects of cadmium exposure in rice

Silicon (Si) has been shown to alleviate Cd stress in rice. Here, we investigated the beneficial effects of foliar Si in an indica rice Huanghuazhan (HHZ). Our results showed that foliar Si increases the dry weight and decreases Cd translocation in Cd-exposed rice at the grain-filling stage only, implying that the filling stage is critical for foliar Si to reduce Cd accumulation. We also investigated the transcriptomics in flag leaves (FLs), spikelets (SPs), and node Is (NIs) of Cd-exposed HHZ after foliar Si application at the filling stage. Importantly, the gene expression profiles associated with the Si-mediated alleviation of Cd stress were tissue specific, while shared pathways were mediated by Si in Cd-exposed rice tissues. Furthermore, after the Si treatment of Cd-exposed rice, the ATP-binding cassette (ABC)-transporters were mostly upregulated in FL and SP, while the bivalent cation transporters were mostly downregulated in FL and NI, possibly helping to reduce Cd accumulation. The genes associated with essential nutrient transporters, carbohydrate and secondary metabolite biosynthesis, and cytochrome oxidase activity were mostly upregulated in Cd-exposed FL and SP, which may help to alleviate oxidative stress and improve plant growth under Cd exposure. Interestingly, genes responsible for signal transduction were negatively regulated in FL, but positively regulated in SP, by foliar Si. Our results provide transcriptomic evidence that foliar Si plays an active role in alleviating the effects of Cd exposure in rice. In particular, foliar Si may alter the expression pattern of genes associated with transport, biosynthesis and metabolism, and oxidation reduction.

Transcriptome Sequencing and Screening of Anthocyanin-Related Genes in the Leaves of Acer truncatum Bunge

Acer truncatum Bunge is generally used as an ornamental tree because of its autumn leaves, although the viewing period is short-approximately 7-15 days. Color improvement of ornamental trees has consistently been an important research topic because color partially determines the value of the commodity; however, a lack of genomic data have limited the progress of molecular breeding research in this area. The purposes of this study were to obtain a transcriptome database for A. truncatum, screen anthocyanin biosynthesis-related genes, and reveal the mechanisms underlying leaf color transformation to provide a basis for increasing the viewing period or breeding cultivars that display red leaves throughout the growing season via gene regulation. In this study, although the use of an Illumina HiSeq 2000 platform and systematic bioinformatics analysis using both young and mature leaves as experimental materials, 233,912,882 clean reads were generated and 121,287 unique transcripts were retrieved. We selected 16 color-related genes (from the transcriptome results) for qRT-PCR to validate the results, and the expression trends of the selected genes were largely consistent with the transcriptome analysis results, with a consistency of 0.875. According to the results of the transcriptome analysis, the validation, and previous studies, we obtained sequences of genes related to anthocyanins, including CHS, CHI, ANS, UFGT, UGT75c1, DFR, BZ1, F3H, F3'H, LAR, ANR, FLS, and those of several transcription factors, including MYB1, BHLH, and WD40. Verifying specific regulation by one or several of these genes in the control of leaf color requires further research. The acquisition of transcriptomic information, especially information concerning anthocyanin biosynthesis-related genes and their base sequences, can provide a theoretical basis for the study of the molecular mechanisms determining changes in leaf color in Acer and is of great importance to the breeding of new cultivars.

Comprehensive Effects of Flowering Locus T-Mediated Stem Growth in Tobacco

In flowering plants, Flowering locus T (FT) encodes a major florigen. It is a key flowering hormone in controlling flowering time and has a wide range of effects on plant development. Although the mechanism by which FT promotes flowering is currently clearly understood, comprehensive effects of the FT gene on plant growth have not been evaluated. Therefore, the effects of FT on vegetative growth need to be explored for a complete understanding of the molecular functions of the FT gene. In this study, the Jatropha curcas L. FT gene was overexpressed in tobacco (JcFT^OE) in order to discover multiple aspects and related mechanisms of how the FT gene affects plant development. In JcFT^OE plants, root, stem, and leaf development was strongly affected. Stem tissues were selected for further transcriptome analysis. In JcFT^OE plants, stem growth was affected because of changes in the nucleus, cytoplasm, and cell wall. In the nucleus of JcFT^OE plants, the primary effect was to weaken all aspects of DNA replication, which ultimately affected the cell cycle and cell division. The number of stem cells decreased significantly in JcFT^OE plants, which decreased the thickness and height of tobacco stems. In the cell wall of JcFT^OE plants, hemicellulose and cellulose contents increased, with the increase in hemicellulose associated with up-regulation of xylan synthase-related genes expression. In the cytoplasm of JcFT^OE plants, the primary effects were on biogenesis of ribonucleoprotein complexes, photosynthesis, carbohydrate biosynthesis, and the cytoskeleton. In addition, in the cytoplasm of JcFT^OE plants, there were changes in certain factors of the core oscillator, expression of many light-harvesting chlorophyll a/b binding proteins was down-regulated, and expression of fructose 1,6-bisphosphatase genes was up-regulated to increase starch content in tobacco stems. Changes in the xylem and phloem of JcFT^OE plants were also identified, and in particular, xylem development was affected by significant increases in expression of irregular xylem genes.

Identification of genes related to tipburn resistance in Chinese cabbage and preliminary exploration of its molecular mechanism

BACKGROUND

Tipburn, also known as leaf tip necrosis, is a severe issue in Chinese cabbage production. One known cause is that plants are unable to provide adequate Ca²⁺ to rapidly expanding leaves. Bacterial infection is also a contributing factor. Different cultivars have varying degrees of tolerance to tipburn. Two inbred lines of Chinese cabbage were employed as resources in this research.

RESULTS

We determined that the inbred line 'J39290' was the tipburn resistant material and the inbred line 'J95822' was the tipburn sensitive material based on the severity of tipburn, and the integrity of cell membrane structure. Ca²⁺ concentration measurements revealed no significant difference in Ca²⁺ concentration between the two materials inner leaves. Transcriptome sequencing technology was also used to find the differentially expressed genes (DEGs) of 'J95822' and 'J39290', and there was no significant difference in the previously reported Ca²⁺ uptake and transport related genes in the two materials. However, it is evident through DEG screening and classification that 23 genes are highly linked to plant-pathogen interactions, and they encode three different types of proteins: CaM/CML, Rboh, and CDPK. These 23 genes mainly function through Ca²⁺-CaM/CML-CDPK signal pathway based on KEGG pathway analysis, protein interaction prediction, and quantitative real-time PCR (qRT-PCR) of key genes.

CONCLUSIONS

By analyzing the Ca²⁺ concentration in the above two materials, the transcription of previously reported genes related to Ca²⁺ uptake and transport, the functional annotation and KEGG pathway of DEGs, it was found that Ca²⁺ deficiency was not the main cause of tipburn in 'J95822', but was probably caused by bacterial infection. This study lays a theoretical foundation for exploring the molecular mechanism of resistance to tipburn in Chinese cabbage, and has important guiding significance for genetics and breeding.

A comparative transcriptome and proteomics study of post-partum ovarian cycle arrest in yaks (Bos grunniens)

Post-partum ovarian cycle arrest is the main factor affecting yak reproductive efficiency. There are few reports regarding the molecular regulatory mechanism of post-partum oestrus at transcriptome and proteome levels in yaks. Our previous studies focussed on the ovaries of yaks with post-partum ovarian cycle arrest and post-partum oestrus yaks. In this study, RNA sequencing transcriptomic study was combined with quantitative proteomic analyses to identify post-partum ovarian cycle-related genes and proteins. Consequently, 1,149 genes and 24 proteins were found to be up- or downregulated during post-partum oestrus. The analysis of differentially regulated genes identified three gene or protein pairs that were synchronously upregulated and no gene or protein pairs that were synchronously downregulated, suggesting that these upregulated genes may regulate the post-partum ovarian cycle. The functional classification of these differentially expressed genes and proteins indicated their connection with the oocyte meiosis, the oestrogen signalling pathway, the progesterone-mediated oocyte maturation and the gonadotrophin-releasing hormone (GnRH) signalling pathway. In this study, a total of six genes and two proteins involved in the oocyte meiosis, the oestrogen signalling pathway, the progesterone-mediated oocyte maturation and the GnRH signalling pathway were identified. The CSNK1A1, M91_09723, M91_11326, M91_21439, M91_19073, SHC2, Atf6b, M91_03062, HSPCA and calmodulin could regulate oestrus, respectively, in the post-partum so as to control the anoestrus status.

VF-4 and DR-8 Derived from Salted Egg White Inhibit Inflammatory Activity via NF-κB/PI3K-Akt/MAPK Signal Transduction Pathways in HT-29 Cells Induced by TNF-α

SCOPE

Inflammation is the pathological basis of many chronic diseases, and persistent intestinal inflammation is a key factor in the further development of colon cancer. Egg-derived peptides have been proven to have anti-intestinal inflammation activity. Egg white treated with salt contains a lot of rich protein, whether its peptides have anti-inflammatory activity and how their mechanism of action is still unclear.

METHODS AND RESULTS

In this study, ELISA is used to determine the anti-inflammatory activity of the peptides (VF-4 and DR-8 from salted egg white), and then RNA-seq is used to explore the mechanism of their anti-inflammatory activity, and then verified by western blotting and inhibitors. The results show that VF-4 and DR-8 significantly inhibit TNF-α-induced IL-8 secretion in HT-29 cells in a concentration-dependent manner, and VF-4 show a more significant anti-inflammatory effect than DR-8. The anti-inflammatory mechanism of VF-4 and DR-8 is through inhibiting the activation of Nuclear factor kappa B (NF-κB), mitogen-activated protein kinase (MAPK), and phosphatidylinositol 3' -kinase(PI3K)-Akt pathways, reducing the production of inflammatory mediators.

CONCLUSION

VF-4 and DR-8 have obvious anti-inflammatory activity, which can reduce intestinal inflammation and inhibit its further development into colon cancer.

Transcriptome analysis reveals differing response and tolerance mechanism of EPSPS and GAT genes among transgenic soybeans

BACKGROUND

Glyphosate is a broad-spectrum, non-selective systemic herbicide. Introduction of glyphosate tolerance genes such as EPSPS or detoxification genes such as GAT can confer glyphosate tolerance on plants. Our previous study revealed that co-expression of EPSPS and GAT genes conferred higher glyphosate tolerance without "yellow flashing". However, the plant response to glyphosate at the transcriptional level was not investigated.

METHODS AND RESULTS

To investigate the glyphosate tolerance mechanism, RNA-seq was conducted using four soybean genotypes, including two non-transgenic (NT) soybeans, ZH10 and MD12, and two GM soybeans, HJ698 and ZH10-6. Differentially expressed genes (DEGs) were identified in these soybeans before and after glyphosate treatment. Similar response to glyphosate in the two NT soybeans and the different effects of glyphosate on the two GM soybeans were identified. As treatment time was prolonged, the expression level of some DEGs involved in shikimate biosynthetic pathway and herbicide targeted cross-pathways was increased or declined continuously in NT soybeans, and altered slightly in HJ698. However, the expression level of some DEGs was altered in ZH10-6 at 12 hpt, while similar expression level of some DEGs involved in shikimate biosynthetic pathway and herbicide targeted cross-pathways was observed in ZH10-6 at 0 hpt and 72 hpt. These observations likely explain the higher glyphosate tolerance in ZH10-6 than in HJ698 and NT soybeans.

CONCLUSIONS

These results suggested that GAT and EPSPS genes together play a crucial role in response to glyphosate, the GAT gene may work at the early stage of glyphosate exposure, whereas the EPSPS gene may be activated after the uptake of glyphosate by plants. These findings will provide valuable insight for the molecular basis underlying glyphosate tolerance or glyphosate detoxication.

Salt Stress Promotes Abscisic Acid Accumulation to Affect Cell Proliferation and Expansion of Primary Roots in Rice

The primary root is the basic component of the root system and plays a key role in early seedling growth in rice. Its growth is easily affected by environmental cues, such as salt stress. Abscisic acid (ABA) plays an essential role in root development, but the molecular mechanism underlying ABA-regulated root growth in response to salt stress remains poorly understood. In this study, we report that salt stress inhibits primary root elongation and promotes primary root swelling. Moreover, salt stress induces the expression of ABA-responsive genes and ABA accumulation in the primary root, revealing that ABA plays an essential role in salt-modulated root growth. Transgenic lines of OsSAPK10-OE and OsABIL2-OE, which constitutively express OsSAPK10 or OsABIL2, with enhanced or attenuated ABA signaling, show increased and decreased sensitivity to salt, correspondingly. Microscopic analysis indicates that salt and ABA inhibits cell proliferation and promotes cell expansion in the root apical meristem. Transcriptome analysis showed that ABA induces the expression of EXPANSIN genes. Further investigations indicate that ABA exerts these effects largely through ABA signaling. Thus, our findings deepen our understanding of the role of ABA in controlling primary root growth in response to salt stress, and this knowledge can be used by breeders to cultivate rice varieties suitable for saline–alkali land.

Multiple responses contribute to the enhanced drought tolerance of the autotetraploid Ziziphus jujuba Mill. var. spinosa

Background

Polyploid plants often exhibit enhanced stress tolerance. The underlying physiological and molecular bases of such mechanisms remain elusive. Here, we characterized the drought tolerance of autotetraploid sour jujube at phenotypic, physiological and molecular levels.

Results

The study findings showed that the autotetraploid sour jujube exhibited a superior drought tolerance and enhanced regrowth potential after dehydration in comparison with the diploid counterpart. Under drought stress, more differentially expressed genes (DEGs) were detected in autotetraploid sour jujube and the physiological responses gradually triggered important functions. Through GO enrichment analysis, many DEGs between the diploid and autotetraploid sour jujube after drought-stress exposure were annotated to the oxidation–reduction process, photosystem, DNA binding transcription factor activity and oxidoreductase activity. Six reactive oxygen species scavenging-related genes were specifically differentially expressed and the larger positive fold-changes of the DEGs involved in glutathione metabolism were detected in autotetraploid. Consistently, the lower O²⁻ level and malonaldehyde (MDA) content and higher antioxidant enzymes activity were detected in the autotetraploid under drought-stress conditions. In addition, DEGs in the autotetraploid after stress exposure were significantly enriched in anthocyanin biosynthesis, DNA replication, photosynthesis and plant hormone, including auxin, abscisic acid and gibberellin signal-transduction pathways. Under osmotic stress conditions, genes associated with the synthesis and transport of osmotic regulators including anthocyanin biosynthesis genes were differentially expressed, and the soluble sugar, soluble protein and proline contents were significantly higher in the autotetraploid. The higher chlorophyll content and DEGs enriched in photosynthesis suggest that the photosynthetic system in the autotetraploid was enhanced compared with diploid during drought stress. Moreover, several genes encoding transcription factors (TFs) including GRAS, Bhlh, MYB, WRKY and NAC were induced specifically or to higher levels in the autotetraploid under drought-stress conditions, and hub genes, LOC107403632, LOC107422279, LOC107434947, LOC107412673 and LOC107432609, related to 18 up-regulated transcription factors in the autotetraploid compared with the diploid were identified.

Conclusion

Taken together, multiple responses contribute to the enhanced drought tolerance of autotetraploid sour jujube. This study could provide an important basis for elucidating the mechanism of tolerance variation after the polyploidization of trees.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-021-00633-1.

RNA-seq analysis of gene expression changes in cuticles during the larval-pupal metamorphosis of Plutella xylostella

The diamondback moth, Plutella xylostella (Lepidoptera: Plutellidae) is a holometabolous insect that its cuticles must undergo the significant changes during the larval-pupal metamorphosis development. To elucidate these changes at molecular levels, RNA-seq analysis of cuticles from LLS (later fourth instar larval stage), PPS (prepupal stage) and PS (pupal stage) were performed in P. xylostella. In this paper, a total of 17,710 transcripts were obtained in the larval-pupal transition of P. xylostella, and out of which 2293 (881 up-regulated and 1412 down-regulated) and 2989 transcripts (2062 up-regulated and 927 down-regulated) were identified to be differentially expressed between LLS and PPS, as well as PPS and PS, respectively. The further GO and KEGG analysis of differentially expressed genes (DEGs) revealed that the 'structural constituent of cuticle', 'chitin metabolic process', 'chitin binding', 'tyrosine metabolism' and 'insect hormone biosynthesis' pathways were significantly enriched, indicating these pathways might be involved in the process of larval pupation in P. xylostella. Then, we found some genes that encoded cuticular proteins, chitinolytic enzymes, chitin synthesis enzymes, and cuticle tanning proteins changed their expression levels remarkably, indicating these genes might play important roles in the restruction (degradation and biosynthesis) of insect cuticles during the larval metamorphosis. Additionally, the significant changes in the mRNA levels of 20-hydroxyecdysone (20E) and juvenile hormone (JH) related genes suggested their crucial roles in regulating cuticle remodeling during the larval metamorphosis of P. xylostella. In conclusion, the present study provide us the comprehensive gene expression profiles to explore the molecular mechanisms of cuticle metamorphosis in P. xylostella, which laid a molecular basis to study roles of specific pathways and genes in insect development.

Molecular mechanism of lateral bud differentiation of Pinus massoniana based on high-throughput sequencing

Knot-free timber cultivation is an important goal of forest breeding, and lateral shoots affect yield and stem shape of tree. The purpose of this study was to analyze the molecular mechanism of lateral bud development by removing the apical dominance of Pinus massoniana young seedlings through transcriptome sequencing and identify key genes involved in lateral bud development. We analyzed hormone contents and transcriptome data for removal of apical dominant of lateral buds as well as apical and lateral buds of normal development ones. Data were analyzed using an comprehensive approach of pathway- and gene-set enrichment analysis, Mapman visualization tool, and gene expression analysis. Our results showed that the contents of auxin (IAA), Zea and strigolactone (SL) in lateral buds significantly increased after removal of apical dominance, while abscisic acid (ABA) decreased. Gibberellin (GA) metabolism, cytokinin (CK), jasmonic acid, zeatin pathway-related genes positively regulated lateral bud development, ABA metabolism-related genes basically negatively regulated lateral bud differentiation, auxin, ethylene, SLs were positive and negative regulation, while only A small number of genes of SA and BRASSINOSTEROID, such as TGA and TCH4, were involved in lateral bud development. In addition, it was speculated that transcription factors such as WRKY, TCP, MYB, HSP, AuxIAA, and AP2 played important roles in the development of lateral buds. In summary, our results provided a better understanding of lateral bud differentiation and lateral shoot formation of P. massoniana from transcriptome level. It provided a basis for molecular characteristics of side branch formation of other timber forests, and contributed to knot-free breeding of forest trees.

Architectural Mediator subunits are differentially essential for global transcription in Saccharomyces cerevisiae

Mediator is a modular coactivator complex involved in the transcription of the majority of RNA polymerase II-regulated genes. However, the degrees to which individual core subunits of Mediator contribute to its activity have been unclear. Here, we investigate the contribution of two essential architectural subunits of Mediator to transcription in Saccharomyces cerevisiae. We show that acute depletion of the main complex scaffold Med14 or the head module nucleator Med17 is lethal and results in global transcriptional downregulation, though Med17 removal has a markedly greater negative effect. Consistent with this, Med17 depletion impairs preinitiation complex (PIC) assembly to a greater extent than Med14 removal. Co-depletion of Med14 and Med17 reduced transcription and TFIIB promoter occupancy similarly to Med17 ablation alone, indicating that the contributions of Med14 and Med17 to Mediator function are not additive. We propose that, while the structural integrity of complete Mediator and the head module are both important for PIC assembly and transcription, the head module plays a greater role in this process and is thus the key functional module of Mediator in this regard.

Hormonal Regulation of Reproductive Diapause That Occurs in the Year-Round Mass Rearing of Bombus terrestris Queens

Adult reproductive diapause is an adaptive strategy under adverse environments for insects and other arthropod species, including bumblebees, which enables queens to survive through a harsh winter and then build new colonies in the following spring. Little research has been done on the molecular regulatory mechanism of reproductive diapause in Bombus terrestris, which is an important pollinator of wild plants and crops. Our previous research identified the conditions that induced reproductive diapause during the year-round mass rearing of B. terrestris. Here, we performed combined transcriptomics and proteomics analyses of reproductive diapause in B. terrestris during and after diapause at three different ecophysiological phases, diapause, postdiapause, and founder postdiapause. The analyses showed that differentially expressed proteins/genes acted in the citrate cycle, insect hormone biosynthesis, insulin and mTOR signaling pathway. To further understand the mechanisms that regulated the reproductive diapause, genes involved in the regulation of JH synthesis, insulin/TOR signal pathway were determined. The BtRheb, BtTOR, BtVg, and BtJHAMT had lower expression levels in diapause queens. The JH III titer levels and the activities of the metabolic enzymes were significantly up-regulated in postdiapause queens. Also, after the microinjection of insulin-like peptides (ILPs) and JH analogue (JHA), hormones, cold-tolerance metabolites, metabolic enzymes, and reproduction showed significant changes. Together with results from other related research, a model of the regulation of reproductive diapause during the year-round mass rearing of B. terrestris was proposed. This study contributes to a comprehensive insight into the molecular regulatory mechanism of reproductive diapause in eusocial insects.

Identification and verification of genes related to pollen development and male sterility induced by high temperature in the thermo-sensitive genic male sterile wheat line

Bioinformatic analysis identified the function of genes regulating wheat fertility. Barley stripe mosaic virus-induced gene silencing verified that the genes TaMut11 and TaSF3 are involved in pollen development and related to fertility conversion. Environment-sensitive genic male sterility is of vital importance to hybrid vigor in crop production and breeding. Therefore, it is meaningful to study the function of the genes related to pollen development and male sterility, which is still not fully understand currently. In this study, YanZhan 4110S, a new thermo-sensitive genic male sterility wheat line, and its near-isogenic line YanZhan 4110 were analyzed. Through comparative transcriptome basic bioinformatics and weighted gene co-expression network to further identify some hub genes, the genes TaMut11 and TaSF3 associated with pollen development and male sterility induced by high-temperature were identified in YanZhan 4110S. Further verification through barley stripe mosaic virus-induced gene silencing elucidated that the silencing of TaMut11 and TaSF3 caused pollen abortion, finally resulting in the declination of fertility. These findings provided data on the abortive mechanism in environment-sensitive genic male sterility wheat.

Cytological Observation and Transcriptome Comparative Analysis of Self-Pollination and Cross-Pollination in Dendrobium Officinale

Dendrobium officinale is a rare and traditional medicinal plant with high pharmacological and nutritional value. The self-incompatibility mechanism of D. officinale reproductive isolation was formed in the long-term evolution process, but intraspecific hybridization of different germplasm resources leads to a large gap in the yield, quality, and medicinal value of D. officinale. To investigate the biological mechanism of self-incompatibility in D. officinale, cytological observation and the transcriptome analysis was carried out on the samples of self-pollination and cross-pollination in D. officinale. Results for self-pollination showed that the pollen tubes could grow in the style at 2 h, but most of pollen tubes stopped growing at 4 h, while a large number of cross-pollinated pollen tubes grew along the placental space to the base of ovary, indicating that the self-incompatibility of D. officinale may be gametophyte self-incompatibility. A total of 63.41 G basesum of D. officinale style samples from non-pollinated, self-pollination, and cross-pollination by RNA-seq were obtained, and a total of 1944, 1758, and 475 differentially expressed genes (DEGs) in the comparison of CK (non-pollinated) vs. HF (cross-pollination sample), CK vs. SF (self-pollination sample) and SF vs. HF were identified, respectively. Forty-one candidate genes related to self-incompatibility were found by function annotation of DEGs, including 6 Ca²⁺ signal genes, 4 armed repeat containing (ARC) related genes, 11 S-locus receptor kinase (SRK) related genes, 2 Exo70 family genes, 9 ubiquitin related genes, 1 fatty acid related gene, 6 amino acid-related genes, 1 pollen-specific leucine-rich repeat extensin-like protein (LRX) related gene and 1 lectin receptor-like kinases (RLKs) related gene, showed that self-incompatibility mechanism of D. officinale involves the interaction of multiple genes and pathways. The results can provide a basis for the study of the self-incompatibility mechanism of D. officinale, and provide ideas for the preservation and utilization of high-quality resources of D. officinale.

Polycystic ovary syndrome is transmitted via a transgenerational epigenetic process

Polycystic ovary syndrome (PCOS) is the most common reproductive and metabolic disorder affecting women of reproductive age. PCOS has a strong heritable component, but its pathogenesis has been unclear. Here, we performed RNA sequencing and genome-wide DNA methylation profiling of ovarian tissue from control and third-generation PCOS-like mice. We found that DNA hypomethylation regulates key genes associated with PCOS and that several of the differentially methylated genes are also altered in blood samples from women with PCOS compared with healthy controls. Based on this insight, we treated the PCOS mouse model with the methyl group donor S-adenosylmethionine and found that it corrected their transcriptomic, neuroendocrine, and metabolic defects. These findings show that the transmission of PCOS traits to future generations occurs via an altered landscape of DNA methylation and propose methylome markers as a possible diagnostic landmark for the condition, while also identifying potential candidates for epigenetic-based therapy.

Transcriptomic Responses of the Honey Bee Brain to Infection with Deformed Wing Virus

Managed colonies of European honey bees (Apis mellifera) are under threat from Varroa destructor mite infestation and infection with viruses vectored by mites. In particular, deformed wing virus (DWV) is a common viral pathogen infecting honey bees worldwide that has been shown to induce behavioral changes including precocious foraging and reduced associative learning. We investigated how DWV infection of bees affects the transcriptomic response of the brain. The transcriptomes of individual brains were analyzed using RNA-Seq after experimental infection of newly emerged adult bees with DWV. Two analytical methods were used to identify differentially expressed genes from the ~15,000 genes in the Apis mellifera genome. The 269 genes that had increased expression in DWV infected brains included genes involved in innate immunity such as antimicrobial peptides (AMPs), Ago2, and Dicer. Single bee brain NMR metabolomics methodology was developed for this work and indicates that proline is strongly elevated in DWV infected brains, consistent with the increased presence of the AMPs abaecin and apidaecin. The 1361 genes with reduced expression levels includes genes involved in cellular communication including G-protein coupled, tyrosine kinase, and ion-channel regulated signaling pathways. The number and function of the downregulated genes suggest that DWV has a major impact on neuron signaling that could explain DWV related behavioral changes.

Transcriptomic analysis reveals somatic embryogenesis-associated signaling pathways and gene expression regulation in maize (Zea mays L.)

Transcriptome analysis of maize embryogenic callus and somatic embryos reveals associated genes reprogramming, hormone signaling pathways and transcriptional regulation involved in somatic embryogenesis in maize. Somatic embryos are widely utilized in propagation and genetic engineering of crop plants. In our laboratory, an elite maize inbred line Y423 that could generate intact somatic embryos was obtained and applied to genetic transformation. To enhance our understanding of regulatory mechanisms during maize somatic embryogenesis, we used RNA-based sequencing (RNA-seq) to characterize the transcriptome of immature embryo (IE), embryogenic callus (EC) and somatic embryo (SE) from maize inbred line Y423. The number of differentially expressed genes (DEGs) in three pairwise comparisons (IE-vs-EC, IE-vs-SE and EC-vs-SE) was 5767, 7084 and 1065, respectively. The expression patterns of DEGs were separated into eight major clusters. Somatic embryogenesis associated genes were mainly grouped into cluster A or B with an expression trend toward up-regulation during dedifferentiation. GO annotation and KEGG pathway analysis revealed that DEGs were implicated in plant hormone signal transduction, stress response and metabolic process. Among the differentially expressed transcription factors, the most frequently represented families were associated with the common stress response or related to cell differentiation, embryogenic patterning and embryonic maturation processes. Genes include hormone response/transduction and stress response, as well as several transcription factors were discussed in this study, which may be potential candidates for further analyses regarding their roles in somatic embryogenesis. Furthermore, the temporal expression patterns of candidate genes were analyzed to reveal their roles in somatic embryogenesis. This transcriptomic data provide insights into future functional studies, which will facilitate further dissections of the molecular mechanisms that control maize somatic embryogenesis.

Gene expression changes in response to low temperatures in embryos of the kelp grouper, Epinephelus moara

The kelp grouper Epinephelus moara has high economic value and is popular in fisheries and aquaculture in China. In the previous study, we treated the embryos at 16-22 somite stage at 4 °C, -25.7 °C, -140 °C and -196 °C, and successfully obtained surviving embryos in each group. To better understand the molecular changes affected by the low temperatures, we conducted a comparative transcriptome analysis among embryos exposed at 4 °C for 30 min, embryos exposed at -25.7 °C for 30 min and the control group. qPCR assays were conducted for the validation. Signal transduction pathways were highly enriched for the differentially expressed genes. c-Fos, c-Jun, JunD, GADD45, involved in MAPK signaling pathway, were upregulated when embryos were treated at low temperatures. As immediate early genes, Egr-1a and b, and IER2, that respond quickly to the environment stress, their expression increased as well. Hsp70 showed similar expression pattern as immediate early genes. Meanwhile, transcription factors Sox, HES, TFIID, muscle movement and protein synthesis-related genes were downregulated. Taken together, our findings suggest that cooling disrupts gene expression patterns in E. moara embryos. The differentially expressed genes may be involved in cellular resistance against low temperatures, possibly through neural activation, apoptosis, proliferation, differentiation, cellular recovery and heat shock regulation. This study also provides transcriptome dataset of E. moara embryos exposed to cold temperatures for future studies focusing on the molecular effects of cryopreservation.

Cytological and molecular characterizations of a novel 2A nullisomic line derived from a widely-grown wheat cultivar Zhoumai 18 conferring male sterility

A dwarf, multi-pistil and male sterile dms mutant was previously reported by us. However, the genetic changes in this dms are unclear. To examine the genetic changes, single nucleotide polymorphism (SNP) association, chromosome counting, and high-resolution chromosome fluorescence in situ hybridization (FISH) techniques were employed. By comparing tall plants (T) with dwarf plants (D) in the offspring of dms mutant plants, SNP association analysis indicated that most SNPs were on chromosome 2A. There were three types in offspring of dms plants, with 42, 41 and 40 chromosomes respectively. High-resolution chromosome painting analysis demonstrated that T plants had all 42 wheat chromosomes; the medium plants (M) had 41 chromosomes, lacking one chromosome 2A; while D plants had 40 wheat chromosomes, and lacked both 2A chromosomes. These data demonstrated that dms resulted from a loss of chromosome 2A. We identified 23 genes on chromosome 2A which might be involved in the development of stamens or pollen grains. These results lay a solid foundation for further analysis of the molecular mechanisms of wheat male sterility. Because D plants can be used as a female parent to cross with other wheat genotypes, dms is a unique germplasm for any functional study of chromosome 2A and wheat breeding specifically targeting genes on 2A.

Metabolic and transcriptome responses of RNAi-mediated AMPKα knockdown in Tribolium castaneum

Background

The AMP-activated protein kinase (AMPK) is an intracellular fuel sensor for lipid and glucose metabolism. In addition to the short-term regulation of metabolic enzymes by phosphorylation, AMPK may also exert long-term effects on the transcription of downstream genes through the regulation of transcription factors and coactivators. In this study, RNA interference (RNAi) was conducted to investigate the effects of knockdown of TcAMPKα on lipid and carbohydrate metabolism in the red flour beetle, Tribolium castaneum, and the transcriptome profiles of dsTcAMPKα-injected and dsEGFP-injected beetles under normal conditions were compared by RNA-sequencing.

Results

RNAi-mediated suppression of TcAMPKα increased whole-body triglyceride (TG) level and the ratio between glucose and trehalose, as was confirmed by in vivo treatment with the AMPK-activating compound, 5-Aminoimidazole-4-carboxamide1-β-D-ribofuranoside (AICAR). A total of 1184 differentially expressed genes (DEGs) were identified between dsTcAMPKα-injected and dsEGFP-injected beetles. These include genes involved in lipid and carbohydrate metabolism as well as insulin/insulin-like growth factor signaling (IIS). Real-time quantitative polymerase chain reaction analysis confirmed the differential expression of selected genes. Interestingly, metabolism-related transcription factors such as sterol regulatory element-binding protein 1 (SREBP1) and carbohydrate response element-binding protein (ChREBP) were also significantly upregulated in dsTcAMPKα-injected beetles.

Conclusions

AMPK plays a critical role in the regulation of beetle metabolism. The findings of DEGs involved in lipid and carbohydrate metabolism provide valuable insight into the role of AMPK signaling in the transcriptional regulation of insect metabolism.

Insights Into the Albinism Mechanism for Two Distinct Color Morphs of Northern Snakehead, Channa argus Through Histological and Transcriptome Analyses

The great northern snakehead (Channa argus) is one of the most important economic and conservational fish in China. In this study, the melanocytes in the skin of two distinct color morphs C. argus were investigated and compared through employment of the microscopic analysis, hematoxylin and eosin (H&E) and Masson Fontana staining. Our results demonstrated the uneven distribution of melanocytes with extremely low density and most of them were in the state of aging or death. Meanwhile, there was no obvious pigment layer and melanocytes distribution pattern found in the albino-type (AT), while the melanocytes were evenly distributed with abundance in the bicolor-type (BT). The transcriptome analysis through Illumina HiSeq sequencing showed that a total of 34.93 Gb Clean Data was obtained, and Q30 base percentage reached 92.66%. The BT and AT northern snakeheads transcriptome data included a total of 56,039,701 and 60,410,063 clean reads (n = 3), respectively. In gene expression analyses, the sample correlation coefficients (r) were ranged between 0.92 and 1.00; the contribution of PC1 and PC2 were 50.25 and 13.73% by using PCA cluster analysis, the total number of DEGs were 1024 (559 up-regulated and 465 down-regulated), and the number of annotated DEGs was 767 (COG 172, KEGG 262, GO 288, SwissProt 548, Pfam 579 and NR 765). Additionally, 46,363 ± 873 and 44,947 ± 392 single nucleotide polymorphisms (SNPs) were compiled via genetic structure analysis, respectively. Ten key pigment-related genes were screened using qRT-PCR. And all of them revealed extremely higher expression levels in the skin of BT than those of AT. This is the first study to analyze the mechanism of albino characteristics of Channa via histology and transcriptomics, and also provide the oretical and practical support for the protection and development of germplasm resources for C. argus.

Comparative transcriptome and histomorphology analysis of testis tissues from mulard and Pekin ducks

Testicular transcriptomes were analyzed to characterize the differentially expressed genes between mulard and Pekin ducks, which will help establish gene expression datasets to assist in further determination of the mechanisms of genetic sterility in mulard ducks. Paraffin sections were made to compare the developmental differences in testis tissue between mulard and Pekin ducks. Comparative transcriptome sequencing of testis tissues was performed, and the expression of candidate genes was verified by quantitative reverse transcription-polymerase chain reaction (qRT-PCR). In mulard ducks, spermatogonia and spermatocytes were arranged in a disordered manner, and no mature sperm were observed in the testis tissue. However, different stages of development of sperm were observed in seminiferous tubules in the testis tissue of Pekin ducks. A total of 43.84 Gb of clean reads were assembled into 193 535 UniGenes. Of these, 2131 transcripts exhibited differential expression (false discover rate <0.001 and fold change ≥2), including 997 upregulated and 1134 downregulated transcripts in mulard ducks as compared to those in Pekin duck testis tissues. Several upregulated genes were related to reproductive functions, including ryanodine receptor 2 (RYR2), calmodulin (CALM), argininosuccinate synthase and delta-1-pyrroline-5-carboxylate synthetase ALDH18A1 (P5CS). Downregulated transcripts included the testis-specific serine/threonine-protein kinase 3, aquaporin-7 (AQP7) and glycerol kinase GlpK (GK). The 10 related transcripts involved in the developmental biological process were identified by GO (Gene Ontology) annotation. The KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways indicated that peroxisome proliferator-activated receptors (PPARs) and calcium signaling pathways were significantly (P<0.001) associated with normal testis physiology. The differential expression of select genes implicated in reproductive processes was verified by qRT-PCR, which was consistent with the expression trend of transcriptome sequencing (RNA-seq). Differentially expressed candidate genes RYR2, CALM, P5CS, AQP7 and GK were identified by transcriptional analysis in mulard and Pekin duck testes. These were important for the normal development of the male duck reproductive system. These data provide a framework for the further exploration of the molecular and genetic mechanisms of sterility in mulard ducks.

Highlights. The mulard duck is an intergeneric sterile hybrid offspring resulting from mating between Muscovy and Pekin ducks. The transcriptomes of testis tissue from mulard and Pekin ducks were systematically characterized, and differentially expressed genes were screened, in order to gain insights into potential gonad gene expression mechanisms contributing to genetic sterility in mulard ducks.

Transcriptional Profiling Suggests T Cells Cluster around Neurons Injected with Toxoplasma gondii Proteins

Like other persistent intracellular pathogens, Toxoplasma gondii, a protozoan parasite, has evolved to evade the immune system and establish a chronic infection in specific cells and organs, including neurons in the CNS. Understanding T. gondii’s persistence in neurons holds the potential to identify novel, curative drug targets. The work presented here offers new insights into the neuron-T. gondii interaction in vivo. By transcriptionally profiling neurons manipulated by T. gondii, we unexpectedly revealed that immune cells, and specifically CD8⁺ T cells, appear to cluster around these neurons, suggesting that CD8⁺ T cells specifically recognize parasite-manipulated neurons. Such a possibility supports evidence from other labs that questions the long-standing dogma that neurons are often persistently infected because they are not directly recognized by immune cells such as CD8⁺ T cells. Collectively, these data suggest we reconsider the broader role of neurons in the context of infection and neuroinflammation.

Toxoplasma gondii’s tropism for and persistence in the central nervous system (CNS) underlies the symptomatic disease that T. gondii causes in humans. Our recent work has shown that neurons are the primary CNS cell with which Toxoplasma interacts and which it infects in vivo. This predilection for neurons suggests that T. gondii’s persistence in the CNS depends specifically upon parasite manipulation of the host neurons. Yet, most work on T. gondii-host cell interactions has been done in vitro and in nonneuronal cells. We address this gap by utilizing our T. gondii-Cre system that allows permanent marking and tracking of neurons injected with parasite effector proteins in vivo. Using laser capture microdissection (LCM) and RNA sequencing using RNA-seq, we isolated and transcriptionally profiled T. gondii-injected neurons (TINs), Bystander neurons (nearby non-T. gondii-injected neurons), and neurons from uninfected mice (controls). These profiles show that TIN transcriptomes significantly differ from the transcriptomes of Bystander and control neurons and that much of this difference is driven by increased levels of transcripts from immune cells, especially CD8⁺ T cells and monocytes. These data suggest that when we used LCM to isolate neurons from infected mice, we also picked up fragments of CD8⁺ T cells and monocytes clustering in extreme proximity around TINs and, to a lesser extent, Bystander neurons. In addition, we found that T. gondii transcripts were primarily found in the TIN transcriptome, not in the Bystander transcriptome. Collectively, these data suggest that, contrary to common perception, neurons that directly interact with or harbor parasites can be recognized by CD8⁺ T cells.

IMPORTANCE Like other persistent intracellular pathogens, Toxoplasma gondii, a protozoan parasite, has evolved to evade the immune system and establish a chronic infection in specific cells and organs, including neurons in the CNS. Understanding T. gondii’s persistence in neurons holds the potential to identify novel, curative drug targets. The work presented here offers new insights into the neuron-T. gondii interaction in vivo. By transcriptionally profiling neurons manipulated by T. gondii, we unexpectedly revealed that immune cells, and specifically CD8⁺ T cells, appear to cluster around these neurons, suggesting that CD8⁺ T cells specifically recognize parasite-manipulated neurons. Such a possibility supports evidence from other labs that questions the long-standing dogma that neurons are often persistently infected because they are not directly recognized by immune cells such as CD8⁺ T cells. Collectively, these data suggest we reconsider the broader role of neurons in the context of infection and neuroinflammation.