Integration of Transcriptomics and Metabolomics Reveals the Responses of Sugar Beet to Continuous Cropping Obstacle

Sugar beet is vulnerable to years of continuous cropping, and allelopathy is one of the important factors leading to continuous cropping disorder. To explore the physiological and molecular mechanisms behind continuous cropping obstacles on sugar beet, this study combined transcriptomics and metabolomics to analyze the effects of different years of continuous cropping on metabolite changes, differential gene expression, and root exudate regulation in sugar beet. We collected sugar beet's root samples from 1-, 3-, and 5-year continuous cropping systems for metabolome and transcriptome analyses. Our data revealed that T3 and T5 had 50 and 33 metabolites significantly different from T1, respectively. The autotoxic substance salicylaldehyde was found to continuously accumulate in root exudates with increasing years of continuous cropping. Sucrose was highly reduced in T3 (4.05-fold decrease) and T5 (2.01-fold decrease) compared to T1. Respectively, 2,660 and 3,515 differentially expressed genes (DEGs) were significantly regulated in T3 and T5 compared to T1. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed that metabolic pathways and biosynthesis of secondary metabolites were perturbed in T3 and T5 vs. T1. Integrated metabolomics analyses identified 73 DEGs involved in enriched metabolic pathways, all of which were the oxidation-reduction process pathways. In conclusion, this study provides evidence that continuous cropping obstacles can change the metabolome and transcriptome of sugar beet, affecting its growth and quality.

A Computational Framework to Identify Biomarkers for Glioma Recurrence and Potential Drugs Targeting Them

Background: Recurrence is still a major obstacle to the successful treatment of gliomas. Understanding the underlying mechanisms of recurrence may help for developing new drugs to combat gliomas recurrence. This study provides a strategy to discover new drugs for recurrent gliomas based on drug perturbation induced gene expression changes. Methods: The RNA-seq data of 511 low grade gliomas primary tumor samples (LGG-P), 18 low grade gliomas recurrent tumor samples (LGG-R), 155 glioblastoma multiforme primary tumor samples (GBM-P), and 13 glioblastoma multiforme recurrent tumor samples (GBM-R) were downloaded from TCGA database. DESeq2, key driver analysis and weighted gene correlation network analysis (WGCNA) were conducted to identify differentially expressed genes (DEGs), key driver genes and coexpression networks between LGG-P vs LGG-R, GBM-P vs GBM-R pairs. Then, the CREEDS database was used to find potential drugs that could reverse the DEGs and key drivers. Results: We identified 75 upregulated and 130 downregulated genes between LGG-P and LGG-R samples, which were mainly enriched in human papillomavirus (HPV) infection, PI3K-Akt signaling pathway, Wnt signaling pathway, and ECM-receptor interaction. A total of 262 key driver genes were obtained with frizzled class receptor 8 (FZD8), guanine nucleotide-binding protein subunit gamma-12 (GNG12), and G protein subunit β2 (GNB2) as the top hub genes. By screening the CREEDS database, we got 4 drugs (Paclitaxel, 6-benzyladenine, Erlotinib, Cidofovir) that could downregulate the expression of up-regulated genes and 5 drugs (Fenofibrate, Oxaliplatin, Bilirubin, Nutlins, Valproic acid) that could upregulate the expression of down-regulated genes. These drugs may have a potential in combating recurrence of gliomas. Conclusion: We proposed a time-saving strategy based on drug perturbation induced gene expression changes to find new drugs that may have a potential to treat recurrent gliomas.

In Silico Identification of Stress-Associated Transposable Elements in Arabidopsis thaliana Using Public Transcriptome Data

Transposable elements (TEs) have been associated with stress response in many plants, making them a key target of study. However, the high variability, genomic repeat-heavy nature, and widely noncoding character of TEs have made them difficult to study using non-specialized methods, whether experimental or computational. In this chapter, we introduce two computational workflows to analyze transposable elements using publicly available transcriptome data. In the first of these methods, we identify TEs, which show differential expression under salt stress using sample transcriptome libraries that includes noncoding transcripts. In the second, we identify protein-coding genes with differential expression under the same conditions, and determine which TEs are enriched in the promoter regions of these stress-related genes.

A Revised Adaptation of the Smart-Seq2 Protocol for Single-Nematode RNA-Seq

The advancement of transcriptomic studies in plant parasitic nematodes will greatly benefit from the development of single-nematode RNA-seq methods. Since many plant parasitic nematodes are obligate parasites, it is often difficult to efficiently obtain sufficient amounts of nematodes for transcriptomic studies. Here we have adapted SMART-Seq2 for single-nematode RNA-seq requiring only an individual nematode for a sample replicate. This protocol provides a detailed step-by-step procedure of the RNA-seq workflow starting from lysis of the nematode to quantification of transcripts using a user-friendly online platform.

Functional Analysis of Wheat NAC Transcription Factor, TaNAC069, in Regulating Resistance of Wheat to Leaf Rust Fungus

NAC transcription factors are one of the largest transcription factor families having functions in a variety of stress responses. Few NACs have been reported for interactions between wheat and the wheat rust fungus Puccinia triticina (Pt). In this study, based on analysis of RNA-seq data from wheat line TcLr19 inoculated by Pt, the NAC transcription factor TaNAC069 was cloned from wheat, and its transcriptional activity and homologous dimer formation were verified. Quantitative real-time PCR analysis showed that the expression of TaNAC069 was induced by Pt and associated signaling molecules. To further characterize the function of the TaNAC069 gene in wheat resistance to Pt, virus-induced gene silencing (VIGS) was utilized, and it revealed that Pt resistance in TaNAC069-silenced plants was significantly reduced. Potential interaction targets of TaNAC069 from wheat and Pt were screened and identified by yeast two-hybrid technology. Eukaryotic elongation factor eEF1A, CBSX3 protein, and cold acclimation protein WCOR410c were screened by yeast one-hybrid technology. The results indicate that the TaNAC069 gene plays a positive regulatory role in wheat resistance to Pt, laying a good foundation to analyze the molecular mechanisms of TaNAC069 and its functional role in wheat resistance to Pt.

Become Competent in Generating RNA-Seq Heat Maps in One Day for Novices Without Prior R Experience

Heat map visualization of RNA-seq data is a commonplace task. However, most laboratories rely on bioinformaticians who are not always available. Biological scientists are afraid to prepare heat maps independently because R is a programming platform. Here, using RNA-seq data for 16 differentially expressed genes in WNT pathway between embryonic stem cells and fibroblasts, I share a tutorial for novices without any prior R experience to master the skills, in one day, required for preparation of heat maps using the pheatmap package. Procedures described include installation of R, RStudio, and the pheatmap package, as well as hands-on practices for some basic R commands, conversion of RNA-seq data frame to a numeric matrix suitable for generation of heat maps, and defining arguments for the pheatmap function to make a desired heat map. More than 20 template scripts are provided to generate heat maps and to control the dimensions and appearances of the heat maps.

Genes That Mediate Starch Metabolism in Developing and Germinated Barley Grain

Starch is synthesized in the endosperm of developing barley grain, where it functions as the primary source of stored carbohydrate. In germinated grain these starch reserves are hydrolyzed to small oligosaccharides and glucose, which are transported to the embryo to support the growth of the developing seedling. Some of the mobilized glucose is transiently stored as starch in the scutellum of germinated grain. These processes are crucial for early seedling vigor, which is a key determinant of crop productivity and global food security. Several starch synthases (SS), starch-branching enzymes (SBEs), and starch debranching enzymes (isoamylases, ISA), together with a limit dextrinase (LD), have been implicated in starch synthesis from nucleotide-sugar precursors. Starch synthesis occurs both in the developing endosperm and in the scutellum of germinated grain. For the complete depolymerization of starch to glucose, α-amylase (Amy), β-amylase (Bmy), isoamylase (ISA), limit dextrinase (LD), and α-glucosidase (AGL) are required. Most of these enzymes are encoded by gene families of up to 10 or more members. Here RNA-seq transcription data from isolated tissues of intact developing and germinated barley grain have allowed us to identify the most important, specific gene family members for each of these processes in vivo and, at the same time, we have defined in detail the spatio-temporal coordination of gene expression in different tissues of the grain. A transcript dataset for 81,280 genes is publicly available as a resource for investigations into other cellular and biochemical processes that occur in the developing grain from 6 days after pollination.

Analysis of Thyroid Hormone Receptor α-Knockout Tadpoles Reveals That the Activation of Cell Cycle Program Is Involved in Thyroid Hormone-Induced Larval Epithelial Cell Death and Adult Intestinal Stem Cell Development During Xenopus tropicalis Metamorphosis

Background: There are two highly conserved thyroid hormone (triiodothyronine [T3]) receptor (TR) genes, TRα and TRβ, in all vertebrates, and the expression of TRα but not TRβ is activated earlier than T3 synthesis during development. In human, high levels of T3 are present during the several months around birth, and T3 deficiency during this period causes severe developmental abnormalities including skeletal and intestinal defects. It is, however, difficult to study this period in mammals as the embryos and neonates depend on maternal supply of nutrients for survival. However, Xenopus tropicalis undergoes a T3-dependent metamorphosis, which drastically changes essentially every organ in a tadpole. Of interest is intestinal remodeling, which involves near complete degeneration of the larval epithelium through apoptosis. Concurrently, adult intestinal stem cells are formed de novo and subsequently give rise to the self-renewing adult epithelial system, resembling intestinal maturation around birth in mammals. We have previously demonstrated that T3 signaling is essential for the formation of adult intestinal stem cells during metamorphosis. Methods: We studied the function of endogenous TRα in the tadpole intestine by using knockout animals and RNA-seq analysis. Results: We observed that removing endogenous TRα caused defects in intestinal remodeling, including drastically reduced larval epithelial cell death and adult intestinal stem cell proliferation. Using RNA-seq on intestinal RNA from premetamorphic wild-type and TRα-knockout tadpoles treated with or without T3 for one day, before any detectable T3-induced cell death and stem cell formation in the tadpole intestine, we identified more than 1500 genes, which were regulated by T3 treatment of the wild-type but not TRα-knockout tadpoles. Gene Ontology and biological pathway analyses revealed that surprisingly, these TRα-regulated genes were highly enriched with cell cycle-related genes, in addition to genes related to stem cells and apoptosis. Conclusions: Our findings suggest that TRα-mediated T3 activation of the cell cycle program is involved in larval epithelial cell death and adult epithelial stem cell development during intestinal remodeling.

Detection of GPCR mRNA Expression in Primary Cells Via qPCR, Microarrays, and RNA-Sequencing

A workflow is described for assaying the expression of G protein-coupled receptors (GPCRs) in cultured cells, using a combination of methods that assess GPCR mRNAs. Beginning from the isolation of cDNA and preparation of mRNA, we provide protocols for designing and testing qPCR primers, assaying mRNA expression using qPCR and high-throughput analysis of GPCR mRNA expression via TaqMan qPCR-based, GPCR-selective arrays. We also provide a workflow for analysis of expression from RNA-sequencing (RNA-seq) assays, which can be queried to yield expression of GPCRs and related genes in samples of interest, as well as to test changes in expression between groups, such as in cells treated with drugs or from healthy and diseased subjects. We place priority on optimized protocols that distinguish signal from noise, as GPCR mRNAs are typically present in low abundance, necessitating techniques that maximize sensitivity while minimizing noise. These methods may also be applicable for assessing the expression of members of families of other low abundance genes via high-throughput analyses of mRNAs, followed by independent confirmation and validation of results via qPCR.

Exogenous Gene Expression and Insect Resistance in Dual Bt Toxin Populus × euramericana 'Neva' Transgenic Plants

Bacillus thuringiensis (Bt) insecticidal protein genes are important tools in efforts to develop insect resistance in poplar. In this study, the Cry1Ac and Cry3A Bt toxin genes were simultaneously transformed into the poplar variety Populus × euramericana 'Neva' by Agrobacterium-mediated transformation to explore the exogenous gene expression and insect resistance, and to examine the effects of Bt toxin on the growth and development of Anoplophora glabripennis larvae after feeding on the transgenic plant. Integration and expression of the transgenes were determined by molecular analyses and the insect resistance of transgenic lines was evaluated in feeding experiments. Sixteen transgenic dual Bt toxin genes Populus × euramericana 'Neva' lines were obtained. The dual Bt toxin genes were expressed at both the transcriptional and translational levels; however, Cry3A protein levels were much higher than those of Cry1Ac. Some of the transgenic lines exhibited high resistance to the first instar larvae of Hyphantria cunea and Micromelalopha troglodyta, and the first and second instar larvae and adults of Plagiodera versicolora. Six transgenic lines inhibited the growth and development of A. glabripennis larvae. The differences in the transcriptomes of A. glabripennis larvae fed transgenic lines or non-transgenic control by RNA-seq analyses were determined to reveal the mechanism by which Bt toxin regulates the growth and development of longicorn beetle larvae. The expression of genes related to Bt prototoxin activation, digestive enzymes, binding receptors, and detoxification and protective enzymes showed significant changes in A. glabripennis larvae fed Bt toxin, indicating that the larvae responded by regulating the expression of genes related to their growth and development. This study lay a theoretical foundation for developing resistance to A. glabripennis in poplar, and provide a foundation for exploring the mechanism of Bt toxin action on Cerambycidae insects.

Integrated Transcriptomics and Metabolomics Analyses Provide Insights Into the Response of Chongyi Wild Mandarin to Candidatus Liberibacter Asiaticus Infection

Candidatus Liberibacter asiaticus (CLas) is the causative agent of Huanglongbing (HLB), which has caused great economic losses to the citrus industry. The molecular mechanism of the host response to CLas in wild citrus germplasm has been reported less. Eighteen weeks after inoculation via grafting, all the CLas-inoculated Chongyi wild mandarin (Citrus reticulata) were positive and showed severe anatomical aberrations, suggesting its susceptibility to HLB. Transcriptomics and metabolomics analyses of leaves, barks, and roots from mock-inoculated (control) and CLas-inoculated seedlings were performed. Comparative transcriptomics identified 3,628, 3,770, and 1,716 differentially expressed genes (DEGs) between CLas-infected and healthy tissues in the leaves, barks, and roots, respectively. The CLas-infected tissues had higher transcripts per kilobase per million values and more genes that reached their maximal expression, suggesting that HLB might cause an overall increase in transcript accumulation. However, HLB-triggered transcriptional alteration showed tissue specificity. In the CLas-infected leaves, many DEGs encoding immune receptors were downregulated. In the CLas-infected barks, nearly all the DEGs involved in signaling and plant-pathogen interaction were upregulated. In the CLas-infected roots, DEGs encoding enzymes or transporters involved in carotenoid biosynthesis and nitrogen metabolism were downregulated. Metabolomics identified 71, 62, and 50 differentially accumulated metabolites (DAMs) in the CLas-infected leaves, barks and roots, respectively. By associating DEGs with DAMs, nitrogen metabolism was the only pathway shared by the three infected tissues and was depressed in the CLas-infected roots. In addition, 26 genes were determined as putative markers of CLas infection, and a hypothesized model for the HLB susceptibility mechanism in Chongyi was proposed. Our study may shed light on investigating the molecular mechanism of the host response to CLas infection in wild citrus germplasm.

Genomic Profiling of Virulence in the Soybean Cyst Nematode Using Single-Nematode Sequencing

Soybean cyst nematode (SCN) is one of the most important diseases in soybean. Currently, the main management strategy relies on planting resistant cultivars. However, the overuse of a single resistance source has led to the selection of virulent SCN populations, although the mechanisms by which the nematode overcomes the resistance genes remain unknown. In this study, we used a nematode-adapted single-cell RNA-seq approach to identify SCN genes potentially involved in resistance breakdown in Peking and PI 88788 parental soybean lines. We established for the first time the full transcriptome of single SCN individuals allowing us to identify a list of putative virulence genes against both major SCN resistance sources. Our analysis identified 48 differentially expressed putative effectors (secreted proteins required for infection) alongside 40 effectors showing evidence of novel structural variants, and 11 effector genes containing phenotype-specific sequence polymorphisms. Additionally, a differential expression analysis revealed an interesting phenomenon of coexpressed gene regions with some containing putative effectors. The selection of virulent SCN individuals on Peking resulted in a profoundly altered transcriptome, especially for genes known to be involved in parasitism. Several sequence polymorphisms were also specific to these virulent nematodes and could potentially play a role in the acquisition of nematode virulence. On the other hand, the transcriptome of virulent individuals on PI 88788 was very similar to avirulent ones with the exception of a few genes, which suggest a distinct virulence strategy to Peking.

Phosphate availability and ectomycorrhizal symbiosis with Pinus sylvestris have independent effects on the Paxillus involutus transcriptome

Many plant species form symbioses with ectomycorrhizal fungi, which help them forage for limiting nutrients in the soil such as inorganic phosphate (Pi). The transcriptional responses to symbiosis and nutrient-limiting conditions in ectomycorrhizal fungal hyphae, however, are largely unknown. An artificial system was developed to study ectomycorrhizal basidiomycete Paxillus involutus growth in symbiosis with its host tree Pinus sylvestris at different Pi concentrations. RNA-seq analysis was performed on P. involutus hyphae growing under Pi-limiting conditions, either in symbiosis or alone. We show that Pi starvation and ectomycorrhizal symbiosis have an independent effect on the P. involutus transcriptome. Notably, low Pi availability induces expression of newly identified putative high-affinity Pi transporter genes, while reducing the expression of putative organic acid transporters. Additionally, low Pi availability induces a close transcriptional interplay between P and N metabolism. GTP-related signalling was found to have a positive effect in the maintenance of ectomycorrhizal symbiosis, whereas multiple putative cytochrome P450 genes were found to be downregulated, unlike arbuscular mycorrhizal fungi. We provide the first evidence of global transcriptional changes induced by low Pi availability and ectomycorrhizal symbiosis in the hyphae of P. involutus, revealing both similarities and differences with better-characterized arbuscular mycorrhizal fungi.

Monocyte Subsets With High Osteoclastogenic Potential and Their Epigenetic Regulation Orchestrated by IRF8

Osteoclasts (OCs) are bone-resorbing cells formed by the serial fusion of monocytes. In mice and humans, three distinct subsets of monocytes exist; however, it is unclear if all of them exhibit osteoclastogenic potential. Here we show that in wild-type (WT) mice, Ly6C^hi and Ly6C^int monocytes are the primary source of OC formation when compared to Ly6C^- monocytes. Their osteoclastogenic potential is dictated by increased expression of signaling receptors and activation of preestablished transcripts, as well as de novo gain in enhancer activity and promoter changes. In the absence of interferon regulatory factor 8 (IRF8), a transcription factor important for myelopoiesis and osteoclastogenesis, all three monocyte subsets are programmed to display higher osteoclastogenic potential. Enhanced NFATc1 nuclear translocation and amplified transcriptomic and epigenetic changes initiated at early developmental stages direct the increased osteoclastogenesis in Irf8-deficient mice. Collectively, our study provides novel insights into the transcription factors and active cis-regulatory elements that regulate OC differentiation. © 2020 American Society for Bone and Mineral Research (ASBMR).

Phenotypic Characteristics and Transcriptome of Cucumber Male Flower Development Under Heat Stress

Cucumber (Cucumis sativus L.) is an important vegetable crop, which is thermophilic not heat resistant. High-temperature stress always results in sterility at reproductive stage. In the present study, we evaluate the male flower developmental changes under normal (CK) and heat stress (HS) condition. After HS, the activities of peroxidase (POD) and superoxide dismutase (SOD) and the contents of malondialdehyde (MDA) were increased. In addition, the pollen fertility was significantly decreased; and abnormal tapetum and microspore were observed by paraffin section. Transcriptome analysis results presented that total of 5828 differentially expressed genes (DEGs) were identified after HS. Among these DEGs, 20 DEGs were found at four stages, including DNA binding transcription factor, glycosyltransferase, and wound-responsive family protein. The gene ontology term of carbohydrate metabolic process was significantly enriched in all anther stages, and many saccharides and starch synthase-related genes, such as invertase, sucrose synthase, and starch branching enzyme, were significantly different expressed in HS compared with CK. Furthermore, co-expression network analysis showed a module (midnightblue) strongly consistent with HS, and two hub genes (CsaV3_6G004180 and CsaV3_5G034860) were found with a high degree of connectivity to other genes. Our results provide comprehensive understandings on male flower development in cucumber under HS.

Extracellular Fragmented Self-DNA Is Involved in Plant Responses to Biotic Stress

A growing body of evidence indicates that extracellular fragmented self-DNA (eDNA), by acting as a signaling molecule, triggers inhibitory effects on conspecific plants and functions as a damage-associated molecular pattern (DAMP). To evaluate early and late events in DAMP-dependent responses to eDNA, we extracted, fragmented, and applied the tomato (Solanum lycopersicum) eDNA to tomato leaves. Non-sonicated, intact self-DNA (intact DNA) was used as control. Early event analyses included the evaluation of plasma transmembrane potentials (Vm), cytosolic calcium variations (Ca²⁺ _{cy t}), the activity and subcellular localization of both voltage-gated and ligand-gated rectified K⁺ channels, and the reactive oxygen species (ROS) subcellular localization and quantification. Late events included RNA-Seq transcriptomic analysis and qPCR validation of gene expression of tomato leaves exposed to tomato eDNA. Application of eDNA induced a concentration-dependent Vm depolarization which was correlated to an increase in (Ca²⁺)_cyt; this event was associated to the opening of K⁺ channels, with particular action on ligand-gated rectified K⁺ channels. Both eDNA-dependent (Ca²⁺)_cyt and K⁺ increases were correlated to ROS production. In contrast, application of intact DNA produced no effects. The plant response to eDNA was the modulation of the expression of genes involved in plant-biotic interactions including pathogenesis-related proteins (PRPs), calcium-dependent protein kinases (CPK1), heat shock transcription factors (Hsf), heat shock proteins (Hsp), receptor-like kinases (RLKs), and ethylene-responsive factors (ERFs). Several genes involved in calcium signaling, ROS scavenging and ion homeostasis were also modulated by application of eDNA. Shared elements among the transcriptional response to eDNA and to biotic stress indicate that eDNA might be a common DAMP that triggers plant responses to pathogens and herbivores, particularly to those that intensive plant cell disruption or cell death. Our results suggest the intriguing hypothesis that some of the plant reactions to pathogens and herbivores might be due to DNA degradation, especially when associated to the plant cell disruption. Fragmented DNA would then become an important and powerful elicitor able to trigger early and late responses to biotic stress.

Genome-Wide Association Analysis Coupled With Transcriptome Analysis Reveals Candidate Genes Related to Salt Stress in Alfalfa (Medicago sativa L.)

Salt stress is the main abiotic factor affecting alfalfa yield and quality. However, knowledge of the genetic basis of the salt stress response in alfalfa is still limited. Here, a genome-wide association study (GWAS) involving 875,023 single-nucleotide polymorphisms (SNPs) was conducted on 220 alfalfa varieties under both normal and salt-stress conditions. Phenotypic analysis showed that breeding status and geographical origin play important roles in the alfalfa salt stress response. For germination ability under salt stress, a total of 15 significant SNPs explaining 9%-14% of the phenotypic variation were identified. For tolerance to salt stress in the seedling stage, a total of 18 significant SNPs explaining 12%-23% of the phenotypic variation were identified. Transcriptome analysis revealed 2,097 and 812 differentially expressed genes (DEGs) that were upregulated and 2,445 and 928 DEGs that were downregulated in the leaves and roots, respectively, under salt stress. Among these DEGs, many encoding transcription factors (TFs) were found, including MYB-, CBF-, NAC-, and bZIP-encoding genes. Combining the results of our GWAS analysis and transcriptome analysis, we identified a total of eight candidate genes (five candidate genes for tolerance to salt stress and three candidate genes for germination ability under salt stress). Two SNPs located within the upstream region of MsAUX28, which encodes an auxin response protein, were significantly associated with tolerance to salt stress. The two significant SNPs within the upstream region of MsAUX28 existed as three different haplotypes in this panel. Hap 1 (G/G, A/A) was under selection in the alfalfa domestication and improvement process.

Concordance of Immune-Related Markers in Lymphocytes and Prefrontal Cortex in Schizophrenia

Schizophrenia is a severe neuropsychiatric disorder associated with a wide array of transcriptomic and neurobiochemical changes. Genome-wide transcriptomic profiling conducted in postmortem brain have provided novel insights into the pathophysiology of this disorder, and identified biological processes including immune/inflammatory-related responses, metabolic, endocrine, and synaptic function. However, few studies have investigated whether similar changes are present in peripheral tissue. Here, we used RNA-sequencing to characterize transcriptomic profiles of lymphocytes in 18 nonpsychotic controls and 19 individuals with schizophrenia. We identified 2819 differentially expressed transcripts (P nominal < .05) in the schizophrenia group when compared to controls. Bioinformatic analyses conducted on a subset of 293 genes (P nominal < .01 and |log2 FC| > 0.5) highlighted immune/inflammatory responses as key biological processes in our dataset. Differentially expressed genes in lymphocytes were highly enriched in gene expression profiles associated with cortex layer 5a and immune cells. Thus, we investigated whether the changes in transcripts levels observed in lymphocytes could also be detected in the prefrontal cortex (PFC, BA10) in a second replication cohort of schizophrenia subjects. Remarkably, mRNA levels detected in the PFC and lymphocytes were in strong agreement, and measurements obtained using RNA-sequencing positively correlated with data obtained by reverse transcriptase-quantitative polymerase chain reaction analysis. Collectively, our work supports a role for immune dysfunction in the pathogenesis of schizophrenia and suggests that peripheral markers can be used as accessible surrogates to investigate putative central nervous system disruptions.

Creation of a Low-Alcohol-Production Yeast by a Mutated SPT15 Transcription Regulator Triggers Transcriptional and Metabolic Changes During Wine Fermentation

There is significant interest in the wine industry to develop methods to reduce the ethanol content of wine. Here the global transcription machinery engineering (gTME) technology was used to engineer a yeast strain with decreased ethanol yield, based on the mutation of the SPT15 gene. We created a strain of Saccharomyces cerevisiae (YS59-409), which possessed ethanol yield reduced by 34.9%; this was accompanied by the increase in CO₂, biomass, and glycerol formation. Five mutation sites were identified in the mutated SPT15 gene of YS59-409. RNA-Seq and metabolome analysis of YS59-409 were conducted compared with control strain, suggesting that ribosome biogenesis, nucleotide metabolism, glycolysis flux, Crabtree effect, NAD⁺/NADH homeostasis and energy metabolism might be regulated by the mutagenesis of SPT15 gene. Furthermore, two genes related to energy metabolism, RGI1 and RGI2, were found to be associated with the weakened ethanol production capacity, although the precise mechanisms involved need to be further elucidated. This study highlighted the importance of applying gTME technology when attempting to reduce ethanol production by yeast, possibly reprogramming yeast’s metabolism at the global level.

Overexpression of ThMYB8 mediates salt stress tolerance by directly activating stress-responsive gene expression

MYB transcription factors are important in abiotic stress responses; however, the detailed mechanisms are unclear. Tamarix hispida contains multiple MYB genes. The present study characterized T. hispida MYB8 (ThMYB8) during salt stress using transgenic T. hispida and Arabidopsis assays. ThMYB8 overexpression and ThMYB8 RNAi analysis demonstrated that ThMYB8 enhanced the salt stress tolerance. Transgenic Arabidopsis ectopic expression of ThMYB8 significantly increased root growth, fresh weight, and seed germination rate compared with that of the wild-type under salt stress. Physiological parameters analysis in T. hispida and Arabidopsis showed that ThMYB8 overexpressing plants had the lowest levels of O², H₂O₂, cell death, malondialdehyde, and electrolyte leakage. Overexpression of ThMYB8 regulated Na⁺ and K⁺ concentrations in plant tissues while maintaining K⁺/Na⁺ homeostasis. Analysis using qRT-PCR and ChIP-PCR identified possible downstream ThMYB8-regulated genes. ThMYB8 regulated the expression of ThCYP450-2 (cytochrome p450-2), Thltk (leucine-rich repeat transmembrane protein kinase), and ThTIP (aquaporin TIP) by binding to the MBSI motif ('CAACTG') in their promoters. The results indicated that ThMYB8 enhanced salt stress tolerance in T. hispida by regulating gene expression related to the activation of stress-associated physiological changes, such as enhanced reactive oxygen species scavenging capability, maintaining K⁺/Na⁺ homeostasis, and decreasing the malondialdehyde content and lipid peroxidation cell membranes.

Transcriptome analysis of sheep follicular development during prerecruitment, dominant, and mature stages after FSH superstimulation

Sheep is usually a monovular animal; superovulation technology is used to increase the number of offspring per individual and shorten generation intervals. To date, mature FSH superstimulatory treatments have been successfully used in sheep breeding, but much remains unknown about genes, pathways, and biological functions involved in follicular development. Therefore, in this study, we performed transcriptome profiling of small follicles (SFs; 2-2.5 mm), medium follicles (MFs; 3.5-4.5 mm), and large follicles (LFs; > 6 mm) in Mongolian ewes after FSH superstimulation. Furthermore, we identified differentially expressed genes and performed Kyoto Encyclopedia of Genes and Genomes pathway and Gene Ontology enrichment analyses in 3 separate pairwise comparisons. We found that ovarian steroidogenesis was significantly enriched in the SFs versus MFs analysis; the associated genes, cytochrome P450 family 19 (CYP19) and Hydroxy-delta-5-steroid dehydrogenase 3 beta- and steroid delta-isomerase 1 (HSD3B1), were significantly upregulated. Moreover, proline metabolism, glutathione metabolism, and PPAR signaling pathways were significantly enriched in the LFs versus SFs analysis; the associated genes, glutamate-cysteine ligase modifier subunit (GCLM) and cystathionine gamma-lyase (CTH), were significantly upregulated, whereas peroxisome proliferator-activated receptor gamma (PPARγ) was significantly downregulated. In summary, our study provides basic data and possible biological direction to further explore the molecular mechanism of sheep follicular development after FSH superstimulation.

Meiocyte Isolation by INTACT and Meiotic Transcriptome Analysis in Arabidopsis

Isolation of nuclei tagged in specific cell types (INTACT) is a method developed to isolate cell-type-specific nuclei that are tagged through in vivo biotin labeling of a nuclear targeting fusion (NTF) protein. In our work, INTACT was used to capture nuclei of meiocytes and to generate a meiotic transcriptome in Arabidopsis. Using the promoter of AtDMC1 recombinase to label meiotic nuclei, we generated transgenic plants carrying AtDMC1:NTF along with biotin ligase enzyme (BirA) under the constitutive ACTIN2 (ACT2) promoter. AtDMC1-driven expression of biotin-labeled NTF allowed us to collect nuclei of meiocytes by streptavidin-coated magnetic beads. The nuclear meiotic transcriptome was obtained by RNA-seq using low-quantity input RNA. Transcripts grouped into different categories according to their expression levels were investigated by gene ontology enrichment analysis (GOEA). The most enriched GO term "DNA demethylation" in mid/high-expression classes suggests that this biological process is particularly relevant to meiosis onset. The majority of genes with established roles in meiosis were distributed in the classes of mid/high and high expression. Meiotic transcriptome was compared with public available transcriptomes from other tissues in Arabidopsis. Bioinformatics analysis by expression network identified a core of more than 1,500 genes related to meiosis landmarks.

Effect of Vitamin E Supplementation on Deposition and Gene Expression Profiling of Abdominal Fat in Broiler Chickens

The aim of this study was to study the regulation of abdominal fat deposition by DL-α-tocopherol acetate (vitamin E) in broilers. Diets supplemented with 50 IU vitamin E significantly diminished abdominal fat deposition in broilers at day 35. Transcriptome sequencing results for abdominal fat tissues of the control (FC) and 50 IU vitamin E-supplemented (FT) groups identified 602 differentially expressed genes (DEGs), which were enriched in cellular process, cell and cell part, and binding Gene Ontology terms. Pathway functional analysis revealed that the DEGs were enriched in 42 metabolic pathways. Notably, the most enriched pathway, fatty acid biosynthesis, was found to play a key role in lipid metabolism. Further, the key regulators of lipid metabolism, including fatty acid synthase, acetyl-CoA carboxylase alpha, and acyl-CoA synthetase long-chain family member 1, demonstrated decreased expression following vitamin E supplementation. Herein, we have identified pathways and genes regulated by vitamin E, thereby providing novel insights into the nutrients regulating abdominal fat deposition in broilers.

Temporal transcriptomes of a marine cyanopodovirus and its Synechococcus host during infection

Marine picocyanobacteria belonging to genera Synechococcus and Prochlorococcus are genetically diverged and distributed into distinct biogeographical patterns, and both are infected by genetically closely related cyanopodoviruses. Previous studies have not fully explored whether the two virus–host systems share similar gene expression patterns during infection. Whole‐genome expression dynamics of T7‐like cyanopodovirus P‐SSP7 and its host Prochlorococcus strain MED4 have already been reported. Here, we conducted genomic and transcriptomic analyses on T7‐like cyanopodovirus S‐SBP1 during its infection on Synechococcus strain WH7803. S‐SBP1 has a latent period of 8 h and phage DNA production of 30 copies per cell. In terms of whole‐genome phylogenetic relationships and average nucleotide identity, S‐SBP1 was most similar to cyanopodovirus S‐RIP2, which also infects Synechococcus WH7803. Three hypervariable genomic islands were identified when comparing the genomes of S‐SBP1 and S‐RIP2. Single nucleotide variants were also observed in three S‐SBP1 genes, which were located within the island regions. Based on RNA‐seq analysis, S‐SBP1 genes clustered into three temporal expression classes, whose gene content was similar to that of P‐SSP7. Thirty‐two host genes were upregulated during phage infection, including those involved in carbon metabolism, ribosome components, and stress response. These upregulated genes were similar to those upregulated by Prochlorococcus MED4 in response to infection by P‐SSP7. Our study demonstrates a programmed temporal expression pattern of cyanopodoviruses and hosts during infection.

Mitochondrial TSPO Deficiency Triggers Retrograde Signaling in MA-10 Mouse Tumor Leydig Cells

The mitochondrial translocator protein (TSPO) has been shown to bind cholesterol with high affinity and is involved in mediating its availability for steroidogenesis. We recently reported that targeted Tspo gene deletion in MA-10 mouse tumor Leydig cells resulted in reduced cAMP-stimulated steroid formation and significant reduction in the mitochondrial membrane potential (ΔΨ_m) compared to control cells. We hypothesized that ΔΨ_m reduction in the absence of TSPO probably reflects the dysregulation and/or maintenance failure of some basic mitochondrial function(s). To explore the consequences of TSPO depletion via CRISPR-Cas9-mediated deletion (indel) mutation in MA-10 cells, we assessed the transcriptome changes in TSPO-mutant versus wild-type (Wt) cells using RNA-seq. Gene expression profiles were validated using real-time PCR. We report herein that there are significant changes in nuclear gene expression in Tspo mutant versus Wt cells. The identified transcriptome changes were mapped to several signaling pathways including the regulation of membrane potential, calcium signaling, extracellular matrix, and phagocytosis. This is a retrograde signaling pathway from the mitochondria to the nucleus and is probably the result of changes in expression of several transcription factors, including key members of the NF-κB pathway. In conclusion, TSPO regulates nuclear gene expression through intracellular signaling. This is the first evidence of a compensatory response to the loss of TSPO with transcriptome changes at the cellular level.