Genomic and Epigenomic Changes in the Progeny of Cold-Stressed Arabidopsis thaliana Plants.

Plants are continuously exposed to various environmental stresses. Because they can not escape stress, they have to develop mechanisms of remembering stress exposures somatically and passing it to the progeny. We studied the Arabidopsis thaliana ecotype Columbia plants exposed to cold stress for 25 continuous generations. Our study revealed that multigenerational exposure to cold stress resulted in the changes in the genome and epigenome (DNA methylation) across generations. Main changes in the progeny were due to the high frequency of genetic mutations rather than epigenetic changes; the difference was primarily in single nucleotide substitutions and deletions. The progeny of cold-stressed plants exhibited the higher rate of missense non-synonymous mutations as compared to the progeny of control plants. At the same time, epigenetic changes were more common in the CHG (C = cytosine, H = cytosine, adenine or thymine, G = guanine) and CHH contexts and favored hypomethylation. There was an increase in the frequency of C to T (thymine) transitions at the CHH positions in the progeny of cold stressed plants; because this type of mutations is often due to the deamination of the methylated cytosines, it can be hypothesized that environment-induced changes in methylation contribute to mutagenesis and may be to microevolution processes and that RNA-dependent DNA methylation plays a crucial role. Our work supports the existence of heritable stress response in plants and demonstrates that genetic changes prevail.

Scatter Irradiation of Rat Brain Triggers Sex- and Brain Region-Specific Changes in the Expression of Non-Coding RNA Fragments.

Non-coding RNA fragments (ncRFs) are small RNA fragments processed from non-coding RNAs (ncRNAs). ncRFs have various functions and are commonly tissue-specific, and their processing is altered by exposure to stress. Information about ncRFs in the brain is scarce. Recently, we reported the brain region-specific and sex-specific expression of ncRNAs and their processing into ncRFs. Here, we analyzed the expression of ncRFs in the frontal cortex (FC), hippocampus (HIP), and cerebellum (CER) of male and female rats exposed to scatter radiation. We found multiple brain region- and sex-specific changes in response to scatter radiation. Specifically, we observed decreased miRNA expression and the increased expression of ra-ncRNA reads in HIP and CER, as well as an increased number of mtR-NA-associated reads in HIP. We also observed the appearance of sense-intronic ncRNAs-in females, in HIP and FC, and in males, in CER. In this work, we also show that tRNA-GlyGCC and tRNA-GlyCCC are most frequently processed to tRFs, in CER in females, as compared to males. An analysis of the targeted pathways revealed that tRFs and snoRFs in scatter radiation samples mapped to genes in several pathways associated with various neuronal functions. While in HIP and CER these pathways were underrepresented, in FC, they were overrepresented. Such changes may play an important role in pathologies that develop in response to scatter radiation, the effect known as "radio-brain", and may in part explain the sex-specific differences observed in animals and humans exposed to radiation and scatter radiation.

Sex-Specific Expression of Non-Coding RNA Fragments in Frontal Cortex, Hippocampus and Cerebellum of Rats.

Non-coding RNA fragments (ncRFs) are processed from various non-coding RNAs (ncRNAs), with the most abundant being those produced from tRNAs. ncRFs were reported in many animal and plant species. Many ncRFs exhibit tissue specificity or/and are affected by stress. There is, however, only a handful of reports that describe differential expression of ncRFs in the brain regions. In this work, we analyzed the abundance of ncRFs processed from four major ncRNAs, including tRNA (tRFs), snoRNA (snoRFs), snRNA (snRFs), and rRNA (rRFs) in the frontal cortex (FC), hippocampus (HIP), and cerebellum (CER) of male and female rats. We found brain-specific and sex-specific differences. Reads mapping to lincRNAs were significantly larger in CER as compared to HIP and CER, while those mapping to snRNAs and tRNA were smaller in HIP than in FC and CER. tRF reads were the most abundant among all ncRF reads, and FC had more reads than HIP and CER. Reads mapping to antisense ncRNAs were significantly larger in females than in males in FC. Additionally, males consistently had more tRF, snRF, and snoRF reads in all brain regions. rRFs were more abundant in males in FC and females in HIP. Several tRFs were significantly underrepresented, including tRF-ValCAC, tRF-ValACC, and tRF-LysCTT in all brain regions. We also found brain- and sex-specific differences in the number of brain function-related mRNA targets. To summarize, we found sex-specific differences in the expression of several ncRNA fragments in various brain regions of healthy rats.

Multigenerational Exposure to Heat Stress Induces Phenotypic Resilience, and Genetic and Epigenetic Variations in Arabidopsis thaliana Offspring.

Plants are sedentary organisms that constantly sense changes in their environment and react to various environmental cues. On a short-time scale, plants respond through alterations in their physiology, and on a long-time scale, plants alter their development and pass on the memory of stress to the progeny. The latter is controlled genetically and epigenetically and allows the progeny to be primed for future stress encounters, thus increasing the likelihood of survival. The current study intended to explore the effects of multigenerational heat stress in Arabidopsis thaliana. Twenty-five generations of Arabidopsis thaliana were propagated in the presence of heat stress. The multigenerational stressed lineage F25H exhibited a higher tolerance to heat stress and elevated frequency of homologous recombination, as compared to the parallel control progeny F25C. A comparison of genomic sequences revealed that the F25H lineage had a three-fold higher number of mutations [single nucleotide polymorphisms (SNPs) and insertions and deletions (INDELs)] as compared control lineages, suggesting that heat stress induced genetic variations in the heat-stressed progeny. The F25H stressed progeny showed a 7-fold higher number of non-synonymous mutations than the F25C line. Methylome analysis revealed that the F25H stressed progeny showed a lower global methylation level in the CHH context than the control progeny. The F25H and F25C lineages were different from the parental control lineage F2C by 66,491 and 80,464 differentially methylated positions (DMPs), respectively. F25H stressed progeny displayed higher frequency of methylation changes in the gene body and lower in the body of transposable elements (TEs). Gene Ontology analysis revealed that CG-DMRs were enriched in processes such as response to abiotic and biotic stimulus, cell organizations and biogenesis, and DNA or RNA metabolism. Hierarchical clustering of these epimutations separated the heat stressed and control parental progenies into distinct groups which revealed the non-random nature of epimutations. We observed an overall higher number of epigenetic variations than genetic variations in all comparison groups, indicating that epigenetic variations are more prevalent than genetic variations. The largest difference in epigenetic and genetic variations was observed between control plants comparison (F25C vs. F2C), which clearly indicated that the spontaneous nature of epigenetic variations and heat-inducible nature of genetic variations. Overall, our study showed that progenies derived from multigenerational heat stress displayed a notable adaption in context of phenotypic, genotypic and epigenotypic resilience.

Role of the ATP-dependent chromatin remodeling enzyme Fun30/Smarcad1 in the regulation of mRNA splicing.

The yeast ATP-dependent chromatin remodeling enzyme Fun30 has been shown to regulate heterochromatin silencing, DNA repair, transcription, and chromatin organization. Although chromatin structure has been proposed to influence splice site recognition and regulation, whether ATP-dependent chromatin remodeling enzyme plays a role in regulating splicing is not known. In this study, we find that pre-mRNA splicing efficiency is impaired and the recruitment of spliceosome is compromised in Fun30-depleted cells. In addition, Fun30 is enriched in the gene body of individual intron-containing genes. Moreover, we show that pre-mRNA splicing efficiency is dependent on the chromatin remodeling activity of Fun30. The function of Fun30 in splicing is further supported by the observation that, Smarcad1, the mammalian homolog of Fun30, regulates alternative splicing. Taken together, these results provide evidence for a novel role of Fun30 in regulating splicing.

System, Method and Software for Calculation of a Cannabis Drug Efficiency Index for the Reduction of Inflammation.

There are many varieties of Cannabis sativa that differ from each other by composition of cannabinoids, terpenes and other molecules. The medicinal properties of these cultivars are often very different, with some being more efficient than others. This report describes the development of a method and software for the analysis of the efficiency of various cannabis extracts to detect the anti-inflammatory properties of the various cannabis extracts. The method uses high-throughput gene expression profiling data but can potentially use other omics data as well. According to the signaling pathway topology, the gene expression profiles are convoluted into the signaling pathway activities using a signaling pathway impact analysis (SPIA) method. The method was tested by inducing inflammation in human 3D epithelial tissues, including intestine, oral and skin, and then exposing these tissues to various extracts and then performing transcriptome analysis. The analysis showed a different efficiency of the various extracts in restoring the transcriptome changes to the pre-inflammation state, thus allowing to calculate a different cannabis drug efficiency index (CDEI).

Transgenerational Response to Heat Stress in the Form of Differential Expression of Noncoding RNA Fragments in Brassica rapa Plants.

Epigenetic regulations in the form of changes in differential expression of noncoding RNAs (ncRNAs) are an essential mechanism of stress response in plants. Previously we showed that heat treatment in L. results in the differential processing and accumulation of ncRNA fragments (ncRFs) stemming from transfer RNAs (tRNAs), ribosomal RNAs (rRNAs), small nuclear RNAs (snRNAs), and small nucleolar RNAs (snoRNAs). In this work, we analyzed whether ncRFs are differentially expressed in the progeny of heat-stressed plants. We found significant changes in the size of tRF reads and a significant decrease in the percentage of tRFs mapping to tRNA-Ala, tRNA-Arg, and tRNA-Tyr and an increase in tRFs mapping to tRNA-Asp. The enrichment analysis showed significant differences in processing of tRFs from tRNA, tRNA, tRNA, tRNA, tRNA, and tRNA isoacceptors. Analysis of potential targets of tRFs showed that they regulate brassinosteroid metabolism, the proton pump ATPase activity, the antiporter activity, the mRNA decay activity as well as nucleosome positioning and the epigenetic regulation of transgenerational response. Gene ontology term analysis of potential targets demonstrated a significant enrichment in tRFs that potentially targeted a cellular component endoplasmic reticulum (ER) and in small nucleolar RNA fragments (snoRFs), the molecular function protein binding. To summarize, our work demonstrated that the progeny of heat-stressed plants exhibit changes in the expression of tRFs and snoRFs but not of small nuclear RNA fragments (snRFs) or ribosomal RNA fragments (rRFs) and these changes likely better prepare the progeny of stressed plants to future stress encounters.

Computational Characterization of ncRNA Fragments in Various Tissues of the Brassica rapa Plant.

Recently, a novel type of non-coding RNA (ncRNA), known as ncRNA fragments or ncRFs, has been characterised in various organisms, including plants. The biogenesis mechanism, function and abundance of ncRFs stemming from various ncRNAs are poorly understood, especially in plants. In this work, we have computationally analysed the composition of ncRNAs and the fragments that derive from them in various tissues of Brassica rapa plants, including leaves, meristem tissue, pollen, unfertilized and fertilized ova, embryo and endosperm. Detailed analysis of transfer RNA (tRNA) fragments (tRFs), ribosomal RNA (rRNA) fragments (rRFs), small nucleolar RNA (snoRNA) fragments (snoRFs) and small nuclear RNA (snRNA) fragments (snRFs) showed a predominance of tRFs, with the 26 nucleotides (nt) fraction being the largest. Mapping ncRF reads to full-length mature ncRNAs showed a strong bias for one or both termini. tRFs mapped predominantly to the 5' end, whereas snRFs mapped to the 3' end, suggesting that there may be specific biogenesis and retention mechanisms. In the case of tRFs, specific isoacceptors were enriched, including tRNAGly(UCC) and tRFAsp(GUC). The analysis showed that the processing of 26-nt tRF5' occurred by cleavage at the last unpaired nucleotide of the loop between the D arm and the anticodon arm. Further support for the functionality of ncRFs comes from the analysis of binding between ncRFs and their potential targets. A higher average percentage of binding at the first half of fragments was observed, with the highest percentage being at 2-6 nt. To summarise, our analysis showed that ncRFs in B. rapa are abundantly produced in a tissue-specific manner, with bias toward a terminus, the bias toward the size of generated fragments and the bias toward the targeting of specific biological processes.

Non-coding RNAs match the deleted genomic regions in humans.

RNA is transcribed from DNA, and therefore, there should be no RNA transcript from the deleted DNA region. Our study attempted to analyse whether any RNA cache that maps the deleted regions is present in human cells. Using data from the 1000 genome project, we selected 41 CEPH (CEU) and 38 Yoruba (YRI) samples that included the data for the entire genome sequence and ncRNA and mRNA sequences. Aligning the ncRNA reads against the genomic DNA in individual samples has revealed that 229 out of 1114 homozygous deletions have ncRNA reads that map to them. Further analysis has revealed that ncRNA reads that map the deleted regions are enriched around the deletion ends and at genic regions of the genome. The read enrichment at deletion ends suggests that these ncRNAs are likely some form of double-strand break induced RNAs. Our analysis suggests that human cells may contain a residual ncRNA cache that is possibly propagated across generations.

Pattern Recognition on Read Positioning in Next Generation Sequencing.

The usefulness and the utility of the next generation sequencing (NGS) technology are based on the assumption that the DNA or cDNA cleavage required to generate short sequence reads is random. Several previous reports suggest the existence of sequencing bias of NGS reads. To address this question in greater detail, we analyze NGS data from four organisms with different GC content, Plasmodium falciparum (19.39%), Arabidopsis thaliana (36.03%), Homo sapiens (40.91%) and Streptomyces coelicolor (72.00%). Using machine learning techniques, we recognize the pattern that the NGS read start is positioned in the local region where the nucleotide distribution is dissimilar from the global nucleotide distribution. We also demonstrate that the mono-nucleotide distribution underestimates sequencing bias, and the recognized pattern is explained largely by the distribution of multi-nucleotides (di-, tri-, and tetra- nucleotides) rather than mono-nucleotides. This implies that the correction of sequencing bias needs to be performed on the basis of the multi-nucleotide distribution. Providing companion software to quantify the effect of the recognized pattern on read positioning, we exemplify that the bias correction based on the mono-nucleotide distribution may not be sufficient to clean sequencing bias.

Sex-specific effects of cytotoxic chemotherapy agents cyclophosphamide and mitomycin C on gene expression, oxidative DNA damage, and epigenetic alterations in the prefrontal cortex and hippocampus - an aging connection.

Recent research shows that chemotherapy agents can be more toxic to healthy brain cells than to the target cancer cells. They cause a range of side effects, including memory loss and cognitive dysfunction that can persist long after the completion of treatment. This condition is known as chemo brain. The molecular and cellular mechanisms of chemo brain remain obscure. Here, we analyzed the effects of two cytotoxic chemotherapy drugs-cyclophosphamide (CPP) and mitomycin C (MMC) - on transcriptomic and epigenetic changes in the murine prefrontal cortex (PFC) and hippocampal regions. We for the first time showed that CPP and MMC treatments led to profound sex- and brain region-specific alterations in gene expression profiles. Gene expression changes were most prominent in the PFC tissues of female mice 3 weeks after MMC treatment, and the gene expression response was much greater for MCC than CPP exposure. MMC exposure resulted in oxidative DNA damage, evidenced by accumulation of 8-oxo-2'-deoxyguanosine (8-oxodG) and a decrease in the level of 8-oxodG repair protein OGG1 in the PFC of female animals 3 weeks after treatment. MMC treatment decreased global DNA methylation and increased DNA hydroxymethylation in the PFC tissues of female mice. The majority of the changes induced by chemotherapy in the PFC tissues of female mice resembled those that occur during the brain's aging processes. Therefore, our study suggests a link between chemotherapy-induced chemo brain and brain aging, and provides an important roadmap for future analysis.

The ATP-dependent chromatin remodeling enzyme Fun30 represses transcription by sliding promoter-proximal nucleosomes.

The evolutionarily conserved ATP-dependent chromatin remodeling enzyme Fun30 has recently been shown to play important roles in heterochromatin silencing and DNA repair. However, how Fun30 remodels nucleosomes is not clear. Here we report a nucleosome sliding activity of Fun30 and its role in transcriptional repression. We observed that Fun30 repressed the expression of genes involved in amino acid and carbohydrate metabolism, the stress response, and meiosis. In addition, Fun30 was localized at the 5' and 3' ends of genes and within the open reading frames of its targets. Consistent with its role in gene repression, we observed that Fun30 target genes lacked histone modifications often associated with gene activation and showed an increased level of ubiquitinated histone H2B. Furthermore, a genome-wide nucleosome mapping analysis revealed that the length of the nucleosome-free region at the 5' end of a subset of genes was changed in Fun30-depleted cells. In addition, the positions of the -1, +2, and +3 nucleosomes at the 5' end of target genes were shifted significantly, whereas the position of the +1 nucleosome remained largely unchanged in the fun30Δ mutant. Finally, we demonstrated that affinity-purified, single-component Fun30 exhibited a nucleosome sliding activity in an ATP-dependent manner. These results define a role for Fun30 in the regulation of transcription and indicate that Fun30 remodels chromatin at the 5' end of genes by sliding promoter-proximal nucleosomes.

NucVoter: A Voting Algorithm for Reliable Nucleosome Prediction Using Next-Generation Sequencing Data.

Nucleosomes, which consist of DNA wrapped around histone octamers, are dynamic, and their structure, including their location, size, and occupancy, can be transformed. Nucleosomes can regulate gene expression by controlling the DNA accessibility of proteins. Using next-generation sequencing techniques along with such laboratory methods as micrococcal nuclease digestion, predicting the genomic locations of nucleosomes is possible. However, the true locations of nucleosomes are unknown, and it is difficult to determine their exact locations using next-generation sequencing data. This paper proposes a novel voting algorithm, NucVoter, for the reliable prediction of nucleosome locations. Multiple models verify the consensus areas in which nucleosomes are placed by the model with the highest priority. NucVoter significantly improves the performance of nucleosome prediction.