ShinyGO: a graphical gene-set enrichment tool for animals and plants

Few Fixed Variants between Trophic Specialist Pupfish Species Reveal Candidate Cis-Regulatory Alleles Underlying Rapid Craniofacial Divergence

Investigating closely related species that rapidly evolved divergent feeding morphology is a powerful approach to identify genetic variation underlying variation in complex traits. This can also lead to the discovery of novel candidate genes influencing natural and clinical variation in human craniofacial phenotypes. We combined whole-genome resequencing of 258 individuals with 50 transcriptomes to identify candidate cis-acting genetic variation underlying rapidly evolving craniofacial phenotypes within an adaptive radiation of Cyprinodon pupfishes. This radiation consists of a dietary generalist species and two derived trophic niche specialists-a molluscivore and a scale-eating species. Despite extensive morphological divergence, these species only diverged 10 kya and produce fertile hybrids in the laboratory. Out of 9.3 million genome-wide SNPs and 80,012 structural variants, we found very few alleles fixed between species-only 157 SNPs and 87 deletions. Comparing gene expression across 38 purebred F1 offspring sampled at three early developmental stages, we identified 17 fixed variants within 10 kb of 12 genes that were highly differentially expressed between species. By measuring allele-specific expression in F1 hybrids from multiple crosses, we found that the majority of expression divergence between species was explained by trans-regulatory mechanisms. We also found strong evidence for two cis-regulatory alleles affecting expression divergence of two genes with putative effects on skeletal development (dync2li1 and pycr3). These results suggest that SNPs and structural variants contribute to the evolution of novel traits and highlight the utility of the San Salvador Island pupfish system as an evolutionary model for craniofacial development.

From Spaceflight to Mars g-Levels: Adaptive Response of A. Thaliana Seedlings in a Reduced Gravity Environment Is Enhanced by Red-Light Photostimulation

The response of plants to the spaceflight environment and microgravity is still not well understood, although research has increased in this area. Even less is known about plants' response to partial or reduced gravity levels. In the absence of the directional cues provided by the gravity vector, the plant is especially perceptive to other cues such as light. Here, we investigate the response of Arabidopsis thaliana 6-day-old seedlings to microgravity and the Mars partial gravity level during spaceflight, as well as the effects of red-light photostimulation by determining meristematic cell growth and proliferation. These experiments involve microscopic techniques together with transcriptomic studies. We demonstrate that microgravity and partial gravity trigger differential responses. The microgravity environment activates hormonal routes responsible for proliferation/growth and upregulates plastid/mitochondrial-encoded transcripts, even in the dark. In contrast, the Mars gravity level inhibits these routes and activates responses to stress factors to restore cell growth parameters only when red photostimulation is provided. This response is accompanied by upregulation of numerous transcription factors such as the environmental acclimation-related WRKY-domain family. In the long term, these discoveries can be applied in the design of bioregenerative life support systems and space farming.

Synergism of Proneurogenic miRNAs Provides a More Effective Strategy to Target Glioma Stem Cells

Simple Summary

miRNAs function as critical regulators of gene expression and have been defined as contributors of cancer phenotypes by acting as oncogenes or tumor suppressors. Based on these findings, miRNA-based therapies have been explored in the treatment of many different malignancies. The use of single miRNAs has faced some challenges and showed limited success. miRNAs cooperate to regulate distinct biological processes and pathways and, therefore, combination of related miRNAs could amplify the repression of oncogenic factors and the effect on cancer relevant pathways. We established that the combination of tumor suppressor miRNAs miR-124, miR-128, and miR-137 is much more effective than single miRNAs in disrupting proliferation and survival of glioma stem cells and neuroblastoma lines and promoting differentiation and response to radiation. Subsequent genomic analyses showed that other combinations of tumor suppressor miRNAs could be equally effective, and its use could provide new routes to target in special cancer-initiating cell populations.

Abstract

Tumor suppressor microRNAs (miRNAs) have been explored as agents to target cancer stem cells. Most strategies use a single miRNA mimic and present many disadvantages, such as the amount of reagent required and the diluted effect on target genes. miRNAs work in a cooperative fashion to regulate distinct biological processes and pathways. Therefore, we propose that miRNA combinations could provide more efficient ways to target cancer stem cells. We have previously shown that miR-124, miR-128, and miR-137 function synergistically to regulate neurogenesis. We used a combination of these three miRNAs to treat glioma stem cells and showed that this treatment was much more effective than single miRNAs in disrupting cell proliferation and survival and promoting differentiation and response to radiation. Transcriptomic analyses indicated that transcription regulation, angiogenesis, metabolism, and neuronal differentiation are among the main biological processes affected by transfection of this miRNA combination. In conclusion, we demonstrated the value of using combinations of neurogenic miRNAs to disrupt cancer phenotypes and glioma stem cell growth. The synergistic effect of these three miRNA amplified the repression of oncogenic factors and the effect on cancer relevant pathways. Future therapeutic approaches would benefit from utilizing miRNA combinations, especially when targeting cancer-initiating cell populations.

Single-cell transcriptomic analysis reveals the immune landscape of lung in steroid-resistant asthma exacerbation

Exaggerated airway hyperresponsiveness and inflammation are hallmarks of asthma, and lipopolysaccharide (LPS) exposure is linked to the severity of the disease and steroid resistance. To investigate the mechanisms underlying asthma exacerbation, we established a mouse model of LPS-induced steroid-resistant exacerbation on the background of house dust mite (HDM)-induced asthma to profile the immune cells in lung by using single-cell RNA deep sequencing. Twenty immune subsets were identified by their molecular and functional properties. Specific cell clusters of basophils, type 2 innate lymphoid cells (ILC2), and CD8+ memory T cells were the predominant sources of interleukin (IL)-4 and IL-13 transcripts whose expressions were dexamethasone resistant. Production of IL-13 by these cells was validated by IL-13-reporter mice. Neutralization of IL-13 abolished HDM/LPS-induced airway hyperresponsiveness, airway inflammation, and decreased mucus hypersecretion. Furthermore, using Ingenuity Pathway Analysis systems, we identified canonical pathways and upstream regulators that regulate the activation of basophils, ILC2, and CD8+ memory T cells. Our study provides mechanistic insights and an important reference resource for further understanding of the immune landscape during asthma exacerbation.

Role of Transportome in the Gills of Chinese Mitten Crabs in Response to Salinity Change: A Meta-Analysis of RNA-Seq Datasets

Chinese mitten crab (CMC) or Eriocheir sinensis is a strong osmoregulator that can keep rigorous cellular homeostasis. CMC can flourish in freshwater, as well as seawater, habitats and represents the most important species for freshwater aquaculture. Salt stress can have direct effects on several stages (e.g., reproduction, molting, growth, etc.) of the CMC life cycle. To get a better overview of the genes involved in the gills of CMC under different salinity conditions, we conducted an RNA-Seq meta-analysis on the transcriptomes of four publicly available datasets. The meta-analysis identified 405 differentially expressed transcripts (DETs), of which 40% were classified into various transporter classes, including accessory factors and primary active transporters as the major transport classes. A network analysis of the DETs revealed that adaptation to salinity is a highly regulated mechanism in which different functional modules play essential roles. To the best of our knowledge, this study is the first to conduct a transcriptome meta-analysis of gills from crab RNA-Seq datasets under salinity. Additionally, this study is also the first to focus on the differential expression of diverse transporters and channels (transportome) in CMC. Our meta-analysis opens new avenues for a better understanding of the osmoregulation mechanism and the selection of potential transporters associated with salinity change.

Arabidopsis bZIP18 and bZIP52 Accumulate in Nuclei Following Heat Stress where They Regulate the Expression of a Similar Set of Genes

Heat stress (HS) is a major abiotic stress that negatively impacts crop yields across the globe. Plants respond to elevated temperatures by changing gene expression, mediated by transcription factors (TFs) functioning to enhance HS tolerance. The involvement of Group I bZIP TFs in the heat stress response (HSR) is not known. In this study, bZIP18 and bZIP52 were investigated for their possible role in the HSR. Localization experiments revealed their nuclear accumulation following heat stress, which was found to be triggered by dephosphorylation. Both TFs were found to possess two motifs containing serine residues that are candidates for phosphorylation. These motifs are recognized by 14–3–3 proteins, and bZIP18 and bZIP52 were found to bind 14–3–3 ε, the interaction of which sequesters them to the cytoplasm. Mutation of both residues abolished 14–3–3 ε interaction and led to a strict nuclear localization for both TFs. RNA-seq analysis revealed coordinated downregulation of several metabolic pathways including energy metabolism and translation, and upregulation of numerous lncRNAs in particular. These results support the idea that bZIP18 and bZIP52 are sequestered to the cytoplasm under control conditions, and that heat stress leads to their re-localization to nuclei, where they jointly regulate gene expression.

Arabidopsis MED18 Interaction With RNA Pol IV and V Subunit NRPD2a in Transcriptional Regulation of Plant Immune Responses

Mediator is a conserved multiprotein complex important for transcription by RNA polymerase II (Pol II). Arabidopsis Mediator subunit MED18 regulates flowering, hormone signaling and plant immunity. Here we report that Arabidopsis MED18 interacted with NUCLEAR RNA POLYMERASE D2a (NRPD2a), the second largest subunit of the nuclear Pol IV and V, which function in RNA-directed DNA methylation and epigenetic regulation of gene expression. Mutants for both MED18 and NRPD2a were compromised in resistance to necrotrophic fungal pathogen Botrytis cinerea. Mutants for NRPD1a, the largest subunit of Pol IV, were also compromised in resistance to Botrytis, supporting a critical role of Pol IV and V in plant defense against Botrytis. Increased Botrytis susceptibility of both the med18 and nrpd2a mutants were associated with reduced accumulation of reactive oxygen species, which are known to promote resistance to Botrytis. Both the basal and pathogen-induced levels of salicylic acid and jasmonic acid were also significantly altered in the med18 and nrpd2a mutants. Transcriptome profiling found that MED18 and NRPD2a affected both unique and overlapping sets of genes in a broad spectrum of biological processes and pathways that influence plant-pathogen interaction. The genes altered in expression in the med18 and nrpd2a mutants include disease resistance proteins, salicylic acid and jasmonic acid signaling and responses, which are known to affect resistance to necrotrophic pathogens. The novel interaction between subunits of Mediator and plant-specific RNA polymerases provides a new mechanism for epigenetic regulation of resistance and expression of defense-related genes in plant immunity.

Chemosensory Neurons Modulate the Response to Oomycete Recognition in Caenorhabditis elegans

Understanding how animals detect and respond to pathogen threats is central to dissecting mechanisms of host immunity. The oomycetes represent a diverse eukaryotic group infecting various hosts from nematodes to humans. We have previously shown that Caenorhabditis elegans mounts a defense response consisting of the induction of chitinase-like (chil) genes in the epidermis to combat infection by its natural oomycete pathogen Myzocytiopsis humicola. We provide here evidence that C. elegans can sense the oomycete by detecting an innocuous extract derived from animals infected with M. humicola. The oomycete recognition response (ORR) leads to changes in the cuticle and reduction in pathogen attachment, thereby increasing animal survival. We also show that TAX-2/TAX-4 function in chemosensory neurons is required for the induction of chil-27 in the epidermis in response to extract exposure. Our findings highlight that neuron-to-epidermis communication may shape responses to oomycete recognition in animal hosts.

•

C. elegans senses its natural oomycete pathogen M. humicola without infection

•

Exposure to a pathogen extract triggers an oomycete recognition response

•

Upon pathogen detection, C. elegans resists infection through changes in the cuticle

•

The response involves signaling between sensory neurons and the epidermis

Bacterial membrane vesicles and phages in blood after consumption of lacticaseibacillus rhamnosus JB-1

Gut microbiota have myriad roles in host physiology, development, and immunity. Though confined to the intestinal lumen by the epithelia, microbes influence distal systems via poorly characterized mechanisms. Recent work has considered the role of extracellular vesicles in interspecies communication, but whether they are involved in systemic microbe-host interaction is unclear. Here, we show that distinctive nanoparticles can be isolated from mouse blood within 2.5 h of consuming Lacticaseibacillus rhamnosus JB-1. In contrast to blood nanoparticles from saline-fed mice, they reproduced lipoteichoic acid-mediated immune functions of the original bacteria, including activation of TLR2 and increased IL-10 expression by dendritic cells. Like the fed bacteria, they also reduced IL-8 induced by TNF in an intestinal epithelial cell line. Though enriched for host neuronal proteins, these isolated nanoparticles also contained proteins and viral (phage) DNA of fed bacterial origin. Our data strongly suggest that oral consumption of live bacteria rapidly leads to circulation of their membrane vesicles and phages and demonstrate a nanoparticulate pathway whereby beneficial bacteria and probiotics may systemically affect their hosts.

Transcriptional and Physiological Analyses to Assess the Effects of a Novel Biostimulant in Tomato

This work aimed to study the effects in tomato (Solanum lycopersicum L.) of foliar applications of a novel calcium-based biostimulant (SOB01) using an omics approach involving transcriptomics and physiological profiling. A calcium-chloride fertilizer (SOB02) was used as a product reference standard. Plants were grown under well-watered (WW) and water stress (WS) conditions in a growth chamber. We firstly compared the transcriptome profile of treated and untreated tomato plants using the software RStudio. Totally, 968 and 1,657 differentially expressed genes (DEGs) (adj-p-value < 0.1 and |log2(fold change)| ≥ 1) were identified after SOB01 and SOB02 leaf treatments, respectively. Expression patterns of 9 DEGs involved in nutrient metabolism and osmotic stress tolerance were validated by real-time quantitative reverse transcription PCR (RT-qPCR) analysis. Principal component analysis (PCA) on RT-qPCR results highlighted that the gene expression profiles after SOB01 treatment in different water regimes were clustering together, suggesting that the expression pattern of the analyzed genes in well water and water stress plants was similar in the presence of SOB01 treatment. Physiological analyses demonstrated that the biostimulant application increased the photosynthetic rate and the chlorophyll content under water deficiency compared to the standard fertilizer and led to a higher yield in terms of fruit dry matter and a reduction in the number of cracked fruits. In conclusion, transcriptome and physiological profiling provided comprehensive information on the biostimulant effects highlighting that SOB01 applications improved the ability of the tomato plants to mitigate the negative effects of water stress.

Towards an understanding of the role of intrinsic protein disorder on plant adaptation to environmental challenges

Intrinsic protein disorder is an interesting structural feature where fully functional proteins lack a three-dimensional structure in solution. In this work, we estimated the relative content of intrinsic protein disorder in 96 plant proteomes including monocots and eudicots. In this analysis, we found variation in the relative abundance of intrinsic protein disorder among these major clades; the relative level of disorder is higher in monocots than eudicots. In turn, there is an inverse relationship between the degree of intrinsic protein disorder and protein length, with smaller proteins being more disordered. The relative abundance of amino acids depends on intrinsic disorder and also varies among clades. Within the nucleus, intrinsically disordered proteins are more abundant than ordered proteins. Intrinsically disordered proteins are specialized in regulatory functions, nucleic acid binding, RNA processing, and in response to environmental stimuli. The implications of this on plants' responses to their environment are discussed.

Low-intensity ultrasound restores long-term potentiation and memory in senescent mice through pleiotropic mechanisms including NMDAR signaling

Advanced physiological aging is associated with impaired cognitive performance and the inability to induce long-term potentiation (LTP), an electrophysiological correlate of memory. Here, we demonstrate in the physiologically aged, senescent mouse brain that scanning ultrasound combined with microbubbles (SUS^+MB), by transiently opening the blood-brain barrier, fully restores LTP induction in the dentate gyrus of the hippocampus. Intriguingly, SUS treatment without microbubbles (SUS^only), i.e., without the uptake of blood-borne factors, proved even more effective, not only restoring LTP, but also ameliorating the spatial learning deficits of the aged mice. This functional improvement is accompanied by an altered milieu of the aged hippocampus, including a lower density of perineuronal nets, increased neurogenesis, and synaptic signaling, which collectively results in improved spatial learning. We therefore conclude that therapeutic ultrasound is a non-invasive, pleiotropic modality that may enhance cognition in elderly humans.

Style transfer with variational autoencoders is a promising approach to RNA-Seq data harmonization and analysis

Motivation

The transcriptomic data are being frequently used in the research of biomarker genes of different diseases and biological states. The most common tasks there are the data harmonization and treatment outcome prediction. Both of them can be addressed via the style transfer approach. Either technical factors or any biological details about the samples which we would like to control (gender, biological state, treatment, etc.) can be used as style components.

Results

The proposed style transfer solution is based on Conditional Variational Autoencoders, Y-Autoencoders and adversarial feature decomposition. To quantitatively measure the quality of the style transfer, neural network classifiers which predict the style and semantics after training on real expression were used. Comparison with several existing style-transfer based approaches shows that proposed model has the highest style prediction accuracy on all considered datasets while having comparable or the best semantics prediction accuracy.

Availability and implementation

https://github.com/NRshka/stvae-source.

Supplementary information

Supplementary data are available at Bioinformatics online.

Comparative Transcriptome Analysis of Iron and Zinc Deficiency in Maize (Zea mays L.)

Globally, one-third of the population is affected by iron (Fe) and zinc (Zn) deficiency, which is severe in developing and underdeveloped countries where cereal-based diets predominate. The genetic biofortification approach is the most sustainable and one of the cost-effective ways to address Fe and Zn malnutrition. Maize is a major source of nutrition in sub-Saharan Africa, South Asia and Latin America. Understanding systems’ biology and the identification of genes involved in Fe and Zn homeostasis facilitate the development of Fe- and Zn-enriched maize. We conducted a genome-wide transcriptome assay in maize inbred SKV616, under –Zn, –Fe and –Fe–Zn stresses. The results revealed the differential expression of several genes related to the mugineic acid pathway, metal transporters, photosynthesis, phytohormone and carbohydrate metabolism. We report here Fe and Zn deficiency-mediated changes in the transcriptome, root length, stomatal conductance, transpiration rate and reduced rate of photosynthesis. Furthermore, the presence of multiple regulatory elements and/or the co-factor nature of Fe and Zn in enzymes indicate their association with the differential expression and opposite regulation of several key gene(s). The differentially expressed candidate genes in the present investigation would help in breeding for Fe and Zn efficient and kernel Fe- and Zn-rich maize cultivars through gene editing, transgenics and molecular breeding.

Nanopore RNA Sequencing Revealed Long Non-Coding and LTR Retrotransposon-Related RNAs Expressed at Early Stages of Triticale SEED Development

The intergenic space of plant genomes encodes many functionally important yet unexplored RNAs. The genomic loci encoding these RNAs are often considered “junk”, DNA as they are frequently associated with repeat-rich regions of the genome. The latter makes the annotations of these loci and the assembly of the corresponding transcripts using short RNAseq reads particularly challenging. Here, using long-read Nanopore direct RNA sequencing, we aimed to identify these “junk” RNA molecules, including long non-coding RNAs (lncRNAs) and transposon-derived transcripts expressed during early stages (10 days post anthesis) of seed development of triticale (AABBRR, 2n = 6x = 42), an interspecific hybrid between wheat and rye. Altogether, we found 796 lncRNAs and 20 LTR retrotransposon-related transcripts (RTE-RNAs) expressed at this stage, with most of them being previously unannotated and located in the intergenic as well as intronic regions. Sequence analysis of the lncRNAs provide evidence for the frequent exonization of Class I (retrotransposons) and class II (DNA transposons) transposon sequences and suggest direct influence of “junk” DNA on the structure and origin of lncRNAs. We show that the expression patterns of lncRNAs and RTE-related transcripts have high stage specificity. In turn, almost half of the lncRNAs located in Genomes A and B have the highest expression levels at 10–30 days post anthesis in wheat. Detailed analysis of the protein-coding potential of the RTE-RNAs showed that 75% of them carry open reading frames (ORFs) for a diverse set of GAG proteins, the main component of virus-like particles of LTR retrotransposons. We further experimentally demonstrated that some RTE-RNAs originate from autonomous LTR retrotransposons with ongoing transposition activity during early stages of triticale seed development. Overall, our results provide a framework for further exploration of the newly discovered lncRNAs and RTE-RNAs in functional and genome-wide association studies in triticale and wheat. Our study also demonstrates that Nanopore direct RNA sequencing is an indispensable tool for the elucidation of lncRNA and retrotransposon transcripts.

The transcriptome of anterior regeneration in earthworm Eudrilus eugeniae

The oligochaete earthworm, Eudrilus eugeniae is capable of regenerating both anterior and posterior segments. The present study focuses on the transcriptome analysis of earthworm E. eugeniae to identify and functionally annotate the key genes supporting the anterior blastema formation and regulating the anterior regeneration of the worm. The Illumina sequencing generated a total of 91,593,182 raw reads which were assembled into 105,193 contigs using CLC genomics workbench. In total, 40,946 contigs were annotated against the NCBI nr and SwissProt database and among them, 15,702 contigs were assigned to 14,575 GO terms. Besides a total of 9389 contigs were mapped to 416 KEGG biological pathways. The RNA-Seq comparison study identified 10,868 differentially expressed genes (DEGs) and of them, 3986 genes were significantly upregulated in the anterior regenerated blastema tissue samples of the worm. The GO enrichment analysis showed angiogenesis and unfolded protein binding as the top enriched functions and the pathway enrichment analysis denoted TCA cycle as the most significantly enriched pathway associated with the upregulated gene dataset of the worm. The identified DEGs and their function and pathway information can be effectively utilized further to interpret the key cellular, genetic and molecular events associated with the regeneration of the worm.

Comparative Transcriptome, Metabolome, and Ionome Analysis of Two Contrasting Common Bean Genotypes in Saline Conditions

Soil salinity is a major abiotic stress factor that limits agricultural productivity worldwide, and this problem is expected to grow in the future. Common bean is an important protein source in developing countries however highly susceptible to salt stress. To understand the underlying mechanism of salt stress responses, transcriptomics, metabolomics, and ion content analysis were performed on both salt-tolerant and susceptible common bean genotypes in saline conditions. Transcriptomics has demonstrated increased photosynthesis in saline conditions for tolerant genotype while the susceptible genotype acted in contrast. Transcriptome also displayed active carbon and amino-acid metabolism for the tolerant genotype. Analysis of metabolites with GC-MS demonstrated the boosted carbohydrate metabolism in the tolerant genotype with increased sugar content as well as better amino-acid metabolism. Accumulation of lysine, valine, and isoleucine in the roots of the susceptible genotype suggested a halted stress response. According to ion content comparison, the tolerant genotype managed to block accumulation of Na⁺ in the leaves while accumulating significantly less Na⁺ in the roots compared to susceptible genotype. K⁺ levels increased in the leaves of both genotype and the roots of the susceptible one but dropped in the roots of the tolerant genotype. Additionally, Zn⁺² and Mn⁺² levels were dropped in the tolerant roots, while Mo⁺² levels were significantly higher in all tissues in both control and saline conditions for tolerant genotype. The results of the presented study have demonstrated the differences in contrasting genotypes and thus provide valuable information on the pivotal molecular mechanisms underlying salt tolerance.

Age related extracellular matrix and interstitial cell phenotype in pulmonary valves

Heart valve disease is a common manifestation of cardiovascular disease and is a significant cause of cardiovascular morbidity and mortality worldwide. The pulmonary valve (PV) is of primary concern because of its involvement in common congenital heart defects, and the PV is usually the site for prosthetic replacement following a Ross operation. Although effects of age on valve matrix components and mechanical properties for aortic and mitral valves have been studied, very little is known about the age-related alterations that occur in the PV. In this study, we isolated PV leaflets from porcine hearts in different age groups (~ 4-6 months, denoted as young versus ~ 2 years, denoted as adult) and studied the effects of age on PV leaflet thickness, extracellular matrix components, and mechanical properties. We also conducted proteomics and RNA sequencing to investigate the global changes of PV leaflets and passage zero PV interstitial cells in their protein and gene levels. We found that the size, thickness, elastic modulus, and ultimate stress in both the radial and circumferential directions and the collagen of PV leaflets increased from young to adult age, while the ultimate strain and amount of glycosaminoglycans decreased when age increased. Young and adult PV had both similar and distinct protein and gene expression patterns that are related to their inherent physiological properties. These findings are important for us to better understand the physiological microenvironments of PV leaflet and valve cells for correctively engineering age-specific heart valve tissues.

Gene expression profiling of corona virus microarray datasets to identify crucial targets in COVID-19 patients

The current outbreak of coronavirus disease (COVID-19) has been affecting millions of people and has caused devastating mortality worldwide. Moreover, it is to be noted that cytokine storm has become an important cause for the rising mortality. However, the efforts for the development of drugs, vaccines and treatment has also been intervened due to poor understanding of host's defense mechanism and also due to the development of cytokine storm against this viral infection. Thus, a deeper understanding of the mechanism behind the immune dysregulation and cytokine storm development might give us clues for the clinical management of the severe cases. Hence, we have implemented differential gene expression analysis together with protein-protein interaction and Gene Ontology (GO) studies with the help of Severe Acute respiratory syndrome coronavirus (SARS-CoV) data sets such as GSE1739 and GSE33267 to give us more knowledge on the host immune response for the pathogenic coronavirus which in turn reduces the mortality. A total of 79 differentially-expressed genes (DEGs) were identified in our data set using the filters such as P-value and log2 fold change values of less than 0.05 and 1.5 respectively. Further, network analysis and GO studies showed that differential expression of two hub genes namely ELANE and LTF which could induce higher levels of pro-inflammatory cytokines in the lungs. We are certain that differential expression of ELANE and LTF results in an excessive inflammatory reaction known as the cytokine storm and ultimately leading to death. Therefore, targeting these key drivers of cytokine storm genes appears to be the potential therapeutic targets for combating the Severe Acute respiratory syndrome coronavirus - 2 (SARS-CoV-2) infection ultimately resulting in reduced mortality. Indeed, this predictive view may open new insights for designing an immune intervention for COVID-19 in the near future resulting in the mitigation of mortality rate.

Evolutionary Dynamics of Transposable Elements Following a Shared Polyploidization Event in the Tribe Andropogoneae

Both polyploidization and transposable element (TE) activity are known to be major drivers of plant genome evolution. Here, we utilize the Zea-Tripsacum clade to investigate TE activity and accumulation after a shared polyploidization event. Comparisons of TE evolutionary dynamics in various Zea and Tripsacum species, in addition to two closely related diploid species, Urelytrum digitatum and Sorghum bicolor, revealed variation in repeat content among all taxa included in the study. The repeat composition of Urelytrum is more similar to that of Zea and Tripsacum compared to Sorghum, despite the similarity in genome size with the latter. Although LTR-retrotransposons were abundant in all species, we observed an expansion of the copia superfamily, specifically in Z. mays and T. dactyloides, species that have adapted to more temperate environments. Additional analyses of the genomic distribution of these retroelements provided evidence of biased insertions near genes involved in various biological processes including plant development, defense, and macromolecule biosynthesis. Specifically, copia insertions in Zea and T. dactyloides were significantly enriched near genes involved in abiotic stress response, suggesting independent evolution post Zea-Tripsacum divergence. The lack of copia insertions near the orthologous genes in S. bicolor suggests that duplicate gene copies generated during polyploidization may offer novel neutral sites for TEs to insert, thereby providing an avenue for subfunctionalization via TE insertional mutagenesis.

ELONGATED HYPOCOTYL5 Negatively Regulates DECREASE WAX BIOSYNTHESIS to Increase Survival during UV-B Stress

Understanding how the distinct cell types of the shoot apical meristem (SAM) withstand ultraviolet radiation (UVR) stress can improve cultivation of plants in high-UVR environments. Here, we show that UV-B irradiation selectively kills epidermal and niche cells in the shoot apex. Plants harboring a mutation in DECREASE WAX BIOSYNTHESIS (DEWAX) are tolerant to UV-B. Our data show that DEWAX negatively regulates genes involved in anthocyanin biosynthesis. ELONGATED HYPOCOTYL5 (HY5) binds to the DEWAX promoter elements and represses its expression to promote the anthocyanin biosynthesis. The HY5-DEWAX regulatory network regulates anthocyanin content in Arabidopsis (Arabidopsis thaliana) and influences the survivability of plants under UV-B irradiation stress. Our cell sorting-based study of the epidermal cell layer transcriptome confirms that core UV-B stress signaling pathway genes are conserved and upregulated in response to UV-B irradiation of the SAM. Furthermore, we show that UV-B induces genes involved in shoot development and organ patterning. We propose that the HY5-DEWAX regulatory relationship is conserved; however, changes in the expression levels of these genes can determine anthocyanin content in planta and, hence, fitness under UV-B irradiation stress.

Contrasting signatures of genomic divergence during sympatric speciation

The transition from 'well-marked varieties' of a single species into 'well-defined species'-especially in the absence of geographic barriers to gene flow (sympatric speciation)-has puzzled evolutionary biologists ever since Darwin1,2. Gene flow counteracts the buildup of genome-wide differentiation, which is a hallmark of speciation and increases the likelihood of the evolution of irreversible reproductive barriers (incompatibilities) that complete the speciation process3. Theory predicts that the genetic architecture of divergently selected traits can influence whether sympatric speciation occurs4, but empirical tests of this theory are scant because comprehensive data are difficult to collect and synthesize across species, owing to their unique biologies and evolutionary histories5. Here, within a young species complex of neotropical cichlid fishes (Amphilophus spp.), we analysed genomic divergence among populations and species. By generating a new genome assembly and re-sequencing 453 genomes, we uncovered the genetic architecture of traits that have been suggested to be important for divergence. Species that differ in monogenic or oligogenic traits that affect ecological performance and/or mate choice show remarkably localized genomic differentiation. By contrast, differentiation among species that have diverged in polygenic traits is genomically widespread and much higher overall, consistent with the evolution of effective and stable genome-wide barriers to gene flow. Thus, we conclude that simple trait architectures are not always as conducive to speciation with gene flow as previously suggested, whereas polygenic architectures can promote rapid and stable speciation in sympatry.

Genome-wide identification of novel long non-coding RNAs and their possible roles in hypoxic zebrafish brain

Long non-coding RNAs (lncRNAs) are the master regulators of numerous biological processes. Hypoxia causes oxidative stress with severe and detrimental effects on brain function and acts as a critical initiating factor in the pathogenesis of Alzheimer's disease (AD). From the RNA-Seq in the forebrain (Fb), midbrain (Mb), and hindbrain (Hb) regions of hypoxic and normoxic zebrafish, we identified novel lncRNAs, whose potential cis targets showed involvement in neuronal development and differentiation pathways. Under hypoxia, several lncRNAs and mRNAs were differentially expressed. Co-expression studies indicated that the Fb and Hb regions' potential lncRNA target genes were involved in the AD pathogenesis. In contrast, those in Mb (cry1b, per1a, cipca) was responsible for regulating circadian rhythm. We identified specific lncRNAs present in the syntenic regions between zebrafish and humans, possibly functionally conserved. We thus identified several conserved lncRNAs as the probable regulators of AD genes (adrb3b, cav1, stat3, bace2, apoeb, psen1, s100b).

Epigenome-wide association scan identifies methylation sites associated with HIV infection

Aim: To investigate the role of epigenetics in HIV pathophysiology. Materials & methods: We conducted an epigenome-wide association scan on HIV infection status among people who inject drugs in the AIDS Linked to the IntraVenous Experience study with primary (n = 397) and validation samples (n = 390). DNA methylation from blood was measured by the Illumina EPIC BeadChip. We controlled for cell type heterogeneity by HIV status. Results: HIV infection status was associated (p < 10^-8) with DNA methylation at 49 CpG sites. Sites were enriched in response to virus, interferon signaling pathway, etc. Among these sites, discovery and validation t-statistics were highly correlated (r = 0.96). Conclusion: In a cohort of people who inject drugs, HIV status was associated with differential DNA methylation at biologically meaningful sites.

Differential proteome analysis of the leaves of lead hyperaccumulator, Rhoeo discolor (L. Her.) Hance

The present study investigated Rhoeo discolor (L. Her.) Hance for its ability to accumulate Pb, which is of relevance to phytoremediation applications. Based on this analysis, plants were found to accumulate greater than 10 mg/g (0.1%) of dry weight Pb in the shoots, which classifies the plant a Pb hyperaccumulator. Further, changes in the leaf proteome profiles in response to Pb stress were investigated. Wild-type plants were subjected to a high concentration of Pb(NO3 )2 , and the levels of Pb that accumulated in different plant tissues were determined using atomic absorption spectrophotometry. Using 2D-difference gel electrophoresis, 181 protein spots were detected to be differentially abundant in response to Pb stress and selected spots exhibiting the strongest differential abundance suggested an impairment of the photosynthetic apparatus of the plant under Pb stress. Subsequently, a more extensive, proteome wide analysis utilizing label-free quantitation further identified a predominant decrease in protein levels and a significant effect on the nuclear proteome, as well as photosynthesis, carbon fixation and metabolism, providing insight into the Pb tolerance of this system in a potential phytoremediation context.

Regulation of Pain Genes-Capsaicin vs Resiniferatoxin: Reassessment of Transcriptomic Data

Emerging evidence has shown a strong association between neuropathic pain and chronic diseases. In recent years, the treatment of neuropathic pain has attracted more attention. Natural products, such as capsaicin and resiniferatoxin, have been well utilized to treat this disease. In this study, we aim to compare the regulatory effects of capsaicin and resiniferatoxin on pain-related genes as well as on genes with no direct association with pain. Public transcriptomic and microarray data on gene expression in the dorsal root ganglia and genes associated with TRPV1 (+) neurons were obtained from the GEO database and then analyzed. Differentially expressed genes were selected for further functional analysis, including pathway enrichment, protein-protein interaction, and regulatory network analysis. Pain-associated genes were extracted with the reference of two pain gene databases and the effects of these two natural drugs on the pain-associated genes were measured. The results of our research indicate that as compared to capsaicin, resiniferatoxin (RTX) regulates more non pain-associated genes and has a negative impact on beneficial genes (off-targets) which are supposed to alleviate nociception and hypersensitivity by themselves. So, based on this study, we may conclude that capsaicin may be less potent when compared to RTX, but it will elicit considerably less adverse effects too. Thereby confirming that capsaicin could be used for the efficient alleviation of neuropathic pain with possibly fewer side effects.

Transcriptome Analysis Identifies SenZfp536, a Sense LncRNA that Suppresses Self-renewal of Cortical Neural Progenitors

Long non-coding RNAs (lncRNAs) regulate transcription to control development and homeostasis in a variety of tissues and organs. However, their roles in the development of the cerebral cortex have not been well elucidated. Here, a bioinformatics pipeline was applied to delineate the dynamic expression and potential cis-regulating effects of mouse lncRNAs using transcriptome data from 8 embryonic time points and sub-regions of the developing cerebral cortex. We further characterized a sense lncRNA, SenZfp536, which is transcribed downstream of and partially overlaps with the protein-coding gene Zfp536. Both SenZfp536 and Zfp536 were predominantly expressed in the proliferative zone of the developing cortex. Zfp536 was cis-regulated by SenZfp536, which facilitates looping between the promoter of Zfp536 and the genomic region that transcribes SenZfp536. Surprisingly, knocking down or activating the expression of SenZfp536 increased or compromised the proliferation of cortical neural progenitor cells (NPCs), respectively. Finally, overexpressing Zfp536 in cortical NPCs reversed the enhanced proliferation of cortical NPCs caused by SenZfp536 knockdown. The study deepens our understanding of how lncRNAs regulate the propagation of cortical NPCs through cis-regulatory mechanisms.

Arabidopsis Response to Inhibitor of Cytokinin Degradation INCYDE: Modulations of Cytokinin Signaling and Plant Proteome

Cytokinins are multifaceted plant hormones that play crucial roles in plant interactions with the environment. Modulations in cytokinin metabolism and signaling have been successfully used for elevating plant tolerance to biotic and abiotic stressors. Here, we analyzed Arabidopsis thaliana response to INhibitor of CYtokinin DEgradation (INCYDE), a potent inhibitor of cytokinin dehydrogenase. We found that at low nanomolar concentration, the effect of INCYCE on seedling growth and development was not significantly different from that of trans-Zeatin treatment. However, an alteration in the spatial distribution of cytokinin signaling was found at low micromolar concentrations, and proteomics analysis revealed a significant impact on the molecular level. An in-depth proteome analysis of an early (24 h) response and a dose-dependent response after 168 h highlighted the effects on primary and secondary metabolism, including alterations in ribosomal subunits, RNA metabolism, modulations of proteins associated with chromatin, and the flavonoid and phenylpropanoid biosynthetic pathway. The observed attenuation in stress-response mechanisms, including abscisic acid signaling and the metabolism of jasmonates, could explain previously reported positive effects of INCYDE under mild stress conditions.

Duchenne muscular dystrophy (DMD) cardiomyocyte-secreted exosomes promote the pathogenesis of DMD-associated cardiomyopathy

Cardiomyopathy is a leading cause of early mortality in Duchenne muscular dystrophy (DMD). There is a need to gain a better understanding of the molecular pathogenesis for the development effective therapies. Exosomes (exo) are secreted vesicles and exert effects via their RNA, lipid and protein cargo. The role of exosomes in disease pathology is unknown. Exosomes derived from stem cells have demonstrated cardioprotection in the murine DMD heart. However, it is unknown how the disease status of the donor cell type influences exosome function. Here, we sought to determine the phenotypic responses of DMD cardiomyocytes (DMD-iCMs) after long-term exposure to DMD cardiac exosomes (DMD-exo). DMD-iCMs were vulnerable to stress, evidenced by production of reactive oxygen species, the mitochondrial membrane potential and cell death levels. Long-term exposure to non-affected exosomes (N-exo) was protective. By contrast, long-term exposure to DMD-exo was not protective, and the response to stress improved with inhibition of DMD-exo secretion in vitro and in vivo The microRNA (miR) cargo, but not exosome surface peptides, was implicated in the pathological effects of DMD-exo. Exosomal surface profiling revealed N-exo peptides associated with PI3K-Akt signaling. Transcriptomic profiling identified unique changes with exposure to either N- or DMD-exo. Furthermore, DMD-exo miR cargo regulated injurious pathways, including p53 and TGF-beta. The findings reveal changes in exosomal cargo between healthy and diseased states, resulting in adverse outcomes. Here, DMD-exo contained miR changes, which promoted the vulnerability of DMD-iCMs to stress. Identification of these molecular changes in exosome cargo and effectual phenotypes might shed new light on processes underlying DMD cardiomyopathy.This article has an associated First Person interview with the first author of the paper.

Genomic landscape and clonal architecture of mouse oral squamous cell carcinomas dictate tumour ecology

To establish whether 4-nitroquinoline N-oxide-induced carcinogenesis mirrors the heterogeneity of human oral squamous cell carcinoma (OSCC), we have performed genomic analysis of mouse tongue lesions. The mutational signatures of human and mouse OSCC overlap extensively. Mutational burden is higher in moderate dysplasias and invasive SCCs than in hyperplasias and mild dysplasias, although mutations in p53, Notch1 and Fat1 occur in early lesions. Laminin-α3 mutations are associated with tumour invasiveness and Notch1 mutant tumours have an increased immune infiltrate. Computational modelling of clonal dynamics indicates that high genetic heterogeneity may be a feature of those mild dysplasias that are likely to progress to more aggressive tumours. These studies provide a foundation for exploring OSCC evolution, heterogeneity and progression.

Deep learning approach for predicting functional Z-DNA regions using omics data

Computational methods to predict Z-DNA regions are in high demand to understand the functional role of Z-DNA. The previous state-of-the-art method Z-Hunt is based on statistical mechanical and energy considerations about B- to Z-DNA transition using sequence information. Z-DNA CHiP-seq experiment results showed little overlap with Z-Hunt predictions implying that sequence information only is not sufficient to explain emergence of Z-DNA at different genomic locations. Adding epigenetic and other functional genomic mark-ups to DNA sequence level can help revealing the functional Z-DNA sites. Here we take advantage of the deep learning approach that can analyze and extract information from large volumes of molecular biology data. We developed a machine learning approach DeepZ that aggregates information from genome-wide maps of epigenetic markers, transcription factor and RNA polymerase binding sites, and chromosome accessibility maps. With the developed model we not only verify the experimental Z-DNA predictions, but also generate the whole-genome annotation, introducing new possible Z-DNA regions, which have not yet been found in experiments and can be of interest to the researchers from various fields.

DNA methylation in canine brains is related to domestication and dog-breed formation

Epigenetic factors such as DNA methylation act as mediators in the interaction between genome and environment. Variation in the epigenome can both affect phenotype and be inherited, and epigenetics has been suggested to be an important factor in the evolutionary process. During domestication, dogs have evolved an unprecedented between-breed variation in morphology and behavior in an evolutionary short period. In the present study, we explore DNA methylation differences in brain, the most relevant tissue with respect to behavior, between wolf and dog breeds. We optimized a combined method of genotype-by-sequencing (GBS) and methylated DNA immunoprecipitation (MeDIP) for its application in canines. Genomic DNA from the frontal cortex of 38 dogs of 8 breeds and three wolves was used. GBS and GBS-MeDIP libraries were prepared and sequenced on Illuma HiSeq2500 platform. The reduced sample represented 1.18 ± 0.4% of the total dog genome (2,4 billion BP), while the GBS-MeDIP covered 11,250,788 ± 4,042,106 unique base pairs. We find substantial DNA methylation differences between wolf and dog and between the dog breeds. The methylation profiles of the different groups imply that epigenetic factors may have been important in the speciation from dog to wolf, but also in the divergence of different dog breeds. Specifically, we highlight methylation differences in genes related to behavior and morphology. We hypothesize that these differences are involved in the phenotypic variation found among dogs, whereas future studies will have to find the specific mechanisms. Our results not only add an intriguing new dimension to dog breeding but are also useful to further understanding of epigenetic involvement.

A Screen for PKN3 Substrates Reveals an Activating Phosphorylation of ARHGAP18

Protein kinase N3 (PKN3) is a serine/threonine kinase implicated in tumor progression of multiple cancer types, however, its substrates and effector proteins still remain largely understudied. In the present work we aimed to identify novel PKN3 substrates in a phosphoproteomic screen using analog sensitive PKN3. Among the identified putative substrates we selected ARHGAP18, a protein from RhoGAP family, for validation of the screen and further study. We confirmed that PKN3 can phosphorylate ARHGAP18 in vitro and we also characterized the interaction of the two proteins, which is mediated via the N-terminal part of ARHGAP18. We present strong evidence that PKN3-ARHGAP18 interaction is increased upon ARHGAP18 phosphorylation and that the phosphorylation of ARHGAP18 by PKN3 enhances its GAP domain activity and contributes to negative regulation of active RhoA. Taken together, we identified new set of potential PKN3 substrates and revealed a new negative feedback regulatory mechanism of Rho signaling mediated by PKN3-induced ARHGAP18 activation.

During bacteremia, Pseudomonas aeruginosa PAO1 adapts by altering the expression of numerous virulence genes including those involved in quorum sensing

Pseudomonas aeruginosa is a Gram-negative opportunistic pathogen that produces numerous virulence factors and causes serious infections in trauma patients and patients with severe burns. We previously showed that the growth of P. aeruginosa in blood from severely burned or trauma patients altered the expression of numerous genes. However, the specific influence of whole blood from healthy volunteers on P. aeruginosa gene expression is not known. Transcriptome analysis of P. aeruginosa grown for 4 h in blood from healthy volunteers compared to that when grown in laboratory medium revealed that the expression of 1085 genes was significantly altered. Quorum sensing (QS), QS-related, and pyochelin synthesis genes were downregulated, while genes of the type III secretion system and those for pyoverdine synthesis were upregulated. The observed effect on the QS and QS-related genes was shown to reside within serum fraction: growth of PAO1 in the presence of 10% human serum from healthy volunteers significantly reduced the expression of QS and QS-regulated genes at 2 and 4 h of growth but significantly enhanced their expression at 8 h. Additionally, the production of QS-regulated virulence factors, including LasA and pyocyanin, was also influenced by the presence of human serum. Serum fractionation experiments revealed that part of the observed effect resides within the serum fraction containing <10-kDa proteins. Growth in serum reduced the production of many PAO1 outer membrane proteins but enhanced the production of others including OprF, a protein previously shown to play a role in the regulation of QS gene expression. These results suggest that factor(s) within human serum: 1) impact P. aeruginosa pathogenesis by influencing the expression of different genes; 2) differentially regulate the expression of QS and QS-related genes in a growth phase- or time-dependent mechanism; and 3) manipulate the production of P. aeruginosa outer membrane proteins.

Host-symbiont coevolution, cryptic structure, and bleaching susceptibility, in a coral species complex (Scleractinia; Poritidae)

The 'species' is a key concept for conservation and evolutionary biology, yet the lines between population and species-level variation are often blurred, especially for corals. The 'Porites lobata species complex' consists of branching and mounding corals that form reefs across the Pacific. We used reduced representation meta-genomic sequencing to examine genetic relationships within this species complex and to identify candidate loci associated with colony morphology, cryptic genetic structure, and apparent bleaching susceptibility. We compared existing Porites data with bleached and unbleached colonies of the branching coral P. compressa collected in Kāne'ohe Bay Hawai'i during the 2015 coral bleaching event. Loci that mapped to coral, symbiont, and microbial references revealed genetic structure consistent with recent host-symbiont co-evolution. Cryptic genetic clades were resolved that previous work has associated with distance from shore, but no genetic structure was associated with bleaching. We identified many candidate loci associated with morphospecies, including candidate host and symbiont loci with fixed differences between branching and mounding corals. We also found many loci associated with cryptic genetic structure, yet relatively few loci associated with bleaching. Recent host-symbiont co-evolution and rapid diversification suggests that variation and therefore the capacity of these corals to adapt may be underappreciated.

Epigenetic Analysis through MSAP-NGS Coupled Technology: The Case Study of White Poplar Monoclonal Populations/Stands

Over the last several decades, several lines of evidence have shown that epigenetic modifications modulate phenotype and mediate an organism’s response to environmental stimuli. Plant DNA is normally highly methylated, although notable differences exist between species. Many biomolecular techniques based on PCR have been developed to analyse DNA methylation status, however a qualitative leap was made with the advent of next-generation sequencing (NGS). In the case of large, repetitive, or not-yet-sequenced genomes characterised by a high level of DNA methylation, the NGS analysis of bisulphite pre-treated DNA is expensive and time consuming, and moreover, in some cases data analysis is a major challenge. Methylation-sensitive amplification polymorphism (MSAP) analysis is a highly effective method to study DNA methylation. The method is based on the comparison of double DNA digestion profiles (EcoRI-HpaII and EcoRI-MspI) to reveal methylation pattern variations. These are often attributable to pedoclimatic and stress conditions which affect all organisms during their lifetime. In our study, five white poplar (Populus alba L.) specimens were collected from different monoclonal stands in the Maltese archipelago, and their DNA was processed by means of an innovative approach where MSAP analysis was followed by NGS. This allowed us to identify genes that were differentially methylated among the different specimens and link them to specific biochemical pathways. Many differentially methylated genes were found to encode transfer RNAs (tRNAs) related to photosynthesis or light reaction pathways. Our results clearly demonstrate that this combinatorial method is suitable for epigenetic studies of unsequenced genomes like P. alba (at the time of study), and to identify epigenetic variations related to stress, probably caused by different and changing pedoclimatic conditions, to which the poplar stands have been exposed.

Sexually dimorphic effects of forkhead box a2 (FOXA2) and uterine glands on decidualization and fetoplacental development

Glands of the uterus are essential for pregnancy establishment. Forkhead box A2 (FOXA2) is expressed specifically in the glands of the uterus and a critical regulator of glandular epithelium (GE) differentiation, development, and function. Mice with a conditional deletion of FOXA2 in the adult uterus, created using the lactotransferrin iCre (Ltf-iCre) model, have a morphologically normal uterus with glands, but lack FOXA2-dependent GE-expressed genes, such as leukemia inhibitory factor (LIF). Adult FOXA2 conditional knockout (cKO; Ltf ^iCre/+ Foxa2 ^f/f ) mice are infertile due to defective embryo implantation arising from a lack of LIF, a critical implantation factor of uterine gland origin. However, intraperitoneal injections of LIF can initiate embryo implantation in the uterus of adult FOXA2 cKO mice with pregnancies maintained to term. Here, we tested the hypothesis that FOXA2-regulated genes in the uterine glands impact development of the decidua, placenta, and fetus. On gestational day 8.5, the antimesometrial and mesometrial decidua transcriptome was noticeably altered in LIF-replaced FOXA2 cKO mice. Viable fetuses were reduced in FOXA2 cKO mice on gestational days 12.5 and 17.5. Sex-dependent differences in fetal weight, placenta histoarchitecture, and the placenta and metrial gland transcriptome were observed between control and FOXA2 cKO mice. The transcriptome of the placenta with a female fetus was considerably more altered than the placenta with a male fetus in FOXA2 cKO dams. These studies reveal previously unrecognized sexually dimorphic effects of FOXA2 and uterine glands on fetoplacental development with potential impacts on offspring health into adulthood.

CeRNA Network Analysis Representing Characteristics of Different Tumor Environments Based on 1p/19q Codeletion in Oligodendrogliomas

Simple Summary

Oligodendroglioma (OD) is a subtype of glioma occurring in the central nervous system. The 1p/19q codeletion is a prognostic marker of OD with an isocitrate dehydrogenase (IDH) mutation and is associated with a clinically favorable overall survival (OS). The long non-coding RNAs (lncRNAs) protects the mRNA from degradation by binding with the same miRNA by acting as a competitive endogenous RNA (ceRNA). Recently, although there is an increasing interest in lncRNAs on glioma studies, however, studies regarding their effects on OD and the 1p/19q codeletion remain limited. In our study, we performed in silico analyses using low-grade gliomas from datasets obtained from The Cancer Genome Atlas to investigate the effects of ceRNA with 1p/19q codeletion on ODs. We constructed 16 coding RNA–miRNA–lncRNA networks and the ceRNA network participated in ion channel activity, insulin secretion, and collagen network and extracellular matrix (ECM) changes. In conclusion, our results can provide insights into the possibility in the different tumor microenvironments and OS following 1p/19q codeletion through changes in the ceRNA network.

Abstract

Oligodendroglioma (OD) is a subtype of glioma occurring in the central nervous system. The 1p/19q codeletion is a prognostic marker of OD with an isocitrate dehydrogenase (IDH) mutation and is associated with a clinically favorable overall survival (OS); however, the exact underlying mechanism remains unclear. Long non-coding RNAs (lncRNAs) have recently been suggested to regulate carcinogenesis and prognosis in cancer patients. Here, we performed in silico analyses using low-grade gliomas from datasets obtained from The Cancer Genome Atlas to investigate the effects of ceRNA with 1p/19q codeletion on ODs. Thus, we selected modules of differentially expressed genes that were closely related to 1p/19q codeletion traits using weighted gene co-expression network analysis and constructed 16 coding RNA–miRNA–lncRNA networks. The ceRNA network participated in ion channel activity, insulin secretion, and collagen network and extracellular matrix (ECM) changes. In conclusion, ceRNAs with a 1p/19q codeletion can create different tumor microenvironments via potassium ion channels and ECM composition changes; furthermore, differences in OS may occur. Moreover, if extrapolated to gliomas, our results can provide insights into the consequences of identical gene expression, indicating the possibility of tracking different biological processes in different subtypes of glioma.

Markers of Angiogenesis, Lymphangiogenesis, and Epithelial-Mesenchymal Transition (Plasticity) in CIN and Early Invasive Carcinoma of the Cervix: Exploring Putative Molecular Mechanisms Involved in Early Tumor Invasion

The establishment of a proangiogenic phenotype and epithelial-to-mesenchymal transition (EMT) are considered as critical events that promote the induction of invasive growth in epithelial tumors, and stimulation of lymphangiogenesis is believed to confer the capacity for early dissemination to cancer cells. Recent research has revealed substantial interdependence between these processes at the molecular level as they rely on common signaling networks. Of great interest are the molecular mechanisms of (lymph-)angiogenesis and EMT associated with the earliest stages of transition from intraepithelial development to invasive growth, as they could provide the source of potentially valuable tools for targeting tumor metastasis. However, in the case of early-stage cervical cancer, the players of (lymph-)angiogenesis and EMT processes still remain substantially uncharacterized. In this study, we used RNA sequencing to compare transcriptomes of HPV(+) preinvasive neoplastic lesions and early-stage invasive carcinoma of the cervix and to identify (lymph-)angiogenesis- and EMT-related genes and pathways that may underlie early acquisition of invasive phenotype and metastatic properties by cervical cancer cells. Second, we applied flow cytometric analysis to evaluate the expression of three key lymphangiogenesis/EMT markers (VEGFR3, MET, and SLUG) in epithelial cells derived from enzymatically treated tissue specimens. Overall, among 201 differentially expressed genes, a considerable number of (lymph-)angiogenesis and EMT regulatory factors were identified, including genes encoding cytokines, growth factor receptors, transcription factors, and adhesion molecules. Pathway analysis confirmed enrichment for angiogenesis, epithelial differentiation, and cell guidance pathways at transition from intraepithelial neoplasia to invasive carcinoma and suggested immune-regulatory/inflammatory pathways to be implicated in initiation of invasive growth of cervical cancer. Flow cytometry showed cell phenotype-specific expression pattern for VEGFR3, MET, and SLUG and revealed correlation with the amount of tumor-infiltrating lymphocytes at the early stages of cervical cancer progression. Taken together, these results extend our understanding of driving forces of angiogenesis and metastasis in HPV-associated cervical cancer and may be useful for developing new treatments.

Sustained Inflammatory Signalling through Stat1/Stat2/IRF9 Is Associated with Amoeboid Phenotype of Melanoma Cells

Simple Summary

Treatment of metastatic cancer is complicated by the ability of cancer cells to utilize various invasion modes when spreading through the body. Here, we studied the transition of melanoma cells between the round, amoeboid and elongated, mesenchymal invasion modes. Our results show that inflammatory signalling, which is commonly upregulated in the tumour microenvironment, is associated with the amoeboid phenotype of cancer cells. Treatment of melanoma cells with interferon beta promotes the amoeboid invasion modes and individual invasion. This suggests that inflammation associated signalling contributes to cancer cell invasion plasticity.

Abstract

The invasive behaviour of cancer cells underlies metastatic dissemination; however, due to the large plasticity of invasion modes, it is challenging to target. It is now widely accepted that various secreted cytokines modulate the tumour microenvironment and pro-inflammatory signalling can promote tumour progression. Here, we report that cells after mesenchymal–amoeboid transition show the increased expression of genes associated with the type I interferon response. Moreover, the sustained activation of type I interferon signalling in response to IFNβ mediated by the Stat1/Stat2/IRF9 complex enhances the round amoeboid phenotype in melanoma cells, whereas its downregulation by various approaches promotes the mesenchymal invasive phenotype. Overall, we demonstrate that interferon signalling is associated with the amoeboid phenotype of cancer cells and suggest a novel role of IFNβ in promoting cancer invasion plasticity, aside from its known role as a tumour suppressor.

The Acidic Stress Response of the Intracellular Pathogen Brucella melitensis: New Insights from a Comparative, Genome-Wide Transcriptome Analysis

The intracellular pathogenic bacteria belonging to the genus Brucella must cope with acidic stress as they penetrate the host via the gastrointestinal route, and again during the initial stages of intracellular infection. A transcription-level regulation has been proposed to explain this but the specific molecular mechanisms are yet to be determined. We recently reported a comparative transcriptomic analysis of the attenuated vaccine Brucella melitensis strain Rev.1 against the virulent strain 16M in cultures grown under either neutral or acidic conditions. Here, we re-analyze the RNA-seq data of 16M from our previous study and compare it to published transcriptomic data of this strain from both an in cellulo and an in vivo model. We identify 588 genes that are exclusively differentially expressed in 16M grown under acidic versus neutral pH conditions, including 286 upregulated genes and 302 downregulated genes that are not differentially expressed in either the in cellulo or the in vivo model. Of these, we highlight 13 key genes that are known to be associated with a bacterial response to acidic stress and, in our study, were highly upregulated under acidic conditions. These genes provide new molecular insights into the mechanisms underlying the acid-resistance of Brucella within its host.

Enzymatic RNA Biotinylation for Affinity Purification and Identification of RNA-Protein Interactions

Throughout their cellular lifetime, RNA transcripts are bound to proteins, playing crucial roles in RNA metabolism, trafficking, and function. Despite the importance of these interactions, identifying the proteins that interact with an RNA of interest in mammalian cells represents a major challenge in RNA biology. Leveraging the ability to site-specifically and covalently label an RNA of interest using E. coli tRNA guanine transglycosylase and an unnatural nucleobase substrate, we establish the identification of RNA-protein interactions and the selective enrichment of cellular RNA in mammalian systems. We demonstrate the utility of this approach through the identification of known binding partners of 7SK snRNA via mass spectrometry. Through a minimal 4-nucleotide mutation of the long noncoding RNA HOTAIR, enzymatic biotinylation enables identification of putative HOTAIR binding partners in MCF7 breast cancer cells that suggest new potential pathways for oncogenic function. Furthermore, using RNA sequencing and qPCR, we establish that an engineered enzyme variant achieves high levels of labeling selectivity against the human transcriptome allowing for 145-fold enrichment of cellular RNA directly from mammalian cell lysates. The flexibility and breadth of this approach suggests that this system could be routinely applied to the functional characterization of RNA, greatly expanding the toolbox available for studying mammalian RNA biology.

Formation of NPR1 Condensates Promotes Cell Survival during the Plant Immune Response

In plants, pathogen effector-triggered immunity (ETI) often leads to programmed cell death, which is restricted by NPR1, an activator of systemic acquired resistance. However, the biochemical activities of NPR1 enabling it to promote defense and restrict cell death remain unclear. Here we show that NPR1 promotes cell survival by targeting substrates for ubiquitination and degradation through formation of salicylic acid-induced NPR1 condensates (SINCs). SINCs are enriched with stress response proteins, including nucleotide-binding leucine-rich repeat immune receptors, oxidative and DNA damage response proteins, and protein quality control machineries. Transition of NPR1 into condensates is required for formation of the NPR1-Cullin 3 E3 ligase complex to ubiquitinate SINC-localized substrates, such as EDS1 and specific WRKY transcription factors, and promote cell survival during ETI. Our analysis of SINCs suggests that NPR1 is centrally integrated into the cell death or survival decisions in plant immunity by modulating multiple stress-responsive processes in this quasi-organelle.

Distinct and overlapping roles of STAG1 and STAG2 in cohesin localization and gene expression in embryonic stem cells

BACKGROUND

The three-dimensional organization of the genome in the nucleus plays an integral role in many biological processes, including gene expression. The genome is folded into DNA loops that bring together distal regulatory elements and genes. Cohesin, a ring-shaped protein complex, is a major player in the formation of DNA loops. Cohesin is composed of a core trimer and one of two variant STAG subunits, STAG1 or STAG2. It is not understood whether variant STAG proteins give rise to cohesin complexes with distinct functions. Recent studies have begun to characterize the roles of STAG1 and STAG2, with partially contradictory results.

RESULTS

Here, we generate stable single-knockout embryonic stem cell lines to investigate the individual contributions of STAG1 and STAG2 in regulating cohesin chromosomal localization and function. We report both overlapping roles for STAG1 and STAG2 in cohesin localization and somewhat distinct roles in gene expression. STAG1 and STAG2 occupy the same sites across the genome, yet do not exist together in a higher order complex. Despite their shared localization, STAG1 and STAG2 have both distinct and redundant effects on gene expression. Loss of both STAG1 and STAG2 causes widespread transcriptome dysregulation, altered cohesin DNA occupancy, and reduced cell proliferation.

CONCLUSIONS

Together, this work reveals the requirement of at least one STAG protein for proper cohesin function. STAG1 and STAG2 have independent roles in cohesin localization and both overlapping and distinct roles in gene expression. The roles of STAG1 and STAG2 in mouse embryonic stem cells may be somewhat different than in other cell types, due to their relative expression levels. These results advance our understanding of the link between mammalian genome organization and gene expression during development and disease contexts.

Sex-Specific Transcriptome Differences in Human Adipose Mesenchymal Stem Cells

In humans, sexual dimorphism can manifest in many ways and it is widely studied in several knowledge fields. It is increasing the evidence that also cells differ according to sex, a correlation still little studied and poorly considered when cells are used in scientific research. Specifically, our interest is on the sex-related dimorphism on the human mesenchymal stem cells (hMSCs) transcriptome. A systematic meta-analysis of hMSC microarrays was performed by using the Transcriptome Mapper (TRAM) software. This bioinformatic tool was used to integrate and normalize datasets from multiple sources and allowed us to highlight chromosomal segments and genes differently expressed in hMSCs derived from adipose tissue (hADSCs) of male and female donors. Chromosomal segments and differentially expressed genes in male and female hADSCs resulted to be related to several processes as inflammation, adipogenic and neurogenic differentiation and cell communication. Obtained results lead us to hypothesize that the donor sex of hADSCs is a variable influencing a wide range of stem cell biologic processes. We believe that it should be considered in biologic research and stem cell therapy.

Identification of SNPs and InDels associated with berry size in table grapes integrating genetic and transcriptomic approaches

BACKGROUND

Berry size is considered as one of the main selection criteria in table grapes breeding programs, due to the consumer preferences. However, berry size is a complex quantitive trait under polygenic control, and its genetic determination of berry weight is not yet fully understood. The aim of this work was to perform marker discovery using a transcriptomic approach, in order to identify and characterize SNP and InDel markers associated with berry size in table grapes. We used an integrative analysis based on RNA-Seq, SNP/InDel search and validation on table grape segregants and varieties with different genetic backgrounds.

RESULTS

Thirty SNPs and eight InDels were identified using a transcriptomic approach (RNA-Seq). These markers were selected from SNP/InDel found among segregants from a Ruby x Sultanina population with contrasting phenotypes for berry size. The set of 38 SNP and InDel markers was distributed in eight chromosomes. Genotype-phenotype association analyses were performed using a set of 13 RxS segregants and 41 table grapes varieties with different genetic backgrounds during three seasons. The results showed several degrees of association of these markers with berry size (10.2 to 30.7%) as other berry-related traits such as length and width. The co-localization of SNP and /or InDel markers and previously reported QTLs and candidate genes associated with berry size were analysed.

CONCLUSIONS

We identified a set of informative and transferable SNP and InDel markers associated with berry size. Our results suggest the suitability of SNPs and InDels as candidate markers for berry weight in seedless table grape breeding. The identification of genomic regions associated with berry weight in chromosomes 8, 15 and 17 was achieved with supporting evidence derived from a transcriptome experiment focused on SNP/InDel search, as well as from a QTL-linkage mapping approach. New regions possibly associated with berry weight in chromosomes 3, 6, 9 and 14 were identified.

Defining the transcriptional responses of Aspergillus nidulans to cation/alkaline pH stress and the role of the transcription factor SltA

Fungi have developed the ability to overcome extreme growth conditions and thrive in hostile environments. The model fungus Aspergillus nidulans tolerates, for example, ambient alkalinity up to pH 10 or molar concentrations of multiple cations. The ability to grow under alkaline pH or saline stress depends on the effective function of at least three regulatory pathways mediated by the zinc-finger transcription factor PacC, which mediates the ambient pH regulatory pathway, the calcineurin-dependent CrzA and the cation homeostasis responsive factor SltA. Using RNA sequencing, we determined the effect of external pH alkalinization or sodium stress on gene expression. The data show that each condition triggers transcriptional responses with a low degree of overlap. By sequencing the transcriptomes of the null mutant, the role of SltA in the above-mentioned homeostasis mechanisms was also studied. The results show that the transcriptional role of SltA is wider than initially expected and implies, for example, the positive control of the PacC-dependent ambient pH regulatory pathway. Overall, our data strongly suggest that the stress response pathways in fungi include some common but mostly exclusive constituents, and that there is a hierarchical relationship among the main regulators of stress response, with SltA controlling pacC expression, at least in A. nidulans.

Mammary tissue transcriptomic analysis for construction of integrated regulatory networks involved in lactogenesis of Ovis aries

The mammary gland experiences vast changes between the onset of lactation and pregnancy. This remodeling involves different functions such as lactation that is controlled by innumerable regulators and various gene networks which are still not completely understood. MicroRNAs (miRNAs) are one of the important non-coding gene regulators which control an extensive range of biological processes. Thus, exploring miRNAs functions is important for solving gene regulation complexity. The main purpose in the present study is to identify the various gene regulative integrated networks involved in lactation progress in mammary gland. We analyzed ovine mammary tissue data sets which included expression profiles of mRNA (genes) and miRNAs related to six ewes in different days of lactation and nutritional treatments. We combined two different types of information: the network that is module inference by mRNAs (RNA-seq data), miRNAs and transcription factors (TFs) expression matrix and prediction of targets via computational methods. To discover the miRNAs regulatory function, 134 modules were predicted by using gene expression data and 14 TFs and 20 miRNAs were allocated to these predicted modules. By applying this integrated computation-based method, 38 miRNA-modules and 35 TF-module interactions were identified from ovine mammary tissue data during lactogenesis. A lot of these modules were involved in lipid and protein metabolism, as well as steroids and vitamin biosynthesis, which would play key roles in mammary tissue and lactation development. These results present new information about the regulatory procedures at the miRNAs and TF levels throughout lactation.

Proteome of Glioblastoma-Derived Exosomes as a Source of Biomarkers

Extracellular vesicles (EV) are involved in important processes of glioblastoma multiforme (GBM), including malignancy and invasion. EV secreted by glioblastoma cells may cross the hematoencephalic barrier and carry molecular cargo derived from the tumor into the peripheral circulation. Therefore, the determination of the molecular composition of exosomes released by glioblastoma cells seems to be a promising approach for the development of non-invasive methods of the detection of the specific exosomal protein markers in the peripheral blood. The present study aimed to determine the common exosomal proteins presented in preparations from different cell lines and search potential glioblastoma biomarkers in exosomes. We have performed proteomics analysis of exosomes obtained from the conditioned culture medium of five glioblastoma cell lines. A list of 133 proteins common for all these samples was generated. Based on the data obtained, virtual two-dimensional electrophoresis (2DE) maps of proteins presented in exosomes of glioblastoma cells were constructed and the gene ontology (GO) analysis of exosome proteins was performed. A correlation between overexpressed in glial cell proteins and their presence in exosomes have been found. Thus, the existence of many potential glioblastoma biomarkers in exosomes was confirmed.

Ecological divergence in sympatry causes gene misexpression in hybrids

Ecological speciation occurs when reproductive isolation evolves as a byproduct of adaptive divergence between populations. Selection favouring gene regulatory divergence between species could result in transgressive levels of gene expression in F1 hybrids that may lower hybrid fitness. We combined 58 resequenced genomes with 124 transcriptomes to identify patterns of hybrid gene misexpression that may be driven by adaptive regulatory divergence within a young radiation of Cyprinodon pupfishes, which consists of a dietary generalist and two trophic specialists-a molluscivore and a scale-eater. We found more differential gene expression between closely related sympatric specialists than between allopatric generalist populations separated by 1,000 km. Intriguingly, 9.6% of genes that were differentially expressed between sympatric species were also misexpressed in F1 hybrids. A subset of these genes were in highly differentiated genomic regions and enriched for functions important for trophic specialization, including head, muscle and brain development. These regions also included genes that showed evidence of hard selective sweeps and were significantly associated with oral jaw length-the most rapidly diversifying skeletal trait in this radiation. Our results indicate that divergent ecological selection in sympatry can contribute to hybrid gene misexpression which may act as a reproductive barrier between nascent species.