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

ideal: an R/Bioconductor package for interactive differential expression analysis

BACKGROUND

RNA sequencing (RNA-seq) is an ever increasingly popular tool for transcriptome profiling. A key point to make the best use of the available data is to provide software tools that are easy to use but still provide flexibility and transparency in the adopted methods. Despite the availability of many packages focused on detecting differential expression, a method to streamline this type of bioinformatics analysis in a comprehensive, accessible, and reproducible way is lacking.

RESULTS

We developed the ideal software package, which serves as a web application for interactive and reproducible RNA-seq analysis, while producing a wealth of visualizations to facilitate data interpretation. ideal is implemented in R using the Shiny framework, and is fully integrated with the existing core structures of the Bioconductor project. Users can perform the essential steps of the differential expression analysis workflow in an assisted way, and generate a broad spectrum of publication-ready outputs, including diagnostic and summary visualizations in each module, all the way down to functional analysis. ideal also offers the possibility to seamlessly generate a full HTML report for storing and sharing results together with code for reproducibility.

CONCLUSION

ideal is distributed as an R package in the Bioconductor project ( http://bioconductor.org/packages/ideal/ ), and provides a solution for performing interactive and reproducible analyses of summarized RNA-seq expression data, empowering researchers with many different profiles (life scientists, clinicians, but also experienced bioinformaticians) to make the ideal use of the data at hand.

Th17-related mammary immunity, but not a high systemic Th1 immune response is associated with protection against E. coli mastitis

Vaccination against bovine mastitis lags behind despite high demand from the dairy industry and margin for efficacy improvement. We previously compared two immunization protocols against E. coli using either only the intramuscular route or a combination of intramuscular and mammary ductal routes, also known as 'prime and pull' strategy. A homologous mammary challenge during the memory phase showed that immunization favorably modified the mastitis course, notably in locally immunized cows in comparison to intramuscular and control adjuvant-only groups. Here, we performed whole-blood profiling through RNA-seq transcriptome and plasma cytokine 15-plex analyses at time points of the E. coli mastitis that showed significant clinical and laboratory differences among the groups. Diminished production of inflammatory cytokines and increased IFNγ were detected in the blood of immunized cows, where a T lymphocyte activation profile was evidenced at 12-h post infection. Acute phase neutropenia was less severe in these cows, and pathways related to neutrophil diapedesis and monocyte activation were also present. Furthermore, three intramammary-immunized cows showing faster healing and shorter mastitis duration had gene profiles that differed from their counterparts, but without any clue for the mastitis susceptibility difference. Inasmuch, when gene expression of CD4 T cells was assessed in mammary tissue, enrichment of IL-17-associated pathways was identified in the quarters of intramammary-immunized cows not only after challenge but also in the control quarters that were not infected. These findings indicate that local immunization mobilizes protective mechanisms that rely on the settlement of type 3 immunity-related CD4 T cells prior to infection.

Manganese Depletion Leads to Multisystem Changes in the Transcriptome of the Opportunistic Pathogen Streptococcus sanguinis

Streptococcus sanguinis is a primary colonizer of teeth and is typically considered beneficial due to its antagonistic relationship with the cariogenic pathogen Streptococcus mutans. However, S. sanguinis can also act as an opportunistic pathogen should it enter the bloodstream and colonize a damaged heart valve, leading to infective endocarditis. Studies have implicated manganese acquisition as an important virulence determinant in streptococcal endocarditis. A knockout mutant lacking the primary manganese import system in S. sanguinis, SsaACB, is severely attenuated for virulence in an in vivo rabbit model. Manganese is a known cofactor for several important enzymes in S. sanguinis, including superoxide dismutase, SodA, and the aerobic ribonucleotide reductase, NrdEF. To determine the effect of manganese depletion on S. sanguinis, we performed transcriptomic analysis on a ΔssaACB mutant grown in aerobic fermentor conditions after the addition of the metal chelator EDTA. Despite the broad specificity of EDTA, analysis of cellular metal content revealed a decrease in manganese, but not in other metals, that coincided with a drop in growth rate. Subsequent supplementation with manganese, but not iron, zinc, or magnesium, restored growth in the fermentor post-EDTA. Reduced activity of Mn-dependent SodA and NrdEF likely contributed to the decreased growth rate post-EDTA, but did not appear entirely responsible. With the exception of the Dps-like peroxide resistance gene, dpr, manganese depletion did not induce stress response systems. By comparing the transcriptome of ΔssaACB cells pre- and post-EDTA, we determined that manganese deprivation led to altered expression of diverse systems. Manganese depletion also led to an apparent induction of carbon catabolite repression in a glucose-independent manner. The combined results suggest that manganese limitation produces effects in S. sanguinis that are diverse and complex, with no single protein or system appearing entirely responsible for the observed growth rate decrease. This study provides further evidence for the importance of this trace element in streptococcal biology. Future studies will focus on determining mechanisms for regulation, as the multitude of changes observed in this study indicate that multiple regulators may respond to manganese levels.

Developing a transcriptomic framework for testing testosterone-mediated handicap hypotheses

Sexually selected traits are hypothesized to be honest signals of individual quality due to the costs associated with their maintenance, development, and/or production. Testosterone, a sex steroid associated with the development and/or production of sexually selected traits, has been proposed to enforce the honesty of sexually selected traits via its immunosuppressive effects (i.e., the Immunocompetence Handicap Hypothesis) and/or by influencing an individual's exposure/susceptibility to oxidative stress (i.e., the Oxidation Handicap Hypothesis). Previous work testing these hypotheses has primarily focused on physiological measurements of immunity or oxidative stress, but little is known about the molecular pathways by which testosterone could influence immunity and/or oxidative stress pathways. To further understand the transcriptomic consequences of experimentally elevated testosterone in the context of handicap hypotheses, we used previously published RNA-seq data from studies that measured the transcriptome of individuals treated with either a testosterone-filled or an empty (i.e., control) implant. Two studies encompassing three species of bird and three tissue types fit our selection criteria and we reanalyzed the data using weighted gene co-expression network analysis. Testosterone-treated individuals exhibited signatures of immunosuppression and our results describe the molecular pathways underlying this effect. We also provide some evidence to suggest that the transcriptomic signature of immunosuppression is evolutionarily conserved between the three species. While our results provide no evidence to suggest testosterone mediates handicaps via pathways associated with oxidative stress, they do support the hypothesis that testosterone enforces the honesty of sexually-selected traits by influencing an individual's immunocompetence. Overall, this study develops a framework for testing testosterone-mediated handicap hypotheses and provides guidelines for future integrative and comparative studies focused on the proximate mechanisms mediating sexually selected traits.

Genetic compensation rather than genetic assimilation drives the evolution of plasticity in response to mild warming across latitudes in a damselfly

Global warming is causing plastic and evolutionary changes in the phenotypes of ectotherms. Yet, we have limited knowledge on how the interplay between plasticity and evolution shapes thermal responses and underlying gene expression patterns. We assessed thermal reaction norm patterns across the transcriptome and identified associated molecular pathways in northern and southern populations of the damselfly Ischnura elegans. Larvae were reared in a common garden experiment at the mean summer water temperatures experienced at the northern (20°C) and southern (24°C) latitudes. This allowed a space-for-time substitution where the current gene expression levels at 24°C in southern larvae are a proxy for the expected responses of northern larvae under gradual thermal evolution to the predicted 4°C warming. Most differentially expressed genes showed fixed differences across temperatures between latitudes, suggesting that thermal genetic adaptation will mainly evolve through changes in constitutive gene expression. Northern populations also frequently showed plastic responses in gene expression to mild warming, while southern populations were much less responsive to temperature. Thermal responsive genes in northern populations showed to a large extent a pattern of genetic compensation, namely gene expression that was induced at 24°C in northern populations remained at a lower constant level in southern populations, and were associated with metabolic and translation pathways. There was instead little evidence for genetic assimilation of an initial plastic response to mild warming. Our data therefore suggest that genetic compensation rather than genetic assimilation may drive the evolution of plasticity in response to mild warming in this damselfly species.

Interleukin-1 promotes autoimmune neuroinflammation by suppressing endothelial heme oxygenase-1 at the blood-brain barrier

The proinflammatory cytokine interleukin 1 (IL-1) is crucially involved in the pathogenesis of multiple sclerosis (MS) and its animal model experimental autoimmune encephalomyelitis (EAE). Herein, we studied the role of IL-1 signaling in blood-brain barrier (BBB) endothelial cells (ECs), astrocytes and microglia for EAE development, using mice with the conditional deletion of its signaling receptor IL-1R1. We found that IL-1 signaling in microglia and astrocytes is redundant for the development of EAE, whereas the IL-1R1 deletion in BBB-ECs markedly ameliorated disease severity. IL-1 signaling in BBB-ECs upregulated the expression of the adhesion molecules Vcam-1, Icam-1 and the chemokine receptor Darc, all of which have been previously shown to promote CNS-specific inflammation. In contrast, IL-1R1 signaling suppressed the expression of the stress-responsive heme catabolizing enzyme heme oxygenase-1 (HO-1) in BBB-ECs, promoting disease progression via a mechanism associated with deregulated expression of the IL-1-responsive genes Vcam1, Icam1 and Ackr1 (Darc). Mechanistically, our data emphasize a functional crosstalk of BBB-EC IL-1 signaling and HO-1, controlling the transcription of downstream proinflammatory genes promoting the pathogenesis of autoimmune neuroinflammation.

Transcriptomic analysis at organ and time scale reveals gene regulatory networks controlling the sulfate starvation response of Solanum lycopersicum

BACKGROUND

Sulfur is a major component of biological molecules and thus an essential element for plants. Deficiency of sulfate, the main source of sulfur in soils, negatively influences plant growth and crop yield. The effect of sulfate deficiency on plants has been well characterized at the physiological, transcriptomic and metabolomic levels in Arabidopsis thaliana and a limited number of crop plants. However, we still lack a thorough understanding of the molecular mechanisms and regulatory networks underlying sulfate deficiency in most plants. In this work we analyzed the impact of sulfate starvation on the transcriptome of tomato plants to identify regulatory networks and key transcriptional regulators at a temporal and organ scale.

RESULTS

Sulfate starvation reduces the growth of roots and leaves which is accompanied by major changes in the organ transcriptome, with the response being temporally earlier in roots than leaves. Comparative analysis showed that a major part of the Arabidopsis and tomato transcriptomic response to sulfate starvation is conserved between these plants and allowed for the identification of processes specifically regulated in tomato at the transcript level, including the control of internal phosphate levels. Integrative gene network analysis uncovered key transcription factors controlling the temporal expression of genes involved in sulfate assimilation, as well as cell cycle, cell division and photosynthesis during sulfate starvation in tomato roots and leaves. Interestingly, one of these transcription factors presents a high identity with SULFUR LIMITATION1, a central component of the sulfate starvation response in Arabidopsis.

CONCLUSIONS

Together, our results provide the first comprehensive catalog of sulfate-responsive genes in tomato, as well as novel regulatory targets for future functional analyses in tomato and other crops.

Bacterial polyphosphates interfere with the innate host defense to infection

Polyphosphates are linear polymers and ubiquitous metabolites. Bacterial polyphosphates are long chains of hundreds of phosphate units. Here, we report that mouse survival of peritoneal Escherichia coli sepsis is compromised by long-chain polyphosphates, and improves with bacterial polyphosphatekinase deficiency or neutralization using recombinant exopolyphosphatase. Polyphosphate activities are chain-length dependent, impair pathogen clearance, antagonize phagocyte recruitment, diminish phagocytosis and decrease production of iNOS and cytokines. Macrophages bind and internalize polyphosphates, in which their effects are independent of P2Y1 and RAGE receptors. The M1 polarization driven by E. coli derived LPS is misdirected by polyphosphates in favor of an M2 resembling phenotype. Long-chain polyphosphates modulate the expression of more than 1800 LPS/TLR4-regulated genes in macrophages. This interference includes suppression of hundreds of type I interferon-regulated genes due to lower interferon production and responsiveness, blunted STAT1 phosphorylation and reduced MHCII expression. In conclusion, prokaryotic polyphosphates disturb multiple macrophage functions for evading host immunity.

Systematic analysis of 1298 RNA-Seq samples and construction of a comprehensive soybean (Glycine max) expression atlas

Soybean (Glycine max [L.] Merr.) is a major crop in animal feed and human nutrition, mainly for its rich protein and oil contents. The remarkable rise in soybean transcriptome studies over the past 5 years generated an enormous amount of RNA-seq data, encompassing various tissues, developmental conditions and genotypes. In this study, we have collected data from 1298 publicly available soybean transcriptome samples, processed the raw sequencing reads and mapped them to the soybean reference genome in a systematic fashion. We found that 94% of the annotated genes (52 737/56 044) had detectable expression in at least one sample. Unsupervised clustering revealed three major groups, comprising samples from aerial, underground and seed/seed-related parts. We found 452 genes with uniform and constant expression levels, supporting their roles as housekeeping genes. On the other hand, 1349 genes showed heavily biased expression patterns towards particular tissues. A transcript-level analysis revealed that 95% (70 963 of 74 490) of the assembled transcripts have intron chains exactly matching those from known transcripts, whereas 3256 assembled transcripts represent potentially novel splicing isoforms. The dataset compiled here constitute a new resource for the community, which can be downloaded or accessed through a user-friendly web interface at http://venanciogroup.uenf.br/resources/. This comprehensive transcriptome atlas will likely accelerate research on soybean genetics and genomics.

Organ of Corti size is governed by Yap/Tead-mediated progenitor self-renewal

Precise control of organ growth and patterning is executed through a balanced regulation of progenitor self-renewal and differentiation. In the auditory sensory epithelium-the organ of Corti-progenitor cells exit the cell cycle in a coordinated wave between E12.5 and E14.5 before the initiation of sensory receptor cell differentiation, making it a unique system for studying the molecular mechanisms controlling the switch between proliferation and differentiation. Here we identify the Yap/Tead complex as a key regulator of the self-renewal gene network in organ of Corti progenitor cells. We show that Tead transcription factors bind directly to the putative regulatory elements of many stemness- and cell cycle-related genes. We also show that the Tead coactivator protein, Yap, is degraded specifically in the Sox2-positive domain of the cochlear duct, resulting in down-regulation of Tead gene targets. Further, conditional loss of the Yap gene in the inner ear results in the formation of significantly smaller auditory and vestibular sensory epithelia, while conditional overexpression of a constitutively active version of Yap, Yap5SA, is sufficient to prevent cell cycle exit and to prolong sensory tissue growth. We also show that viral gene delivery of Yap5SA in the postnatal inner ear sensory epithelia in vivo drives cell cycle reentry after hair cell loss. Taken together, these data highlight the key role of the Yap/Tead transcription factor complex in maintaining inner ear progenitors during development, and suggest new strategies to induce sensory cell regeneration.

Cell-Specific Transcriptional Responses to Heat Shock in the Mouse Utricle Epithelium

Sensory epithelia of the inner ear contain mechanosensory hair cells (HCs) and glia-like supporting cells (SCs), both of which are required for hearing and balance functions. Each of these cell types has unique responses to ototoxic and cytoprotective stimuli. Non-lethal heat stress in the mammalian utricle induces heat shock proteins (HSPs) and protects against ototoxic drug-induced hair cell death. Induction of HSPs in the utricle demonstrates cell-type specificity at the protein level, with HSP70 induction occurring primarily in SCs, while HSP32 (also known as heme oxygenase 1, HMOX1) is induced primarily in resident macrophages. Neither of these HSPs are robustly induced in HCs, suggesting that HCs may have little capacity for induction of stress-induced protective responses. To determine the transcriptional responses to heat shock of these different cell types, we performed cell-type-specific transcriptional profiling using the RiboTag method, which allows for immunoprecipitation (IP) of actively translating mRNAs from specific cell types. RNA-Seq differential gene expression analyses demonstrated that the RiboTag method identified known cell type-specific markers as well as new markers for HCs and SCs. Gene expression differences suggest that HCs and SCs exhibit differential transcriptional heat shock responses. The chaperonin family member Cct8 was significantly enriched only in heat-shocked HCs, while Hspa1l (HSP70 family), and Hspb1 and Cryab (HSP27 and HSP20 families, respectively) were enriched only in SCs. Together our data indicate that HCs exhibit a limited but unique heat shock response, and SCs exhibit a broader and more robust transcriptional response to protective heat stress.

Circulating RNA Transcriptome of Pregnant Women with TSH Just Above the Trimester-Specific Reference and its Correlation with the Hypertensive Phenotype

During gestation, a woman’s body undergoes physiological changes that alter thyroid function. Pregnant women with hypothyroidism may exhibit gestational complications, including hypertension and preeclampsia. We investigated differentially expressed genes (DEGs) in circulating RNAs from pregnant women with TSH levels just above the normal range to determine the impact of a mild elevation of TSH in pregnancy. We selected three women with healthy thyroid pregnancy (HTP), three pregnant women with gestational hypothyroidism (GHT), and three nonpregnant women (NPG) to construct transcriptome libraries. We also compared our results with data from the GEO dataset and DisGeNET. We identified 1500 DEG in GHT and 1656 DEG in HTP. From GEO dataset, we recognized 453 DEGs in trimester-specific plasma RNA, 1263 DEGs in placental tissues from healthy women, 1031 DEGs from preeclamptic uteroplacental tissues and 1657 DEGs from placental tissues from severely preeclamptic women. In this scenario, 12.26% and 12.86% genes were shared between these datasets in GHT and HTP, respectively. We stablished 62 genes in GHT DEGs related to hypertensive phenotype hallmarks. In conclusion, even in women with a mild TSH increment, we were able to detect some DEGs that could be associated with a hypertensive phenotype.

Tissue Determinants of Human NK Cell Development, Function, and Residence

Immune responses in diverse tissue sites are critical for protective immunity and homeostasis. Here, we investigate how tissue localization regulates the development and function of human natural killer (NK) cells, innate lymphocytes important for anti-viral and tumor immunity. Integrating high-dimensional analysis of NK cells from blood, lymphoid organs, and mucosal tissue sites from 60 individuals, we identify tissue-specific patterns of NK cell subset distribution, maturation, and function maintained across age and between individuals. Mature and terminally differentiated NK cells with enhanced effector function predominate in blood, bone marrow, spleen, and lungs and exhibit shared transcriptional programs across sites. By contrast, precursor and immature NK cells with reduced effector capacity populate lymph nodes and intestines and exhibit tissue-resident signatures and site-specific adaptations. Together, our results reveal anatomic control of NK cell development and maintenance as tissue-resident populations, whereas mature, terminally differentiated subsets mediate immunosurveillance through diverse peripheral sites. VIDEO ABSTRACT.

IL-17+ CD8+ T cell suppression by dimethyl fumarate associates with clinical response in multiple sclerosis

IL-17-producing CD8⁺ (Tc17) cells are enriched in active lesions of patients with multiple sclerosis (MS), suggesting a role in the pathogenesis of autoimmunity. Here we show that amelioration of MS by dimethyl fumarate (DMF), a mechanistically elusive drug, associates with suppression of Tc17 cells. DMF treatment results in reduced frequency of Tc17, contrary to Th17 cells, and in a decreased ratio of the regulators RORC-to-TBX21, along with a shift towards cytotoxic T lymphocyte gene expression signature in CD8⁺ T cells from MS patients. Mechanistically, DMF potentiates the PI3K-AKT-FOXO1-T-BET pathway, thereby limiting IL-17 and RORγt expression as well as STAT5-signaling in a glutathione-dependent manner. This results in chromatin remodeling at the Il17 locus. Consequently, T-BET-deficiency in mice or inhibition of PI3K-AKT, STAT5 or reactive oxygen species prevents DMF-mediated Tc17 suppression. Overall, our data disclose a DMF-AKT-T-BET driven immune modulation and suggest putative therapy targets in MS and beyond.

Dimethyl fumarate (DMF) is a therapy for multiple sclerosis (MS) with undetermined mechanism of action. Here the authors find that clinical response to DMF associates with decrease in IL-17-producing CD8⁺ T cells (Tc17), delineate molecular pathways involved, and show that DMF suppresses Tc17 pathogenicity in a mouse model of MS.

Gene expression profile of Aedes aegypti females in courtship and mating

Aedes aegypti is the most synanthropic and anthropophilic mosquito of Culicidae. This species always cohabits with humans and is extremely opportunistic. Vector dispersal is directly related to the ability of the females on successfully finding a mate in a generally patchy urban scenario. In the present work, we investigate transcriptional changes in Ae. aegypti females during the courtship process and after mating. We observe a substantial alteration in gene expression triggered just upon contact with Ae. aegypti males, which in turn was not fully correlated to the changes triggered by the contact. After analysing shared significant differentially regulated genes between conspecific contact and insemination, the major part of the observed transcriptomic change triggered by contact is reversed after mating, indicating an intermediary situation between naive and mating conditions that we hypothesize to be crucial for mating success. Upon contact, several chemosensory related genes are repressed, especially odorant binding proteins. Most of these genes return to higher expression rates after mating. None of these genes are significantly regulated by the encounter of a different species, Aedes albopictus. The results presented here might be applied to an innovative control approach focusing on the semiochemical systems of mosquitoes in an effort to disrupt undesirable host–insect interaction to reduce the risk of pathogen transmission to humans.