pcaExplorer: an R/Bioconductor package for interacting with RNA-seq principal components

Differences in gene expression profiles in early and late stage rhodesiense HAT individuals in Malawi

T. b. rhodesiense is the causative agent of Rhodesian human African trypanosomiasis (r-HAT) in Malawi. Clinical presentation of r-HAT in Malawi varies between foci and differs from East African HAT clinical phenotypes. The purpose of this study was to gain more insights into the transcriptomic profiles of patients with early stage 1 and late stage 2 HAT disease in Malawi. Whole blood from individuals infected with T. b. rhodesiense was used for RNA-Seq. Control samples were from healthy trypanosome negative individuals matched on sex, age range, and disease foci. Illumina sequence FASTQ reads were aligned to the GRCh38 release 84 human genome sequence using HiSat2 and differential analysis was done in R Studio using the DESeq2 package. XGR, ExpressAnalyst and InnateDB algorithms were used for functional annotation and gene enrichment analysis of significant differentially expressed genes. RNA-seq was done on 23 r-HAT case samples and 28 healthy controls with 7 controls excluded for downstream analysis as outliers. A total of 4519 genes were significant differentially expressed (p adjusted <0.05) in individuals with early stage 1 r-HAT disease (n = 12) and 1824 genes in individuals with late stage 2 r-HAT disease (n = 11) compared to controls. Enrichment of innate immune response genes through neutrophil activation was identified in individuals with both early and late stages of the disease. Additionally, lipid metabolism genes were enriched in late stage 2 disease. We further identified uniquely upregulated genes (log2 Fold Change 1.4-2.0) in stage 1 (ZNF354C) and stage 2 (TCN1 and MAGI3) blood. Our data add to the current understanding of the human transcriptome profiles during T. b. rhodesiense infection. We further identified biological pathways and transcripts enriched than were enriched during stage 1 and stage 2 r-HAT. Lastly, we have identified transcripts which should be explored in future research whether they have potential of being used in combination with other markers for staging or r-HAT.

Complementary HLH susceptibility factors converge on CD8 T-cell hyperactivation

Hemophagocytic lymphohistiocytosis (HLH) and macrophage activation syndrome (MAS) are life-threatening hyperinflammatory syndromes. Familial HLH is caused by genetic impairment of granule-mediated cytotoxicity (eg, perforin deficiency). MAS is linked to excess activity of the inflammasome-activated cytokine interleukin-18 (IL-18). Though individually tolerated, mice with dual susceptibility (Prf1⁻/⁻Il18tg; DS) succumb to spontaneous, lethal hyperinflammation. We hypothesized that understanding how these susceptibility factors synergize would uncover key pathomechanisms in the activation, function, and persistence of hyperactivated CD8 T cells. In IL-18 transgenic (Il18tg) mice, IL-18 effects on CD8 T cells drove MAS after a viral (lymphocytic choriomeningitis virus), but not innate (toll like receptor 9), trigger. In vitro, CD8 T cells also required T-cell receptor (TCR) stimulation to fully respond to IL-18. IL-18 induced but perforin deficiency impaired immunoregulatory restimulation-induced cell death (RICD). Paralleling hyperinflammation, DS mice displayed massive postthymic oligoclonal CD8 T-cell hyperactivation in their spleens, livers, and bone marrow as early as 3 weeks. These cells increased proliferation and interferon gamma production, which contrasted with increased expression of receptors and transcription factors associated with exhaustion. Broad-spectrum antibiotics and antiretrovirals failed to ameliorate the disease. Attempting to genetically "fix" TCR antigen-specificity instead demonstrated the persistence of spontaneous HLH and hyperactivation, chiefly on T cells that had evaded TCR fixation. Thus, drivers of HLH may preferentially act on CD8 T cells: IL-18 amplifies activation and demand for RICD, whereas perforin supplies critical immunoregulation. Together, these factors promote a terminal CD8 T-cell activation state, combining features of exhaustion and effector function. Therefore, susceptibility to hyperinflammation may converge on a unique, unrelenting, and antigen-dependent state of CD8 T-cell hyperactivation.

Exploring the Potential of Olfactory Receptor Circulating RNA Measurement for Preeclampsia Prediction and Its Linkage to Mild Gestational Hypothyroidism

Gestational hypothyroidism may lead to preeclampsia development. However, this pathophysiological is unknown. We expect to find a shared mechanism by comparing hypothyroidism and preeclampsia. From our transcriptome data, we recognized olfactory receptors as that fingerprint. The reduction of taste and smell in hypothyroid patients has been known for a long time. Therefore, we decided to look to the olfactory receptors and aimed to identify genes capable of predicting preeclampsia (PEC). Methods: An Ion Proton Sequencer (Thermo Fisher Scientific, Waltham, MA, USA) was used to construct the transcriptome databases. RStudio with packages Limma v.3.50.0, GEOquery v.2.62.2, and umap v.0.2.8.8 were used to analyze the differentially expressed genes in GSE149440 from the Gene Expression Omnibus (GEO). The 7500 Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) was used for RT-qPCR amplification of OR6X1 and OR4E2. Results: Our transcriptomic datasets analysis revealed 25.08% and 26.75% downregulated olfactory receptor (ORs) in mild nontreated gestational hypothyroidism (GHT) and PEC, respectively. In the GSE149440 GEO dataset, we found OR5H1, OR5T3, OR51A7, OR51B6, OR10J5, OR6C6, and OR2AG2 as predictors of early-onset PEC. We also evaluate two chosen biomarkers' responses to levothyroxine. The RT-qPCR demonstrated a difference in OR6X1 and OR4E2 expression between GHT and healthy pregnancy (p < 0.05). Those genes presented a negative correlation with TSH (r: -0.51, p < 0.05; and r: -0.44, p < 0.05), a strong positive correlation with each other (r: 0.89; p < 0.01) and the levothyroxine-treated group had no difference from the healthy one. We conclude that ORs could be used as biomarkers at the beginning of gestation, and the downregulated ORs found in GHT may be improved with levothyroxine treatment.

Vaginal microbes alter epithelial transcriptomic and epigenomic modifications providing insight into the molecular mechanisms for susceptibility to adverse reproductive outcomes

The cervicovaginal microbiome is highly associated with women's health with microbial communities dominated by Lactobacillus spp. being considered optimal. Conversely, a lack of lactobacilli and a high abundance of strict and facultative anaerobes including Gardnerella vaginalis , have been associated with adverse reproductive outcomes. However, the molecular pathways modulated by microbe interactions with the cervicovaginal epithelia remain unclear. Using RNA-sequencing, we characterize the in vitro cervicovaginal epithelial transcriptional response to different vaginal bacteria and their culture supernatants. We showed that G. vaginalis upregulated genes were associated with an activated innate immune response including anti-microbial peptides and inflammasome pathways, represented by NLRP3-mediated increases in caspase-1, IL-1β and cell death. Cervicovaginal epithelial cells exposed to L. crispatus showed limited transcriptomic changes, while exposure to L. crispatus culture supernatants resulted in a shift in the epigenomic landscape of cervical epithelial cells. ATAC-sequencing confirmed epigenetic changes with reduced chromatin accessibility. This study reveals new insight into host-microbe interactions in the lower reproductive tract and suggest potential therapeutic strategies leveraging the vaginal microbiome to improve reproductive health.

Comparative analysis of syngeneic mouse models of high-grade serous ovarian cancer

Ovarian cancers exhibit high rates of recurrence and poor treatment response. Preclinical models that recapitulate human disease are critical to develop new therapeutic approaches. Syngeneic mouse models allow for the generation of tumours comprising the full repertoire of non-malignant cell types but have expanded in number, varying in the cell type of origin, method for transformation, and ultimately, the properties of the tumours they produce. Here we have performed a comparative analysis of high-grade serous ovarian cancer models based on transcriptomic profiling of 22 cell line models, and intrabursal and intraperitoneal tumours from 12. Among cell lines, we identify distinct signalling activity, such as elevated inflammatory signalling in STOSE and OVE16 models, and MAPK/ERK signalling in ID8 and OVE4 models; metabolic differences, such as reduced glycolysis-associated expression in several engineered ID8 subclones; and relevant functional properties, including differences in EMT activation, PD-L1 and MHC class I expression, and predicted chemosensitivity. Among tumour samples, we observe increased variability and stromal content among intrabursal tumours. Finally, we predict differences in the microenvironment of ID8 models engineered with clinically relevant mutations. We anticipate that this work will serve as a valuable resource, providing new insight to help select models for specific experimental objectives.

Transcriptome analysis of peripheral blood of Schistosoma mansoni infected children from the Albert Nile region in Uganda reveals genes implicated in fibrosis pathology

Over 290 million people are infected by schistosomes worldwide. Schistosomiasis control efforts focus on mass drug treatment with praziquantel (PZQ), a drug that kills the adult worm of all Schistosoma species. Nonetheless, re-infections have continued to be detected in endemic areas with individuals living in the same area presenting with varying infection intensities. Our objective was to characterize the transcriptome profiles in peripheral blood of children between 10-15 years with varying intensities of Schistosoma mansoni infection living along the Albert Nile in Uganda. RNA extracted from peripheral blood collected from 44 S. mansoni infected (34 high and 10 low by circulating anodic antigen [CAA] level) and 20 uninfected children was sequenced using Illumina NovaSeq S4 and the reads aligned to the GRCh38 human genome. Differential gene expression analysis was done using DESeq2. Principal component analysis revealed clustering of gene expression by gender when S. mansoni infected children were compared with uninfected children. In addition, we identified 14 DEGs between S. mansoni infected and uninfected individuals, 56 DEGs between children with high infection intensity and uninfected individuals, 33 DEGs between those with high infection intensity and low infection intensity and no DEGs between those with low infection and uninfected individuals. We also observed upregulation and downregulation of some DEGs that are associated with fibrosis and its regulation. These data suggest expression of fibrosis associated genes as well as genes that regulate fibrosis in S. mansoni infection. The relatively few significant DEGS observed in children with schistosomiasis suggests that chronic S. mansoni infection is a stealth infection that does not stimulate a strong immune response.

PlateletSeq: A novel method for discovery of blood-based biomarkers

Platelets are small circulating fragments of cells that play important roles in thrombosis, haemostasis, immune response, inflammation and cancer growth. Although anucleate, they contain a rich RNA repertoire which offers an opportunity to characterise changes in platelet gene expression in health and disease. Whilst this can be achieved with conventional RNA sequencing, a large input of high-quality RNA, and hence blood volume, is required (unless a pre-amplification step is added), along with specialist bioinformatic skills for data analysis and interpretation. We have developed a transcriptomics next-generation sequencing-based approach that overcomes these limitations. Termed PlateletSeq, this method requires very low levels of RNA input and does not require specialist bioinformatic analytical skills. Here we describe the methodology, from sample collection to processing and data analysis. Specifically, blood samples can be stored for up to 8 days at 4 °C prior to analysis. Platelets are isolated using multi-step centrifugation and a purity of ≤ 1 leucocyte per 0.26x106 platelets is optimal for gene expression analysis. We have applied PlateletSeq to normal adult blood samples and show there are no age-associated variations and only minor gender-associated differences. In contrast, platelets from patients with myeloproliferative neoplasms show differences in platelet transcript profiles from normal and between disease subtypes. This illustrates the potential applicability of PlateletSeq for biomarker discovery and studying platelet biology in patient samples. It also opens avenues for assessing platelet quality in other fields such as transfusion research.

Genomic and phenotypic comparison of two variants of multidrug-resistant Salmonella enterica serovar Heidelberg isolated during the 2015-2017 multi-state outbreak in cattle

Salmonella enterica subspecies enterica serovar Heidelberg (Salmonella Heidelberg) has caused several multistate foodborne outbreaks in the United States, largely associated with the consumption of poultry. However, a 2015-2017 multidrug-resistant (MDR) Salmonella Heidelberg outbreak was linked to contact with dairy beef calves. Traceback investigations revealed calves infected with outbreak strains of Salmonella Heidelberg exhibited symptoms of disease frequently followed by death from septicemia. To investigate virulence characteristics of Salmonella Heidelberg as a pathogen in bovine, two variants with distinct pulse-field gel electrophoresis (PFGE) patterns that differed in morbidity and mortality during the multistate outbreak were genotypically and phenotypically characterized and compared. Strain SX 245 with PFGE pattern JF6X01.0523 was identified as a dominant and highly pathogenic variant causing high morbidity and mortality in affected calves, whereas strain SX 244 with PFGE pattern JF6X01.0590 was classified as a low pathogenic variant causing less morbidity and mortality. Comparison of whole-genome sequences determined that SX 245 lacked ~200 genes present in SX 244, including genes associated with the IncI1 plasmid and phages; SX 244 lacked eight genes present in SX 245 including a second YdiV Anti-FlhC(2)FlhD(4) factor, a lysin motif domain containing protein, and a pentapeptide repeat protein. RNA-sequencing revealed fimbriae-related, flagella-related, and chemotaxis genes had increased expression in SX 245 compared to SX 244. Furthermore, SX 245 displayed higher invasion of human and bovine epithelial cells than SX 244. These data suggest that the presence and up-regulation of genes involved in type 1 fimbriae production, flagellar regulation and biogenesis, and chemotaxis may play a role in the increased pathogenicity and host range expansion of the Salmonella Heidelberg isolates involved in the bovine-related outbreak.

A transcriptional and regulatory map of mouse somite maturation

The mammalian body plan is shaped by rhythmic segmentation of mesoderm into somites, which are transient embryonic structures that form down each side of the neural tube. We have analyzed the genome-wide transcriptional and chromatin dynamics occurring within nascent somites, from early inception of somitogenesis to the latest stages of body plan establishment. We created matched gene expression and open chromatin maps for the three leading pairs of somites at six time points during mouse embryonic development. We show that the rate of somite differentiation accelerates as development progresses. We identified a conserved maturation program followed by all somites, but somites from more developed embryos concomitantly switch on differentiation programs from derivative cell lineages soon after segmentation. Integrated analysis of the somitic transcriptional and chromatin activities identified opposing regulatory modules controlling the onset of differentiation. Our results provide a powerful, high-resolution view of the molecular genetics underlying somitic development in mammals.

Characterization of 3D organotypic epithelial tissues reveals tonsil-specific differences in tonic interferon signaling

Three-dimensional (3D) culturing techniques can recapitulate the stratified nature of multicellular epithelial tissues. Organotypic 3D epithelial tissue culture methods have several applications, including the study of tissue development and function, drug discovery and toxicity testing, host-pathogen interactions, and the development of tissue-engineered constructs for use in regenerative medicine. We grew 3D organotypic epithelial tissues from foreskin, cervix, and tonsil-derived primary cells and characterized the transcriptome of these in vitro tissue equivalents. Using the same 3D culturing method, all three tissues yielded stratified squamous epithelium, validated histologically using basal and superficial epithelial cell markers. The goal of this study was to use RNA-seq to compare gene expression patterns in these three types of epithelial tissues to gain a better understanding of the molecular mechanisms underlying their function and identify potential therapeutic targets for various diseases. Functional profiling by over-representation and gene set enrichment analysis revealed tissue-specific differences: i.e., cutaneous homeostasis and lipid metabolism in foreskin, extracellular matrix remodeling in cervix, and baseline innate immune differences in tonsil. Specifically, tonsillar epithelia may play an active role in shaping the immune microenvironment of the tonsil balancing inflammation and immune responses in the face of constant exposure to microbial insults. Overall, these data serve as a resource, with gene sets made available for the research community to explore, and as a foundation for understanding the epithelial heterogeneity and how it may impact their in vitro use. An online resource is available to investigate these data (https://viz.datascience.arizona.edu/3DEpiEx/).

ADAD1 is required for normal translation of nuclear pore and transport protein transcripts in spermatids of Mus musculus†

ADAD1 is a testis-specific RNA-binding protein expressed in post-meiotic spermatids whose loss leads to defective sperm and male infertility. However, the drivers of the Adad1 phenotype remain unclear. Morphological and functional analysis of Adad1 mutant sperm showed defective DNA compaction, abnormal head shaping, and reduced motility. Mutant testes demonstrated minimal transcriptome changes; however, ribosome association of many transcripts was reduced, suggesting ADAD1 may be required for their translational activation. Further, immunofluorescence of proteins encoded by select transcripts showed delayed protein accumulation. Additional analyses demonstrated impaired subcellular localization of multiple proteins, suggesting protein transport is also abnormal in Adad1 mutants. To clarify the mechanism giving rise to this, the manchette, a protein transport microtubule network, and the LINC (linker of nucleoskeleton and cytoskeleton) complex, which connects the manchette to the nuclear lamin, were assessed across spermatid development. Proteins of both displayed delayed translation and/or localization in mutant spermatids implicating ADAD1 in their regulation, even in the absence of altered ribosome association. Finally, ADAD1's impact on the NPC (nuclear pore complex), a regulator of both the manchette and the LINC complex, was examined. Reduced ribosome association of NPC encoding transcripts and reduced NPC protein abundance along with abnormal localization in Adad1 mutants confirmed ADAD1's impact on translation is required for a NPC in post-meiotic germ cells. Together, these studies lead to a model whereby ADAD1's influence on nuclear transport leads to deregulation of the LINC complex and the manchette, ultimately generating the range of physiological defects observed in the Adad1 phenotype.

Investigating the role of the ROS/CncC signaling pathway in the response to xenobiotics in Spodoptera frugiperda using Sf9 cells

Spodoptera frugiperda (fall armyworm, FAW) is an invasive polyphagous lepidopteran pest that has developed sophisticated resistance mechanisms involving detoxification enzymes to eliminate toxic compounds it encounters in its diet including insecticides. Although its inventory of detoxification enzymes is known, the mechanisms that enable an adapted response depending on the toxic compound remain largely unexplored. Sf9 cells were used to investigate the role of the transcription factors, Cap n' collar isoform C (CncC) and musculoaponeurotic fibrosarcoma (Maf) in the regulation of the detoxification response. We overexpressed CncC, Maf or both genes, and knocked out (KO) CncC or its repressor Kelch-like ECH associated protein 1 (Keap1). Joint overexpression of CncC and Maf is required to confer increased tolerance to indole 3-carbinol (I3C), a plant secondary metabolite, and to methoprene, an insecticide. Both molecules induce reactive oxygen species (ROS) pulses in the different cell lines. The use of an antioxidant reversed ROS pulses and restored the tolerance to I3C and methoprene. The activity of detoxification enzymes varied according to the cell line. Suppression of Keap1 significantly increased the activity of cytochrome P450s, carboxylesterases and glutathione S-transferases. RNAseq experiments showed that CncC mainly regulates the expression of detoxification genes but is also at the crossroads of several signaling pathways (reproduction and immunity) maintaining homeostasis. We present new data in Sf9 cell lines suggesting that the CncC:Maf pathway plays a central role in FAW response to natural and synthetic xenobiotics. This knowledge helps to better understand detoxification gene expression and may help to design next-generation pest insect control measures.

Analysis of subcellular RNA fractions demonstrates significant genetic regulation of gene expression in human brain post-transcriptionally

Gaining insight into the genetic regulation of gene expression in human brain is key to the interpretation of genome-wide association studies for major neurological and neuropsychiatric diseases. Expression quantitative trait loci (eQTL) analyses have largely been used to achieve this, providing valuable insights into the genetic regulation of steady-state RNA in human brain, but not distinguishing between molecular processes regulating transcription and stability. RNA quantification within cellular fractions can disentangle these processes in cell types and tissues which are challenging to model in vitro. We investigated the underlying molecular processes driving the genetic regulation of gene expression specific to a cellular fraction using allele-specific expression (ASE). Applying ASE analysis to genomic and transcriptomic data from paired nuclear and cytoplasmic fractions of anterior prefrontal cortex, cerebellar cortex and putamen tissues from 4 post-mortem neuropathologically-confirmed control human brains, we demonstrate that a significant proportion of genetic regulation of gene expression occurs post-transcriptionally in the cytoplasm, with genes undergoing this form of regulation more likely to be synaptic. These findings have implications for understanding the structure of gene expression regulation in human brain, and importantly the interpretation of rapidly growing single-nucleus brain RNA-sequencing and eQTL datasets, where cytoplasm-specific regulatory events could be missed.

The Effect of Plastic-Related Compounds on Transcriptome-Wide Gene Expression on CYP2C19-Overexpressing HepG2 Cells

In recent years, plastic and especially microplastic in the oceans have caused huge problems to marine flora and fauna. Recently, such particles have also been detected in blood, breast milk, and placenta, underlining their ability to enter the human body, presumably via the food chain and other yet-unknown mechanisms. In addition, plastic contains plasticizers, antioxidants, or lubricants, whose impact on human health is also under investigation. At the cellular level, the most important enzymes involved in the metabolism of xenobiotic compounds are the cytochrome P450 monooxygenases (CYPs). Despite their extensive characterization in the maintenance of cellular balance, their interactions with plastic and related products are unexplored. In this study, the possible interactions between several plastic-related compounds and one of the most important cytochromes, CYP2C19, were analyzed. By applying virtual compound screening and molecular docking to more than 1000 commercially available plastic-related compounds, we identified candidates that are likely to interact with this protein. A growth inhibition assay confirmed their cytotoxic activity on a CYP2C19-transfected hepatic cell line. Subsequently, we studied the effect of the selected compounds on the transcriptome-wide gene expression level by conducting RNA sequencing. Three candidate molecules were identified, i.e., 2,2'-methylene bis(6-tert-butyl-4-methylphenol), 1,1-bis(3,5-di-tert-butyl-2-hydroxyphenyl) ethane, and 2,2'-methylene bis(6-cyclohexyl-4-methylphenol)), which bound with a high affinity to CYP2C19 in silico. They exerted a profound cytotoxicity in vitro and interacted with several metabolic pathways, of which the 'cholesterol biosynthesis process' was the most affected. In addition, other affected pathways involved mitosis, DNA replication, and inflammation, suggesting an increase in hepatotoxicity. These results indicate that plastic-related compounds could damage the liver by affecting several molecular pathways.

Integrative analysis of transcriptome dynamics during human craniofacial development identifies candidate disease genes

Craniofacial disorders arise in early pregnancy and are one of the most common congenital defects. To fully understand how craniofacial disorders arise, it is essential to characterize gene expression during the patterning of the craniofacial region. To address this, we performed bulk and single-cell RNA-seq on human craniofacial tissue from 4-8 weeks post conception. Comparisons to dozens of other human tissues revealed 239 genes most strongly expressed during craniofacial development. Craniofacial-biased developmental enhancers were enriched +/- 400 kb surrounding these craniofacial-biased genes. Gene co-expression analysis revealed that regulatory hubs are enriched for known disease causing genes and are resistant to mutation in the normal healthy population. Combining transcriptomic and epigenomic data we identified 539 genes likely to contribute to craniofacial disorders. While most have not been previously implicated in craniofacial disorders, we demonstrate this set of genes has increased levels of de novo mutations in orofacial clefting patients warranting further study.

Transcriptional analysis in multiple barley varieties identifies signatures of waterlogging response

Waterlogging leads to major crop losses globally, particularly for waterlogging-sensitive crops such as barley. Waterlogging reduces oxygen availability and results in additional stresses, leading to the activation of hypoxia and stress response pathways that promote plant survival. Although certain barley varieties have been shown to be more tolerant to waterlogging than others and some tolerance-related quantitative trait loci have been identified, the molecular mechanisms underlying this trait are mostly unknown. Transcriptomics approaches can provide very valuable information for our understanding of waterlogging tolerance. Here, we surveyed 21 barley varieties for the differential transcriptional activation of conserved hypoxia-response genes under waterlogging and selected five varieties with different levels of induction of core hypoxia-response genes. We further characterized their phenotypic response to waterlogging in terms of shoot and root traits. RNA sequencing to evaluate the genome-wide transcriptional responses to waterlogging of these selected varieties led to the identification of a set of 98 waterlogging-response genes common to the different datasets. Many of these genes are orthologs of the so-called "core hypoxia response genes," thus highlighting the conservation of plant responses to waterlogging. Hierarchical clustering analysis also identified groups of genes with intrinsic differential expression between varieties prior to waterlogging stress. These genes could constitute interesting candidates to study "predisposition" to waterlogging tolerance or sensitivity in barley.

Interaction of high-fat diet and brain trauma alters adipose tissue macrophages and brain microglia associated with exacerbated cognitive dysfunction

Obesity increases the morbidity and mortality of traumatic brain injury (TBI). We performed a detailed analysis of transcriptomic changes in the brain and adipose tissue to examine the interactive effects between high-fat diet-induced obesity (DIO) and TBI in relation to central and peripheral inflammatory pathways, as well as neurological function. Adult male mice were fed a high-fat diet (HFD) for 12 weeks prior to experimental TBI and continuing after injury. Combined TBI and HFD resulted in additive dysfunction in the Y-Maze, novel object recognition (NOR), and Morris water maze (MWM) cognitive function tests. We also performed high-throughput transcriptomic analysis using Nanostring panels of cellular compartments in the brain and total visceral adipose tissue (VAT), followed by unsupervised clustering, principal component analysis, and IPA pathway analysis to determine shifts in gene expression programs and molecular pathway activity. Analysis of cellular populations in the cortex and hippocampus as well as in visceral adipose tissue during the chronic phase after combined TBI-HFD showed amplification of central and peripheral microglia/macrophage responses, including superadditive changes in select gene expression signatures and pathways. These data suggest that HFD-induced obesity and TBI can independently prime and support the development of altered states in brain microglia and visceral adipose tissue macrophages, including the disease-associated microglia/macrophage (DAM) phenotype observed in neurodegenerative disorders. The interaction between HFD and TBI promotes a shift toward chronic reactive microglia/macrophage transcriptomic signatures and associated pro-inflammatory disease-altered states that may, in part, underlie the exacerbation of cognitive deficits. Targeting of HFD-induced reactive cellular phenotypes, including in peripheral adipose tissue macrophages, may serve to reduce microglial maladaptive states after TBI, attenuating post-traumatic neurodegeneration and neurological dysfunction.

Inhibiting NLRP3 signaling in aging podocytes improves their life- and health-span

The decrease in the podocyte's lifespan and health-span that typify healthy kidney aging cause a decrease in their normal structure, physiology and function. The ability to halt and even reverse these changes becomes clinically relevant when disease is superimposed on an aged kidney. RNA-sequencing of podocytes from middle-aged mice showed an inflammatory phenotype with increases in the NLRP3 inflammasome, signaling for IL2/Stat5, IL6 and TNF, interferon gamma response, allograft rejection and complement, consistent with inflammaging. Furthermore, injury-induced NLRP3 signaling in podocytes was further augmented in aged mice compared to young ones. The NLRP3 inflammasome (NLRP3, Caspase-1, IL1β IL-18) was also increased in podocytes of middle-aged humans. Higher transcript expression for NLRP3 in human glomeruli was accompanied by reduced podocyte density and increased global glomerulosclerosis and glomerular volume. Pharmacological inhibition of NLRP3 with MCC950, or gene deletion, reduced podocyte senescence and the genes typifying aging in middle-aged mice, which was accompanied by an improved podocyte lifespan and health-span. Moreover, modeling the injury-dependent increase in NLRP3 signaling in human kidney organoids confirmed the anti-senescence effect of MC9950. Finally, NLRP3 also impacted liver aging. Together, these results suggest a critical role for the NLRP3 inflammasome in podocyte and liver aging.

Screening autism-associated environmental factors in differentiating human neural progenitors with fractional factorial design-based transcriptomics

Research continues to identify genetic variation, environmental exposures, and their mixtures underlying different diseases and conditions. There is a need for screening methods to understand the molecular outcomes of such factors. Here, we investigate a highly efficient and multiplexable, fractional factorial experimental design (FFED) to study six environmental factors (lead, valproic acid, bisphenol A, ethanol, fluoxetine hydrochloride and zinc deficiency) and four human induced pluripotent stem cell line derived differentiating human neural progenitors. We showcase the FFED coupled with RNA-sequencing to identify the effects of low-grade exposures to these environmental factors and analyse the results in the context of autism spectrum disorder (ASD). We performed this after 5-day exposures on differentiating human neural progenitors accompanied by a layered analytical approach and detected several convergent and divergent, gene and pathway level responses. We revealed significant upregulation of pathways related to synaptic function and lipid metabolism following lead and fluoxetine exposure, respectively. Moreover, fluoxetine exposure elevated several fatty acids when validated using mass spectrometry-based metabolomics. Our study demonstrates that the FFED can be used for multiplexed transcriptomic analyses to detect relevant pathway-level changes in human neural development caused by low-grade environmental risk factors. Future studies will require multiple cell lines with different genetic backgrounds for characterising the effects of environmental exposures in ASD.

Comparative transcriptomic analysis of male and females in the dioecious weeds Amaranthus palmeri and Amaranthus tuberculatus

BACKGROUND

Waterhemp (Amaranthus tuberculatus (Moq.) Sauer) and Palmer amaranth (Amaranthus palmeri S. Wats.) are two dioecious and important weed species in the world that can rapidly evolve herbicide-resistance traits. Understanding these two species' dioecious and sex-determination mechanisms could open opportunities for new tools to control them. This study aims to identify the differential expression patterns between males and females in A. tuberculatus and A. palmeri. Multiple analyses, including differential expression, co-expression, and promoter analyses, used RNA-seq data from multiple tissue types to identify putative essential genes for sex determination in both dioecious species.

RESULTS

Genes were identified as potential key players for sex determination in A. palmeri. Genes PPR247, WEX, and ACD6 were differentially expressed between the sexes and located at scaffold 20 within or near the male-specific Y (MSY) region. Multiple genes involved with flower development were co-expressed with these three genes. For A. tuberculatus, no differentially expressed gene was identified within the MSY region; however, multiple autosomal class B and C genes were identified as differentially expressed and possible candidate genes.

CONCLUSIONS

This is the first study comparing the global expression profile between males and females in dioecious weedy Amaranthus species. Results narrow down putative essential genes for sex-determination in A. palmeri and A. tuberculatus and also strengthen the hypothesis of two different evolutionary events for dioecy within the genus.

Twist1-IRF9 Interaction Is Necessary for IFN-Stimulated Gene Anti-Zika Viral Infection

An efficient immune defense against pathogens requires sufficient basal sensing mechanisms that can deliver prompt responses. Type I IFNs are protective against acute viral infections and respond to viral and bacterial infections, but their efficacy depends on constitutive basal activity that promotes the expression of downstream genes known as IFN-stimulated genes (ISGs). Type I IFNs and ISGs are constitutively produced at low quantities and yet exert profound effects essential for numerous physiological processes beyond antiviral and antimicrobial defense, including immunomodulation, cell cycle regulation, cell survival, and cell differentiation. Although the canonical response pathway for type I IFNs has been extensively characterized, less is known regarding the transcriptional regulation of constitutive ISG expression. Zika virus (ZIKV) infection is a major risk for human pregnancy complications and fetal development and depends on an appropriate IFN-β response. However, it is poorly understood how ZIKV, despite an IFN-β response, causes miscarriages. We have uncovered a mechanism for this function specifically in the context of the early antiviral response. Our results demonstrate that IFN regulatory factor (IRF9) is critical in the early response to ZIKV infection in human trophoblast. This function is contingent on IRF9 binding to Twist1. In this signaling cascade, Twist1 was not only a required partner that promotes IRF9 binding to the IFN-stimulated response element but also an upstream regulator that controls basal levels of IRF9. The absence of Twist1 renders human trophoblast cells susceptible to ZIKV infection.

Glucosylceramide Synthase Inhibitors Induce Ceramide Accumulation and Sensitize H3K27 Mutant Diffuse Midline Glioma to Irradiation

H3K27M mutant (mut) diffuse midline glioma (DMG) is a lethal cancer with no effective cure. The glycosphingolipids (GSL) metabolism is altered in these tumors and could be exploited to develop new therapies. We tested the effect of the glucosylceramide synthase inhibitors (GSI) miglustat and eliglustat on cell proliferation, alone or in combination with temozolomide or ionizing radiation. Miglustat was included in the therapy protocol of two pediatric patients. The effect of H3.3K27 trimethylation on GSL composition was analyzed in ependymoma. GSI reduced the expression of the ganglioside GD2 in a concentration and time-dependent manner and increased the expression of ceramide, ceramide 1-phosphate, sphingosine, and sphingomyelin but not of sphingosine 1-phosphate. Miglustat significantly increased the efficacy of irradiation. Treatment with miglustat according to dose recommendations for patients with Niemann-Pick disease was well tolerated with manageable toxicities. One patient showed a mixed response. In ependymoma, a high concentration of GD2 was found only in the presence of the loss of H3.3K27 trimethylation. In conclusion, treatment with miglustat and, in general, targeting GSL metabolism may offer a new therapeutic opportunity and can be administered in close proximity to radiation therapy. Alterations in H3K27 could be useful to identify patients with a deregulated GSL metabolism.

Cell Type Specific CAG Repeat Expansions and Toxicity of Mutant Huntingtin in Human Striatum and Cerebellum

Brain region-specific degeneration and somatic expansions of the mutant Huntingtin (mHTT) CAG tract are key features of Huntington's disease (HD). However, the relationships between CAG expansions, death of specific cell types, and molecular events associated with these processes are not established. Here we employed fluorescence-activated nuclear sorting (FANS) and deep molecular profiling to gain insight into the properties of cell types of the human striatum and cerebellum in HD and control donors. CAG expansions arise in striatal medium spiny neurons (MSNs) and cholinergic interneurons, in cerebellar Purkinje neurons, and at mATXN3 in MSNs from SCA3 donors. CAG expansions in MSNs are associated with higher levels of MSH2 and MSH3 (forming MutSβ), which can inhibit nucleolytic excision of CAG slip-outs by FAN1 in a concentration-dependent manner. Our data indicate that ongoing CAG expansions are not sufficient for cell death, and identify transcriptional changes associated with somatic CAG expansions and striatal toxicity.

Gene coexpression network during ontogeny in the yellow fever mosquito, Aedes aegypti

BACKGROUND

The behaviors and ontogeny of Aedes aegypti are closely related to the spread of diseases caused by dengue (DENV), chikungunya (CHIKV), Zika (ZIKV), and yellow fever (YFV) viruses. During the life cycle, Ae. aegypti undergoes drastic morphological, metabolic, and functional changes triggered by gene regulation and other molecular mechanisms. Some essential regulatory factors that regulate insect ontogeny have been revealed in other species, but their roles are still poorly investigated in the mosquito.

RESULTS

Our study identified 6 gene modules and their intramodular hub genes that were highly associated with the ontogeny of Ae. aegypti in the constructed network. Those modules were found to be enriched in functional roles related to cuticle development, ATP generation, digestion, immunity, pupation control, lectins, and spermatogenesis. Additionally, digestion-related pathways were activated in the larvae and adult females but suppressed in the pupae. The integrated protein‒protein network also identified cilium-related genes. In addition, we verified that the 6 intramodular hub genes encoding proteins such as EcKinase regulating larval molt were only expressed in the larval stage. Quantitative RT‒PCR of the intramodular hub genes gave similar results as the RNA-Seq expression profile, and most hub genes were ontogeny-specifically expressed.

CONCLUSIONS

The constructed gene coexpression network provides a useful resource for network-based data mining to identify candidate genes for functional studies. Ultimately, these findings will be key in identifying potential molecular targets for disease control.

Low-Protein Diets Composed of Protein Recovered from Food Processing Supported Growth, but Induced Mild Hepatic Steatosis Compared with a No-Protein Diet in Young Female Rats

BACKGROUND

Living in low-income countries often restricts the consumption of adequate protein and animal protein.

OBJECTIVES

This study aimed to investigate the effects of feeding low-protein diets on growth and liver health using proteins recovered from animal processing.

METHODS

Female Sprague-Dawley rats (aged 28 d) were randomly assigned (n = 8 rats/group) to be fed standard purified diets with 0% or 10% kcal protein that was comprised of either carp, whey, or casein.

RESULTS

Rats that were fed low-protein diets showed higher growth but developed mild hepatic steatosis compared to rats that were fed a no-protein diet, regardless of the protein source. Real-time quantitative polymerase chain reactions targeting the expression of genes involved in liver lipid homeostasis were not significantly different among groups. Global RNA-sequencing technology identified 9 differentially expressed genes linked to folate-mediated 1-carbon metabolism, endoplasmic reticulum (ER) stress, and metabolic diseases. Canonical pathway analysis revealed that mechanisms differed depending on the protein source. ER stress and dysregulated energy metabolism were implicated in hepatic steatosis in carp- and whey-fed rats. In contrast, impaired liver one-carbon methylations, lipoprotein assembly, and lipid export were implicated in casein-fed rats.

CONCLUSIONS

Carp sarcoplasmic protein showed comparable results to commercially available casein and whey protein. A better understanding of the molecular mechanisms in hepatic steatosis development can assist formulation of proteins recovered from food processing into a sustainable source of high-quality protein.

Profiling of Dermanyssus gallinae genes involved in acaricide resistance

The poultry red mite (PRM), Dermanyssus gallinae, is a major threat for the poultry industry worldwide. Chemical compounds have been extensively used for PRM control, leading to selection of resistant mites. Molecular mechanisms of resistance have been investigated in arthropods, showing the role of target-site insensitivity and enhanced detoxification. Few studies are available about those mechanisms in D. gallinae, and none have yet focused on the expression levels of detoxification enzymes and other defense-related genes through RNA-seq. We tested PRM populations from Italy for their susceptibility to the acaricidal compounds phoxim and cypermethrin. Mutations in the voltage-gated sodium channel (vgsc) and in acetylcholinesterase (AChE) were investigated, detecting mutations known to be associated with acaricide/insecticide resistance in arthropods, including M827I and M918L/T in the vgsc and G119S in the AChE. RNA-seq analysis was performed to characterize metabolic resistance in fully susceptible PRM and in cypermethrin-resistant PRM exposed and unexposed to cypermethrin as well as phoxim resistant PRM exposed and unexposed to phoxim. Detoxification enzymes (including P450 monooxygenases and glutathione-S-transferases), ABC transporters and cuticular proteins were constitutively overexpressed in phoxim and cypermethrin resistant mites. In addition, heat shock proteins were found constitutively and inductively upregulated in phoxim resistant mites, while in cypermethrin resistant mites esterases and an aryl hydrocarbon receptor were constitutively highly expressed. The findings suggest that acaricide resistance in D. gallinae is due to both target-site insensitivity and overexpression of detoxification enzymes and other xenobiotic defense-related genes, which is mostly constitutive and not induced by treatment. Understanding the molecular basis of resistance could be useful to screen or test PRM populations in order to select targeted acaricides and to avoid the abuse/misuse of the few available compounds.

c-Met Mediated Cytokine Network Promotes Brain Metastasis of Breast Cancer by Remodeling Neutrophil Activities

The brain is one of the most common metastatic sites among breast cancer patients, especially in those who have Her2-positive or triple-negative tumors. The brain microenvironment has been considered immune privileged, and the exact mechanisms of how immune cells in the brain microenvironment contribute to brain metastasis remain elusive. In this study, we found that neutrophils are recruited and influenced by c-Met high brain metastatic cells in the metastatic sites, and depletion of neutrophils significantly suppressed brain metastasis in animal models. Overexpression of c-Met in tumor cells enhances the secretion of a group of cytokines, including CXCL1/2, G-CSF, and GM-CSF, which play critical roles in neutrophil attraction, granulopoiesis, and homeostasis. Meanwhile, our transcriptomic analysis demonstrated that conditioned media from c-Met high cells significantly induced the secretion of lipocalin 2 (LCN2) from neutrophils, which in turn promotes the self-renewal of cancer stem cells. Our study unveiled the molecular and pathogenic mechanisms of how crosstalk between innate immune cells and tumor cells facilitates tumor progression in the brain, which provides novel therapeutic targets for treating brain metastasis.

Acute frataxin knockdown in induced pluripotent stem cell-derived cardiomyocytes activates a type I interferon response

Friedreich ataxia, the most common hereditary ataxia, is a neuro- and cardio-degenerative disorder caused, in most cases, by decreased expression of the mitochondrial protein frataxin. Cardiomyopathy is the leading cause of premature death. Frataxin functions in the biogenesis of iron-sulfur clusters, which are prosthetic groups that are found in proteins involved in many biological processes. To study the changes associated with decreased frataxin in human cardiomyocytes, we developed a novel isogenic model by acutely knocking down frataxin, post-differentiation, in cardiomyocytes derived from induced pluripotent stem cells (iPSCs). Transcriptome analysis of four biological replicates identified severe mitochondrial dysfunction and a type I interferon response as the pathways most affected by frataxin knockdown. We confirmed that, in iPSC-derived cardiomyocytes, loss of frataxin leads to mitochondrial dysfunction. The type I interferon response was activated in multiple cell types following acute frataxin knockdown and was caused, at least in part, by release of mitochondrial DNA into the cytosol, activating the cGAS-STING sensor pathway.

Intrinsic sexual dimorphism in the placenta determines the differential response to benzene exposure

Maternal immune activation (MIA) by environmental challenges is linked to severe developmental complications, such as neurocognitive disorders, autism, and even fetal/maternal death. Benzene is a major toxic compound in air pollution that affects the mother as well as the fetus and has been associated with reproductive complications. Our objective was to elucidate whether benzene exposure during gestation triggers MIA and its impact on fetal development. We report that benzene exposure during pregnancy leads MIA associated with increased fetal resorptions, fetal growth, and abnormal placenta development. Furthermore, we demonstrate the existence of a sexual dimorphic response to benzene exposure in male and female placentas. The sexual dimorphic response is a consequence of inherent differences between male and female placenta. These data provide crucial information on the origins or sexual dimorphism and how exposure to environmental factors can have a differential impact on the development of male and female offspring.

RNA-sequencing revealed apple pomace ameliorates expression of genes in the hypothalamus associated with neurodegeneration in female rats fed a Western diet during adolescence to adulthood

OBJECTIVES

Apple pomace, a waste byproduct of apple processing, is rich in nutrients (e.g. polyphenols and soluble fiber) with the potential to be neuroprotective. The aim of this study was to employ RNA-sequencing (RNASeq) technology to investigate diet-gene interactions in the hypothalamus of rats after feeding a Western diet calorically substituted with apple pomace.

METHODS

Adolescent (age 21-29 days) female Sprague-Dawley rats were randomly assigned (n = 8 rats/group) to consume either a purified standard diet, Western (WE) diet, or Western diet calorically substituted with 10% apple pomace (WE/AP) for 8 weeks. RNA-seq was performed (n = 5 rats/group) to determine global differentially expressed genes in the hypothalamus.

RESULTS

RNA-seq results comparing rats fed WE to WE/AP revealed 15 differentially expressed genes in the hypothalamus. Caloric substitution of WE diet with 10% apple pomace downregulated (q < 0.06) five genes implicated in brain aging and neurodegenerative disorders: synuclein alpha, phospholipase D family member 5, NADH dehydrogenase Fe-S protein 6, choline O-acetyltransferase, and frizzled class receptor 6.

DISCUSSION

Altered gene expression of these five genes suggests that apple pomace ameliorated synthesis of the neurotransmitter, acetylcholine, in rats fed a WE diet. Apple pomace, a rich source of antioxidant polyphenols and soluble fiber, has been shown to reverse nonalcoholic fatty liver disease (NAFLD). Diet-induced NAFLD decreases hepatic de novo synthesis of choline, a precursor to acetylcholine. Based on preclinical evidence, apple pomace has the potential to be a sustainable functional food for maintaining brain function and for reducing the risk of neurodegeneration.

High-throughput phenotyping of infection by diverse microsporidia species reveals a wild C. elegans strain with opposing resistance and susceptibility traits

Animals are under constant selective pressure from a myriad of diverse pathogens. Microsporidia are ubiquitous animal parasites, but the influence they exert on shaping animal genomes is mostly unknown. Using multiplexed competition assays, we measured the impact of four different species of microsporidia on 22 wild isolates of Caenorhabditis elegans. This resulted in the identification and confirmation of 13 strains with significantly altered population fitness profiles under infection conditions. One of these identified strains, JU1400, is sensitive to an epidermal-infecting species by lacking tolerance to infection. JU1400 is also resistant to an intestinal-infecting species and can specifically recognize and destroy this pathogen. Genetic mapping of JU1400 demonstrates that these two opposing phenotypes are caused by separate loci. Transcriptional analysis reveals the JU1400 sensitivity to epidermal microsporidia infection results in a response pattern that shares similarity to toxin-induced responses. In contrast, we do not observe JU1400 intestinal resistance being regulated at the transcriptional level. The transcriptional response to these four microsporidia species is conserved, with C. elegans strain-specific differences in potential immune genes. Together, our results show that phenotypic differences to microsporidia infection amongst C. elegans are common and that animals can evolve species-specific genetic interactions.

Renal control of life-threatening malarial anemia

Iron recycling prevents the development of anemia under homeostatic conditions. Whether iron recycling was co-opted as a defense strategy to prevent the development of anemia in response to infection is unclear. We find that in severe Plasmodium falciparum malaria, the onset of life-threatening anemia is associated with acute kidney injury (AKI), irrespective of parasite load. Using a well-established experimental rodent model of malaria anemia, we identify a transcriptional response that endows renal proximal tubule epithelial cells (RPTECs) with the capacity to store and recycle iron during P. chabaudi chabaudi (Pcc) infection. This response encompasses the induction of ferroportin 1/SLC40A1, which exports iron from RPTECs and counteracts AKI while supporting compensatory erythropoiesis and preventing the onset of life-threatening malarial anemia. Iron recycling by myeloid cells is dispensable to this protective response, suggesting that RPTECs provide an iron-recycling salvage pathway that prevents the pathogenesis of life-threatening malarial anemia.

Endothelial Differentiation of CCM1 Knockout iPSCs Triggers the Establishment of a Specific Gene Expression Signature

Cerebral cavernous malformation (CCM) is a neurovascular disease that can lead to seizures and stroke-like symptoms. The familial form is caused by a heterozygous germline mutation in either the CCM1, CCM2, or CCM3 gene. While the importance of a second-hit mechanism in CCM development is well established, it is still unclear whether it immediately triggers CCM development or whether additional external factors are required. We here used RNA sequencing to study differential gene expression in CCM1 knockout induced pluripotent stem cells (CCM1^-/- iPSCs), early mesoderm progenitor cells (eMPCs), and endothelial-like cells (ECs). Notably, CRISPR/Cas9-mediated inactivation of CCM1 led to hardly any gene expression differences in iPSCs and eMPCs. However, after differentiation into ECs, we found the significant deregulation of signaling pathways well known to be involved in CCM pathogenesis. These data suggest that a microenvironment of proangiogenic cytokines and growth factors can trigger the establishment of a characteristic gene expression signature upon CCM1 inactivation. Consequently, CCM1^-/- precursor cells may exist that remain silent until entering the endothelial lineage. Collectively, not only downstream consequences of CCM1 ablation but also supporting factors must be addressed in CCM therapy development.

Transcription factors underlying photoreceptor diversity

During development, retinal progenitors navigate a complex landscape of fate decisions to generate the major cell classes necessary for proper vision. Transcriptional regulation is critical to generate diversity within these major cell classes. Here, we aim to provide the resources and techniques required to identify transcription factors necessary to generate and maintain diversity in photoreceptor subtypes, which are critical for vision. First, we generate a key resource: a high-quality and deep transcriptomic profile of each photoreceptor subtype in adult zebrafish. We make this resource openly accessible, easy to explore, and have integrated it with other currently available photoreceptor transcriptomic datasets. Second, using our transcriptomic profiles, we derive an in-depth map of expression of transcription factors in photoreceptors. Third, we use efficient CRISPR-Cas9 based mutagenesis to screen for null phenotypes in F0 larvae (F0 screening) as a fast, efficient, and versatile technique to assess the involvement of candidate transcription factors in the generation of photoreceptor subtypes. We first show that known phenotypes can be easily replicated using this method: loss of S cones in foxq2 mutants and loss of rods in nr2e3 mutants. We then identify novel functions for the transcription factor Tbx2, demonstrating that it plays distinct roles in controlling the generation of all photoreceptor subtypes within the retina. Our study provides a roadmap to discover additional factors involved in this process. Additionally, we explore four transcription factors of unknown function (Skor1a, Sall1a, Lrrfip1a, and Xbp1), and find no evidence for their involvement in the generation of photoreceptor subtypes. This dataset and screening method will be a valuable way to explore the genes involved in many other essential aspects of photoreceptor biology.

Dosage differences in 12-OXOPHYTODIENOATE REDUCTASE genes modulate wheat root growth

Wheat, an essential crop for global food security, is well adapted to a wide variety of soils. However, the gene networks shaping different root architectures remain poorly understood. We report here that dosage differences in a cluster of monocot-specific 12-OXOPHYTODIENOATE REDUCTASE genes from subfamily III (OPRIII) modulate key differences in wheat root architecture, which are associated with grain yield under water-limited conditions. Wheat plants with loss-of-function mutations in OPRIII show longer seminal roots, whereas increased OPRIII dosage or transgenic over-expression result in reduced seminal root growth, precocious development of lateral roots and increased jasmonic acid (JA and JA-Ile). Pharmacological inhibition of JA-biosynthesis abolishes root length differences, consistent with a JA-mediated mechanism. Transcriptome analyses of transgenic and wild-type lines show significant enriched JA-biosynthetic and reactive oxygen species (ROS) pathways, which parallel changes in ROS distribution. OPRIII genes provide a useful entry point to engineer root architecture in wheat and other cereals.

Vitamin D3 and deconvoluting a rash

BACKGROUNDAdverse drug reactions are unpredictable immunologic events presenting frequent challenges to clinical management. Systemically administered cholecalciferol (vitamin D3) has immunomodulatory properties. In this randomized, double-blinded, placebo-controlled interventional trial of healthy human adults, we investigated the clinical and molecular immunomodulatory effects of a single high dose of oral vitamin D3 on an experimentally induced chemical rash.METHODSSkin inflammation was induced with topical nitrogen mustard (NM) in 28 participants. Participant-specific inflammatory responses to NM alone were characterized using clinical measures, serum studies, and skin tissue analysis over the next week. All participants underwent repeat NM exposure to the opposite arm and then received placebo or 200,000 IU cholecalciferol intervention. The complete rash reaction was followed by multi-omic analysis, clinical measures, and serum studies over 6 weeks.RESULTSCholecalciferol mitigated acute inflammation in all participants and achieved 6 weeks of durable responses. Integrative analysis of skin and blood identified an unexpected divergence in response severity to NM, corroborated by systemic neutrophilia and significant histopathologic and clinical differences. Multi-omic and pathway analyses revealed a 3-biomarker signature (CCL20, CCL2, CXCL8) unique to exaggerated responders that is suppressed by cholecalciferol and implicates IL-17 signaling involvement.CONCLUSIONHigh-dose systemic cholecalciferol may be an effective treatment for severe reactions to topical chemotherapy. Our findings have broad implications for cholecalciferol as an antiinflammatory intervention against the development of exaggerated immune responses.TRIAL REGISTRATIONclinicaltrials.gov (NCT02968446).FUNDINGNIH and National Institute of Arthritis and Musculoskeletal and Skin Diseases (NIAMS; grants U01AR064144, U01AR071168, P30 AR075049, U54 AR079795, and P30 AR039750 (CWRU)).

Characterization of 3D organotypic epithelial tissues reveals tonsil-specific differences in tonic interferon signaling

Three-dimensional (3D) culturing techniques can recapitulate the stratified nature of multicellular epithelial tissues. Organotypic 3D epithelial tissue culture methods have several applications, including the study of tissue development and function, drug discovery and toxicity testing, host-pathogen interactions, and the development of tissue-engineered constructs for use in regenerative medicine. We grew 3D organotypic epithelial tissues from foreskin, cervix, and tonsil-derived primary cells and characterized the transcriptome of these in vitro tissue equivalents. Using the same 3D culturing method, all three tissues yielded stratified squamous epithelium, validated histologically using basal and superficial epithelial cell markers. The goal of this study was to use RNA-seq to compare gene expression patterns in these three types of epithelial tissues to gain a better understanding of the molecular mechanisms underlying their function and identify potential therapeutic targets for various diseases. Functional profiling by over-representation and gene set enrichment analysis revealed tissue-specific differences: i.e. , cutaneous homeostasis and lipid metabolism in foreskin, extracellular matrix remodeling in cervix, and baseline innate immune differences in tonsil. Specifically, tonsillar epithelia may play an active role in shaping the immune microenvironment of the tonsil balancing inflammation and immune responses in the face of constant exposure to microbial insults. Overall, these data serve as a resource, with gene sets made available for the research community to explore, and as a foundation for understanding the epithelial heterogeneity and how it may impact their in vitro use. An online resource is available to investigate these data ( https://viz.datascience.arizona.edu/3DEpiEx/ ).

Establishing Physalis as a Solanaceae model system enables genetic reevaluation of the inflated calyx syndrome

The highly diverse Solanaceae family contains several widely studied models and crop species. Fully exploring, appreciating, and exploiting this diversity requires additional model systems. Particularly promising are orphan fruit crops in the genus Physalis, which occupy a key evolutionary position in the Solanaceae and capture understudied variation in traits such as inflorescence complexity, fruit ripening and metabolites, disease and insect resistance, self-compatibility, and most notable, the striking inflated calyx syndrome (ICS), an evolutionary novelty found across angiosperms where sepals grow exceptionally large to encapsulate fruits in a protective husk. We recently developed transformation and genome editing in Physalis grisea (groundcherry). However, to systematically explore and unlock the potential of this and related Physalis as genetic systems, high-quality genome assemblies are needed. Here, we present chromosome-scale references for P. grisea and its close relative Physalis pruinosa and use these resources to study natural and engineered variations in floral traits. We first rapidly identified a natural structural variant in a bHLH gene that causes petal color variation. Further, and against expectations, we found that CRISPR-Cas9-targeted mutagenesis of 11 MADS-box genes, including purported essential regulators of ICS, had no effect on inflation. In a forward genetics screen, we identified huskless, which lacks ICS due to mutation of an AP2-like gene that causes sepals and petals to merge into a single whorl of mixed identity. These resources and findings elevate Physalis to a new Solanaceae model system and establish a paradigm in the search for factors driving ICS.

Creation of Philadelphia chromosome by CRISPR/Cas9-mediated double cleavages on BCR and ABL1 genes as a model for initial event in leukemogenesis

The Philadelphia (Ph) chromosome was the first translocation identified in leukemia. It is supposed to be generated by aberrant ligation between two DNA double-strand breaks (DSBs) at the BCR gene located on chromosome 9q34 and the ABL1 gene located on chromosome 22q11. Thus, mimicking the initiation process of translocation, we induced CRISPR/Cas9-mediated DSBs simultaneously at the breakpoints of the BCR and ABL1 genes in a granulocyte-macrophage colony-stimulating factor (GM-CSF) dependent human leukemia cell line. After transfection of two single guide RNAs (sgRNAs) targeting intron 13 of the BCR gene and intron 1 of the ABL1 gene, a factor-independent subline was obtained. In the subline, p210 BCR::ABL1 and its reciprocal ABL1::BCR fusions were generated as a result of balanced translocation corresponding to the Ph chromosome. Another set of sgRNAs targeting intron 1 of the BCR gene and intron 1 of the ABL1 gene induced a factor-independent subline expressing p190 BCR::ABL1. Both p210 and p190 BCR::ABL1 induced factor-independent growth by constitutively activating intracellular signaling pathways for transcriptional regulation of cell cycle progression and cell survival that are usually regulated by GM-CSF. These observations suggested that simultaneous DSBs at the BCR and ABL1 gene breakpoints are initiation events for oncogenesis in Ph+ leukemia. (200/200 words).

Conservation of Nematocida microsporidia gene expression and host response in Caenorhabditis nematodes

Microsporidia are obligate intracellular parasites that are known to infect most types of animals. Many species of microsporidia can infect multiple related hosts, but it is not known if microsporidia express different genes depending upon which host species is infected or if the host response to infection is specific to each microsporidia species. To address these questions, we took advantage of two species of Nematocida microsporidia, N. parisii and N. ausubeli, that infect two species of Caenorhabditis nematodes, C. elegans and C. briggsae. We performed RNA-seq at several time points for each host infected with either microsporidia species. We observed that Nematocida transcription was largely independent of its host. We also observed that the host transcriptional response was similar when infected with either microsporidia species. Finally, we analyzed if the host response to microsporidia infection was conserved across host species. We observed that although many of the genes upregulated in response to infection are not direct orthologs, the same expanded gene families are upregulated in both Caenorhabditis hosts. Together our results describe the transcriptional interactions of Nematocida infection in Caenorhabditis hosts and demonstrate that these responses are evolutionarily conserved.

TRPV4 Activation Increases the Expression of CD207 (Langerin) of Monocyte-Derived Langerhans Cells without Affecting their Maturation

Langerhans cells (LCs) are the sole professional antigen-presenting cell normally found in the human epidermal compartment. Research into their physiological role is hindered by the fact that they are invariably activated during isolation from the skin. To overcome this challenge, we turned to a monocyte-derived LC (moLC) model, which we characterized with RNA sequencing, and compared the transcriptome of moLCs with that of donor-matched immature dendritic cells. We found that moLCs express markers characteristic of LC2 cells as well as TRPV4. TRPV4 is especially important in the skin because it has been linked to the conservation of the skin barrier, immunological responses, as well as acute and chronic itch, but we know little about its function on LCs. Our results show that TRPV4 activation increased the expression of Langerin and led to increased intracellular calcium concentration in moLCs. Regarding the functionality of moLCs, we found that TRPV4 agonism had a mitigating effect on their inflammatory responses because it decreased their cytokine production and T-cell activating capability. Because TRPV4 has emerged as a potential therapeutic target in dermatological conditions, it is important to highlight LCs as, to our knowledge, a previously unreported target of these therapies.

Magnetique: an interactive web application to explore transcriptome signatures of heart failure

Background

Despite a recent increase in the number of RNA-seq datasets investigating heart failure (HF), accessibility and usability remain critical issues for medical researchers. We address the need for an intuitive and interactive web application to explore the transcriptional signatures of heart failure with this work.

Methods

We reanalysed the Myocardial Applied Genomics Network RNA-seq dataset, one of the largest publicly available datasets of left ventricular RNA-seq samples from patients with dilated (DCM) or hypertrophic (HCM) cardiomyopathy, as well as unmatched non-failing hearts (NFD) from organ donors and patient characteristics that allowed us to model confounding factors. We analyse differential gene expression, associated pathway signatures and reconstruct signaling networks based on inferred transcription factor activities through integer linear programming. We additionally focus, for the first time, on differential RNA transcript isoform usage (DTU) changes and predict RNA-binding protein (RBP) to target transcript interactions using a Global test approach. We report results for all pairwise comparisons (DCM, HCM, NFD).

Results

Focusing on the DCM versus HCM contrast (DCMvsHCM), we identified 201 differentially expressed genes, some of which can be clearly associated with changes in ERK1 and ERK2 signaling. Interestingly, the signs of the predicted activity for these two kinases have been inferred to be opposite to each other: In the DCMvsHCM contrast, we predict ERK1 to be consistently less activated in DCM while ERK2 was more activated in DCM. In the DCMvsHCM contrast, we identified 149 differently used transcripts. One of the top candidates is the O-linked N-acetylglucosamine (GlcNAc) transferase (OGT), which catalyzes a common post-translational modification known for its role in heart arrhythmias and heart hypertrophy. Moreover, we reconstruct RBP – target interaction networks and showcase the examples of CPEB1, which is differentially expressed in the DCMvsHCM contrast.

Conclusion

Magnetique (https://shiny.dieterichlab.org/app/magnetique) is the first online application to provide an interactive view of the HF transcriptome at the RNA isoform level and to include transcription factor signaling and RBP:RNA interaction networks. The source code for both the analyses (https://github.com/dieterich-lab/magnetiqueCode2022) and the web application (https://github.com/AnnekathrinSilvia/magnetique) is available to the public. We hope that our application will help users to uncover the molecular basis of heart failure.

ScRNA-seq expression of IFI27 and APOC2 identifies four alveolar macrophage superclusters in healthy BALF

Alveolar macrophages (AMs) reside on the luminal surface of the airways and alveoli, ensuring proper gas exchange by ingesting cellular debris and pathogens, and regulating inflammatory responses. Therefore, understanding the heterogeneity and diverse roles played by AMs, interstitial macrophages, and recruited monocytes is critical for treating airway diseases. We performed single-cell RNA sequencing on 113,213 bronchoalveolar lavage cells from four healthy and three uninflamed cystic fibrosis subjects and identified two MARCKS+LGMN+IMs, FOLR2+SELENOP+ and SPP1+PLA2G7+ IMs, monocyte subtypes, DC1, DC2, migDCs, plasmacytoid DCs, lymphocytes, epithelial cells, and four AM superclusters (families) based on the gene expression of IFI27 and APOC2 These four AM families have at least eight distinct functional members (subclusters) named after their differentially expressed gene(s): IGF1, CCL18, CXCL5, cholesterol, chemokine, metallothionein, interferon, and small-cluster AMs. Interestingly, the chemokine cluster further divides with each subcluster selectively expressing a unique combination of chemokines. One of the most striking observations, besides the heterogeneity, is the conservation of AM family members in relatively equal ratio across all AM superclusters and individuals. Transcriptional data and TotalSeq technology were used to investigate cell surface markers that distinguish resident AMs from recruited monocytes. Last, other AM datasets were projected onto our dataset. Similar AM superclusters and functional subclusters were observed, along with a significant increase in chemokine and IFN AM subclusters in individuals infected with COVID-19. Overall, functional specializations of the AM subclusters suggest that there are highly regulated AM niches with defined programming states, highlighting a clear division of labor.

Rapid and robust directed differentiation of mouse epiblast stem cells into definitive endoderm and forebrain organoids

Directed differentiation of pluripotent stem cells (PSCs) is a powerful model system for deconstructing embryonic development. Although mice are the most advanced mammalian model system for genetic studies of embryonic development, state-of-the-art protocols for directed differentiation of mouse PSCs into defined lineages require additional steps and generates target cell types with lower purity than analogous protocols for human PSCs, limiting their application as models for mechanistic studies of development. Here, we examine the potential of mouse epiblast stem cells cultured in media containing Wnt pathway inhibitors as a starting point for directed differentiation. As a proof of concept, we focused our efforts on two specific cell/tissue types that have proven difficult to generate efficiently and reproducibly from mouse embryonic stem cells: definitive endoderm and neural organoids. We present new protocols for rapid generation of nearly pure definitive endoderm and forebrain-patterned neural organoids that model the development of prethalamic and hippocampal neurons. These differentiation models present new possibilities for combining mouse genetic tools with in vitro differentiation to characterize molecular and cellular mechanisms of embryonic development.

Dithranol as novel co-adjuvant for non-invasive dermal vaccination

Transcutaneous immunization (TCI) utilizing the TLR7 agonist imiquimod (IMQ-TCI) induces T cell-driven protective immunity upon application onto intact skin. In our present work, we combine the anti-psoriatic agent dithranol with IMQ-TCI to boost vaccination efficacy (Dithranol/IMQ-based transcutaneous vaccination (DIVA)). Using ovalbumin-derived peptides as model antigens in mice, DIVA induced superior cytolytic CD8⁺ T cells and CD4⁺ T cells with a T_H1 cytokine profile in the priming as well as in the memory phase. Regarding the underlying mechanisms, dithranol induced an oxidant-dependent, monocyte-attracting inflammatory milieu in the skin boosting TLR7-dependent activation of dendritic cells and macrophages leading to superior T cell priming and protective immunity in vaccinia virus infection. In conclusion, we introduce the non-invasive vaccination method DIVA to induce strong primary and memory T cell responses upon a single local treatment. This work provides relevant insights in cutaneous vaccination approaches, paving the way for clinical development in humans.

easyMF: A Web Platform for Matrix Factorization-Based Gene Discovery from Large-scale Transcriptome Data

With the development of high-throughput experimental technologies, large-scale RNA sequencing (RNA-Seq) data have been and continue to be produced, but have led to challenges in extracting relevant biological knowledge hidden in the produced high-dimensional gene expression matrices. Here, we develop easyMF ( https://github.com/cma2015/easyMF ), a web platform that can facilitate functional gene discovery from large-scale transcriptome data using matrix factorization (MF) algorithms. Compared with existing MF-based software packages, easyMF exhibits several promising features, such as greater functionality, flexibility and ease of use. The easyMF platform is equipped using the Big-Data-supported Galaxy system with user-friendly graphic user interfaces, allowing users with little programming experience to streamline transcriptome analysis from raw reads to gene expression, carry out multiple-scenario MF analysis, and perform multiple-way MF-based gene discovery. easyMF is also powered with the advanced packing technology to enhance ease of use under different operating systems and computational environments. We illustrated the application of easyMF for seed gene discovery from temporal, spatial, and integrated RNA-Seq datasets of maize (Zea mays L.), resulting in the identification of 3,167 seed stage-specific, 1,849 seed compartment-specific, and 774 seed-specific genes, respectively. The present results also indicated that easyMF can prioritize seed-related genes with superior prediction performance over the state-of-art network-based gene prioritization system MaizeNet. As a modular, containerized and open-source platform, easyMF can be further customized to satisfy users' specific demands of functional gene discovery and deployed as a web service for broad applications.

Specialist root herbivore modulates plant transcriptome and downregulates defensive secondary metabolites in a brassicaceous plant

Plants face attackers aboveground and belowground. Insect root herbivores can lead to severe crop losses, yet the underlying transcriptomic responses have rarely been studied. We studied the dynamics of the transcriptomic response of Brussels sprouts (Brassica oleracea var. gemmifera) primary roots to feeding damage by cabbage root fly larvae (Delia radicum), alone or in combination with aboveground herbivory by cabbage aphids (Brevicoryne brassicae) or diamondback moth caterpillars (Plutella xylostella). This was supplemented with analyses of phytohormones and the main classes of secondary metabolites; aromatic, indole and aliphatic glucosinolates. Root herbivory leads to major transcriptomic rearrangement that is modulated by aboveground feeding caterpillars, but not aphids, through priming soon after root feeding starts. The root herbivore downregulates aliphatic glucosinolates. Knocking out aliphatic glucosinolate biosynthesis with CRISPR-Cas9 results in enhanced performance of the specialist root herbivore, indicating that the herbivore downregulates an effective defence. This study advances our understanding of how plants cope with root herbivory and highlights several novel aspects of insect-plant interactions for future research. Further, our findings may help breeders develop a sustainable solution to a devastating root pest.

Innate immune signaling in trophoblast and decidua organoids defines differential antiviral defenses at the maternal-fetal interface

Infections at the maternal-fetal interface can directly harm the fetus and induce complications that adversely impact pregnancy outcomes. Innate immune signaling by both fetal-derived placental trophoblasts and the maternal decidua must provide antimicrobial defenses at this critical interface without compromising its integrity. Here, we developed matched trophoblast (TO) and decidua organoids (DO) from human placentas to define the relative contributions of these cells to antiviral defenses at the maternal-fetal interface. We demonstrate that TO and DO basally secrete distinct immunomodulatory factors, including the constitutive release of the antiviral type III interferon IFN-λ2 from TOs, and differentially respond to viral infections through the induction of organoid-specific factors. Finally, we define the differential susceptibility and innate immune signaling of TO and DO to human cytomegalovirus (HCMV) and develop a co-culture model of TO and DO which showed that trophoblast-derived factors protect decidual cells from HCMV infection. Our findings establish matched TO and DO as ex vivo models to study vertically transmitted infections and highlight differences in innate immune signaling by fetal-derived trophoblasts and the maternal decidua.

Priming of Arabidopsis resistance to herbivory by insect egg deposition depends on the plant's developmental stage

While traits of plant resistance to herbivory often change during ontogeny, it is unknown whether the primability of this resistance depends on the plant's developmental stage. Resistance in non-flowering Arabidopsis thaliana against Pieris brassicae larvae is known to be primable by prior egg deposition on leaves. We investigated whether this priming effect is maintained in plants at the flowering stage. Larval performance assays revealed that flowering plants' resistance to herbivory was not primable by egg deposition. Accordingly, transcriptomes of flowering plants showed almost no response to eggs. In contrast, egg deposition on non-flowering plants enhanced the expression of genes induced by subsequent larval feeding. Strikingly, flowering plants showed constitutively high expression levels of these genes. Larvae performed generally worse on flowering than on non-flowering plants, indicating that flowering plants constitutively resist herbivory. Furthermore, we determined the seed weight in regrown plants that had been exposed to eggs and larvae during the non-flowering or flowering stage. Non-flowering plants benefitted from egg priming with a smaller loss in seed yield. The seed yield of flowering plants was unaffected by the treatments, indicating tolerance towards the larvae. Our results show that the primability of anti-herbivore defences in Arabidopsis depends on the plant's developmental stage.