Stochastic Effects in Retrotransposon Dynamics Revealed by Modeling under Competition for Cellular Resources

Transposons are genomic elements that can relocate within a host genome using a ‘cut’- or ‘copy-and-paste’ mechanism. They make up a significant part of many genomes, serve as a driving force for genome evolution, and are linked with Mendelian diseases and cancers. Interactions between two specific retrotransposon types, autonomous (e.g., LINE1/L1) and nonautonomous (e.g., Alu), may lead to fluctuations in the number of these transposons in the genome over multiple cell generations. We developed and examined a simple model of retrotransposon dynamics under conditions where transposon replication machinery competed for cellular resources: namely, free ribosomes and available energy (i.e., ATP molecules). Such competition is likely to occur in stress conditions that a malfunctioning cell may experience as a result of a malignant transformation. The modeling revealed that the number of actively replicating LINE1 and Alu elements in a cell decreases with the increasing competition for resources; however, stochastic effects interfere with this simple trend. We stochastically simulated the transposon dynamics in a cell population and showed that the population splits into pools with drastically different transposon behaviors. The early extinction of active Alu elements resulted in a larger number of LINE1 copies occurring in the first pool, as there was no competition between the two types of transposons in this pool. In the other pool, the competition process remained and the number of L1 copies was kept small. As the level of available resources reached a critical value, both types of dynamics demonstrated an increase in noise levels, and both the period and the amplitude of predator–prey oscillations rose in one of the cell pools. We hypothesized that the presented dynamical effects associated with the impact of the competition for cellular resources inflicted on the dynamics of retrotransposable elements could be used as a characteristic feature to assess a cell state, or to control the transposon activity.

Tissue-Specific Knockdown of Genes of the Argonaute Family Modulates Lifespan and Radioresistance in Drosophila Melanogaster

Small RNAs are essential for the coordination of many cellular processes, including the regulation of gene expression patterns, the prevention of genomic instability, and the suppression of mutagenic transposon activity. These processes determine aging, longevity, and sensitivity of cells and an organism to stress factors (particularly, ionizing radiation). The biogenesis and activity of small RNAs are provided by proteins of the Argonaute family. These proteins participate in the processing of small RNA precursors and the formation of an RNA-induced silencing complex. However, the role of Argonaute proteins in the regulation of lifespan and radioresistance remains poorly explored. We studied the effect of knockdown of Argonaute genes (AGO1, AGO2, AGO3, piwi) in various tissues on the Drosophila melanogaster lifespan and survival after the γ-irradiation at a dose of 700 Gy. In most cases, these parameters were reduced or did not change significantly in flies with tissue-specific RNA interference. Surprisingly, piwi knockdown in both the fat body and the nervous system caused a lifespan increase. But changes in radioresistance depended on the tissue in which the gene was knocked out. In addition, analysis of changes in retrotransposon levels and expression of stress response genes allowed us to determine associated molecular mechanisms.

DNA METHYLATION IN EMBRYOS AND SEEDLINGS OF PLANTS

A complex system of genetic and epigenetic networks controls growth and development of plants. DNA methylation is one of the epigenetic mechanisms involved in suppression of transposon activity, chromatin reorganization, genomic imprinting, and regulation of gene expression. Modulation of the degree of genome methylation is observed during the implementation of morphogenetic development programs and in response to external influences. The change in the degree and pattern of methylation at the embryo-seedling stage in this mini-review is considered. The issues of molecular mechanisms of methylation and its role in the processes of embryo formation, germination and seedling development are discussed.

Involvement of DNA Repair Genes and System of Radiation-Induced Activation of Transposons in Formation of Transgenerational Effects

The study of the genetic basis of the manifestation of radiation-induced effects and their transgenerational inheritance makes it possible to identify the mechanisms of adaptation and possible effective strategies for the survival of organisms in response to chronic radioactive stress. One persistent hypothesis is that the activation of certain genes involved in cellular defense is a specific response of the cell to irradiation. There is also data indicating the important role of transposable elements in the formation of radiosensitivity/radioresistance of biological systems. In this work, we studied the interaction of the systems of hobo transposon activity and DNA repair in the cell under conditions of chronic low-dose irradiation and its participation in the inheritance of radiation-induced transgenerational instability in Drosophila. Our results showed a significant increase of sterility and locus-specific mutability, a decrease of survival, fertility and genome stability (an increase the frequency of dominant lethal mutations and DNA damage) in non-irradiated F1/F2 offspring of irradiated parents with dysfunction of the mus304 gene which is responsible for excision and post-replicative recombination repair and repair of double-stranded DNA breaks. The combined action of dysfunction of the mus309 gene and transpositional activity of hobo elements also led to the transgenerational effects of irradiation but only in the F1 offspring. Dysfunction of the genes of other DNA repair systems (mus101 and mus210) showed no visible effects inherited from irradiated parents subjected to hobo transpositions. The mei-41 gene showed specificity in this type of interaction, which consists in its higher efficiency in sensing events induced by transpositional activity rather than irradiation.

THE ROLE OF TRANSPOSABLE ELEMENTS IN ENDOCRINE CHANGES DURING AGING

Одним из ключевых механизмов старения является изменение гормональной регуляции, для эффективного воздействия на которую с целью продления жизни необходимо определение первопричины данных процессов. В качестве молекулярных драйверов, управляющих динамикой уровня гормонов, могут служить транспозоны. Это связано с их использованием в качестве источников нуклеотидных последовательностей, воспринимающих специфические сигналы рибозимов, транскрипционных факторов, гормонов и их мессенджеров. В то же время, в эволюции транспозоны являются источниками рибозимов и белков, обладающих ДНК-связывающими доменами. Начиная с деления зиготы, видоспецифический состав и распределение транспозонов в геноме могут использоваться как биологическая кодировка, необходимая для последовательной и специфической для типов клеток экспрессии генов. Сделано предположение, что гормональная регуляция является одним из компонентов сложной системы управления онтогенезом под влиянием мобильных элементов. В качестве подтверждения приведены работы о роли транспозонов в управлении генами эндокринной системы, а также о влиянии гормонов на активность транспозонов. Исследование этих взаимосвязей может иметь перспективы для разработки методов продления жизни, так как эпигенетические изменения под влиянием транспозонов носят обратимый характер. Species-specific changes in the endocrine system are key factors in aging. Therefore, to prolong life, it is necessary to find regulators of the highest level, the changes of which lead to physiological aging. The molecular drivers that control dynamics of hormone levels can be transposons. This is due to the use of nucleotide sequences of transposons as binding sites that perceive specific signals of ribozymes, transcription factors, hormones and their messengers. At the same time, transposons are evolutionary sources of ribozymes and proteins that have DNA-binding domains. Starting from zygote division, the species-specific composition and distribution of transposons in the genome serves as a biological coding, which is necessary for the sequential expression of genes specific to cell types and stage of development. We suggest that hormonal regulation is one of the components of this complex system of regulation of ontogenesis under the control of transposons. To confirm our hypothesis, this review contains articles that prove the importance of transposons for species-specific control of endocrine system genes, as well as the effect of hormones on transposon activity. The research of these relationships is promising for the development of methods for the effective prolongation of life, since epigenetic changes under the influence of transposons are reversible.

A Maize Practical Haplotype Graph Leverages Diverse NAM Assemblies

As a result of millions of years of transposon activity, multiple rounds of ancient polyploidization, and large populations that preserve diversity, maize has an extremely structurally diverse genome, evidenced by high-quality genome assemblies that capture substantial levels of both tropical and temperate diversity. We generated a pangenome representation (the Practical Haplotype Graph, PHG) of these assemblies in a database, representing the pangenome haplotype diversity and providing an initial estimate of structural diversity. We leveraged the pangenome to accurately impute haplotypes and genotypes of taxa using various kinds of sequence data, ranging from WGS to extremely-low coverage GBS. We imputed the genotypes of the recombinant inbred lines of the NAM population with over 99% mean accuracy, while unrelated germplasm attained a mean imputation accuracy of 92 or 95% when using GBS or WGS data, respectively. Most of the imputation errors occur in haplotypes within European or tropical germplasm, which have yet to be represented in the maize PHG database. Also, the PHG stores the imputation data in a 30,000-fold more space-efficient manner than a standard genotype file, which is a key improvement when dealing with large scale data.

Potential silencing of gene expression by PIWI-interacting RNAs (piRNAs) in somatic tissues in mollusk

PIWI/piRNA suppress transposon activity in animals, thereby safeguarding the genome from detrimental insertion mutagenesis. Recently, evidence revealed additional piRNA targets and functions in various animals. Although piRNAs are ubiquitously expressed in somatic tissues of the pearl oyster Pinctada fucata, their role is not well-characterized. Here, we report a PIWI/piRNA pathway, including piRNA biogenesis and piRNA-mediated gene regulation in P. fucata. A locked-nucleic-acid modified oligonucleotide (LNA-antagonist) was used to silence a single piRNA (piRNA0001) expression in P. fucata, which resulted in the differential expression of hundreds of endogenous genes. Target prediction analysis revealed that, following silencing, tens of endogenous genes were targeted by piRNA0001, including twelve up-regulated and nine down-regulated genes. Bioinformatic analyses suggested that different piRNA populations participate in the ping-pong amplification loop in a tissue-specific manner. These findings have improved our knowledge of the role of piRNA in mollusks, and provided evidence to understand the regulatory function of the PIWI/piRNA pathway on protein-coding genes outside of germline cells.

Identification and characterization of retro-DNAs, a new type of retrotransposons originated from DNA transposons, in primate genomes

Mobile elements (MEs) can be divided into two major classes based on their transposition mechanisms as retrotransposons and DNA transposons. DNA transposons move in the genomes directly in the form of DNA in a cut-and-paste style, while retrotransposons utilize an RNA-intermediate to transpose in a “copy-and-paste” fashion. In addition to the target site duplications (TSDs), a hallmark of transposition shared by both classes, the DNA transposons also carry terminal inverted repeats (TIRs). DNA transposons constitute ~3% of primate genomes and they are thought to be inactive in the recent primate genomes since ~37My ago despite their success during early primate evolution. Retrotransposons can be further divided into Long Terminal Repeat retrotransposons (LTRs), which are characterized by the presence of LTRs at the two ends, and non-LTRs, which lack LTRs. In the primate genomes, LTRs constitute ~9% of genomes and have a low level of ongoing activity, while non-LTR retrotransposons represent the major types of MEs, contributing to ~37% of the genomes with some members being very young and currently active in retrotransposition. The four known types of non-LTR retrotransposons include LINEs, SINEs, SVAs, and processed pseudogenes, all characterized by the presence of a polyA tail and TSDs, which mostly range from 8 to 15 bp in length. All non-LTR retrotransposons are known to utilize the L1-based target-primed reverse transcription (TPRT) machineries for retrotransposition. In this study, we report a new type of non-LTR retrotransposon, which we named as retro-DNAs, to represent DNA transposons by sequence but non-LTR retrotransposons by the transposition mechanism in the recent primate genomes. By using a bioinformatics comparative genomics approach, we identified a total of 1,750 retro-DNAs, which represent 748 unique insertion events in the human genome and nine non-human primate genomes from the ape and monkey groups. These retro-DNAs, mostly as fragments of full-length DNA transposons, carry no TIRs but longer TSDs with ~23.5% also carrying a polyA tail and with their insertion site motifs and TSD length pattern characteristic of non-LTR retrotransposons. These features suggest that these retro-DNAs are DNA transposon sequences likely mobilized by the TPRT mechanism. Further, at least 40% of these retro-DNAs locate to genic regions, presenting significant potentials for impacting gene function. More interestingly, some retro-DNAs, as well as their parent sites, show certain levels of current transcriptional expression, suggesting that they have the potential to create more retro-DNAs in the current primate genomes. The identification of retro-DNAs, despite small in number, reveals a new mechanism in propagating the DNA transposons sequences in the primate genomes with the absence of canonical DNA transposon activity. It also suggests that the L1 TPRT machinery may have the ability to retrotranspose a wider variety of DNA sequences than what we currently know.

Tissue-specific transposon-associated small RNAs in the gymnosperm tree, Norway spruce

Background Small RNAs (sRNAs) are regulatory molecules impacting on gene expression and transposon activity. MicroRNAs (miRNAs) are responsible for tissue-specific and environmentally-induced gene repression. Short interfering RNAs (siRNA) are constitutively involved in transposon silencing across different type of tissues. The male gametophyte in angiosperms has a unique set of sRNAs compared to vegetative tissues, including phased siRNAs from intergenic or genic regions, or epigenetically activated siRNAs. This is contrasted by a lack of knowledge about the sRNA profile of the male gametophyte of gymnosperms. Results Here, we isolated mature pollen from male cones of Norway spruce and investigated its sRNA profiles. While 21-nt sRNAs is the major size class of sRNAs in needles, in pollen 21-nt and 24-nt sRNAs are the most abundant size classes. Although the 24-nt sRNAs were exclusively derived from TEs in pollen, both 21-nt and 24-nt sRNAs were associated with TEs. We also investigated sRNAs from somatic embryonic callus, which has been reported to contain 24-nt sRNAs. Our data show that the 24-nt sRNA profiles are tissue-specific and differ between pollen and cell culture. Conclusion Our data reveal that gymnosperm pollen, like angiosperm pollen, has a unique sRNA profile, differing from vegetative leaf tissue. Thus, our results reveal that angiosperm and gymnosperm pollen produce new size classes not present in vegetative tissues; while in angiosperm pollen 21-nt sRNAs are generated, in the gymnosperm Norway spruce 24-nt sRNAs are generated. The tissue-specific production of distinct TE-derived sRNAs in angiosperms and gymnosperms provides insights into the diversification process of sRNAs in TE silencing pathways between the two groups of seed plants.

Efficient depletion of ribosomal RNA for RNA sequencing in planarians

Background The astounding regenerative abilities of planarian flatworms prompt steadily growing interest in examining their molecular foundation. Planarian regeneration was found to require hundreds of genes and is hence a complex process. Thus, RNA interference followed by transcriptome-wide gene expression analysis by RNA-seq is a popular technique to study the impact of any particular planarian gene on regeneration. Typically, the removal of ribosomal RNA (rRNA) is the first step of all RNA-seq library preparation protocols. To date, rRNA removal in planarians was primarily achieved by the enrichment of polyadenylated (poly(A)) transcripts. However, to better reflect transcriptome dynamics and to cover also non-poly(A) transcripts, a procedure for the targeted removal of rRNA in planarians is needed. Results In this study, we describe a workflow for the efficient depletion of rRNA in the planarian model species S. mediterranea. Our protocol is based on subtractive hybridization using organism-specific probes. Importantly, the designed probes also deplete rRNA of other freshwater triclad families, a fact that considerably broadens the applicability of our protocol. We tested our approach on total RNA isolated from stem cells (termed neoblasts) of S. mediterranea and compared ribodepleted libraries with publicly available poly(A)-enriched ones. Overall, mRNA levels after ribodepletion were consistent with poly(A) libraries. However, ribodepleted libraries revealed higher transcript levels for transposable elements and histone mRNAs that remained underrepresented in poly(A) libraries. As neoblasts experience high transposon activity this suggests that ribodepleted libraries better reflect the transcriptional dynamics of planarian stem cells. Furthermore, the presented ribodepletion procedure was successfully expanded to the removal of ribosomal RNA from the gram-negative bacterium Salmonella typhimurium. Conclusions The ribodepletion protocol presented here ensures the efficient rRNA removal from low input total planarian RNA, which can be further processed for RNA-seq applications. Resulting libraries contain less than 2% rRNA. Moreover, for a cost-effective and efficient removal of rRNA prior to sequencing applications our procedure might be adapted to any prokaryotic or eukaryotic species of choice.