CirGO: an alternative circular way of visualising gene ontology terms

Decoding bee cleptoparasitism through comparative transcriptomics of Coelioxoides waltheriae and its host Tetrapedia diversipes

Cleptoparasitism, also known as brood parasitism, is a widespread strategy among bee species in which the parasite lays eggs into the nests of the host species. Even though this behavior has significant ecological implications for the dynamics of several species, little is known about the molecular pathways associated with cleptoparasitism. To shed some light on this issue, we used gene expression data to perform a comparative analysis between two solitary neotropical bees: Coelioxoides waltheriae, an obligate parasite, and their specific host Tetrapedia diversipes. We found that ortholog genes involved in signal transduction, sensory perception, learning, and memory formation were differentially expressed between the cleptoparasite and the host. We hypothesize that these genes and their associated molecular pathways are engaged in cleptoparasitism-related processes and, hence, are appealing subjects for further investigation into functional and evolutionary aspects of cleptoparasitism in bees.

Transcriptional reprogramming in sound-treated Micro-Tom plants inoculated with Pseudomonas syringae pv. tomato DC3000

Sound vibrations (SV) are known to influence molecular and physiological processes that can improve crop performance and yield. In this study, the effects of three audible frequencies (100, 500 and 1000 Hz) at constant amplitude (90 dB) on tomato Micro-Tom physiological responses were evaluated 1 and 3 days post-treatment. Moreover, the potential use of SV treatment as priming agent for improved Micro-Tom resistance to Pseudomonas syringae pv. tomato DC3000 was tested by microarray. Results showed that the SV-induced physiological changes were frequency- and time-dependent, with the largest changes registered at 1000 Hz at day 3. SV treatments tended to alter the foliar content of photosynthetic pigments, soluble proteins, sugars, phenolic composition, and the enzymatic activity of polyphenol oxidase, peroxidase, superoxide dismutase and catalase. Microarray data revealed that 1000 Hz treatment is effective in eliciting transcriptional reprogramming in tomato plants grown under normal conditions, but particularly after the infection with Pst DC3000. Broadly, in plants challenged with Pst DC3000, the 1000 Hz pretreatment provoked the up-regulation of unique differentially expressed genes (DEGs) involved in cell wall reinforcement, phenylpropanoid pathway and defensive proteins. In addition, in those plants, DEGs associated with enhancing plant basal immunity, such as proteinase inhibitors, pathogenesis-related proteins, and carbonic anhydrase 3, were notably up-regulated in comparison with non-SV pretreated, infected plants. These findings provide new insights into the modulation of Pst DC3000-tomato interaction by sound and open up prospects for further development of strategies for plant disease management through the reinforcement of defense mechanisms in Micro-Tom plants.

Detailing Early Shoot Growth Arrest in Kro-0 x BG-5 Hybrids of Arabidopsis thaliana

Shoot growth directly impacts plant productivity. Plants adjust their shoot growth in response to varying environments to maximize resource capture and stress resilience. While several factors controlling shoot growth are known, the complexity of the regulation and the input of the environment are not fully understood. We have investigated shoot growth repression induced by low ambient temperatures in hybrids of Arabidopsis thaliana Kro-0 and BG-5 accessions. To continue our previous studies, we confirmed that the Kro-0 allele of DYNAMIN-RELATED PROTEIN 3B causes stunted shoot growth in the BG-5 background. We also found that shoot growth repression was most pronounced near the apex at a lower temperature and that the cells in the hybrid stem failed to elongate correctly. Furthermore, we observed that shoot growth repression in hybrids depended on light availability. Global gene expression analysis indicated the involvement of hormones, especially strigolactone, associated with the dwarf phenotype. Altogether, this study enhances our knowledge on the genetic, physiological and environmental factors associated with shoot growth regulation.

Assaying Effector Cell-to-Cell Mobility in Plant Tissues Identifies Hypermobility and Indirect Manipulation of Plasmodesmata

In plants, plasmodesmata establish cytoplasmic continuity between cells to allow for communication and resource exchange across the cell wall. While plant pathogens use plasmodesmata as a pathway for both molecular and physical invasion, the benefits of molecular invasion (cell-to-cell movement of pathogen effectors) are poorly understood. To establish a methodology for identification and characterization of the cell-to-cell mobility of effectors, we performed a quantitative live imaging-based screen of candidate effectors of the fungal pathogen Colletotrichum higginsianum. We predicted C. higginsianum effectors by their expression profiles, the presence of a secretion signal, and their predicted and in planta localization when fused to green fluorescent protein. We assayed for cell-to-cell mobility of nucleocytosolic effectors and identified 14 that are cell-to-cell mobile. We identified that three of these effectors are "hypermobile," showing cell-to-cell mobility greater than expected for a protein of that size. To explore the mechanism of hypermobility, we chose two hypermobile effectors and measured their impact on plasmodesmata function and found that even though they show no direct association with plasmodesmata, each increases the transport capacity of plasmodesmata. Thus, our methods for quantitative analysis of cell-to-cell mobility of candidate microbe-derived effectors, or any suite of host proteins, can identify cell-to-cell hypermobility and offer greater understanding of how proteins affect plasmodesmal function and intercellular connectivity. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Quantitative subcellular reconstruction reveals a lipid mediated inter-organelle biogenesis network

The structures and functions of organelles in cells depend on each other but have not been systematically explored. We established stable knockout cell lines of peroxisomal, Golgi and endoplasmic reticulum genes identified in a whole-genome CRISPR knockout screen for inducers of mitochondrial biogenesis stress, showing that defects in peroxisome, Golgi and endoplasmic reticulum metabolism disrupt mitochondrial structure and function. Our quantitative total-organelle profiling approach for focussed ion beam scanning electron microscopy revealed in unprecedented detail that specific organelle dysfunctions precipitate multi-organelle biogenesis defects, impair mitochondrial morphology and reduce respiration. Multi-omics profiling showed a unified proteome response and global shifts in lipid and glycoprotein homeostasis that are elicited when organelle biogenesis is compromised, and that the resulting mitochondrial dysfunction can be rescued with precursors for ether-glycerophospholipid metabolic pathways. This work defines metabolic and morphological interactions between organelles and how their perturbation can cause disease.

Machine learning extracts marks of thiamine's role in cold acclimation in the transcriptome of Vitis vinifera

Introduction

The escalating challenge of climate change has underscored the critical need to understand cold defense mechanisms in cultivated grapevine Vitis vinifera. Temperature variations can affect the growth and overall health of vine.

Methods

We used Self Organizing Maps machine learning method to analyze gene expression data from leaves of five Vitis vinifera cultivars each treated by four different temperature conditions. The algorithm generated sample-specific "portraits" of the normalized gene expression data, revealing distinct patterns related to the temperature conditions applied.

Results

Our analysis unveiled a connection with vitamin B1 (thiamine) biosynthesis, suggesting a link between temperature regulation and thiamine metabolism, in agreement with thiamine related stress response established in Arabidopsis before. Furthermore, we found that epigenetic mechanisms play a crucial role in regulating the expression of stress-responsive genes at low temperatures in grapevines.

Discussion

Application of Self Organizing Maps portrayal to vine transcriptomics identified modules of coregulated genes triggered under cold stress. Our machine learning approach provides a promising option for transcriptomics studies in plants.

De-novo whole genome assembly of the orange jewelweed, Impatiens capensis Meerb. (Balsaminaceae) using nanopore long-read sequencing

The plant family Balsaminaceae comprises only two genera, and they are a study in contrasts. While Impatiens is an impressively prolific genus, with over 1,000 species and more being discovered each year, its sister genus, Hydrocera, has one solitary species, H. triflora. The two genera also differ in geographic distribution and habitat type (Impatiens species are widely distributed in much of the Old World and N. America, while H. triflora is confined to wetlands specific to S. India, Sri Lanka, and SE Asia). Other contrasting features include plant habit, habitat, floral architecture, mode of seed dispersal, and a host of other traits. The family Balsaminaceae is therefore an excellent model for studying speciation and character evolution as well as understanding the proximal and evolutionary forces that have driven the two genera to adopt such contrasting evolutionary paths. Various species of the Impatiens genus are also commercially important in the ornamental flower industry and as sources of phytochemicals that are of medicinal and other commercial value. As a preliminary step towards studying the genomic basis of the contrasting features of the two genera, we have sequenced and assembled, de novo, the genome of an iconic Impatiens species from N. America, namely I. capensis, and report our findings here.

Proteomic interrogation of the meninges reveals the molecular identities of structural components and regional distinctions along the CNS axis

The meninges surround the brain and spinal cord, affording physical protection while also serving as a niche of neuroimmune activity. Though possessing stromal qualities, its complex cellular and extracellular makeup has yet to be elaborated, and it remains unclear whether the meninges vary along the neuroaxis. Hence, studies were carried-out to elucidate the protein composition and structural organization of brain and spinal cord meninges in normal, adult Biozzi ABH mice. First, shotgun, bottom-up proteomics was carried-out. Prominent proteins at both brain and spinal levels included Type II collagen and Type II keratins, representing extracellular matrix (ECM) and cytoskeletal categories, respectively. While the vast majority of total proteins detected was shared between both meningeal locales, more were uniquely detected in brain than in spine. This pattern was also seen when total proteins were subdivided by cellular compartment, except in the case of the ECM category where brain and spinal meninges each had near equal number of unique proteins, and Type V and type III collagen registered exclusively in the spine. Quantitative analysis revealed differential expression of several collagens and cytoskeletal proteins between brain and spinal meninges. High-resolution immunofluorescence and immunogold-scanning electronmicroscopy on sections from whole brain and spinal cord - still encased within bone -identified major proteins detected by proteomics, and highlighted their association with cellular and extracellular elements of variously shaped arachnoid trabeculae. Western blotting aligned with the proteomic and immunohistological analyses, reinforcing differential appearance of proteins in brain vs spinal meninges. Results could reflect regional distinctions in meninges that govern protective and/or neuroimmune functions.

A comparative analysis of gene and protein expression in chronic and acute models of photoreceptor degeneration in adult zebrafish

Background: Adult zebrafish are capable of photoreceptor (PR) regeneration following acute phototoxic lesion (AL). We developed a chronic low light (CLL) exposure model that more accurately reflects chronic PR degeneration observed in many human retinal diseases. Methods: Here, we characterize the morphological and transcriptomic changes associated with acute and chronic models of PR degeneration at 8 time-points over a 28-day window using immunohistochemistry and 3'mRNA-seq. Results: We first observed a differential sensitivity of rod and cone PRs to CLL. Next, we found no evidence for Müller glia (MG) gliosis or regenerative cell-cycle re-entry in the CLL model, which is in contrast to the robust gliosis and proliferative response from resident MG in the AL model. Differential responses of microglia between the models was also observed. Transcriptomic comparisons between the models revealed gene-specific networks of PR regeneration and degeneration, including genes that are activated under conditions of chronic PR stress. Finally, we showed that CLL is at least partially reversible, allowing for rod and cone outer segment outgrowth and replacement of rod cell nuclei via an apparent upregulation of the existing rod neurogenesis mechanism. Discussion: Collectively, these data provide a direct comparison of the morphological and transcriptomic PR degeneration and regeneration models in zebrafish.

Benzo[a]pyrene-induced multigenerational changes in gene expression, behavior, and DNA methylation are primarily influenced by paternal exposure

Benzo[a]pyrene (BaP), a polycyclic aromatic hydrocarbon (PAH), is implicated in many developmental and behavioral adverse outcomes in offspring of exposed parents. The objective of this study was to investigate sex-dependent multigenerational effects of preconceptional effects of BaP exposure. Adult wild-type (5D) zebrafish were fed 708 μg BaP/g diet (measured) at a rate of 1% body weight twice/day (14 μg BaP/g fish/day) for 21 days. Fish were spawned using a crossover design, and parental (F0) behavior and reproductive indexes were measured. In offspring, behavioral effects were measured at 96 h post fertilization (hpf) in F1 & F2 larvae, and again when F1s were adults. Compared to controls, there was no significant effect on F0 adult behavior immediately following exposure, but locomotor activity was significantly increased in F1 adults of both sexes. Larval behavior (96 hpf, photomotor response assay) was significantly altered in both the F1 and F2 generations. To assess molecular changes associated with BaP exposure, we conducted transcriptome and DNA methylation profiling in F0 gametes (sperm and eggs) and F1 embryos (10 hpf) from all four crosses. Embryos resulting from the BaP male and control female cross had the most differentially expressed genes (DEGs) and differentially methylated regions (DMRs). Some DMRs were associated with genes encoding chromatin modifying enzymes suggesting regulation of chromatin conformation by DNA methylation. Overall, these results suggest that parental dietary BaP exposure significantly contributes to the multigenerational adverse outcomes.

Methadone alters transcriptional programs associated with synapse formation in human cortical organoids

Opioid use disorder (OUD) among pregnant women has become an epidemic in the United States. Pharmacological interventions for maternal OUD most commonly involve methadone, a synthetic opioid analgesic that attenuates withdrawal symptoms and behaviors linked with drug addiction. However, evidence of methadone's ability to readily accumulate in neural tissue, and cause long-term neurocognitive sequelae, has led to concerns regarding its effect on prenatal brain development. We utilized human cortical organoid (hCO) technology to probe how this drug impacts the earliest mechanisms of cortico-genesis. Bulk mRNA sequencing of 2-month-old hCOs chronically treated with a clinically relevant dose of 1 μM methadone for 50 days revealed a robust transcriptional response to methadone associated with functional components of the synapse, the underlying extracellular matrix (ECM), and cilia. Co-expression network and predictive protein-protein interaction analyses demonstrated that these changes occurred in concert, centered around a regulatory axis of growth factors, developmental signaling pathways, and matricellular proteins (MCPs). TGFβ1 was identified as an upstream regulator of this network and appeared as part of a highly interconnected cluster of MCPs, of which thrombospondin 1 (TSP1) was most prominently downregulated and exhibited dose-dependent reductions in protein levels. These results demonstrate that methadone exposure during early cortical development alters transcriptional programs associated with synaptogenesis, and that these changes arise by functionally modulating extra-synaptic molecular mechanisms in the ECM and cilia. Our findings provide novel insight into the molecular underpinnings of methadone's putative effect on cognitive and behavioral development and a basis for improving interventions for maternal opioid addiction.

Copy number variation in tRNA isodecoder genes impairs mammalian development and balanced translation

The number of tRNA isodecoders has increased dramatically in mammals, but the specific molecular and physiological reasons for this expansion remain elusive. To address this fundamental question we used CRISPR editing to knockout the seven-membered phenylalanine tRNA gene family in mice, both individually and combinatorially. Using ATAC-Seq, RNA-seq, ribo-profiling and proteomics we observed distinct molecular consequences of single tRNA deletions. We show that tRNA-Phe-1-1 is required for neuronal function and its loss is partially compensated by increased expression of other tRNAs but results in mistranslation. In contrast, the other tRNA-Phe isodecoder genes buffer the loss of each of the remaining six tRNA-Phe genes. In the tRNA-Phe gene family, the expression of at least six tRNA-Phe alleles is required for embryonic viability and tRNA-Phe-1-1 is most important for development and survival. Our results reveal that the multi-copy configuration of tRNA genes is required to buffer translation and viability in mammals.

Dehydration Alters Transcript Levels in the Mosquito Midgut, Likely Facilitating Rapid Rehydration following a Bloodmeal

The mosquito midgut is an important site for bloodmeal regulation while also acting as a primary site for pathogen exposure within the mosquito. Recent studies show that exposure to dehydrating conditions alters mosquito bloodfeeding behaviors as well as post-feeding regulation, likely altering how pathogens interact with the mosquito. Unfortunately, few studies have explored the underlying dynamics between dehydration and bloodmeal utilization, and the overall impact on disease transmission dynamics remains veiled. In this study, we find that dehydration-based feeding in the yellow fever mosquito, Aedes aegypti, prompts alterations to midgut gene expression, as well as subsequent physiological factors involving water control and post-bloodfeeding (pbf) regulation. Altered expression of ion transporter genes and aquaporin 2 (AQP2) in the midgut of dehydrated mosquitoes as well as the rapid reequilibration of hemolymph osmolality after a bloodmeal indicate an ability to expedite fluid and ion processing. These alterations ultimately indicate that female A. aegypti employ mechanisms to ameliorate the detriments of dehydration by imbibing a bloodmeal, providing an effective avenue for rehydration. Continued research into bloodmeal utilization and the resulting effects on arthropod-borne transmission dynamics becomes increasingly important as drought prevalence is increased by climate change.

Global transcriptome and gene co-expression network analyses reveal regulatory and non-additive effects of drought and heat stress in grapevine

Despite frequent co-occurrence of drought and heat stress, the molecular mechanisms governing plant responses to these stresses in combination have not often been studied. This is particularly evident in non-model, perennial plants. We conducted large scale physiological and transcriptome analyses to identify genes and pathways associated with grapevine response to drought and/or heat stress during stress progression and recovery. We identified gene clusters with expression correlated to leaf temperature and water stress and five hub genes for the combined stress co-expression network. Several differentially expressed genes were common to the individual and combined stresses, but the majority were unique to the individual or combined stress treatments. These included heat-shock proteins, mitogen-activated kinases, sugar metabolizing enzymes, and transcription factors, while phenylpropanoid biosynthesis and histone modifying genes were unique to the combined stress treatment. Following physiological recovery, differentially expressed genes were found only in plants under heat stress, both alone and combined with drought. Taken collectively, our results suggest that the effect of the combined stress on physiology and gene expression is more severe than that of individual stresses, but not simply additive, and that epigenetic chromatin modifications may play an important role in grapevine responses to combined drought and heat stress.

Coordinated control of neuronal differentiation and wiring by sustained transcription factors

The large diversity of cell types in nervous systems presents a challenge in identifying the genetic mechanisms that encode it. Here, we report that nearly 200 distinct neurons in the Drosophila visual system can each be defined by unique combinations of on average 10 continuously expressed transcription factors. We show that targeted modifications of this terminal selector code induce predictable conversions of neuronal fates that appear morphologically and transcriptionally complete. Cis-regulatory analysis of open chromatin links one of these genes to an upstream patterning factor that specifies neuronal fates in stem cells. Experimentally validated network models describe the synergistic regulation of downstream effectors by terminal selectors and ecdysone signaling during brain wiring. Our results provide a generalizable framework of how specific fates are implemented in postmitotic neurons.

Common Assays in Mammalian Golgi Studies

The Golgi is a complex structure characterized by stacks of tightly aligned flat cisternae. In mammalian cells, Golgi stacks often concentrate in the perinuclear region and link together to form a ribbon. This structure is dynamic to accommodate continuous cargo flow in and out of the Golgi in both directions and undergoes morphological changes under physiological and pathological conditions. The fine, stacked Golgi structure makes it difficult to study by conventional light or even super-resolution microscopy. Furthermore, efforts to understand how Golgi structural dynamics impact cellular processes have been slow because of the knowledge gap in the protein machinery that maintains the complex and dynamic Golgi structure. In this method article, we list the common assays used in our research to help new and established researchers select the most appropriate method to properly address their questions.

Early prediction and longitudinal modeling of preeclampsia from multiomics

Preeclampsia is a complex disease of pregnancy whose physiopathology remains unclear. We developed machine-learning models for early prediction of preeclampsia (first 16 weeks of pregnancy) and over gestation by analyzing six omics datasets from a longitudinal cohort of pregnant women. For early pregnancy, a prediction model using nine urine metabolites had the highest accuracy and was validated on an independent cohort (area under the receiver-operating characteristic curve [AUC] = 0.88, 95% confidence interval [CI] [0.76, 0.99] cross-validated; AUC = 0.83, 95% CI [0.62,1] validated). Univariate analysis demonstrated statistical significance of identified metabolites. An integrated multiomics model further improved accuracy (AUC = 0.94). Several biological pathways were identified including tryptophan, caffeine, and arachidonic acid metabolisms. Integration with immune cytometry data suggested novel associations between immune and proteomic dynamics. While further validation in a larger population is necessary, these encouraging results can serve as a basis for a simple, early diagnostic test for preeclampsia.

Phylogenomic diversity of Vibrio species and other Gammaproteobacteria isolated from Pacific oysters (Crassostrea gigas) during a summer mortality outbreak

The Pacific oyster (PO), Crassostrea gigas, is an important commercial marine species but periodically experiences large stock losses due to disease events known as summer mortality. Summer mortality has been linked to environmental perturbations and numerous viral and bacterial agents, indicating this disease is multifactorial in nature. In 2013 and 2014, several summer mortality events occurred within the Port Stephens estuary (NSW, Australia). Extensive culture and molecular-based investigations were undertaken and several potentially pathogenic Vibrio species were identified. To improve species identification and genomically characterise isolates obtained from this outbreak, whole-genome sequencing (WGS) and subsequent genomic analyses were performed on 48 bacterial isolates, as well as a further nine isolates from other summer mortality studies using the same batch of juveniles. Average nucleotide identity (ANI) identified most isolates to the species level and included members of the Photobacterium, Pseudoalteromonas, Shewanella and Vibrio genera, with Vibrio species making up more than two-thirds of all species identified. Construction of a phylogenomic tree, ANI analysis, and pan-genome analysis of the 57 isolates represents the most comprehensive culture-based phylogenomic survey of Vibrios during a PO summer mortality event in Australian waters and revealed large genomic diversity in many of the identified species. Our analysis revealed limited and inconsistent associations between isolate species and their geographical origins, or host health status. Together with ANI and pan-genome results, these inconsistencies suggest that to determine the role that microbes may have in Pacific oyster summer mortality events, isolate identification must be at the taxonomic level of strain. Our WGS data (specifically, the accessory genomes) differentiated bacterial strains, and coupled with associated metadata, highlight the possibility of predicting a strain's environmental niche and level of pathogenicity.

Chromosome-scale genome assemblies and annotations for Poales species Carex cristatella, Carex scoparia, Juncus effusus, and Juncus inflexus

The majority of sequenced genomes in the monocots are from species belonging to Poaceae, which include many commercially important crops. Here, we expand the number of sequenced genomes from the monocots to include the genomes of 4 related cyperids: Carex cristatella and Carex scoparia from Cyperaceae and Juncus effusus and Juncus inflexus from Juncaceae. The high-quality, chromosome-scale genome sequences from these 4 cyperids were assembled by combining whole-genome shotgun sequencing of Nanopore long reads, Illumina short reads, and Hi-C sequencing data. Some members of the Cyperaceae and Juncaceae are known to possess holocentric chromosomes. We examined the repeat landscapes in our sequenced genomes to search for potential repeats associated with centromeres. Several large satellite repeat families, comprising 3.2–9.5% of our sequenced genomes, showed dispersed distribution of large satellite repeat clusters across all Carex chromosomes, with few instances of these repeats clustering in the same chromosomal regions. In contrast, most large Juncus satellite repeats were clustered in a single location on each chromosome, with sporadic instances of large satellite repeats throughout the Juncus genomes. Recognizable transposable elements account for about 20% of each of the 4 genome assemblies, with the Carex genomes containing more DNA transposons than retrotransposons while the converse is true for the Juncus genomes. These genome sequences and annotations will facilitate better comparative analysis within monocots.

Nine quick tips for pathway enrichment analysis

Pathway enrichment analysis (PEA) is a computational biology method that identifies biological functions that are overrepresented in a group of genes more than would be expected by chance and ranks these functions by relevance. The relative abundance of genes pertinent to specific pathways is measured through statistical methods, and associated functional pathways are retrieved from online bioinformatics databases. In the last decade, along with the spread of the internet, higher availability of computational resources made PEA software tools easy to access and to use for bioinformatics practitioners worldwide. Although it became easier to use these tools, it also became easier to make mistakes that could generate inflated or misleading results, especially for beginners and inexperienced computational biologists. With this article, we propose nine quick tips to avoid common mistakes and to out a complete, sound, thorough PEA, which can produce relevant and robust results. We describe our nine guidelines in a simple way, so that they can be understood and used by anyone, including students and beginners. Some tips explain what to do before starting a PEA, others are suggestions of how to correctly generate meaningful results, and some final guidelines indicate some useful steps to properly interpret PEA results. Our nine tips can help users perform better pathway enrichment analyses and eventually contribute to a better understanding of current biology.

A Path-Based Analysis of Infected Cell Line and COVID-19 Patient Transcriptome Reveals Novel Potential Targets and Drugs Against SARS-CoV-2

Most transcriptomic studies of SARS-CoV-2 infection have focused on differentially expressed genes, which do not necessarily reveal the genes mediating the transcriptomic changes. In contrast, exploiting curated biological network, our PathExt tool identifies central genes from the differentially active paths mediating global transcriptomic response. Here we apply PathExt to multiple cell line infection models of SARS-CoV-2 and other viruses, as well as to COVID-19 patient-derived PBMCs. The central genes mediating SARS-CoV-2 response in cell lines were uniquely enriched for ATP metabolic process, G1/S transition, leukocyte activation and migration. In contrast, PBMC response reveals dysregulated cell-cycle processes. In PBMC, the most frequently central genes are associated with COVID-19 severity. Importantly, relative to differential genes, PathExt-identified genes show greater concordance with several benchmark anti-COVID-19 target gene sets. We propose six novel anti-SARS-CoV-2 targets ADCY2, ADSL, OCRL, TIAM1, PBK, and BUB1, and potential drugs targeting these genes, such as Bemcentinib, Phthalocyanine, and Conivaptan.

Transcriptional Basis for Haustorium Formation and Host Establishment in Hemiparasitic Psittacanthus schiedeanus Mistletoes

The mistletoe Psittacanthus schiedeanus, a keystone species in interaction networks between plants, pollinators, and seed dispersers, infects a wide range of native and non-native tree species of commercial interest. Here, using RNA-seq methodology we assembled the whole circularized quadripartite structure of P. schiedeanus chloroplast genome and described changes in the gene expression of the nuclear genomes across time of experimentally inoculated seeds. Of the 140,467 assembled and annotated uniGenes, 2,000 were identified as differentially expressed (DEGs) and were classified in six distinct clusters according to their expression profiles. DEGs were also classified in enriched functional categories related to synthesis, signaling, homoeostasis, and response to auxin and jasmonic acid. Since many orthologs are involved in lateral or adventitious root formation in other plant species, we propose that in P. schiedeanus (and perhaps in other rootless mistletoe species), these genes participate in haustorium formation by complex regulatory networks here described. Lastly, and according to the structural similarities of P. schiedeanus enzymes with those that are involved in host cell wall degradation in fungi, we suggest that a similar enzymatic arsenal is secreted extracellularly and used by mistletoes species to easily parasitize and break through tissues of the host.

Common effects of bipolar disorder medications on expression quantitative trait loci genes

The molecular mechanism(s) underpinning the clinical efficacy of the current drugs for bipolar disorder (BD) are largely unknown. This study evaluated the transcriptional perturbations potentially playing roles in the therapeutic efficacy of four commonly prescribed psychotropic drugs used to treat BD. NT2-N cells were treated with lamotrigine, lithium, quetiapine, valproate or vehicle control for 24 h. Genome-wide mRNA expression was quantified by RNA-sequencing. Incorporating drug-induced gene expression profiles with BD-associated transcriptional changes from post-mortem brains, we identified potential therapeutic-relevant genes associated with both drug treatments and BD pathophysiology and focused on expression quantitative trait loci (eQTL) genes with genome-wide association with BD. Each eQTL gene was ranked based on its potential role in the therapeutic effect across multiple drugs. The expression of highest-ranked eQTL genes were measured by RT-qPCR to confirm their transcriptional changes observed in RNA-seq. We found 775 genes for which at least 2 drugs reversed expression levels relative to the differential expression in post-mortem brains. Pathway analysis identified enriched biological processes highlighting mitochondrial and endoplasmic reticulum function. Differential expression of SRPK2 and CHDH was confirmed by RT-qPCR following multiple-dose treatments. We pinpointed potential genes involved in the beneficial effects of drugs used for BD and their main associated biological pathways. CHDH, which encodes a mitochondrial protein, had a significant dose-responsive downregulation following treatment with increasing doses of quetiapine and lamotrigine, which in combination with the enriched mitochondrial pathways suggests potential therapeutic roles and demand more studies on mitochondrial involvement in BD to identify novel treatment targets.

Cross-species transcriptomes reveal species-specific and shared molecular adaptations for plants development on iron-rich rocky outcrops soils

Background

Canga is the Brazilian term for the savanna-like vegetation harboring several endemic species on iron-rich rocky outcrops, usually considered for mining activities. Parkia platycephala Benth. and Stryphnodendron pulcherrimum (Willd.) Hochr. naturally occur in the cangas of Serra dos Carajás (eastern Amazonia, Brazil) and the surrounding forest, indicating high phenotypic plasticity. The morphological and physiological mechanisms of the plants’ establishment in the canga environment are well studied, but the molecular adaptative responses are still unknown. To understand these adaptative responses, we aimed to identify molecular mechanisms that allow the establishment of these plants in the canga environment.

Results

Plants were grown in canga and forest substrates collected in the Carajás Mineral Province. RNA was extracted from pooled leaf tissue, and RNA-seq paired-end reads were assembled into representative transcriptomes for P. platycephala and S. pulcherrimum containing 31,728 and 31,311 primary transcripts, respectively. We identified both species-specific and core molecular responses in plants grown in the canga substrate using differential expression analyses. In the species-specific analysis, we identified 1,112 and 838 differentially expressed genes for P. platycephala and S. pulcherrimum, respectively. Enrichment analyses showed that unique biological processes and metabolic pathways were affected for each species. Comparative differential expression analysis was based on shared single-copy orthologs. The overall pattern of ortholog expression was species-specific. Even so, we identified almost 300 altered genes between plants in canga and forest substrates with conserved responses in the two species. The genes were functionally associated with the response to light stimulus and the circadian rhythm pathway.

Conclusions

Plants possess species-specific adaptative responses to cope with the substrates. Our results also suggest that plants adapted to both canga and forest environments can adjust the circadian rhythm in a substrate-dependent manner. The circadian clock gene modulation might be a central mechanism regulating the plants’ development in the canga substrate in the studied legume species. The mechanism may be shared as a common mechanism to abiotic stress compensation in other native species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08449-0.

High Resistance to Quinclorac in Multiple-Resistant Echinochloa colona Associated with Elevated Stress Tolerance Gene Expression and Enriched Xenobiotic Detoxification Pathway

Echinochloa colona and other species in this genus are a threat to global rice production and food security. Quinclorac, an auxin mimic, is a common herbicide for grass weed control in rice, and Echinochloa spp. have evolved resistance to it. The complete mode of quinclorac action and subsequent evolution of resistance is not fully understood. We analyzed the de novo transcriptome of multiple-herbicide-resistant (ECO-R) and herbicide-susceptible genotypes in response to quinclorac. Several biological processes were constitutively upregulated in ECO-R, including carbon metabolism, photosynthesis, and ureide metabolism, indicating improved metabolic efficiency. The transcriptional change in ECO-R following quinclorac treatment indicates an efficient response, with upregulation of trehalose biosynthesis, which is also known for abiotic stress mitigation. Detoxification-related genes were induced in ECO-R, mainly the UDP-glycosyltransferase (UGT) family, most likely enhancing quinclorac metabolism. The transcriptome data also revealed that many antioxidant defense elements were uniquely elevated in ECO-R to protect against the auxin-mediated oxidative stress. We propose that upon quinclorac treatment, ECO-R detoxifies quinclorac utilizing UGT genes, which modify quinclorac using the sufficient supply of UDP-glucose from the elevated trehalose pathway. Thus, we present the first report of upregulation of trehalose synthesis and its association with the herbicide detoxification pathway as an adaptive mechanism to herbicide stress in Echinochloa, resulting in high resistance.

Nuclear genome of a pedinophyte pinpoints genomic innovation and streamlining in the green algae

The genomic diversity underpinning high ecological and species diversity in the green algae (Chlorophyta) remains little known. Here, we aimed to track genome evolution in the Chlorophyta, focusing on loss and gain of homologous genes, and lineage-specific innovations of the core Chlorophyta. We generated a high-quality nuclear genome for pedinophyte YPF701, a sister lineage to others in the core Chlorophyta and incorporated this genome in a comparative analysis with 25 other genomes from diverse Viridiplantae taxa. The nuclear genome of pedinophyte YPF701 has an intermediate size and gene number between those of most prasinophytes and the remainder of the core Chlorophyta. Our results suggest positive selection for genome streamlining in the Pedinophyceae, independent from genome minimisation observed among prasinophyte lineages. Genome expansion was predicted along the branch leading to the UTC clade (classes Ulvophyceae, Trebouxiophyceae and Chlorophyceae) after divergence from their last common ancestor with pedinophytes, with genomic novelty implicated in a range of basic biological functions. Results emphasise multiple independent signals of genome minimisation within the Chlorophyta, as well as the genomic novelty arising before diversification in the UTC clade, which may underpin the success of this species-rich clade in a diversity of habitats.

The genome and preliminary single-nuclei transcriptome of Lemna minuta reveals mechanisms of invasiveness

The ability to trace every cell in some model organisms has led to the fundamental understanding of development and cellular function. However, in plants the complexity of cell number, organ size, and developmental time makes this a challenge even in the diminutive model plant Arabidopsis (Arabidopsis thaliana). Duckweed, basal nongrass aquatic monocots, provide an opportunity to follow every cell of an entire plant due to their small size, reduced body plan, and fast clonal growth habit. Here we present a chromosome-resolved genome for the highly invasive Lesser Duckweed (Lemna minuta) and generate a preliminary cell atlas leveraging low cell coverage single-nuclei sequencing. We resolved the 360 megabase genome into 21 chromosomes, revealing a core nonredundant gene set with only the ancient tau whole-genome duplication shared with all monocots, and paralog expansion as a result of tandem duplications related to phytoremediation. Leveraging SMARTseq2 single-nuclei sequencing, which provided higher gene coverage yet lower cell count, we profiled 269 nuclei covering 36.9% (8,457) of the L. minuta transcriptome. Since molecular validation was not possible in this nonmodel plant, we leveraged gene orthology with model organism single-cell expression datasets, gene ontology, and cell trajectory analysis to define putative cell types. We found that the tissue that we computationally defined as mesophyll expressed high levels of elemental transport genes consistent with this tissue playing a role in L. minuta wastewater detoxification. The L. minuta genome and preliminary cell map provide a paradigm to decipher developmental genes and pathways for an entire plant.

Differential translational regulation of host exosomal proteins play key role in immunomodulation in antimony resistance in Visceral Leishmaniasis: A proteomic profiling study

In host-pathogen interactions, exosomal secretions are crucial for cell to cell communication and have an established role in immunomodulation. Protozoans, including Leishmania, modulates their host vesicular secretions for better survival; although the role of exosomal secretions in unresponsive against sodium antimony gluconate (SAG) has never been documented. In this study, the exosomal proteome of RAW macrophages infected with either SAG responsive (SAG^S) or SAG unresponsive (SAG^R) L. donovani parasites has been compared with uninfected RAW macrophages. Proteins isolated from exosomes were labelled with iTRAQ reagents; followed by subsequent LC-TOF/-MS analysis. In total, 394 proteins (p < 0.05) were identified which were shared common among all sets. Highly differentially expressed proteins were sorted by log2 value -1 and +1 as down regulated and up regulated respectively which yielded 58 proteins in SAG^R and 41 proteins during SAG^S infection. Out of the 58 proteins identified during SAG^R infection, 17 proteins were of immune modulatory function. Network visualization model and pathway analysis revealed the interactions among these proteins via different immunological pathways with reported involvement of some proteins in SAG resistance and host immune modulation. Hence, the differential abundance of immune pathway related proteins in exosomes of infected host during SAG^R infection supports the immune modulatory strategy adopted by SAG resistant parasites for enhanced survival .

Genome resequencing reveals rapid, repeated evolution in the Colorado potato beetle

Insecticide resistance and rapid pest evolution threatens food security and the development of sustainable agricultural practices, yet the evolutionary mechanisms that allow pests to rapidly adapt to control tactics remains unclear. Here we examine how a global super-pest, the Colorado potato beetle (CPB), Leptinotarsa decemlineata, rapidly evolves resistance to insecticides. Using whole genome resequencing and transcriptomic data focused on its ancestral and pest range in North America, we assess evidence for three, non-mutually exclusive models of rapid evolution: pervasive selection on novel mutations, rapid regulatory evolution, and repeated selection on standing genetic variation. Population genomic analysis demonstrates that CPB is geographically structured, even among recently established pest populations. Pest populations exhibit similar levels of nucleotide diversity, relative to non-pest populations, and show evidence of recent expansion. Genome scans provide clear signatures of repeated adaptation across CPB populations, with especially strong evidence of selection on insecticide resistance genes in different populations. Analyses of gene expression show that constitutive upregulation of candidate insecticide resistance genes drives distinctive population patterns. CPB evolves insecticide resistance repeatedly across agricultural regions, leveraging similar genetic pathways but different genes, demonstrating a polygenic trait architecture for insecticide resistance that can evolve from standing genetic variation. Despite expectations, we do not find support for strong selection on novel mutations, or rapid evolution from selection on regulatory genes. These results suggest that integrated pest management practices must mitigate the evolution of polygenic resistance phenotypes among local pest populations, in order to maintain the efficacy and sustainability of novel control techniques.

OUP accepted manuscript

Exposure to environmental toxicants during preconception has been shown to affect offspring health and epigenetic mechanisms such as DNA methylation are hypothesized to be involved in adverse outcomes. However, studies addressing the effects of exposure to environmental toxicants during preconception on epigenetic changes in gametes are limited. The objective of this study is to determine the effect of preconceptional exposure to a dioxin-like polychlorinated biphenyl (3,3',4,4',5-pentachlorobiphenyl [PCB126]) on DNA methylation and gene expression in testis. Adult zebrafish were exposed to 3 and 10 nM PCB126 for 24 h and testis tissue was sampled at 7 days postexposure for histology, DNA methylation, and gene expression profiling. Reduced representation bisulfite sequencing revealed 37 and 92 differentially methylated regions (DMRs) in response to 3 and 10 nM PCB126 exposures, respectively. Among them, 19 DMRs were found to be common between both PCB126 treatment groups. Gene ontology (GO) analysis of DMRs revealed that enrichment of terms such as RNA processing, iron-sulfur cluster assembly, and gluconeogenesis. Gene expression profiling showed differential expression of 40 and 1621 genes in response to 3 and 10 nM PCB126 exposures, respectively. GO analysis of differentially expressed genes revealed enrichment of terms related to xenobiotic metabolism, oxidative stress, and immune function. There is no overlap in the GO terms or individual genes between DNA methylation and RNA sequencing results, but functionally many of the altered pathways have been shown to cause spermatogenic defects.

GeneTonic: an R/Bioconductor package for streamlining the interpretation of RNA-seq data

BACKGROUND

The interpretation of results from transcriptome profiling experiments via RNA sequencing (RNA-seq) can be a complex task, where the essential information is distributed among different tabular and list formats-normalized expression values, results from differential expression analysis, and results from functional enrichment analyses. A number of tools and databases are widely used for the purpose of identification of relevant functional patterns, yet often their contextualization within the data and results at hand is not straightforward, especially if these analytic components are not combined together efficiently.

RESULTS

We developed the GeneTonic software package, which serves as a comprehensive toolkit for streamlining the interpretation of functional enrichment analyses, by fully leveraging the information of expression values in a differential expression context. GeneTonic is implemented in R and Shiny, leveraging packages that enable HTML-based interactive visualizations for executing drilldown tasks seamlessly, viewing the data at a level of increased detail. GeneTonic is integrated with the core classes of existing Bioconductor workflows, and can accept the output of many widely used tools for pathway analysis, making this approach applicable to a wide range of use cases. Users can effectively navigate interlinked components (otherwise available as flat text or spreadsheet tables), bookmark features of interest during the exploration sessions, and obtain at the end a tailored HTML report, thus combining the benefits of both interactivity and reproducibility.

CONCLUSION

GeneTonic is distributed as an R package in the Bioconductor project ( https://bioconductor.org/packages/GeneTonic/ ) under the MIT license. Offering both bird's-eye views of the components of transcriptome data analysis and the detailed inspection of single genes, individual signatures, and their relationships, GeneTonic aims at simplifying the process of interpretation of complex and compelling RNA-seq datasets for many researchers with different expertise profiles.

Shared brain transcriptomic signature in TDP-43 type A FTLD patients with or without GRN mutations

Frontotemporal lobar degeneration with TDP-43 inclusions (FTLD-TDP) is a complex heterogeneous neurodegenerative disorder for which mechanisms are poorly understood.

To explore transcriptional changes underlying FTLD-TDP, we performed RNA-sequencing on 66 genetically unexplained FTLD-TDP patients, 24 FTLD-TDP patients with GRN mutations and 24 control participants. Using principal component analysis, hierarchical clustering, differential expression and coexpression network analyses, we showed that GRN mutation carriers and FTLD-TDP-A patients without a known mutation shared a common transcriptional signature that is independent of GRN loss-of-function. After combining both groups, differential expression as compared to the control group and coexpression analyses revealed alteration of processes related to immune response, synaptic transmission, RNA metabolism, angiogenesis and vesicle-mediated transport. Deconvolution of the data highlighted strong cellular alterations that were similar in FTLD-TDP-A and GRN mutation carriers with NSF as a potentially important player in both groups. We propose several potentially druggable pathways such as the GABAergic, GDNF and sphingolipid pathways.

Our findings underline new disease mechanisms and strongly suggest that affected pathways in GRN mutation carriers extend beyond GRN and contribute to genetically unexplained forms of FTLD-TDP-A.

Towards a Cardoon (Cynara cardunculus var. altilis)-Based Biorefinery: A Case Study of Improved Cell Cultures via Genetic Modulation of the Phenylpropanoid Pathway

Cultivated cardoon (Cynara cardunculus var. altilis L.) is a promising candidate species for the development of plant cell cultures suitable for large-scale biomass production and recovery of nutraceuticals. We set up a protocol for Agrobacterium tumefaciens-mediated transformation, which can be used for the improvement of cardoon cell cultures in a frame of biorefinery. As high lignin content determines lower saccharification yields for the biomass, we opted for a biotechnological approach, with the purpose of reducing lignin content; we generated transgenic lines overexpressing the Arabidopsis thaliana MYB4 transcription factor, a known repressor of lignin/flavonoid biosynthesis. Here, we report a comprehensive characterization, including metabolic and transcriptomic analyses of AtMYB4 overexpression cardoon lines, in comparison to wild type, underlining favorable traits for their use in biorefinery. Among these, the improved accessibility of the lignocellulosic biomass to degrading enzymes due to depletion of lignin content, the unexpected increased growth rates, and the valuable nutraceutical profiles, in particular for hydroxycinnamic/caffeoylquinic and fatty acids profiles.

The Plasmodium NOT1-G paralogue is an essential regulator of sexual stage maturation and parasite transmission

Productive transmission of malaria parasites hinges upon the execution of key transcriptional and posttranscriptional regulatory events. While much is now known about how specific transcription factors activate or repress sexual commitment programs, far less is known about the production of a preferred mRNA homeostasis following commitment and through the host-to-vector transmission event. Here, we show that in Plasmodium parasites, the NOT1 scaffold protein of the CAF1/CCR4/Not complex is duplicated, and one paralogue is dedicated for essential transmission functions. Moreover, this NOT1-G paralogue is central to the sex-specific functions previously associated with its interacting partners, as deletion of not1-g in Plasmodium yoelii leads to a comparable or complete arrest phenotype for both male and female parasites. We show that, consistent with its role in other eukaryotes, PyNOT1-G localizes to cytosolic puncta throughout much of the Plasmodium life cycle. PyNOT1-G is essential to both the complete maturation of male gametes and to the continued development of the fertilized zygote originating from female parasites. Comparative transcriptomics of wild-type and pynot1-g⁻ parasites shows that loss of PyNOT1-G leads to transcript dysregulation preceding and during gametocytogenesis and shows that PyNOT1-G acts to preserve mRNAs that are critical to sexual and early mosquito stage development. Finally, we demonstrate that the tristetraprolin (TTP)-binding domain, which acts as the typical organization platform for RNA decay (TTP) and RNA preservation (ELAV/HuR) factors is dispensable for PyNOT1-G’s essential blood stage functions but impacts host-to-vector transmission. Together, we conclude that a NOT1-G paralogue in Plasmodium fulfills the complex transmission requirements of both male and female parasites.

Malaria parasites face two bottlenecks in their life cycle: their two transmission events. This study shows that Plasmodium has taken the unorthodox approach of duplicating the gene for the NOT1 RNA regulatory scaffold protein, allowing it to dedicate one paralog to functions that are essential for transmission from mammalian hosts to the mosquito vector.

Identification and Elucidation of the Protective isomiRs in Lung Cancer Patient Prognosis

MicroRNAs (miRNAs) are approximately 20–22 nucleotides in length, which are well known to participate in the post-transcriptional modification. The mature miRNAs were observed to be varied on 5′ or 3′ that raise another term—the isoforms of mature miRNAs (isomiRs), which have been proven not the artifacts and discussed widely recently. In our research, we focused on studying the 5′ isomiRs in lung adenocarcinoma (LUAD) in The Cancer Genome Atlas (TCGA). We characterized 75 isomiRs significantly associated with better prognosis and 43 isomiRs with poor prognosis. The 75 protective isomiRs can successfully distinguish tumors from normal samples and are expressed differently between patients of early and late stages. We also found that most of the protective isomiRs tend to be with downstream shift and upregulated compared with those with upstream shift, implying that a possible selection occurs during cancer development. Among these protective isomiRs, we observed a highly positive and significant correlation, as well as in harmful isomiRs, suggesting cooperation within the group. However, between protective and harmful, there is no such a concordance but conversely more negative correlation, suggesting the possible antagonistic effect between protective and harmful isomiRs. We also identified that two isomiRs miR-181a-3p|-3 and miR-181a-3p|2, respectively, belong to the harmful and protective groups, suggesting a bidirectional regulation of their originated archetype—miR-181a-3p. Additionally, we found that the protective isomiRs of miR-21-5p, which is an oncomiR, may be evolved as the tumor suppressors through producing isomiRs to hinder metastasis. In summary, these results displayed the characteristics of the protective isomiRs and their potential for developing the treatment of lung cancer.

Gelatin methacrylate hydrogels culture model for glioblastoma cells enriches for mesenchymal-like state and models interactions with immune cells

Glioblastoma is the most lethal primary malignant brain tumor in adults. Simplified two-dimensional (2D) cell culture and neurospheres in vitro models fail to recapitulate the complexity of the tumor microenvironment, limiting its ability to predict therapeutic response. Three-dimensional (3D) scaffold-based models have emerged as a promising alternative for addressing these concerns. One such 3D system is gelatin methacrylate (GelMA) hydrogels, and we aimed to understand the suitability of using this system to mimic treatment-resistant glioblastoma cells that reside in specific niches. We characterized the phenotype of patient-derived glioma cells cultured in GelMA hydrogels (3D-GMH) for their tumorigenic properties using invasion and chemoresponse assays. In addition, we used integrated single-cell and spatial transcriptome analysis to compare cells cultured in 3D-GMH to neoplastic cells in vivo. Finally, we assessed tumor-immune cell interactions with a macrophage infiltration assay and a cytokine array. We show that the 3D-GMH system enriches treatment-resistant mesenchymal cells that are not represented in neurosphere cultures. Cells cultured in 3D-GMH resemble a mesenchymal-like cellular phenotype found in perivascular and hypoxic regions and recruit macrophages by secreting cytokines, a hallmark of the mesenchymal phenotype. Our 3D-GMH model effectively mimics the phenotype of glioma cells that are found in the perivascular and hypoxic niches of the glioblastoma core in situ, in contrast to the neurosphere cultures that enrich cells of the infiltrative edge of the tumor. This contrast highlights the need for due diligence in selecting an appropriate model when designing a study's objectives.

Examining the Developmental Trajectory of an in Vitro Model of Mouse Primordial Germ Cells following Exposure to Environmentally Relevant Bisphenol A Levels

BACKGROUND

Animal-based studies indicate that bisphenol A (BPA) exposure is detrimental to reproductive health, but its impact on the earliest stages of germ cell development remains poorly defined.

OBJECTIVES

Using a murine in vitro model of early germ cell specification and differentiation, we sought to assess whether exposure to low levels of BPA prior to formation of primordial germ cells (PGCs) alters their differentiation trajectory and unique molecular program.

METHODS

We used an established method of in vitro differentiation of mouse embryonic stem cells (ESCs) into epiblast-like cells (EpiLCs) followed by PGC-like cells (PGCLCs), which together recapitulate defined stages of early germ cell development. Cellular consequences were determined using hemocytometer-based cell counting, fixation, and intracellular staining, followed by flow cytometry/fluorescence-activated cell sorting (FACS) of cells exposed to increasing concentrations (range: 1 nM-10 μM) of BPA. To interrogate and characterize gene expression differences resulting from BPA exposure, we also generated RNA-seq libraries from RNA extracted from FACS-purified PGCLCs and performed transcriptome analysis using bioinformatics-based approaches.

RESULTS

Exposure of EpiLCs to BPA resulted in higher numbers of cells that were associated with a higher proportion of cells in S-phase as well as a lower proportion undergoing apoptosis; this difference occurred in a concentration-dependent manner. Exposure also resulted in a greater fraction of EpiLCs showing signs of DNA damage. Remarkably, EpiLC exposure did not negatively affect PGC specification and resulted in a concentration-dependent effect on PGCLC proliferation in XX but not XY cells. PGCLC transcriptome analysis revealed an aberrant program with significant deregulation of X-linked genes and retrotransposon expression. Differential gene expression analysis also revealed the deregulation of genes associated with lipid metabolism as well as deregulated expression of genes associated with later stages of gametogenesis.

CONCLUSIONS

To the best of our knowledge our findings represent the first characterization of the consequences of early BPA exposure on a model of mammalian PGC development, highlighting altered cell behavior, altered underlying pathways, and altered molecular processes. https://doi.org/10.1289/EHP8196.

Effects of chronic exposure to low levels of IR on Medaka (Oryzias latipes): a proteomic and bioinformatic approach

PURPOSE

Chronic exposure to ionizing radiation (IR) at low doses (<100 mGy) has been insufficiently studied to understand fully the risk to health. Relatively little knowledge exists regarding how species and healthy tissues respond at the protein level to chronic exposure to low doses of IR, and mass spectrometric-based profiling of protein expression is a powerful tool for studying changes in protein abundance.

MATERIALS AND METHODS

SDS gel electrophoresis, LC-MS/MS mass spectrometry-based approaches and bioinformatic data analytics were used to detect proteomic changes following chronic exposure to moderate/low doses of radiation in adults and normally developed Medaka fish (Oryzias latipes).

RESULTS

Significant variations in the abundance of proteins involved in thyroid hormone signaling and lipid metabolism were detected, which could be related to the gonadal regression phenotype observed after 21.04 mGy and 204.3 mGy/day exposure. The global proteomic change was towards overexpression of proteins in muscle and skin, while the opposite effect was observed in internal organs.

CONCLUSION

The present study provides information on the impacts of biologically relevant low doses of IR, which will be useful in future research for the identification of potential biomarkers of IR exposure and allow for a better assessment of radiation biosafety regulations.

Functional characterization of the cnidarian antiviral immune response reveals ancestral complexity

Animals evolved a broad repertoire of innate immune sensors and downstream effector cascades for defense against RNA viruses. Yet, this system varies greatly among different bilaterian animals, masking its ancestral state. In this study, we aimed to characterize the antiviral immune response of the cnidarian Nematostella vectensis and decipher the function of the retinoic acid-inducible gene I (RIG-I)-like receptors (RLRs) known to detect viral double-stranded RNA (dsRNA) in bilaterians but activate different antiviral pathways in vertebrates and nematodes. We show that polyinosinic:polycytidylic acid (poly(I:C)), a mimic of long viral dsRNA and a primary ligand for the vertebrate RLR melanoma differentiation-associated protein 5 (MDA5), triggers a complex antiviral immune response bearing features distinctive for both vertebrate and invertebrate systems. Importantly, a well-characterized agonist of the vertebrate RIG-I receptor does not induce a significant transcriptomic response that bears signature of the antiviral immune response, which experimentally supports the results of a phylogenetic analysis indicating clustering of the two N. vectensis RLR paralogs (NveRLRa and NveRLRb) with MDA5. Furthermore, the results of affinity assays reveal that NveRLRb binds poly(I:C) and long dsRNA and its knockdown impairs the expression of putative downstream effector genes including RNA interference components. Our study provides for the first time the functional evidence for the conserved role of RLRs in initiating immune response to dsRNA that originated before the cnidarian–bilaterian split and lay a strong foundation for future research on the evolution of the immune responses to RNA viruses.

PIGNON: a protein-protein interaction-guided functional enrichment analysis for quantitative proteomics

Background

Quantitative proteomics studies are often used to detect proteins that are differentially expressed across different experimental conditions. Functional enrichment analyses are then typically used to detect annotations, such as biological processes that are significantly enriched among such differentially expressed proteins to provide insights into the molecular impacts of the studied conditions. While common, this analytical pipeline often heavily relies on arbitrary thresholds of significance. However, a functional annotation may be dysregulated in a given experimental condition, while none, or very few of its proteins may be individually considered to be significantly differentially expressed. Such an annotation would therefore be missed by standard approaches.

Results

Herein, we propose a novel graph theory-based method, PIGNON, for the detection of differentially expressed functional annotations in different conditions. PIGNON does not assess the statistical significance of the differential expression of individual proteins, but rather maps protein differential expression levels onto a protein–protein interaction network and measures the clustering of proteins from a given functional annotation within the network. This process allows the detection of functional annotations for which the proteins are differentially expressed and grouped in the network. A Monte-Carlo sampling approach is used to assess the clustering significance of proteins in an expression-weighted network. When applied to a quantitative proteomics analysis of different molecular subtypes of breast cancer, PIGNON detects Gene Ontology terms that are both significantly clustered in a protein–protein interaction network and differentially expressed across different breast cancer subtypes. PIGNON identified functional annotations that are dysregulated and clustered within the network between the HER2+, triple negative and hormone receptor positive subtypes. We show that PIGNON’s results are complementary to those of state-of-the-art functional enrichment analyses and that it highlights functional annotations missed by standard approaches. Furthermore, PIGNON detects functional annotations that have been previously associated with specific breast cancer subtypes.

Conclusion

PIGNON provides an alternative to functional enrichment analyses and a more comprehensive characterization of quantitative datasets. Hence, it contributes to yielding a better understanding of dysregulated functions and processes in biological samples under different experimental conditions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04042-6.

A high-quality carabid genome assembly provides insights into beetle genome evolution and cold adaptation

The hyperdiverse order Coleoptera comprises a staggering ~25% of known species on Earth. Despite recent breakthroughs in next generation sequencing, there remains a limited representation of beetle diversity in assembled genomes. Most notably, the ground beetle family Carabidae, comprising more than 40,000 described species, has not been studied in a comparative genomics framework using whole genome data. Here we generate a high-quality genome assembly for Nebria riversi, to examine sources of novelty in the genome evolution of beetles, as well as genetic changes associated with specialization to high-elevation alpine habitats. In particular, this genome resource provides a foundation for expanding comparative molecular research into mechanisms of insect cold adaptation. Comparison to other beetles shows a strong signature of genome compaction, with N. riversi possessing a relatively small genome (~147 Mb) compared to other beetles, with associated reductions in repeat element content and intron length. Small genome size is not, however, associated with fewer protein-coding genes, and an analysis of gene family diversity shows significant expansions of genes associated with cellular membranes and membrane transport, as well as protein phosphorylation and muscle filament structure. Finally, our genomic analyses show that these high-elevation beetles have endosymbiotic Spiroplasma, with several metabolic pathways (e.g., propanoate biosynthesis) that might complement N. riversi, although its role as a beneficial symbiont or as a reproductive parasite remains equivocal.

Membrane-Enriched Proteomics Link Ribosome Accumulation and Proteome Reprogramming With Cold Acclimation in Barley Root Meristems

Due to their sessile nature, plants rely on root systems to mediate many biotic and abiotic cues. To overcome these challenges, the root proteome is shaped to specific responses. Proteome-wide reprogramming events are magnified in meristems due to their active protein production. Using meristems as a test system, here, we study the major rewiring that plants undergo during cold acclimation. We performed tandem mass tag-based bottom-up quantitative proteomics of two consecutive segments of barley seminal root apexes subjected to suboptimal temperatures. After comparing changes in total and ribosomal protein (RP) fraction-enriched contents with shifts in individual protein abundances, we report ribosome accumulation accompanied by an intricate translational reprogramming in the distal apex zone. Reprogramming ranges from increases in ribosome biogenesis to protein folding factors and suggests roles for cold-specific RP paralogs. Ribosome biogenesis is the largest cellular investment; thus, the vast accumulation of ribosomes and specific translation-related proteins during cold acclimation could imply a divergent ribosomal population that would lead to a proteome shift across the root. Consequently, beyond the translational reprogramming, we report a proteome rewiring. First, triggered protein accumulation includes spliceosome activity in the root tip and a ubiquitous upregulation of glutathione production and S-glutathionylation (S-GSH) assemblage machineries in both root zones. Second, triggered protein depletion includes intrinsically enriched proteins in the tip-adjacent zone, which comprise the plant immune system. In summary, ribosome and translation-related protein accumulation happens concomitantly to a proteome reprogramming in barley root meristems during cold acclimation. The cold-accumulated proteome is functionally implicated in feedbacking transcript to protein translation at both ends and could guide cold acclimation.

Characterization of the Survival Influential Genes in Carcinogenesis

The genes influencing cancer patient mortality have been studied by survival analysis for many years. However, most studies utilized them only to support their findings associated with patient prognosis: their roles in carcinogenesis have not yet been revealed. Herein, we applied an in silico approach, integrating the Cox regression model with effect size estimated by the Monte Carlo algorithm, to screen survival-influential genes in more than 6000 tumor samples across 16 cancer types. We observed that the survival-influential genes had cancer-dependent properties. Moreover, the functional modules formed by the harmful genes were consistently associated with cell cycle in 12 out of the 16 cancer types and pan-cancer, showing that dysregulation of the cell cycle could harm patient prognosis in cancer. The functional modules formed by the protective genes are more diverse in cancers; the most prevalent functions are relevant for immune response, implying that patients with different cancer types might develop different mechanisms against carcinogenesis. We also identified a harmful set of 10 genes, with potential as prognostic biomarkers in pan-cancer. Briefly, our results demonstrated that the survival-influential genes could reveal underlying mechanisms in carcinogenesis and might provide clues for developing therapeutic targets for cancers.

Transcriptomic Diversity in the Livers of South African Sardines Participating in the Annual Sardine Run

During austral winter, the southern and eastern coastlines of South Africa witness one of the largest animal migrations on the planet, the KwaZulu-Natal sardine run. Hundreds of millions of temperate sardines, Sardinops sagax, form large shoals that migrate north-east towards the subtropical Indian Ocean. Recent studies have highlighted the role that genetic and environmental factors play in sardine run formation. In the present study, we used massively parallel sequencing to assemble and annotate the first reference transcriptome from the liver cells of South African sardines, and to investigate the functional content and transcriptomic diversity. A total of 1,310,530 transcripts with an N50 of 1578 bp were assembled de novo. Several genes and core biochemical pathways that modulate energy production, energy storage, digestion, secretory processes, immune responses, signaling, regulatory processes, and detoxification were identified. The functional content of the liver transcriptome from six individuals that participated in the 2019 sardine run demonstrated heterogeneous levels of variation. Data presented in the current study provide new insights into the complex function of the liver transcriptome in South African sardines.

Temporal modulation of the NF-κB RelA network in response to different types of DNA damage

Different types of DNA damage can initiate phosphorylation-mediated signalling cascades that result in stimulus specific pro- or anti-apoptotic cellular responses. Amongst its many roles, the NF-κB transcription factor RelA is central to these DNA damage response pathways. However, we still lack understanding of the co-ordinated signalling mechanisms that permit different DNA damaging agents to induce distinct cellular outcomes through RelA. Here, we use label-free quantitative phosphoproteomics to examine the temporal effects of exposure of U2OS cells to either etoposide (ETO) or hydroxyurea (HU) by monitoring the phosphorylation status of RelA and its protein binding partners. Although few stimulus-specific differences were identified in the constituents of phosphorylated RelA interactome after exposure to these DNA damaging agents, we observed subtle, but significant, changes in their phosphorylation states, as a function of both type and duration of treatment. The DNA double strand break (DSB)-inducing ETO invoked more rapid, sustained responses than HU, with regulated targets primarily involved in transcription, cell division and canonical DSB repair. Kinase substrate prediction of ETO-regulated phosphosites suggest abrogation of CDK and ERK1 signalling, in addition to the known induction of ATM/ATR. In contrast, HU-induced replicative stress mediated temporally dynamic regulation, with phosphorylated RelA binding partners having roles in rRNA/mRNA processing and translational initiation, many of which contained a 14-3-3ε binding motif, and were putative substrates of the dual specificity kinase CLK1. Our data thus point to differential regulation of key cellular processes and the involvement of distinct signalling pathways in modulating DNA damage-specific functions of RelA.

Ouabain-Induced Gene Expression Changes in Human iPSC-Derived Neuron Culture Expressing Dopamine and cAMP-Regulated Phosphoprotein 32 and GABA Receptors

Cardiotonic steroids (CTS) are specific inhibitors and endogenous ligands of a key enzyme in the CNS-the Na⁺, K⁺-ATPase, which maintains and creates an ion gradient on the plasma membrane of neurons. CTS cause the activation of various signaling cascades and changes in gene expression in neurons and other cell types. It is known that intracerebroventricular injection of cardiotonic steroid ouabain causes mania-like behavior in rodents, in part due to activation of dopamine-related signaling cascades in the dopamine and cAMP-regulated phosphoprotein 32 (DARPP-32) expressing medium spiny neurons in the striatum. Dopaminergic projections in the striatum innervate these GABAergic medium spiny neurons. The objective of this study was to assess changes in the expression of all genes in human iPSC-derived expressing DARPP-32 and GABA receptors neurons under the influence of ouabain. We noted a large number of statistically significant upregulated and downregulated genes after a 16-h incubation with non-toxic concentration (30 nM) of ouabain. These changes in the transcriptional activity were accomplished with activation of MAP-kinase ERK1/2 and transcriptional factor cAMP response element-binding protein (CREB). Thus, it can be concluded that 30 nM ouabain incubated for 16 h with human iPSC-derived expressing DARPP-32 and GABA receptors neurons activates genes associated with neuronal maturation and synapse formation, by increasing the expression of genes associated with translation, vesicular transport, and increased electron transport chain function. At the same time, the expression of genes associated with proliferation, migration, and early development of neurons decreases. These data indicate that non-toxic concentrations of ouabain may induce neuronal maturation, neurite growth, and increased synaptogenesis in dopamine-receptive GABAergic neurons, suggesting formation of plasticity and the establishment of new neuronal junctions.

Generation of an Interactome for the Tetratricopeptide Repeat Domain of O-GlcNAc Transferase Indicates a Role for the Enzyme in Intellectual Disability

The O-GlcNAc transferase (OGT) modifies nuclear and cytoplasmic proteins with β-N-acetyl-glucosamine (O-GlcNAc). With thousands of O-GlcNAc-modified proteins but only one OGT encoded in the mammalian genome, a prevailing question is how OGT selects its substrates. Prior work has indicated that the tetratricopeptide repeat (TPR) domain of OGT is involved in substrate selection. Furthermore, several variants of OGT causal for X-linked intellectual disability (XLID) occur in the TPR domain. Therefore, we adapted the BioID labeling method to identify interactors of a TPR-BirA* fusion protein in HeLa cells. We identified 115 interactors representing known and novel O-GlcNAc-modified proteins and OGT interactors (raw data deposited in MassIVE, Dataset ID MSV000085626). The interactors are enriched in known OGT processes (e.g., chromatin remodeling) as well as processes in which OGT has yet to be implicated (e.g., pre-mRNA processing). Importantly, the identified TPR interactors are linked to several disease states but most notably are enriched in pathologies featuring intellectual disability that may underlie the mechanism by which mutations in OGT lead to XLID. This interactome for the TPR domain of OGT serves as a jumping-off point for future research exploring the role of OGT, the TPR domain, and its protein interactors in multiple cellular processes and disease mechanisms, including intellectual disability.

Shotgun proteomics of peach fruit reveals major metabolic pathways associated to ripening

BACKGROUND

Fruit ripening in Prunus persica melting varieties involves several physiological changes that have a direct impact on the fruit organoleptic quality and storage potential. By studying the proteomic differences between the mesocarp of mature and ripe fruit, it would be possible to highlight critical molecular processes involved in the fruit ripening.

RESULTS

To accomplish this goal, the proteome from mature and ripe fruit was assessed from the variety O'Henry through shotgun proteomics using 1D-gel (PAGE-SDS) as fractionation method followed by LC/MS-MS analysis. Data from the 131,435 spectra could be matched to 2740 proteins, using the peach genome reference v1. After data pre-treatment, 1663 proteins could be used for comparison with datasets assessed using transcriptomic approaches and for quantitative protein accumulation analysis. Close to 26% of the genes that code for the proteins assessed displayed higher expression at ripe fruit compared to other fruit developmental stages, based on published transcriptomic data. Differential accumulation analysis between mature and ripe fruit revealed that 15% of the proteins identified were modulated by the ripening process, with glycogen and isocitrate metabolism, and protein localization overrepresented in mature fruit, as well as cell wall modification in ripe fruit. Potential biomarkers for the ripening process, due to their differential accumulation and gene expression pattern, included a pectin methylesterase inhibitor, a gibbellerin 2-beta-dioxygenase, an omega-6 fatty acid desaturase, a homeobox-leucine zipper protein and an ACC oxidase. Transcription factors enriched in NAC and Myb protein domains would target preferentially the genes encoding proteins more abundant in mature and ripe fruit, respectively.

CONCLUSIONS

Shotgun proteomics is an unbiased approach to get deeper into the proteome allowing to detect differences in protein abundance between samples. This technique provided a resolution so that individual gene products could be identified. Many proteins likely involved in cell wall and sugar metabolism, aroma and color, change their abundance during the transition from mature to ripe fruit.

Urbanization is associated with modifications in DNA methylation in a small passerine bird

Urbanization represents a fierce driver of phenotypic change, yet the molecular mechanisms underlying observed phenotypic patterns are poorly understood. Epigenetic changes are expected to facilitate more rapid adaption to changing or novel environments, such as our towns and cities, compared with slow changes in gene sequence. A comparison of liver and blood tissue from great tits Parus major originating from an urban and a forest site demonstrated that urbanization is associated with variation in genome-wide patterns of DNA methylation. Combining reduced representation bisulphite sequencing with transcriptome data, we revealed habitat differences in DNA methylation patterns that suggest a regulated and coordinated response to the urban environment. In the liver, genomic sites that were differentially methylated between urban- and forest-dwelling birds were over-represented in regulatory regions of the genome and more likely to occur in expressed genes. DNA methylation levels were also inversely correlated with gene expression at transcription start sites. Furthermore, differentially methylated CpG sites, in liver, were over-represented in pathways involved in (i) steroid biosynthesis, (ii) superoxide metabolism, (iii) secondary alcohol metabolism, (iv) chylomicron remodelling, (v) cholesterol transport, (vi) reactive oxygen species (ROS) metabolic process and (vii) epithelial cell proliferation. This corresponds with earlier studies identifying diet and exposure to ROS as two of the main drivers of divergence between organisms in urban and nonurban environments. Conversely, in blood, sites that were differentially methylated between urban- and forest-dwelling birds were under-represented in regulatory regions, more likely to occur in nonexpressed genes and not over-represented in specific biological pathways. It remains to be determined whether diverging patterns of DNA methylation represent adaptive evolutionary responses and whether the conclusions can be more widely attributed to urbanization.