Transposon variants and their effects on gene expression in Arabidopsis

The evolution of transposable elements in Brachypodium distachyon is governed by purifying selection, while neutral and adaptive processes play a minor role

Understanding how plants adapt to changing environments and the potential contribution of transposable elements (TEs) to this process is a key question in evolutionary genomics. While TEs have recently been put forward as active players in the context of adaptation, few studies have thoroughly investigated their precise role in plant evolution. Here, we used the wild Mediterranean grass Brachypodium distachyon as a model species to identify and quantify the forces acting on TEs during the adaptation of this species to various conditions, across its entire geographic range. Using sequencing data from more than 320 natural B. distachyon accessions and a suite of population genomics approaches, we reveal that putatively adaptive TE polymorphisms are rare in wild B. distachyon populations. After accounting for changes in past TE activity, we show that only a small proportion of TE polymorphisms evolved neutrally (<10%), while the vast majority of them are under moderate purifying selection regardless of their distance to genes. TE polymorphisms should not be ignored when conducting evolutionary studies, as they can be linked to adaptation. However, our study clearly shows that while they have a large potential to cause phenotypic variation in B. distachyon, they are not favored during evolution and adaptation over other types of mutations (such as point mutations) in this species.

Whole genome doubling-induced the enrichment of H3K27me3 in genes carrying specific TEs in Aegilops tauschii

Polyploidization plays important roles in the evolution and breeding of the common wheat. Aegilops tauschii, the D-genome progenitor of the common wheat, provides a valuable pool of resistance genes to multiple diseases. Extensive studies focus on the exploration of these genes for wheat improvement. However, few studies have unveiled alternations on genome-wide expression pattern and histone modifications induced by whole-genome doubling (WGD) process. In this study, we conducted transcriptome analysis for the diploid and tetraploid Ae. taushcii lines using the leaf and root tissues. Both lines tend to display similar tissue-specific pattern. Interestingly, we found that TEs located in genic regions were depleted of the repressive histone mark H3K27me3, whereas their adjacent chromatin was enriched with H3K27me3. The tetraploid line exhibited higher levels of H3K27me3 in those regions than the diploid line, particularly for genic regions associated with TEs of the long interspersed nuclear elements (LINEs), CACTA, PIF/Harbinger, Tc1/Mariner and unclassed DNA transposon. Surprisingly, the expression levels of these TEs cognate genes were negatively associated with the levels of H3K27me3 between the tetraploid and diploid lines, suggesting the five types of TEs located within genic regions might be involved in the regulation of the ploidy-related gene expression, possibly through differential enrichment of H3K27me3 in the genic regions. These findings will help to understand the potential role of specific types of TEs on transcription in response to WGD.

Methylome, transcriptome, and phenotype changes induced by temperature conditions experienced during sexual reproduction in Fragaria vesca

Temperature conditions experienced during embryogenesis and seed development may induce epigenetic changes that increase phenotypic variation in plants. Here we investigate if embryogenesis and seed development at two different temperatures (28 vs. 18°C) result in lasting phenotypic effects and DNA methylation changes in woodland strawberry (Fragaria vesca). Using five European ecotypes from Spain (ES12), Iceland (ICE2), Italy (IT4), and Norway (NOR2 and NOR29), we found statistically significant differences between plants from seeds produced at 18 or 28°C in three of four phenotypic features investigated under common garden conditions. This indicates the establishment of a temperature-induced epigenetic memory-like response during embryogenesis and seed development. The memory effect was significant in two ecotypes: in NOR2 flowering time, number of growth points and petiole length were affected, and in ES12 number of growth points was affected. This indicates that genetic differences between ecotypes in their epigenetic machinery, or other allelic differences, impact this type of plasticity. We observed statistically significant differences between ecotypes in DNA methylation marks in repetitive elements, pseudogenes, and genic elements. Leaf transcriptomes were also affected by embryonic temperature in an ecotype-specific manner. Although we observed significant and lasting phenotypic change in at least some ecotypes, there was considerable variation in DNA methylation between individual plants within each temperature treatment. This within-treatment variability in DNA methylation marks in F. vesca progeny may partly be a result of allelic redistribution from recombination during meiosis and subsequent epigenetic reprogramming during embryogenesis.

Chromosome-level genome assembly of Murraya paniculata sheds light on biosynthesis of floral volatiles

BACKGROUND

Murraya paniculata (L.) Jack, commonly called orange jessamine in the family Rutaceae, is an important ornamental plant in tropical and subtropical regions which is famous for its strong fragrance. Although genome assemblies have been reported for many Rutaceae species, mainly in the genus Citrus, full genomic information has not been reported for M. paniculata, which is a prerequisite for in-depth genetic studies on Murraya and manipulation using genetic engineering techniques. Here, we report a high-quality chromosome-level genome assembly of M. paniculata and aim to provide insights on the molecular mechanisms of flower volatile biosynthesis.

RESULTS

The genome assembly with a contig N50 of 18.25 Mb consists of 9 pseudomolecules and has a total length of 216.86 Mb. Phylogenetic analysis revealed that M. paniculata diverged from the common ancestor approximately 25 million years ago and has not undergone any species-specific whole genome duplication events. Genome structural annotation and comparative genomics analysis revealed that there are obvious differences in transposon contents among the genomes of M. paniculata and Citrus species, especially in the upstream regions of genes. Research on the flower volatiles of M. paniculata and C. maxima at three flowering stages revealed significant differences in volatile composition with the flowers of C. maxima lacking benzaldehyde and phenylacetaldehyde. Notably, there are transposons inserted in the upstream region of the phenylacetaldehyde synthase (PAAS) genes Cg1g029630 and Cg1g029640 in C. maxima, but not in the upstream region of three PAAS genes Me2G_2379, Me2G_2381, and Me2G_2382 in M. paniculata. Our results indicated that compared to the low expression levels of PAAS genes in C. maxima, the higher expression levels of the three PAAS genes in M. paniculata are the main factor affecting the phenylacetaldehyde biosynthesis and causing the content difference of phenylacetaldehyde. The phenylacetaldehyde synthetic activities of the enzymes encoded by M. paniculata PAAS genes were validated by in vitro analyses.

CONCLUSIONS

Our study provides useful genomic resources of M. paniculata for further research on Rutaceae plants, identifies new PAAS genes, and provides insights into how transposons contribute to variations in flower volatiles among Murraya and Citrus plants.

Off-target piRNA gene silencing in Drosophila melanogaster rescued by a transposable element insertion

Transposable elements (TE) are selfish genetic elements that can cause harmful mutations. In Drosophila, it has been estimated that half of all spontaneous visible marker phenotypes are mutations caused by TE insertions. Several factors likely limit the accumulation of exponentially amplifying TEs within genomes. First, synergistic interactions between TEs that amplify their harm with increasing copy number are proposed to limit TE copy number. However, the nature of this synergy is poorly understood. Second, because of the harm posed by TEs, eukaryotes have evolved systems of small RNA-based genome defense to limit transposition. However, as in all immune systems, there is a cost of autoimmunity and small RNA-based systems that silence TEs can inadvertently silence genes flanking TE insertions. In a screen for essential meiotic genes in Drosophila melanogaster, a truncated Doc retrotransposon within a neighboring gene was found to trigger the germline silencing of ald, the Drosophila Mps1 homolog, a gene essential for proper chromosome segregation in meiosis. A subsequent screen for suppressors of this silencing identified a new insertion of a Hobo DNA transposon in the same neighboring gene. Here we describe how the original Doc insertion triggers flanking piRNA biogenesis and local gene silencing. We show that this local gene silencing occurs in cis and is dependent on deadlock, a component of the Rhino-Deadlock-Cutoff (RDC) complex, to trigger dual-strand piRNA biogenesis at TE insertions. We further show how the additional Hobo insertion leads to de-silencing by reducing flanking piRNA biogenesis triggered by the original Doc insertion. These results support a model of TE-mediated gene silencing by piRNA biogenesis in cis that depends on local determinants of transcription. This may explain complex patterns of off-target gene silencing triggered by TEs within populations and in the laboratory. It also provides a mechanism of sign epistasis among TE insertions, illuminates the complex nature of their interactions and supports a model in which off-target gene silencing shapes the evolution of the RDC complex.

Megakaryopoiesis impairment through acute innate immune signaling activation by azacitidine

Thrombocytopenia, prevalent in the majority of patients with myeloid malignancies, such as myelodysplastic syndrome (MDS) or acute myeloid leukemia (AML), is an independent adverse prognostic factor. Azacitidine (AZA), a mainstay therapeutic agent for stem cell transplant-ineligible patients with MDS/AML, often transiently induces or further aggravates disease-associated thrombocytopenia by an unknown mechanism. Here, we uncover the critical role of an acute type-I interferon (IFN-I) signaling activation in suppressing megakaryopoiesis in AZA-mediated thrombocytopenia. We demonstrate that megakaryocytic lineage-primed progenitors present IFN-I receptors and, upon AZA exposure, engage STAT1/SOCS1-dependent downstream signaling prematurely attenuating thrombopoietin receptor (TPO-R) signaling and constraining megakaryocytic progenitor cell growth and differentiation following TPO-R stimulation. Our findings directly implicate RNA demethylation and IFN-I signal activation as a root cause for AZA-mediated thrombocytopenia and suggest mitigation of TPO-R inhibitory innate immune signaling as a suitable therapeutic strategy to support platelet production, particularly during the early phases of AZA therapy.

The evolution and function of transposons in epigenetic regulation in response to the environment

Transposable elements (TEs) make up a major proportion of plant genomes. Despite their prevalence genome-wide, TEs are often tossed aside as "junk DNA" since they rarely cause phenotypes, and epigenetic mechanisms silence TEs to prevent them from causing deleterious mutations through movement. While this bleak picture of TEs in genomes is true on average, a growing number of examples across many plant species point to TEs as drivers of phenotypic diversity and novel stress responses. Examples of TE-influenced phenotypes illustrate the many ways that novel transposition events can alter local gene expression and how this relates to potential variation in plant responses to environmental stress. Since TE families and insertions at the locus level lack evolutionary conservation, advancements in the field will require TE experts across diverse species to identify and utilize TE variation in their own systems as a means of crop improvement.

Transposable Elements in the Revealing of Polymorphism-Based Differences in the Seeds of Flax Varieties Grown in Remediated Chernobyl Area

The nuclear reactor accident in Chernobyl, Ukraine, resulted in effects both locally and farther away. Most of the contaminated areas were the agricultural fields and forests. Experimental fields were established near Chernobyl—radioactively contaminated fields localized 5 km from Chernobyl Nuclear Power Plant as well as the remediated soil that is localized directly in the Chernobyl town. Two flax varieties growing under chronic exposition to ionizing radiation were used for this study—the local Ukrainian variety Kyivskyi and a commercial variety Bethune. The screening of the length polymorphism generated by transposable elements insertions were performed. All known types of common flax transposon, retrotransposons and iPBS approach were used. In the iPBS multiplex analyze, for the Kyivskyi variety, a unique addition was found in the seeds from the radioactive contaminated field and for the Bethune variety, a total of five amplicon additions were obtained and one deletion. For the TRIM Cassandra fingerprints, two amplicon additions were generated in the seeds from radioactive contaminated fields for the Bethune variety. In summary, the obtained data represent the genetic diversity between control and irradiated subgroups of flax seeds from Chernobyl area and the presence of activated transposable elements due to the irradiation stress.

DNA methylation in transposable elements buffers the connection between three-dimensional chromatin organization and gene transcription upon rice genome duplication

INTRODUCTION

Polyploidy is a major force in plant evolution and the domestication of cultivated crops.

OBJECTIVES

The study aimed to explore the relationship and underlying mechanism between three-dimensional (3D) chromatin organization and gene transcription upon rice genome duplication.

METHODS

The 3D chromatin structures between diploid (2C) and autotetraploid (4C) rice were compared using high-throughput chromosome conformation capture (Hi-C) analysis. The study combined genetics, transcriptomics, whole-genome bisulfite sequencing (WGBS-seq) and 3D genomics approaches to uncover the mechanism for DNA methylation in modulating gene transcription through 3D chromatin architectures upon rice genome duplication.

RESULTS

We found that 4C rice presents weakened intra-chromosomal interactions compared to its 2C progenitor in some chromosomes. In addition, we found that changes of 3D chromatin organizations including chromatin compartments, topologically associating domains (TADs), and loops, are uncorrelated with gene transcription. Moreover, DNA methylations in the regulatory sequences of genes in compartment A/B switched regions and TAD boundaries are unrelated to their expression. Importantly, although there was no significant difference in the methylation levels in transposable elements (TEs) in differentially expressed gene (DEG) and non-DEG promoters between 2C and 4C rice, we found that the hypermethylated TEs across genes in compartment A/B switched regions and TAD boundaries may suppress the expression of these genes.

CONCLUSION

The study proposed that the rice genome doubling might modulate TE methylation to buffer the effects of chromatin architecture on gene transcription in compartment A/B switched regions and TAD boundaries, resulting in the disconnection between 3D chromatin structure alteration and gene transcription upon rice genome duplication.

The evolutionary history of small RNAs in Solanaceae

The Solanaceae or "nightshade" family is an economically important group with remarkable diversity. To gain a better understanding of how the unique biology of the Solanaceae relates to the family's small RNA (sRNA) genomic landscape, we downloaded over 255 publicly available sRNA data sets that comprise over 2.6 billion reads of sequence data. We applied a suite of computational tools to predict and annotate two major sRNA classes: (1) microRNAs (miRNAs), typically 20- to 22-nucleotide (nt) RNAs generated from a hairpin precursor and functioning in gene silencing and (2) short interfering RNAs (siRNAs), including 24-nt heterochromatic siRNAs typically functioning to repress repetitive regions of the genome via RNA-directed DNA methylation, as well as secondary phased siRNAs and trans-acting siRNAs generated via miRNA-directed cleavage of a polymerase II-derived RNA precursor. Our analyses described thousands of sRNA loci, including poorly understood clusters of 22-nt siRNAs that accumulate during viral infection. The birth, death, expansion, and contraction of these sRNA loci are dynamic evolutionary processes that characterize the Solanaceae family. These analyses indicate that individuals within the same genus share similar sRNA landscapes, whereas comparisons between distinct genera within the Solanaceae reveal relatively few commonalities.

TE Density: a tool to investigate the biology of transposable elements

BACKGROUND

Transposable elements (TEs) are powerful creators of genotypic and phenotypic diversity due to their inherent mutagenic capabilities and in this way they serve as a deep reservoir of sequences for genomic variation. As agents of genetic disruption, a TE's potential to impact phenotype is partially a factor of its location in the genome. Previous research has shown TEs' ability to impact the expression of neighboring genes, however our understanding of this trend is hampered by the exceptional amount of diversity in the TE world, and a lack of publicly available computational methods that quantify the presence of TEs relative to genes.

RESULTS

Here, we have developed a tool to more easily quantify TE presence relative to genes through the use of only a gene and TE annotation, yielding a new metric we call TE Density. Briefly defined as the proportion of TE-occupied base-pairs relative to a window-size of the genome. This new pipeline reports TE density for each gene in the genome, for each type descriptor of TE (order and superfamily), and for multiple positions and distances relative to the gene (upstream, intragenic, and downstream) over sliding, user-defined windows. In this way, we overcome previous limitations to the study of TE-gene relationships by focusing on all TE types present in the genome, utilizing flexible genomic distances for measurement, and reporting a TE presence metric for every gene in the genome.

CONCLUSIONS

Together, this new tool opens up new avenues for studying TE-gene relationships, genome architecture, comparative genomics, and the tremendous diversity present of the TE world. TE Density is open-source and freely available at: https://github.com/sjteresi/TE_Density .

DNA Methylome and LncRNAome Analysis Provide Insights Into Mechanisms of Genome-Dosage Effects in Autotetraploid Cassava

Whole genome duplication (WGD) increases the dosage of all coding and non-coding genes, yet the molecular implications of genome-dosage effects remain elusive. In this study, we generated integrated maps of the methylomes and lncRNAomes for diploid and artificially generated autotetraploid cassava (Manihot esculenta Crantz). We found that transposable elements (TEs) suppressed adjacent protein coding gene (PCG)-expression levels, while TEs activated the expression of nearby long non-coding RNAs (lncRNAs) in the cassava genome. The hypermethylation of DNA transposons in mCG and mCHH sites may be an effective way to suppress the expression of nearby PCGs in autotetraploid cassava, resulting in similar expression levels for most of PCGs between autotetraploid and diploid cassava. In the autotetraploid, decreased methylation levels of retrotransposons at mCHG and mCHH sites contributed to reduced methylation of Gypsy-neighboring long intergenic non-coding RNAs, potentially preserving diploid-like expression patterns in the major of lncRNAs. Collectively, our study highlighted that WGD-induced DNA methylation variation in DNA transposons and retrotransposons may be as direct adaptive responses to dosage of all coding-genes and lncRNAs, respectively.

Intergenic spaces: a new frontier to improving plant health

To more sustainably mitigate the impact of crop diseases on plant health and productivity, there is a need for broader spectrum, long-lasting resistance traits. Defense response (DR) genes, located throughout the genome, participate in cellular and system-wide defense mechanisms to stave off infection by diverse pathogens. This multigenic resistance avoids rapid evolution of a pathogen to overcome host resistance. DR genes reside within resistance-associated quantitative trait loci (QTL), and alleles of DR genes in resistant varieties are more active during pathogen attack relative to susceptible haplotypes. Differential expression of DR genes results from polymorphisms in their regulatory regions, that includes cis-regulatory elements such as transcription factor binding sites as well as features that influence epigenetic structural changes to modulate chromatin accessibility during infection. Many of these elements are found in clusters, known as cis-regulatory modules (CRMs), which are distributed throughout the host genome. Regulatory regions involved in plant-pathogen interactions may also contain pathogen effector binding elements that regulate DR gene expression, and that, when mutated, result in a change in the plants' response. We posit that CRMs and the multiple regulatory elements that comprise them are potential targets for marker-assisted breeding for broad-spectrum, durable disease resistance.

The Dynamism of Transposon Methylation for Plant Development and Stress Adaptation

Plant development processes are regulated by epigenetic alterations that shape nuclear structure, gene expression, and phenotypic plasticity; these alterations can provide the plant with protection from environmental stresses. During plant growth and development, these processes play a significant role in regulating gene expression to remodel chromatin structure. These epigenetic alterations are mainly regulated by transposable elements (TEs) whose abundance in plant genomes results in their interaction with genomes. Thus, TEs are the main source of epigenetic changes and form a substantial part of the plant genome. Furthermore, TEs can be activated under stress conditions, and activated elements cause mutagenic effects and substantial genetic variability. This introduces novel gene functions and structural variation in the insertion sites and primarily contributes to epigenetic modifications. Altogether, these modifications indirectly or directly provide the ability to withstand environmental stresses. In recent years, many studies have shown that TE methylation plays a major role in the evolution of the plant genome through epigenetic process that regulate gene imprinting, thereby upholding genome stability. The induced genetic rearrangements and insertions of mobile genetic elements in regions of active euchromatin contribute to genome alteration, leading to genomic stress. These TE-mediated epigenetic modifications lead to phenotypic diversity, genetic variation, and environmental stress tolerance. Thus, TE methylation is essential for plant evolution and stress adaptation, and TEs hold a relevant military position in the plant genome. High-throughput techniques have greatly advanced the understanding of TE-mediated gene expression and its associations with genome methylation and suggest that controlled mobilization of TEs could be used for crop breeding. However, development application in this area has been limited, and an integrated view of TE function and subsequent processes is lacking. In this review, we explore the enormous diversity and likely functions of the TE repertoire in adaptive evolution and discuss some recent examples of how TEs impact gene expression in plant development and stress adaptation.

Genetic toolbox and regulatory circuits of plant-nematode associations

Plant-nematode associations are the most imperative area of study that forms the basis to understand their regulatory networks and coordinated functional aspects. Nematodes are highly parasitic organisms known so far, to cause relentless damage towards agricultural crops on a global scale. They pierce the roots of host plants and form neo-plastic feeding structures to extract out resources for their functional development. Moreover, they undergo re-differentiation within plant cells to form giant multi-nucleate feeding structures or syncytium. All these processes are facilitated by numerous transcriptomic, proteomic, metabolomic and epigenetic modifications, that regulate different biological attractions among plants and nematodes. Nevertheless, these mechanisms are quite remarkable and have been explored in the present review. Here, we have shed light on genomic as well as genetic approaches to acquire an effective understanding regarding plant-nematode associations. Transcriptomics have revealed an extensive network to unravel feeding mechanism of nematodes through gene-expression programming of target genes. Also, the regulatory circuits of epigenetic alterations through DNA-methylation, non-coding RNAs and histone modifications very well explain epigenetic profiling within plants. Since decades, research have observed many intricacies to elucidate the dynamic nature of epigenetic modulations in plant-nematode attractions. By this review, we have highlighted the functional aspects of small RNAs in inducing plant-nematode parasitism along with the putative role of miRNAs. These RNAs act as chief genetic elements to mediate the expressional changes in plants through post-transcriptional silencing of various effector proteins as well as transcriptional factors. A pragmatic role of miRNAs in modulating gene expression in nematode infection and feeding site development have also been reviewed. Hence, they have been considered master regulators for functional reprogramming the expression during establishment of feeding sites. We have also encapsulated the advancement of genome-broadened DNA-methylation and untangled the nematode mediated dynamic alterations within plant methylome along with assessing transcriptional activities of various genes and transposons. In particular, we have highlighted the role of effector proteins in stimulating epigenetic changes. Finally, we have emerged towards a molecular-based core understanding about plant-nematode associations.

Aphid feeding induces the relaxation of epigenetic control and the associated regulation of the defense response in Arabidopsis

Environmentally induced changes in the epigenome help individuals to quickly adapt to fluctuations in the conditions of their habitats. We explored those changes in Arabidopsis thaliana plants subjected to multiple biotic and abiotic stresses, and identified transposable element (TE) activation in plants infested with the green peach aphid, Myzus persicae. We performed a genome-wide analysis mRNA expression, small RNA accumulation and DNA methylation Our results demonstrate that aphid feeding induces loss of methylation of hundreds of loci, mainly TEs. This loss of methylation has the potential to regulate gene expression and we found evidence that it is involved in the control of plant immunity genes. Accordingly, mutant plants deficient in DNA and H3K9 methylation (kyp) showed increased resistance to M. persicae infestation. Collectively, our results show that changes in DNA methylation play a significant role in the regulation of the plant transcriptional response and induction of defense response against aphid feeding.

Non-perfectly matching small RNAs can induce stable and heritable epigenetic modifications and can be used as molecular markers to trace the origin and fate of silencing RNAs

Short non-coding RNA molecules (sRNAs) play a fundamental role in gene regulation and development in higher organisms. They act as molecular postcodes and guide AGO proteins to target nucleic acids. In plants, sRNA-targeted mRNAs are degraded, reducing gene expression. In contrast, sRNA-targeted DNA sequences undergo cytosine methylation referred to as RNA-directed DNA methylation (RdDM). Cytosine methylation can suppress transcription, thus sRNAs are potent regulators of gene expression. sRNA-mediated RdDM is involved in genome stability through transposon silencing, mobile signalling for epigenetic gene control and hybrid vigour. Since cytosine methylation can be passed on to subsequent generations, RdDM contributes to transgenerational inheritance of the epigenome. Using a novel approach, which can differentiate between primary (inducer) and secondary (amplified) sRNAs, we show that initiation of heritable RdDM does not require complete sequence complementarity between the sRNAs and their nuclear target sequences. sRNAs with up to four regularly interspaced mismatches are potent inducers of RdDM, however, the number and disruptive nature of nucleotide polymorphisms negatively correlate with their efficacy. Our findings contribute to understanding how sRNA can directly shape the epigenome and may be used in designing the next generation of RNA silencing constructs.

Hybridization History and Repetitive Element Content in the Genome of a Homoploid Hybrid, Yucca gloriosa (Asparagaceae)

Hybridization in plants results in phenotypic and genotypic perturbations that can have dramatic effects on hybrid physiology, ecology, and overall fitness. Hybridization can also perturb epigenetic control of transposable elements, resulting in their proliferation. Understanding the mechanisms that maintain genomic integrity after hybridization is often confounded by changes in ploidy that occur in hybrid plant species. Homoploid hybrid species, which have no change in chromosome number relative to their parents, offer an opportunity to study the genomic consequences of hybridization in the absence of change in ploidy. Yucca gloriosa (Asparagaceae) is a young homoploid hybrid species, resulting from a cross between Yucca aloifolia and Yucca filamentosa. Previous analyses of ∼11 kb of the chloroplast genome and nuclear-encoded microsatellites implicated a single Y. aloifolia genotype as the maternal parent of Y. gloriosa. Using whole genome resequencing, we assembled chloroplast genomes from 41 accessions of all three species to re-assess the hybrid origins of Y. gloriosa. We further used re-sequencing data to annotate transposon abundance in the three species and mRNA-seq to analyze transcription of transposons. The chloroplast phylogeny and haplotype analysis suggest multiple hybridization events contributing to the origin of Y. gloriosa, with both parental species acting as the maternal donor. Transposon abundance at the superfamily level was significantly different between the three species; the hybrid was frequently intermediate to the parental species in TE superfamily abundance or appeared more similar to one or the other parent. In only one case-Copia LTR transposons-did Y. gloriosa have a significantly higher abundance relative to either parent. Expression patterns across the three species showed little increased transcriptional activity of transposons, suggesting that either no transposon release occurred in Y. gloriosa upon hybridization, or that any transposons that were activated via hybridization were rapidly silenced. The identification and quantification of transposon families paired with expression evidence paves the way for additional work seeking to link epigenetics with the important trait variation seen in this homoploid hybrid system.

Elucidating the Role of Virulence Traits in the Survival of Pathogenic E. coli PI-7 Following Disinfection

Reuse and discharge of treated wastewater can result in dissemination of microorganisms into the environment. Deployment of disinfection strategies is typically proposed as a last stage remediation effort to further inactivate viable microorganisms. In this study, we hypothesize that virulence traits, including biofilm formation, motility, siderophore, and curli production along with the capability to internalize into mammalian cells play a role in survival against disinfectants. Pathogenic E. coli PI-7 strain was used as a model bacterium that was exposed to diverse disinfection strategies such as chlorination, UV and solar irradiation. To this end, we used a random transposon mutagenesis library screening approach to generate 14 mutants that exhibited varying levels of virulence traits. In these 14 isolated mutants, we observed that an increase in virulence traits such as biofilm formation, motility, curli production, and internalization capability, increased the inactivation half-lives of mutants compared to wild-type E. coli PI-7. In addition, oxidative stress response and EPS production contributed to lengthening the lag phase duration (defined as the time required for exposure to disinfectant prior to decay). However, traits related to siderophore production did not help with survival against the tested disinfection strategies. Taken together, the findings suggested that selected virulence traits facilitate survival of pathogenic E. coli PI-7, which in turn could account for the selective enrichment of pathogens over the non-pathogenic ones after wastewater treatment. Further, the study also reflected on the effectiveness of UV as a more viable disinfection strategy for inactivation of pathogens.

Impact of Transposable Elements on Methylation and Gene Expression across Natural Accessions of Brachypodium distachyon

Transposable elements (TEs) constitute a large fraction of plant genomes and are mostly present in a transcriptionally silent state through repressive epigenetic modifications, such as DNA methylation. TE silencing is believed to influence the regulation of adjacent genes, possibly as DNA methylation spreads away from the TE. Whether this is a general principle or a context-dependent phenomenon is still under debate, pressing for studying the relationship between TEs, DNA methylation, and nearby gene expression in additional plant species. Here, we used the grass Brachypodium distachyon as a model and produced DNA methylation and transcriptome profiles for 11 natural accessions. In contrast to what is observed in Arabidopsis thaliana, we found that TEs have a limited impact on methylation spreading and that only few TE families are associated with a low expression of their adjacent genes. Interestingly, we found that a subset of TE insertion polymorphisms is associated with differential gene expression across accessions. Thus, although not having a global impact on gene expression, distinct TE insertions may contribute to specific gene expression patterns in B. distachyon.

Identification of Differentially Methylated miRNA Genes During Compatible and Incompatible Interactions Between Soybean and Soybean Cyst Nematode

DNA methylation is a widespread epigenetic mark that affects gene expression and transposon mobility during plant development and stress responses. However, the role of DNA methylation in regulating the expression of microRNA (miRNA) genes remains largely unexplored. Here, we analyzed DNA methylation changes of miRNA genes using a pair of soybean (Glycine max) near-isogenic lines (NILs) differing in their response to soybean cyst nematode (SCN; Heterodera glycines). Differences in global DNA methylation levels over miRNA genes in response to SCN infection were observed between the isogenic lines. miRNA genes with significant changes in DNA methylation levels in the promoter and primary transcript-coding regions were detected in both lines. In the susceptible isogenic line (NIL-S), 82 differentially methylated miRNAs were identified in response to SCN infection whereas, in the resistant isogenic line (NIL-R), only 16 differentially methylated miRNAs were identified. Interestingly, gma-miR5032, gma-miR5043, gma-miR1520b, and gma-2107-ch16 showed opposite methylation patterns in the isogenic lines. In addition, the miRNA paralogs gma-miR5770a and gma-miR5770b showed hypermethylation and hypomethylation in NIL-S and NIL-R, respectively. Gene expression quantification of gma-miR5032, gma-miR5043, gma-miR1520b, and gma-miR5770a/b and their confirmed targets indicated a role of DNA methylation in regulating miRNA expression and, thus, their targets upon SCN infection. Furthermore, overexpression of these four miRNAs in NIL-S using transgenic hairy root system enhanced plant resistance to SCN to various degrees with a key role observed for miR5032. Together, our results provide new insights into the role of epigenetic mechanisms in controlling miRNA regulatory function during SCN-soybean interactions.[Formula: see text] Copyright © 2020 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

MSH1-induced heritable enhanced growth vigor through grafting is associated with the RdDM pathway in plants

Plants transmit signals long distances, as evidenced in grafting experiments that create distinct rootstock-scion junctions. Noncoding small RNA is a signaling molecule that is graft transmissible, participating in RNA-directed DNA methylation; but the meiotic transmissibility of graft-mediated epigenetic changes remains unclear. Here, we exploit the MSH1 system in Arabidopsis and tomato to introduce rootstock epigenetic variation to grafting experiments. Introducing mutations dcl2, dcl3 and dcl4 to the msh1 rootstock disrupts siRNA production and reveals RdDM targets of methylation repatterning. Progeny from grafting experiments show enhanced growth vigor relative to controls. This heritable enhancement-through-grafting phenotype is RdDM-dependent, involving 1380 differentially methylated genes, many within auxin-related gene pathways. Growth vigor is associated with robust root growth of msh1 graft progeny, a phenotype associated with auxin transport based on inhibitor assays. Large-scale field experiments show msh1 grafting effects on tomato plant performance, heritable over five generations, demonstrating the agricultural potential of epigenetic variation.

The meiotic transmissibility and progeny phenotypic influence of graft-mediated epigenetic changes remain unclear. Here, the authors use the msh1 mutant in the rootstock to trigger heritable enhanced growth vigor in Arabidopsis and tomato, and show it is associated with the RNA-directed DNA methylation pathway.

Genomic re-assessment of the transposable element landscape of the potato genome

We provide a comprehensive and reliable potato TE landscape, based on a wide variety of identification tools and integrative approaches, producing clear and ready-to-use outputs for the scientific community. Transposable elements (TEs) are DNA sequences with the ability to autoreplicate and move throughout the host genome. TEs are major drivers in stress response and genome evolution. Given their significance, the development of clear and efficient TE annotation pipelines has become essential for many species. The latest de novo TE discovery tools, along with available TEs from Repbase and sRNA-seq data, allowed us to perform a reliable potato TEs detection, classification and annotation through an open-source and freely available pipeline ( https://github.com/DiegoZavallo/TE_Discovery ). Using a variety of tools, approaches and rules, we were able to provide a clearly annotated of characterized TEs landscape. Additionally, we described the distribution of the different types of TEs across the genome, where LTRs and MITEs present a clear clustering pattern in pericentromeric and subtelomeric/telomeric regions respectively. Finally, we analyzed the insertion age and distribution of LTR retrotransposon families which display a distinct pattern between the two major superfamilies. While older Gypsy elements concentrated around heterochromatic regions, younger Copia elements located predominantly on euchromatic regions. Overall, we delivered not only a reliable, ready-to-use potato TE annotation files, but also all the necessary steps to perform de novo detection for other species.

Epigenetic Mechanisms in Nematode-Plant Interactions

Epigenetic mechanisms play fundamental roles in regulating numerous biological processes in various developmental and environmental contexts. Three highly interconnected epigenetic control mechanisms, including small noncoding RNAs, DNA methylation, and histone modifications, contribute to the establishment of plant epigenetic profiles. During the past decade, a growing body of experimental work has revealed the intricate, diverse, and dynamic roles that epigenetic modifications play in plant-nematode interactions. In this review, I summarize recent progress regarding the functions of small RNAs in mediating plant responses to infection by cyst and root-knot nematodes, with a focus on the functions of microRNAs. I also recapitulate recent advances in genome-wide DNA methylation analysis and discuss how cyst nematodes induce extensive and dynamic changes in the plant methylome that impact the transcriptional activity of genes and transposable elements. Finally, the potential role of nematode effector proteins in triggering such epigenome changes is discussed.

DNA methylation reprogramming during seed development and its functional relevance in seed size/weight determination in chickpea

Seed development is orchestrated via complex gene regulatory networks and pathways. Epigenetic factors may also govern seed development and seed size/weight. Here, we analyzed DNA methylation in a large-seeded chickpea cultivar (JGK 3) during seed development stages. Progressive gain of CHH context DNA methylation in transposable elements (TEs) and higher frequency of small RNAs in hypermethylated TEs during seed development suggested a role of the RNA-dependent DNA methylation pathway. Frequency of intragenic TEs was higher in CHH context differentially methylated region (DMR) associated differentially expressed genes (DEGs). CG context hyper/hypomethylation within the gene body was observed for most of DMR-associated DEGs in JGK 3 as compared to small-seeded chickpea cultivar (Himchana 1). We identified candidate genes involved in seed size/weight determination exhibiting CG context hypermethylation within the gene body and higher expression in JGK 3. This study provides insights into the role of DNA methylation in seed development and seed size/weight determination in chickpea.

Rajkumar et al. report progressive gain of CHH context DNA methylation in transposable elements during seed development in chickpea, of which hypermethylation is associated with small RNAs. The candidate genes that determine seed size/weight in chickpea show CG context hypermethylation in the gene body and higher expression in large-seeded cultivar.

AthCNV: A Map of DNA Copy Number Variations in the Arabidopsis Genome

Copy number variations (CNVs) greatly contribute to intraspecies genetic polymorphism and phenotypic diversity. Recent analyses of sequencing data for >1000 Arabidopsis (Arabidopsis thaliana) accessions focused on small variations and did not include CNVs. Here, we performed genome-wide analysis and identified large indels (50 to 499 bp) and CNVs (500 bp and larger) in these accessions. The CNVs fully overlap with 18.3% of protein-coding genes, with enrichment for evolutionarily young genes and genes involved in stress and defense. By combining analysis of both genes and transposable elements (TEs) affected by CNVs, we revealed that the variation statuses of genes and TEs are tightly linked and jointly contribute to the unequal distribution of these elements in the genome. We also determined the gene copy numbers in a set of 1060 accessions and experimentally validated the accuracy of our predictions by multiplex ligation-dependent probe amplification assays. We then successfully used the CNVs as markers to analyze population structure and migration patterns. Finally, we examined the impact of gene dosage variation triggered by a CNV spanning the SEC10 gene on SEC10 expression at both the transcript and protein levels. The catalog of CNVs, CNV-overlapping genes, and their genotypes in a top model dicot will stimulate the exploration of the genetic basis of phenotypic variation.

Inter-retrotransposon amplified polymorphism markers revealed long terminal repeat retrotransposon insertion polymorphism in flax cultivated on the experimental fields around Chernobyl

Ionizing radiation in environment comes from various natural and anthropogenic sources. The effect of radioactivity released after the CNPP (Chernobyl Nuclear Power Plant) on plant systems remains of great interest. Even now, more than three decades after the nuclear accident, the long-lived radionuclides represent a strong stress factor. Herein, the emphasis has been placed on analysis of genetic variability represented by activation of LTR (Long Terminal Repeat)-retrotransposons. Polymorphism in LTR-retrotransposon insertions has been investigated throughout the genome of two flax varieties, Kyivskyi and Bethune. For this purpose, two retrotransposon-based marker techniques, IRAP (Inter-Retrotransposon Amplified Polymorphism) and iPBS (inter-Primer Binding Site), have been employed. The hypothesis that chronic radioactive stress may induce mechanism of retransposition has been supported by the activation of FL9, FL11 and FL12 LTR-retrotransposons in flax seeds harvested from radioactive environment. Out of two retrotransposon-based approaches, IRAP appears to be more suitable for identification of LTR-retrotransposon polymorphism. Even though the LTR-retrotransposon polymorphism was identified, the results suggest the high level of plant adaptation in the radioactive Chernobyl area. However, it is not really surprising that plants developed an effective strategy to survive in radio-contaminated environment over the past 30 years.

Transposable Elements Are Important Contributors to Standing Variation in Gene Expression in Capsella Grandiflora

Transposable elements (TEs) make up a significant portion of eukaryotic genomes and are important drivers of genome evolution. However, the extent to which TEs affect gene expression variation on a genome-wide scale in comparison with other types of variants is still unclear. We characterized TE insertion polymorphisms and their association with gene expression in 124 whole-genome sequences from a single population of Capsella grandiflora, and contrasted this with the effects of single nucleotide polymorphisms (SNPs). Population frequency of insertions was negatively correlated with distance to genes, as well as density of conserved noncoding elements, suggesting that the negative effects of TEs on gene regulation are important in limiting their abundance. Rare TE variants strongly influence gene expression variation, predominantly through downregulation. In contrast, rare SNPs contribute equally to up- and down-regulation, but have a weaker individual effect than TEs. An expression quantitative trait loci (eQTL) analysis shows that a greater proportion of common TEs are eQTLs as opposed to common SNPs, and a third of the genes with TE eQTLs do not have SNP eQTLs. In contrast with rare TE insertions, common insertions are more likely to increase expression, consistent with recent models of cis-regulatory evolution favoring enhancer alleles. Taken together, these results imply that TEs are a significant contributor to gene expression variation and are individually more likely than rare SNPs to cause extreme changes in gene expression.

DNA Methylation Readers in Plants

In plants, DNA methylation occurs in distinct sequence contexts, including CG, CHG, and CHH. Thus, plants have developed a surprisingly diverse set of DNA methylation readers to cope with an extended repertoire of methylated sites. The Arabidopsis genome contains twelve Methyl-Binding Domain proteins (MBD), and nine SET and RING finger-associated (SRA) domain containing proteins belonging to the SUVH clade, in addition to three homologs of UHRF1, namely VIM1-3, all containing SRA domains. In this review, we will highlight several research questions that remain unresolved with respect to the function of plant DNA methylation readers, which can have both de novo demethylase and maintenance activity. We argue that maintenance of CG methylation in plants likely involved actors not found in their mammalian counterparts, and that new evidence suggests significant reprogramming of DNA methylation during plant reproduction as an important new development in the field.

SINE Retrotransposon variation drives Ecotypic disparity in natural populations of Coilia nasus

Background

SINEs are a type of nonautonomous retrotransposon that can transpose from one site to be integrated elsewhere in an organism genome. SINE insertion can give rise to genetic variants and regulate gene expression, allowing organisms to acquire new adaptive capacity. Studies on this subject have focused on the impacts of SINEs on genes. However, ecological disparities in fish have not yet been explained by SINEs.

Results

New SINEs were isolated from Coilia nasus, which has two ecotypes—migratory and resident—that differ in their spawning and migration behaviors. The SINEs possess two structures that resemble a tRNA gene and a LINE retrotransposon tail. Comparison of olfactory tissue transcriptomes, intact SINE transcript copies were detected in only the migratory fish at the initial retrotransposition stage. The SINE DNA copy numbers were higher in the resident type than in the migratory type, while the frequency of SINE insertion was higher in the migratory type than in the resident type. Furthermore, SINE insertions can lead to new repeats of short DNA fragments in the genome, along with target site duplications. SINEs in the resident type have undergone excision via a mechanism in which predicted cleavage sites are formed by mutations, resulting in gaps that are then filled by microsatellites via microhomology-induced replication.

Conclusions

Notably, SINEs in the resident type have undergone strong natural selection, causing genomic heteroplasmy and driving ecological diversity of C. nasus. Our results reveal possible evolutionary mechanisms underlying the ecological diversity at the interface between SINE mobilization and organism defense.

Population Genomics of Transitions to Selfing in Brassicaceae Model Systems

Many plants harbor complex mechanisms that promote outcrossing and efficient pollen transfer. These include floral adaptations as well as genetic mechanisms, such as molecular self-incompatibility (SI) systems. The maintenance of such systems over long evolutionary timescales suggests that outcrossing is favorable over a broad range of conditions. Conversely, SI has repeatedly been lost, often in association with transitions to self-fertilization (selfing). This transition is favored when the short-term advantages of selfing outweigh the costs, primarily inbreeding depression. The transition to selfing is expected to have major effects on population genetic variation and adaptive potential, as well as on genome evolution. In the Brassicaceae, many studies on the population genetic, gene regulatory, and genomic effects of selfing have centered on the model plant Arabidopsis thaliana and the crucifer genus Capsella. The accumulation of population genomics datasets have allowed detailed investigation of where, when and how the transition to selfing occurred. Future studies will take advantage of the development of population genetics theory on the impact of selfing, especially regarding positive selection. Furthermore, investigation of systems including recent transitions to selfing, mixed mating populations and/or multiple independent replicates of the same transition will facilitate dissecting the effects of mating system variation from processes driven by demography.

Elimination of a Retrotransposon for Quenching Genome Instability in Modern Rice

Transposable elements (TEs) constitute the most abundant portions of plant genomes and can dramatically shape host genomes during plant evolution. They also play important roles in crop domestication. However, whether TEs themselves are also selected during crop domestication has remained unknown. Here, we identify an active long terminal repeat (LTR) retrotransposon, HUO, as a potential target of selection during rice domestication and breeding. HUO is a low-copy-number LTR retrotransposon, and is active under natural growth conditions and transmitted through male gametogenesis, preferentially inserting into genomic regions capable of transcription. HUO exists in all wild rice accessions and about half of the archaeological rice grains (1200-7000 years ago) and landraces surveyed, but is absent in almost all modern varieties, indicating its gradual elimination during rice domestication and breeding. Further analyses showed that HUO is subjected to strict gene silencing through the RNA-directed DNA methylation pathway. Our results also suggest that multiple HUO copies may trigger genomic instability through altering genome-wide DNA methylation and small RNA biogenesis and changing global gene expression, resulting in decreased disease resistance and yield, coinciding with its elimination during rice breeding. Together, our study suggests that negative selection of an active retrotransposon might be important for genome stability during crop domestication and breeding.

Identification of an active miniature inverted-repeat transposable element mJing in rice

Miniature inverted-repeat transposable elements (MITEs) are structurally homogeneous non-autonomous DNA transposons with high copy numbers that play important roles in genome evolution and diversification. Here, we analyzed the rice high-tillering dwarf (htd) mutant in an advanced backcross population between cultivated and wild rice, and identified an active MITE named miniature Jing (mJing). The mJing element belongs to the PIF/Harbinger superfamily. japonica rice var. Nipponbare and indica var. 93-11 harbor 72 and 79 mJing family members, respectively, have undergone multiple rounds of amplification bursts during the evolution of Asian cultivated rice (Oryza sativa L.). A heterologous transposition experiment in Arabidopsis thaliana indicated that the autonomous element Jing is likely to have provides the transposase needed for mJing mobilization. We identified 297 mJing insertion sites and their presence/absence polymorphism among 71 rice samples through targeted high-throughput sequencing. The results showed that the copy number of mJing varies dramatically among Asian cultivated rice (O. sativa), its wild ancestor (O. rufipogon), and African cultivated rice (O. glaberrima) and that some mJing insertions are subject to directional selection. These findings suggest that the amplification and removal of mJing elements have played an important role in rice genome evolution and species diversification.

Evolutionary Epigenomics of Retrotransposon-Mediated Methylation Spreading in Rice

Plant genomes contain numerous transposable elements (TEs), and many hypotheses on the evolutionary drivers that restrict TE activity have been postulated. Few models, however, have focused on the evolutionary epigenomic interaction between the plant host and its TE. The host genome recruits epigenetic factors, such as methylation, to silence TEs but methylation can spread beyond the TE sequence and influence the expression of nearby host genes. In this study, we investigated this epigenetic trade-off between TE and proximal host gene silencing by studying the epigenomic regulation of repressing long terminal repeat (LTR) retrotransposons (RTs) in Oryza sativa. Results showed significant evidence of methylation spreading originating from the LTR-RT sequences, and the extent of spreading was dependent on five factors: 1) LTR-RT family, 2) time since the LTR-RT insertion, 3) recombination rate of the LTR-RT region, 4) level of LTR-RT sequence methylation, and 5) chromosomal location. Methylation spreading had negative effects by reducing host gene expression, but only on host genes with LTR-RT inserted in its introns. Our results also suggested high levels of LTR-RT methylation might have a role in suppressing TE-mediated deleterious ectopic recombination. In the end, despite the methylation spreading, no strong epigenetic trade-off was detected and majority of LTR-RT may have only minor epigenetic effects on nearby host genes.

A DNA methylation reader complex that enhances gene transcription

DNA methylation generally functions as a repressive transcriptional signal, but it is also known to activate gene expression. In either case, the downstream factors remain largely unknown. By using comparative interactomics, we isolated proteins in Arabidopsis thaliana that associate with methylated DNA. Two SU(VAR)3-9 homologs, the transcriptional antisilencing factor SUVH1, and SUVH3, were among the methyl reader candidates. SUVH1 and SUVH3 bound methylated DNA in vitro, were associated with euchromatic methylation in vivo, and formed a complex with two DNAJ domain-containing homologs, DNAJ1 and DNAJ2. Ectopic recruitment of DNAJ1 enhanced gene transcription in plants, yeast, and mammals. Thus, the SUVH proteins bind to methylated DNA and recruit the DNAJ proteins to enhance proximal gene expression, thereby counteracting the repressive effects of transposon insertion near genes.

Root-Specific Expression of a Jacalin Lectin Family Protein Gene Requires a Transposable Element Sequence in the Promoter

Transposable elements (TEs) are widespread in the plant genome and can impact on the expression of neighbouring genes. Our previous studies have identified a number of DNA demethylase-regulated defence-related genes that contain TE sequences in the promoter and show tissue-specific expression in Arabidopsis. In this study we investigated the role of the promoter TE insertions in the root-specific expression of a DNA demethylase-regulated gene, AT5G38550, encoding a Jacalin lectin family protein. Using a promoter:GUS fusion reporter gene approach, we first demonstrated that the full-length promoter fragment, carrying four TE sequences, contained the essential regulatory information required for root-specific expression and DNA demethylase regulation in Arabidopsis. By successive deletion of the four TE sequences, we showed that one of the four TE insertions, a 201-bp TE fragment of the hAT DNA transposon family, was required for root-specific expression: Deletion of this TE, but not the first two TE sequences, converted the root-specific expression pattern to a constitutive expression pattern in Arabidopsis plants. Our study provides an example indicating an important role of TE insertions in tissue-specific expression of plant defence-related genes.

Maintenance of grafting-induced epigenetic variations in the asexual progeny of Brassica oleracea and B. juncea chimera

Grafting-induced variations have been observed in many plant species, but the heritability of variation in progeny is not well understood. In our study, adventitious shoots from the C cell lineage of shoot apical meristem (SAM) grafting chimera TCC (where the origin of the outmost, middle and innermost cell layers, respectively, of SAM is designated by 'T' for tuber mustard and 'C' for red cabbage) were induced and identified as r-CCC (r = regenerated). To investigate the maintenance of grafting variations during cell propagation and regeneration, different generations of asexual progeny (r-CCCn, n = generation) were established through successive regeneration of axillary shoots from r-CCC. The fourth generation of r-CCC (r-CCC4) was selected to perform whole genome bisulfite sequencing for comparative analysis of hetero-grafting-induced global methylation changes relative to r-s-CCC4 (s = self-grafting). Increased CHH methylation levels and proportions were observed in r-CCC4, with substantial changes occurring in the repeat elements. Small RNA sequencing revealed 1135 specific small interfering RNA (siRNA) tags that were typically expressed in r-CCC, r-CCC2 and r-CCC4. Notably, 65% of these specific siRNAs were associated with repeat elements, termed RE siRNAs. Subsequent analysis revealed that the CHH methylation of RE siRNA-overlapping regions was mainly hypermethylation in r-CCC4, indicating that they were responsible for directing and maintaining grafting-induced CHH methylation. Moreover, the expression of 13 differentially methylated genes (DMGs) correlated with the phenotypic variation, showing differential expression levels between r-CCC4 and r-s-CCC4. These DMGs were predominantly CG hypermethylated, their methylation modifications corresponded to the transcription of relative methyltransferase.

Shaping Plant Adaptability, Genome Structure and Gene Expression through Transposable Element Epigenetic Control: Focus on Methylation

In plants, transposable elements (TEs) represent a large fraction of the genome, with potential to alter gene expression and produce genomic rearrangements. Epigenetic control of TEs is often used to stop unrestricted movement of TEs that would result in detrimental effects due to insertion in essential genes. The current review focuses on the effects of methylation on TEs and their genomic context, and how this type of epigenetic control affects plant adaptability when plants are faced with different stresses and changes. TEs mobilize in response to stress elicitors, including biotic and abiotic cues, but also developmental transitions and ‘genome shock’ events like polyploidization. These events transitionally lift TE repression, allowing TEs to move to new genomic locations. When TEs fall close to genes, silencing through methylation can spread to nearby genes, resulting in lower gene expression. The presence of TEs in gene promoter regions can also confer stress inducibility modulated through alternative methylation and demethylation of the TE. Bursts of transposition triggered by events of genomic shock can increase genome size and account for differences seen during polyploidization or species divergence. Finally, TEs have evolved several mechanisms to suppress their own repression, including the use of microRNAs to control genes that promote methylation. The interplay between silencing, transient TE activation, and purifying selection allows the genome to use TEs as a reservoir of potential beneficial modifications but also keeps TEs under control to stop uncontrolled detrimental transposition.

Stress affects the epigenetic marks added by natural transposable element insertions in Drosophila melanogaster

Transposable elements are emerging as an important source of cis-acting regulatory sequences and epigenetic marks that could influence gene expression. However, few studies have dissected the role of specific transposable element insertions on epigenetic gene regulation. Bari-Jheh is a natural transposon that mediates resistance to oxidative stress by adding cis-regulatory sequences that affect expression of nearby genes. In this work, we integrated publicly available ChIP-seq and piRNA data with chromatin immunoprecipitation experiments to get a more comprehensive picture of Bari-Jheh molecular effects. We showed that Bari-Jheh was enriched for H3K9me3 in nonstress conditions, and for H3K9me3, H3K4me3 and H3K27me3 in oxidative stress conditions, which is consistent with expression changes in adjacent genes. We further showed that under oxidative stress conditions, H3K4me3 and H3K9me3 spread to the promoter region of Jheh1 gene. Finally, another insertion of the Bari1 family was associated with increased H3K27me3 in oxidative stress conditions suggesting that Bari1 histone marks are copy-specific. We concluded that besides adding cis-regulatory sequences, Bari-Jheh influences gene expression by affecting the local chromatin state.

Transcriptional control and exploitation of an immune-responsive family of plant retrotransposons

Mobilization of transposable elements (TEs) in plants has been recognized as a driving force of evolution and adaptation, in particular by providing genes with regulatory modules that impact their transcription. In this study, we employed an ATCOPIA93 long-terminal repeat (LTR) promoter-GUS fusion to show that this retrotransposon behaves like an immune-responsive gene during pathogen defense in Arabidopsis We also showed that the endogenous ATCOPIA93 copy "EVD", which is activated in the presence of bacterial stress, is negatively regulated by both DNA methylation and polycomb-mediated silencing, a mode of repression typically found at protein-coding and microRNA genes. Interestingly, an ATCOPIA93-derived soloLTR is located upstream of the disease resistance gene RPP4 and is devoid of DNA methylation and H3K27m3 marks. Through loss-of-function experiments, we demonstrate that this soloLTR is required for the proper expression of RPP4 during plant defense, thus linking the responsiveness of ATCOPIA93 to biotic stress and the co-option of its LTR for plant immunity.

Dispersion Profiles and Gene Associations of Repetitive DNAs in the Euchromatin of the Beetle Tribolium castaneum

Satellite DNAs are tandemly repeated sequences clustered within heterochromatin. However, in some cases, such as the major TCAST1 satellite DNA from the beetle Tribolium castaneum, they are found partially dispersed within euchromatin. Such organization together with transcriptional activity enables TCAST1 to modulate the activity of neighboring genes. In order to explore if other T. castaneum repetitive families have features that could provide them with a possible gene-modulatory role, we compare here the structure, organization, dispersion profiles, and transcription activity of 10 distinct TCAST repetitive families including TCAST1. The genome organization of TCAST families exhibit either satellite-like or transposon-like characteristics. In addition to heterochromatin localization, bioinformatic searches of the assembled genome have revealed dispersion of all families within euchromatin, preferentially in the form of single repeats. Dispersed TCAST repeats are mutually correlated in distribution and are grouped in distinct regions of euchromatin. The repeats are associated with genes, are enriched in introns relative to intergenic regions, and very rarely overlap exons. In spite of the different mechanisms of repeat proliferation, such as transposition and homologous recombination, all TCAST families share a similar frequency of spreading as well as dispersion and gene association profiles. Additionally, TCAST families are transcribed and their transcription is significantly activated by heat stress. A possibility that such common features of TCAST families might be related to their potential gene-modulatory role is discussed.

Small RNAs Reflect Grandparental Environments in Apomictic Dandelion

Plants can show long-term effects of environmental stresses and in some cases a stress “memory” has been reported to persist across generations, potentially mediated by epigenetic mechanisms. However, few documented cases exist of transgenerational effects that persist for multiple generations and it remains unclear if or how epigenetic mechanisms are involved. Here, we show that the composition of small regulatory RNAs in apomictic dandelion lineages reveals a footprint of drought stress and salicylic acid treatment experienced two generations ago. Overall proportions of 21 and 24 nt RNA pools were shifted due to grandparental treatments. While individual genes did not show strong up- or downregulation of associated sRNAs, the subset of genes that showed the strongest shifts in sRNA abundance was significantly enriched for several GO terms including stress-specific functions. This suggests that a stress-induced signal was transmitted across multiple unexposed generations leading to persistent changes in epigenetic gene regulation.

Stress response regulation by epigenetic mechanisms: changing of the guards

Plants are sessile organisms that lack a specialized immune system to cope with biotic and abiotic stress. Instead, plants have complex regulatory networks that determine the appropriate distribution of resources between the developmental and the defense programs. In the last years, epigenetic regulation of repeats and gene expression has evolved as an important player in the transcriptional regulation of stress-related genes. Here, we review the current knowledge about how different stresses interact with different levels of epigenetic control of the genome. Moreover, we analyze the different examples of transgenerational epigenetic inheritance and connect them with the known features of genome epigenetic regulation. Although yet to be explored, the interplay between epigenetics and stress resistance seems to be a relevant and dynamic player of the interaction of plants with their environments.

Analysis of Transposable Elements in Coccidioides Species

Coccidioides immitis and C. posadasii are primary pathogenic fungi that cause disease in immunologically-normal animals and people. The organism is found exclusively in arid regions of the Southwestern United States, Mexico, and South America, but not in other parts of the world. This study is a detailed analysis of the transposable elements (TE) in Coccidioides spp. As is common in most fungi, Class I and Class II transposons were identified and the LTR Gypsy superfamily is the most common. The minority of Coccidioides Gypsy transposons contained regions highly homologous to polyprotein domains. Phylogenetic analysis of the integrase and reverse transcriptase sequences revealed that many, but not all, of the Gypsy reverse transcriptase and integrase domains clustered by species suggesting extensive transposition after speciation of the two Coccidiodies spp. The TEs were clustered and the distribution is enriched for the ends on contigs. Analysis of gene expression data from C. immitis found that protein-coding genes within 1 kB of hAT or Gypsy TEs were poorly expressed. The expression of C. posadasii genes within 1 kB of Gypsy TEs was also significantly lower compared to all genes but the difference in expression was smaller than C. immitis. C. posadasii orthologs of C. immitis Gyspsy-associated genes were also likely to be TE-associated. In both C. immitis and C. posadasii the TEs were preferentially associated with genes annotated with protein kinase gene ontology terms. These observations suggest that TE may play a role in influencing gene expression in Coccidioides spp. Our hope is that these bioinformatic studies of the potential TE influence on expression and evolution of Coccidioides will prompt the development of testable hypotheses to better understand the role of TEs in the biology and gene regulation of Coccidioides spp.

Convergent adaptive evolution in marginal environments: unloading transposable elements as a common strategy among mangrove genomes

Several clades of mangrove trees independently invade the interface between land and sea at the margin of woody plant distribution. As phenotypic convergence among mangroves is common, the possibility of convergent adaptation in their genomes is quite intriguing. To study this molecular convergence, we sequenced multiple mangrove genomes. In this study, we focused on the evolution of transposable elements (TEs) in relation to the genome size evolution. TEs, generally considered genomic parasites, are the most common components of woody plant genomes. Analyzing the long terminal repeat-retrotransposon (LTR-RT) type of TE, we estimated their death rates by counting solo-LTRs and truncated elements. We found that all lineages of mangroves massively and convergently reduce TE loads in comparison to their nonmangrove relatives; as a consequence, genome size reduction happens independently in all six mangrove lineages; TE load reduction in mangroves can be attributed to the paucity of young elements; the rarity of young LTR-RTs is a consequence of fewer births rather than access death. In conclusion, mangrove genomes employ a convergent strategy of TE load reduction by suppressing element origination in their independent adaptation to a new environment.

Maize RNA PolIV affects the expression of genes with nearby TE insertions and has a genome-wide repressive impact on transcription

BACKGROUND

RNA-directed DNA methylation (RdDM) is a plant-specific epigenetic process that relies on the RNA polymerase IV (Pol IV) for the production of 24 nucleotide small interfering RNAs (siRNA) that guide the cytosine methylation and silencing of genes and transposons. Zea mays RPD1/RMR6 gene encodes the largest subunit of Pol IV and is required for normal plant development, paramutation, transcriptional repression of certain transposable elements (TEs) and transcriptional regulation of specific alleles.

RESULTS

In this study we applied a total RNA-Seq approach to compare the B73 and rpd1/rmr6 leaf transcriptomes. Although previous studies indicated that loss of siRNAs production in RdDM mutants provokes a strong loss of CHH DNA methylation but not massive gene or TEs transcriptional activation in both Arabidopsis and maize, our total RNA-Seq analysis of rpd1/rmr6 transcriptome reveals that loss of Pol IV activity causes a global increase in the transcribed fraction of the maize genome. Our results point to the genes with nearby TE insertions as being the most strongly affected by Pol IV-mediated gene silencing. TEs modulation of nearby gene expression is linked to alternative methylation profiles on gene flanking regions, and these profiles are strictly dependent on specific characteristics of the TE member inserted. Although Pol IV is essential for the biogenesis of siRNAs, the genes with associated siRNA loci are less affected by the pol IV mutation.

CONCLUSIONS

This deep and integrated analysis of gene expression, TEs distribution, smallRNA targeting and DNA methylation levels, reveals that loss of Pol IV activity globally affects genome regulation, pointing at TEs as modulator of nearby gene expression and indicating the existence of multiple level epigenetic silencing mechanisms. Our results also suggest a predominant role of the Pol IV-mediated RdDM pathway in genome dominance regulation, and subgenome stability and evolution in maize.

A protein complex regulates RNA processing of intronic heterochromatin-containing genes in Arabidopsis

In several eukaryotic organisms, heterochromatin (HC) in the introns of genes can regulate RNA processing, including polyadenylation, but the mechanism underlying this regulation is poorly understood. By promoting distal polyadenylation, the bromo-adjacent homology (BAH) domain-containing and RNA recognition motif-containing protein ASI1 and the H3K9me2-binding protein EDM2 are required for the expression of functional full-length transcripts of intronic HC-containing genes in Arabidopsis Here we report that ASI1 and EDM2 form a protein complex in vivo via a bridge protein, ASI1-Immunoprecipitated Protein 1 (AIPP1), which is another RNA recognition motif-containing protein. The complex also may contain the Pol II CTD phosphatase CPL2, the plant homeodomain-containing protein AIPP2, and another BAH domain protein, AIPP3. As is the case with dysfunction of ASI1 and EDM2, dysfunction of AIPP1 impedes the use of distal polyadenylation sites at tested intronic HC-containing genes, such as the histone demethylase gene IBM1, resulting in a lack of functional full-length transcripts. A mutation in AIPP1 causes silencing of the 35S-SUC2 transgene and genome-wide CHG hypermethylation at gene body regions, consistent with the lack of full-length functional IBM1 transcripts in the mutant. Interestingly, compared with asi1, edm2, and aipp1 mutations, mutations in CPL2, AIPP2, and AIPP3 cause the opposite effects on the expression of intronic HC-containing genes and other genes, suggesting that CPL2, AIPP2, and AIPP3 may form a distinct subcomplex. These results advance our understanding of the interplay between heterochromatic epigenetic modifications and RNA processing in higher eukaryotes.

Ectopic application of the repressive histone modification H3K9me2 establishes post-zygotic reproductive isolation in Arabidopsis thaliana

Hybrid seed lethality as a consequence of interspecies or interploidy hybridizations is a major mechanism of reproductive isolation in plants. This mechanism is manifested in the endosperm, a dosage-sensitive tissue supporting embryo growth. Deregulated expression of imprinted genes such as ADMETOS (ADM) underpin the interploidy hybridization barrier in Arabidopsis thaliana; however, the mechanisms of their action remained unknown. In this study, we show that ADM interacts with the AT hook domain protein AHL10 and the SET domain-containing SU(VAR)3-9 homolog SUVH9 and ectopically recruits the heterochromatic mark H3K9me2 to AT-rich transposable elements (TEs), causing deregulated expression of neighboring genes. Several hybrid incompatibility genes identified in Drosophila encode for dosage-sensitive heterochromatin-interacting proteins, which has led to the suggestion that hybrid incompatibilities evolve as a consequence of interspecies divergence of selfish DNA elements and their regulation. Our data show that imbalance of dosage-sensitive chromatin regulators underpins the barrier to interploidy hybridization in Arabidopsis, suggesting that reproductive isolation as a consequence of epigenetic regulation of TEs is a conserved feature in animals and plants.

LTR-retrotransposons in plants: Engines of evolution

LTR retrotransposons are the most abundant group of transposable elements (TEs) in plants. These elements can fall inside or close to genes, and therefore influence their expression and evolution. This review aims to examine how LTR retrotransposons, especially Ty1-copia elements, mediate gene regulation and evolution. Various stimuli, including polyploidization and biotic and abiotic elicitors, result in the transcription and movement of these retrotransposons, and can facilitate adaptation. The presence of cis-regulatory motifs in the LTRs are central to their stress-mediated responses and are shared with host stress-responsive genes, showing a complex evolutionary history in which TEs provide new regulatory units to genes. The presence of retrotransposon remnants in genes that are necessary for normal gene function, demonstrates the importance of exaptation and co-option, and is also a consequence of the abundance of these elements in plant genomes. Furthermore, insertions of LTR retrotransposons in and around genes provide potential for alternative splicing, epigenetic control, transduction, duplication and recombination. These characteristics can become an active part of the evolution of gene families as in the case of resistance genes (R-genes). The character of TEs as exclusively selfish is now being re-evaluated. Since genome-wide reprogramming via TEs is a long evolutionary process, the changes we can examine are case-specific and their fitness advantage may not be evident until TE-derived motifs and domains have been completely co-opted and fixed. Nevertheless, the presence of LTR retrotransposons inside genes and as part of gene promoter regions is consistent with their roles as engines of plant genome evolution.