Transposable element origins of epigenetic gene regulation

cis-Regulatory Elements in Plant Development, Adaptation, and Evolution

cis-Regulatory elements encode the genomic blueprints that ensure the proper spatiotemporal patterning of gene expression necessary for appropriate development and responses to the environment. Accumulating evidence implicates changes to gene expression as a major source of phenotypic novelty in eukaryotes, including acute phenotypes such as disease and cancer in mammals. Moreover, genetic and epigenetic variation affecting cis-regulatory sequences over longer evolutionary timescales has become a recurring theme in studies of morphological divergence and local adaptation. Here, we discuss the functions of and methods used to identify various classes of cis-regulatory elements, as well as their role in plant development and response to the environment. We highlight opportunities to exploit cis-regulatory variants underlying plant development and environmental responses for crop improvement efforts. Although a comprehensive understanding of cis-regulatory mechanisms in plants has lagged behind that in animals, we showcase several breakthrough findings that have profoundly influenced plant biology and shaped the overall understanding of transcriptional regulation in eukaryotes.

Recent transposable element bursts are associated with the proximity to genes in a fungal plant pathogen

The activity of transposable elements (TEs) contributes significantly to pathogen genome evolution. TEs often destabilize genome integrity but may also confer adaptive variation in pathogenicity or resistance traits. De-repression of epigenetically silenced TEs often initiates bursts of transposition activity that may be counteracted by purifying selection and genome defenses. However, how these forces interact to determine the expansion routes of TEs within a pathogen species remains largely unknown. Here, we analyzed a set of 19 telomere-to-telomere genomes of the fungal wheat pathogen Zymoseptoria tritici. Phylogenetic reconstruction and ancestral state estimates of individual TE families revealed that TEs have undergone distinct activation and repression periods resulting in highly uneven copy numbers between genomes of the same species. Most TEs are clustered in gene poor niches, indicating strong purifying selection against insertions near coding sequences, or as a consequence of insertion site preferences. TE families with high copy numbers have low sequence divergence and strong signatures of defense mechanisms (i.e., RIP). In contrast, small non-autonomous TEs (i.e., MITEs) are less impacted by defense mechanisms and are often located in close proximity to genes. Individual TE families have experienced multiple distinct burst events that generated many nearly identical copies. We found that a Copia element burst was initiated from recent copies inserted substantially closer to genes compared to older copies. Overall, TE bursts tended to initiate from copies in GC-rich niches that escaped inactivation by genomic defenses. Our work shows how specific genomic environments features provide triggers for TE proliferation in pathogen genomes.

Epigenetic targeting of transposon relics: beating the dead horses of the genome?

Transposable elements (TEs) have been seen as selfish genetic elements that can propagate in a host genome. Their propagation success is however hindered by a combination of mechanisms such as mutations, selection, and their epigenetic silencing by the host genome. As a result, most copies of TEs in a given genome are dead relics: their sequence is too degenerated to allow any transposition. Nevertheless, these TE relics often, but not always, remain epigenetically silenced, and if not to prevent transposition anymore, one can wonder the reason for this phenomenon. The mere self-perpetuating loop inherent to epigenetic silencing could alone explain that even when inactive, TE copies remain silenced. Beyond this process, nevertheless, antagonistic selective forces are likely to act on TE relic silencing. Especially, without the benefit of preventing transposition, TE relic silencing may prove deleterious to the host fitness, suggesting that the maintenance of TE relic silencing is the result of a fine, and perhaps case-by-case, evolutionary trade-off between beneficial and deleterious effects. Ultimately, the release of TE relics silencing may provide a 'safe' ground for adaptive epimutations to arise. In this review, we provide an overview of these questions in both plants and animals.

Characterization of Transposon-Derived Accessible Chromatin Regions in Rice (Oryza Sativa)

Growing evidence indicates that transposons or transposable elements (TEs)-derived accessible chromatin regions (ACRs) play essential roles in multiple biological processes by interacting with trans-acting factors. However, the function of TE-derived ACRs in the regulation of gene expression in the rice genome has not been well characterized. In this study, we examined the chromatin dynamics in six types of rice tissues and found that ~8% of ACRs were derived from TEs and exhibited distinct levels of accessibility and conservation as compared to those without TEs. TEs exhibited a TE subtype-dependent impact on ACR formation, which can be mediated by changes in the underlying DNA methylation levels. Moreover, we found that tissue-specific TE-derived ACRs might function in the tissue development through the modulation of nearby gene expression. Interestingly, many genes in domestication sweeps were found to overlap with TE-derived ACRs, suggesting their potential functions in the rice domestication. In addition, we found that the expression divergence of 1070 duplicate gene pairs were associated with TE-derived ACRs and had distinct distributions of TEs and ACRs around the transcription start sites (TSSs), which may experience different selection pressures. Thus, our study provides some insights into the biological implications of TE-derived ACRs in the rice genome. Our results imply that these ACRs are likely involved in the regulation of tissue development, rice domestication and functional divergence of duplicated genes.

An Eruption of LTR Retrotransposons in the Autopolyploid Genomes of Chrysanthemum nankingense (Asteraceae)

Whole genome duplication, associated with the induction of widespread genetic changes, has played an important role in the evolution of many plant taxa. All extant angiosperm species have undergone at least one polyploidization event, forming either an auto- or allopolyploid organism. Compared with allopolyploidization, however, few studies have examined autopolyploidization, and few studies have focused on the response of genetic changes to autopolyploidy. In the present study, newly synthesized C. nankingense autotetraploids (Asteraceae) were employed to characterize the genome shock following autopolyploidization. Available evidence suggested that the genetic changes primarily involved the loss of old fragments and the gain of novel fragments, and some novel sequences were potential long terminal repeat (LTR) retrotransposons. As Ty1-copia and Ty3-gypsy elements represent the two main superfamilies of LTR retrotransposons, the dynamics of Ty1-copia and Ty3-gypsy were evaluated using RT-PCR, transcriptome sequencing, and LTR retrotransposon-based molecular marker techniques. Additionally, fluorescence in situ hybridization(FISH)results suggest that autopolyploidization might also be accompanied by perturbations of LTR retrotransposons, and emergence retrotransposon insertions might show more rapid divergence, resulting in diploid-like behaviour, potentially accelerating the evolutionary process among progenies. Our results strongly suggest a need to expand the current evolutionary framework to include a genetic dimension when seeking to understand genomic shock following autopolyploidization in Asteraceae.

Silencing and anti-silencing mechanisms that shape the epigenome in plants

Epigenome information mediates genome function and maintenance by regulating gene expression and chromatin organization. Because the epigenome pattern can change in response to internal and external environments, it may underlie an adaptive genome response that modulates phenotypes during development and in changing environments. Here I summarize recent progress in our understanding of how epigenome patterns are shaped and modulated by concerted actions of silencing and anti-silencing factors mainly in Arabidopsis thaliana. I discuss the dynamic nature of epigenome regulation, which is realized by cooperation and counteraction among silencing and anti-silencing factors, and how the dynamic epigenome mediates robust and plastic responses of plants to fluctuating environments.

Patterns of stability and change in the maize genome: a case study of small RNA transcriptomes in two recombinant inbred lines and their progenitors

Small RNAs (sRNAs) are epigenetic regulators of eukaryotic genes and transposable elements (TEs). Diverse sRNA expression patterns exist within a species, but how this diversity arises is not well understood. To provide a window into the dynamics of maize sRNA patterning, sRNA and mRNA transcriptomes were examined in two related Zea mays recombinant inbred lines (RILs) and their inbred parents. Analysis of these RILs revealed that most clusters of sRNA expression retained the parental sRNA expression level. However, expression states that differ from the parental allele were also observed, predominantly reflecting decreases in sRNA expression. When RIL sRNA expression differed from the parental allele, the new state was frequently similar between the two RILs, and similar to the expression state found at the allele in the other parent. Novel sRNA expression patterns, distinct from those of either parent, were rare. Additionally, examination of sRNA expression over TEs revealed one TE family, Gyma, which showed consistent enrichment for RIL sRNA expression differences compared to those found in parental alleles. These findings provide insights into how sRNA silencing might evolve over generations and suggest that further investigation into the molecular nature of sRNA trans regulators is warranted.

Integrated Genomic Analyses From Low-Depth Sequencing Help Resolve Phylogenetic Incongruence in the Bamboos (Poaceae: Bambusoideae)

The bamboos (Bambusoideae, Poaceae) comprise a major grass lineage with a complex evolutionary history involving ancient hybridization and allopolyploidy. About 1700 described species are classified into three tribes, Olyreae (herbaceous bamboos), Bambuseae (tropical woody bamboos), and Arundinarieae (temperate woody bamboos). Nuclear analyses strongly support monophyly of the woody tribes, whereas plastome analyses strongly support paraphyly, with Bambuseae sister to Olyreae. Our objectives were to clarify the origin(s) of the woody bamboo tribes and resolve the nuclear vs. plastid conflict using genomic tools. For the first time, plastid and nuclear genomic information from the same bamboo species were combined in a single study. We sampled 51 species of bamboos representing the three tribes, estimated their genome sizes and generated low-depth sample sequence data, from which plastomes were assembled and nuclear repeats were analyzed. The distribution of repeat families was found to agree with nuclear gene phylogenies, but also provides novel insights into nuclear evolutionary history. We infer two early, independent hybridization events, one between an Olyreae ancestor and a woody ancestor giving rise to the two Bambuseae lineages, and another between two woody ancestors giving rise to the Arundinarieae. Retention of the Olyreae plastome associated with differential dominance of nuclear genomes and subsequent diploidization in some lineages explains the paraphyly observed in plastome phylogenetic estimations. We confirm ancient hybridization and allopolyploidy in the origins of the extant woody bamboo lineages and propose biased fractionation and diploidization as important factors in their evolution.

Transposon dynamics and the epigenetic switch hypothesis

The recent explosion of interest in epigenetics is often portrayed as the dawning of a scientific revolution that promises to transform biomedical science along with developmental and evolutionary biology. Much of this enthusiasm surrounds what we call the epigenetic switch hypothesis, which regards certain examples of epigenetic inheritance as an adaptive organismal response to environmental change. This interpretation overlooks an alternative explanation in terms of coevolutionary dynamics between parasitic transposons and the host genome. This raises a question about whether epigenetics researchers tend to overlook transposon dynamics more generally. To address this question, we surveyed a large sample of scientific publications on the topics of epigenetics and transposons over the past fifty years. We found that enthusiasm for epigenetics is often inversely related to interest in transposon dynamics across the four disciplines we examined. Most surprising was a declining interest in transposons within biomedical science and cellular and molecular biology over the past two decades. Also notable was a delayed and relatively muted enthusiasm for epigenetics within evolutionary biology. An analysis of scientific abstracts from the past twenty-five years further reveals systematic differences among disciplines in their uses of the term epigenetic, especially with respect to heritability commitments and functional interpretations. Taken together, these results paint a nuanced picture of the rise of epigenetics and the possible neglect of transposon dynamics, especially among biomedical scientists.

Identification of Active Transposable Elements in Plants: The Mobilome-Seq Approach

Transposable elements (TEs) are the main component of eukaryotic genomes. Besides their impact on genome size, TEs are also functionally important as they can alter gene expression and influence phenotypic variation. In plants, most top-down studies focus on extremely clear phenotypes such as the shape or the color of individuals and do not explore fully the role of TEs in evolution. Assessing the impact of TEs in a more systematic manner, however, requires identifying active TEs to further study their impact on phenotypes. In this chapter, we describe an in planta approach that consists in activating TEs by interfering with pathways involved in their silencing. It enables to directly investigate the functional impact of single TE families at low cost.

Chromatin Dynamics Contribute to the Spatiotemporal Expression Pattern of Virulence Genes in a Fungal Plant Pathogen

Dynamic changes in transcription profiles are key for the success of pathogens in colonizing their hosts. In many pathogens, genes associated with virulence, such as effector genes, are located in regions of the genome that are rich in transposable elements and heterochromatin. The contribution of chromatin modifications to gene expression in pathogens remains largely unknown. Using a combination of a reporter gene-based approach and chromatin immunoprecipitation, we show that the heterochromatic environment of effector genes in the fungal plant pathogen Zymoseptoria tritici is a key regulator of their specific spatiotemporal expression patterns. Enrichment in trimethylated lysine 27 of histone H3 dictates the repression of effector genes in the absence of the host. Chromatin decondensation during host colonization, featuring a reduction in this repressive modification, indicates a major role for epigenetics in effector gene induction. Our results illustrate that chromatin modifications triggered during host colonization determine the specific expression profile of effector genes at the cellular level and, hence, provide new insights into the regulation of virulence in fungal plant pathogens.IMPORTANCE Fungal plant pathogens possess a large repertoire of genes encoding putative effectors, which are crucial for infection. Many of these genes are expressed at low levels in the absence of the host but are strongly induced at specific stages of the infection. The mechanisms underlying this transcriptional reprogramming remain largely unknown. We investigated the role of the genomic environment and associated chromatin modifications of effector genes in controlling their expression pattern in the fungal wheat pathogen Zymoseptoria tritici Depending on their genomic location, effector genes are epigenetically repressed in the absence of the host and during the initial stages of infection. Derepression of effector genes occurs mainly during and after penetration of plant leaves and is associated with changes in histone modifications. Our work demonstrates the role of chromatin in shaping the expression of virulence components and, thereby, the interaction between fungal pathogens and their plant hosts.

Distinct Expression and Methylation Patterns for Genes with Different Fates following a Single Whole-Genome Duplication in Flowering Plants

For most sequenced flowering plants, multiple whole-genome duplications (WGDs) are found. Duplicated genes following WGD often have different fates that can quickly disappear again, be retained for long(er) periods, or subsequently undergo small-scale duplications. However, how different expression, epigenetic regulation, and functional constraints are associated with these different gene fates following a WGD still requires further investigation due to successive WGDs in angiosperms complicating the gene trajectories. In this study, we investigate lotus (Nelumbo nucifera), an angiosperm with a single WGD during the K-pg boundary. Based on improved intraspecific-synteny identification by a chromosome-level assembly, transcriptome, and bisulfite sequencing, we explore not only the fundamental distinctions in genomic features, expression, and methylation patterns of genes with different fates after a WGD but also the factors that shape post-WGD expression divergence and expression bias between duplicates. We found that after a WGD genes that returned to single copies show the highest levels and breadth of expression, gene body methylation, and intron numbers, whereas the long-retained duplicates exhibit the highest degrees of protein-protein interactions and protein lengths and the lowest methylation in gene flanking regions. For those long-retained duplicate pairs, the degree of expression divergence correlates with their sequence divergence, degree in protein-protein interactions, and expression level, whereas their biases in expression level reflecting subgenome dominance are associated with the bias of subgenome fractionation. Overall, our study on the paleopolyploid nature of lotus highlights the impact of different functional constraints on gene fate and duplicate divergence following a single WGD in plant.

Bioinformatics Analysis of Evolution and Human Disease Related Transposable Element-Derived microRNAs

Transposable element (TE) has the ability to insert into certain parts of the genome, and due to this event, it is possible for TEs to generate new factors and one of these factors are microRNAs (miRNA). miRNAs are non-coding RNAs made up of 19 to 24 nucleotides and numerous miRNAs are derived from TE. In this study, to support general knowledge on TE and miRNAs derived from TE, several bioinformatics tools and databases were used to analyze miRNAs derived from TE in two aspects: evolution and human disease. The distribution of TEs in diverse species presents that almost half of the genome is covered with TE in mammalians and less than a half in other vertebrates and invertebrates. Based on selected evolution-related miRNAs studies, a total of 51 miRNAs derived from TE were found and analyzed. For the human disease-related miRNAs, total of 34 miRNAs derived from TE were organized from the previous studies. In summary, abundant miRNAs derived from TE are found, however, the function of miRNAs derived from TE is not informed either. Therefore, this study provides theoretical understanding of miRNAs derived from TE by using various bioinformatics tools.

Screening diversity and distribution of Copia retrotransposons reveals a specific amplification of BARE1 elements in genomes of the polyploid Hordeum murinum complex

We explored diversity, distribution and evolutionary dynamics of Ty1-Copia retrotransposons in the genomes of the Hordeum murinum polyploid complex and related taxa. Phylogenetic and fluorescent in situ hybridization (FISH) analyses of reverse transcriptase sequences identified four Copia families in these genomes: the predominant BARE1 (including three groups or subfamilies, A, B and C), and the less represented RIRE1, IKYA and TAR-1. Within the BARE1 family, BARE1-A elements and a subgroup of BARE1-B elements (named B1) have proliferated in the allopolyploid members of the H. murinum complex (H. murinum and H. leporinum), and in their extant diploid progenitor, subsp. glaucum. Moreover, we found a specific amplification of BARE1-B elements within each Hordeum species surveyed. The low occurrence of RIRE1, IKYA and TAR-1 elements in the allopolyploid cytotypes suggests that they are either weakly represented or highly degenerated in their diploid progenitors. The results demonstrate that BARE1-A and BARE1-B1 Copia elements are particularly well represented in the genomes of the H. murinum complex and constitute its genomic hallmark. No BARE1-A and -B1 homologs were detected in the reference barley genome. The similar distribution of RT-Copia probes across chromosomes of diploid, tetraploid and hexaploid taxa of the murinum complex shows no evidence of proliferation following polyploidization.

Copy Number Variation of Transposable Elements in Thinopyrum intermedium and Its Diploid Relative Species

Diploid and polyploid wild species of Triticeae have complex relationships, and the understanding of their evolution and speciation could help to increase the usability of them in wheat breeding as a source of genetic diversity. The diploid species Pseudoroegneria spicata (St), Thinopyrum bessarabicum (Jb), Dasypyrum villosum (V) derived from a hypothetical common ancestor are considered to be possible subgenome donors in hexaploid species Th. intermedium (JrJvsSt, where indices r, v, and s stand for the partial relation to the genomes of Secale, Dasypyrum, and Pseudoroegneria, respectively). We quantified 10 families of transposable elements (TEs) in P. spicata, Th. bessarabicum, D. villosum (per one genome), and Th. intermedium (per one average subgenome) using the quantitative real time PCR assay and compared their abundance within the studied genomes as well as between them. Sabrina was the most abundant among all studied elements in P. spicata, D. villosum, and Th. intermedium, and among Ty3/Gypsy elements in all studied species. Among Ty1/Copia elements, Angela-A and WIS-A showed the highest and close abundance with the exception of D. villosum, and comprised the majority of all studied elements in Th. bessarabicum. Sabrina, BAGY2, and Angela-A showed similar abundance among diploids and in Th. intermedium hexaploid; Latidu and Barbara demonstrated sharp differences between diploid genomes. The relationships between genomes of Triticeae species based on the studied TE abundance and the role of TEs in speciation and polyploidization in the light of the current phylogenetic models is discussed.

Dynamic Patterns of Transcript Abundance of Transposable Element Families in Maize

Transposable Elements (TEs) are mobile elements that contribute the majority of DNA sequences in the maize genome. Due to their repetitive nature, genomic studies of TEs are complicated by the difficulty of properly attributing multi-mapped short reads to specific genomic loci. Here, we utilize a method to attribute RNA-seq reads to TE families rather than particular loci in order to characterize transcript abundance for TE families in the maize genome. We applied this method to assess per-family expression of transposable elements in >800 published RNA-seq libraries representing a range of maize development, genotypes, and hybrids. While a relatively small proportion of TE families are transcribed, expression is highly dynamic with most families exhibiting tissue-specific expression. A large number of TE families were specifically detected in pollen and endosperm, consistent with reproductive dynamics that maintain silencing of TEs in the germ line. We find that B73 transcript abundance is a poor predictor of TE expression in other genotypes and that transcript levels can differ even for shared TEs. Finally, by assessing recombinant inbred line and hybrid transcriptomes, complex patterns of TE transcript abundance across genotypes emerged. Taken together, this study reveals a dynamic contribution of TEs to maize transcriptomes.

Fine mapping and candidate gene analysis of a QTL associated with leaf rolling index on chromosome 4 of maize (Zea mays L.)

One QTL qLRI4 controlling leaf rolling index on chromosome 4 was finely mapped, and ZmOCL5, a member of the HD-Zip class IV genes, is likely a candidate. Leaf rolling is an important agronomic trait related to plant architecture that can change the light condition and photosynthetic efficiency of the population. Here, we isolated one EMS-induced mutant in Chang7-2 background with extreme abaxial rolling leaf, named abrl1. Histological analysis showed that the increased number and area of bulliform cells may contribute to abaxial rolling leaf in abrl1. The F₂ and F_2:3 populations derived from Wu9086 with flat leaves and abrl1 were developed to map abrl1. Non-Mendelian segregation of phenotypic variation was observed in these populations and five genomic regions controlling the leaf rolling index (LRI) were identified, which could be due to the phenotypic difference between Chang7-2 and Wu9086. Moreover, one major QTL qLRI4 on chromosome 4 was further validated and finely mapped to a genetic interval between InDel13 and InDel10, with a physical distance of approximately 277 kb using NIL populations, among which one 602-bp insertion was identified in the promoter region of HD-Zip class IV gene Zm00001d049443 (named as ZmOCL5) of abrl1 compared with wild-type Chang7-2. Remarkably, the 602-bp InDel was associated with LRI in an F₂ population developed by crossing abrl1 mutant and its wild-type. In addition, the 602-bp insertion increased ZmOCL5 promoter activity and expression. Haplotype analysis demonstrated that the 602-bp insertion was a rare mutation event. Taken together, we propose that the rolled leaf in the abrl1 mutant may be partially attributed to the 602-bp insertion, which may be an attractive target for the genetic improvement of LRI in maize.

Monitoring the interplay between transposable element families and DNA methylation in maize

DNA methylation and epigenetic silencing play important roles in the regulation of transposable elements (TEs) in many eukaryotic genomes. A majority of the maize genome is derived from TEs that can be classified into different orders and families based on their mechanism of transposition and sequence similarity, respectively. TEs themselves are highly methylated and it can be tempting to view them as a single uniform group. However, the analysis of DNA methylation profiles in flanking regions provides evidence for distinct groups of chromatin properties at different TE families. These differences among TE families are reproducible in different tissues and different inbred lines. TE families with varying levels of DNA methylation in flanking regions also show distinct patterns of chromatin accessibility and modifications within the TEs. The differences in the patterns of DNA methylation flanking TE families arise from a combination of non-random insertion preferences of TE families, changes in DNA methylation triggered by the insertion of the TE and subsequent selection pressure. A set of nearly 70,000 TE polymorphisms among four assembled maize genomes were used to monitor the level of DNA methylation at haplotypes with and without the TE insertions. In many cases, TE families with high levels of DNA methylation in flanking sequence are enriched for insertions into highly methylated regions. The majority of the >2,500 TE insertions into unmethylated regions result in changes in DNA methylation in haplotypes with the TE, suggesting the widespread potential for TE insertions to condition altered methylation in conserved regions of the genome. This study highlights the interplay between TEs and the methylome of a major crop species.

The rhg1-a (Rhg1 low-copy) nematode resistance source harbors a copia-family retrotransposon within the Rhg1-encoded α-SNAP gene

Soybean growers widely use the Resistance to Heterodera glycines 1 (Rhg1) locus to reduce yield losses caused by soybean cyst nematode (SCN). Rhg1 is a tandemly repeated four gene block. Two classes of SCN resistance-conferring Rhg1 haplotypes are recognized: rhg1-a ("Peking-type," low-copy number, three or fewer Rhg1 repeats) and rhg1-b ("PI 88788-type," high-copy number, four or more Rhg1 repeats). The rhg1-a and rhg1-b haplotypes encode α-SNAP (alpha-Soluble NSF Attachment Protein) variants α-SNAP Rhg1 LC and α-SNAP Rhg1 HC, respectively, with differing atypical C-terminal domains, that contribute to SCN resistance. Here we report that rhg1-a soybean accessions harbor a copia retrotransposon within their Rhg1 Glyma.18G022500 (α-SNAP-encoding) gene. We termed this retrotransposon "RAC," for Rhg1 alpha-SNAP copia. Soybean carries multiple RAC-like retrotransposon sequences. The Rhg1 RAC insertion is in the Glyma.18G022500 genes of all true rhg1-a haplotypes we tested and was not detected in any examined rhg1-b or Rhg1WT (single-copy) soybeans. RAC is an intact element residing within intron 1, anti-sense to the rhg1-a α-SNAP open reading frame. RAC has intrinsic promoter activities, but overt impacts of RAC on transgenic α-SNAP Rhg1 LC mRNA and protein abundance were not detected. From the native rhg1-a RAC+ genomic context, elevated α-SNAP Rhg1 LC protein abundance was observed in syncytium cells, as was previously observed for α-SNAP Rhg1 HC (whose rhg1-b does not carry RAC). Using a SoySNP50K SNP corresponding with RAC presence, just ~42% of USDA accessions bearing previously identified rhg1-a SoySNP50K SNP signatures harbor the RAC insertion. Subsequent analysis of several of these putative rhg1-a accessions lacking RAC revealed that none encoded α-SNAPRhg1LC, and thus, they are not rhg1-a. rhg1-a haplotypes are of rising interest, with Rhg4, for combating SCN populations that exhibit increased virulence against the widely used rhg1-b resistance. The present study reveals another unexpected structural feature of many Rhg1 loci, and a selectable feature that is predictive of rhg1-a haplotypes.

Evidence of developmental escape from transcriptional gene silencing in MESSI retrotransposons

Transposable elements (TEs) are ubiquitous genomic features. 'Copy-and-paste' long-terminal-repeat (LTR) retrotransposons have been particularly successful during evolution of the plant kingdom, representing a substantial proportion of genomes. For survival in copious numbers, these TEs may have evolved replicative mobilization strategies that circumvented hosts' epigenetic silencing. Stressful circumstances are known to trigger the majority of known mobilizing plant retrotransposons, leading to the idea that most are activated by environmental signals. However, previous research revealed that plant developmental programs include steps of silencing relaxation, suggesting that developmental signals may also be of importance for thriving parasitic elements. Here, we uncover an unusual family of giant LTR retrotransposons from the Solanum clade, named MESSI, with transcriptional competence in shoot apical meristems of tomato. Despite being recognized and targeted by the host epigenetic surveillance, this family is activated in specific meristematic areas fundamental for plant shoot development, which are involved in meristem formation and maintenance. Our work provides initial evidence that some retrotransposons may evolve developmentally associated escape strategies to overcome transcriptional gene silencing in vegetative tissues contributing to the host's next generation. This implies that not only environmental but also developmental signals could be exploited by selfish elements for survival within the plant kingdom.

Genome-wide variation in nucleotides and retrotransposons in alpine populations of Arabis alpina (Brassicaceae)

Advances in high-throughput sequencing have promoted the collection of reference genomes and genome-wide diversity. However, the assessment of genomic variation among populations has hitherto mainly been surveyed through single-nucleotide polymorphisms (SNPs) and largely ignored the often major fraction of genomes represented by transposable elements (TEs). Despite accumulating evidence supporting the evolutionary significance of TEs, comprehensive surveys remain scarce. Here, we sequenced the full genomes of 304 individuals of Arabis alpina sampled from four nearby natural populations to genotype SNPs as well as polymorphic long terminal repeat retrotransposons (polymorphic TEs; i.e., presence/absence of TE insertions at specific loci). We identified 291,396 SNPs and 20,548 polymorphic TEs, comparing their contributions to genomic diversity and divergence across populations. Few SNPs were shared among populations and overall showed high population-specific variation, whereas most polymorphic TEs segregated among populations. The genomic context of these two classes of variants further highlighted candidate adaptive loci having a putative impact on functional genes. In particular, 4.96% of the SNPs were identified as nonsynonymous or affecting start/stop codons. In contrast, 43% of the polymorphic TEs were present next to Arabis genes enriched in functional categories related to the regulation of reproduction and responses to biotic as well as abiotic stresses. This unprecedented data set, mapping variation gained from SNPs and complementary polymorphic TEs within and among populations, will serve as a rich resource for addressing microevolutionary processes shaping genome variation.

LtrDetector: A tool-suite for detecting long terminal repeat retrotransposons de-novo

BACKGROUND

Long terminal repeat retrotransposons are the most abundant transposons in plants. They play important roles in alternative splicing, recombination, gene regulation, and defense mechanisms. Large-scale sequencing projects for plant genomes are currently underway. Software tools are important for annotating long terminal repeat retrotransposons in these newly available genomes. However, the available tools are not very sensitive to known elements and perform inconsistently on different genomes. Some are hard to install or obsolete. They may struggle to process large plant genomes. None can be executed in parallel out of the box and very few have features to support visual review of new elements. To overcome these limitations, we developed LtrDetector, which uses techniques inspired by signal-processing.

RESULTS

We compared LtrDetector to LTR_Finder and LTRharvest, the two most successful predecessor tools, on six plant genomes. For each organism, we constructed a ground truth data set based on queries from a consensus sequence database. According to this evaluation, LtrDetector was the most sensitive tool, achieving 16-23% improvement in sensitivity over LTRharvest and 21% improvement over LTR_Finder. All three tools had low false positive rates, with LtrDetector achieving 98.2% precision, in between its two competitors. Overall, LtrDetector provides the best compromise between high sensitivity and low false positive rate while requiring moderate time and utilizing memory available on personal computers.

CONCLUSIONS

LtrDetector uses a novel methodology revolving around k-mer distributions, which allows it to produce high-quality results using relatively lightweight procedures. It is easy to install and use. It is not species specific, performing well using its default parameters on genomes of varying size and repeat content. It is automatically configured for parallel execution and runs efficiently on an ordinary personal computer. It includes a k-mer scores visualization tool to facilitate manual review of the identified elements. These features make LtrDetector an attractive tool for future annotation projects involving long terminal repeat retrotransposons.

LINE-1 hypomethylation in human hepatocellular carcinomas correlates with shorter overall survival and CIMP phenotype

Reactivation of interspersed repetitive sequences due to loss of methylation is associated with genomic instability, one of the hallmarks of cancer cells. LINE-1 hypomethylation is a surrogate marker for global methylation loss and is potentially a new diagnostic and prognostic biomarker in tumors. However, the correlation of LINE-1 hypomethylation with clinicopathological parameters and the CpG island methylator phenotype (CIMP) in patients with liver tumors is not yet well defined, particularly in Caucasians who show quite low rates of HCV/HBV infection and a higher incidence of liver steatosis. Therefore, quantitative DNA methylation analysis of LINE-1, RASSF1A, and CCND2 using pyrosequencing was performed in human hepatocellular carcinomas (HCC, n = 40), hepatocellular adenoma (HCA, n = 10), focal nodular hyperplasia (FNH, n = 5), and corresponding peritumoral liver tissues as well as healthy liver tissues (n = 5) from Caucasian patients. Methylation results were correlated with histopathological findings and clinical data. We found loss of LINE-1 DNA methylation only in HCC. It correlated significantly with poor survival (log rank test, p = 0.007). An inverse correlation was found for LINE-1 and RASSF1A DNA methylation levels (r² = -0.47, p = 0.002). LINE-1 hypomethylation correlated with concurrent RASSF1/CCND2 hypermethylation (Fisher’s exact test, p = 0.02). Both LINE-1 hypomethylation and RASSF1A/CCND2 hypermethylation were not found in benign hepatocellular tumors (HCA and FNH). Our results show that LINE-1 hypomethylation and RASSF1A/CCND2 hypermethylation are epigenetic aberrations specific for the process of malignant liver transformation. In addition, LINE-1 hypomethylation might serve as a future predictive biomarker to identify HCC patients with unfavorable overall survival.

Transposable elements are regulated by context-specific patterns of chromatin marks in mouse embryonic stem cells

The majority of mammalian genomes are devoted to transposable elements (TEs). Whilst TEs are increasingly recognized for their important biological functions, they are a potential danger to genomic stability and are carefully regulated by the epigenetic system. However, the full complexity of this regulatory system is not understood. Here, using mouse embryonic stem cells, we show that TEs are suppressed by heterochromatic marks like H3K9me3, and are also labelled by all major types of chromatin modification in complex patterns, including bivalent activatory and repressive marks. We identified 29 epigenetic modifiers that significantly deregulated at least one type of TE. The loss of Setdb1, Ncor2, Rnf2, Kat5, Prmt5, Uhrf1, and Rrp8 caused widespread changes in TE expression and chromatin accessibility. These effects were context-specific, with different chromatin modifiers regulating the expression and chromatin accessibility of specific subsets of TEs. Our work reveals the complex patterns of epigenetic regulation of TEs.

Transposable elements (TEs) fulfill essential but poorly understood roles in genome organization and gene expression control. Here the authors show that the regulation of TEs occurs through overlapping epigenetic mechanisms that control the expression and chromatin signatures at TEs.

The role of transposable elements in the evolution of aluminium resistance in plants

Aluminium (Al) toxicity can severely reduce root growth and consequently affect plant development and yield. A mechanism by which many species resist the toxic effects of Al relies on the efflux of organic anions (OAs) from the root apices via OA transporters. Several of the genes encoding these OA transporters contain transposable elements (TEs) in the coding sequences or in flanking regions. Some of the TE-induced mutations impact Al resistance by modifying the level and/or location of gene expression so that OA efflux from the roots is increased. The importance of genomic modifications for improving the adaptation of plants to acid soils has been raised previously, but the growing number of examples linking TEs with these changes requires highlighting. Here, we review the role of TEs in creating genetic modifications that enhance the adaptation of plants to acid soils by increasing the release of OAs from the root apices. We argue that TEs have been an important source of beneficial mutations that have co-opted OA transporter proteins with other functions to perform this role. These changes have occurred relatively recently in the evolution of many species and likely facilitated their expansion into regions with acidic soils.

Engineering modular diterpene biosynthetic pathways in Physcomitrella patens

Modular assembly and heterologous expression in the moss Physcomitrella patens of pairs of diterpene synthases results in accumulation of modern land plant diterpenoids. Physcomitrella patens is a representative of the ancient bryophyte plant lineage with a genome size of 511 Mb, dominant haploid life cycle and limited chemical and metabolic complexity. For these plants, exceptional capacity for genome editing through homologous recombination is met with recently demonstrated in vivo assembly of multiple heterologous DNA fragments. These traits earlier made P. patens an attractive choice as a biotechnological chassis for photosynthesis-driven production of recombinant peptides. The lack of diterpene gibberellic acid phytohormones in P. patens combined with the recent targeted disruption of the single bifunctional diterpene synthase yielded lines devoid of endogenous diterpenoid metabolites and well-suited for engineering of terpenoid production. Here, we mimicked the modular nature of diterpene biosynthetic pathways found in modern land plants by developing a flexible pipeline to install three combinations of class II and class I diterpene synthases in P. patens to access industrially relevant diterpene biomaterials. In addition to a well-established neutral locus for targeted integration, we also explored loci created by a class of Long Terminal Repeat Retrotransposon present at moderate number in the genome of P. patens. Assembly of the pathways and production of the enzymes from the neutral locus led to accumulation of diterpenes matching the reported activities in the angiosperm sources. In contrast, insights gained with the retrotransposon loci indicate their suitability for targeting, but reveal potentially inherent complications which may require adaptation of the experimental design.

Discovery and annotation of a novel transposable element family in Gossypium

BACKGROUND

Fluorescence in situ hybridization (FISH) is an efficient cytogenetic technology to study chromosome structure. Transposable element (TE) is an important component in eukaryotic genomes and can provide insights in the structure and evolution of eukaryotic genomes.

RESULTS

A FISH probe derived from bacterial artificial chromosome (BAC) clone 299N22 generated striking signals on all 26 chromosomes of the cotton diploid A genome (AA, 2x=26) but very few on the diploid D genome (DD, 2x=26). All 26 chromosomes of the A sub genome (At) of tetraploid cotton (AADD, 2n=4x=52) also gave positive signals with this FISH probe, whereas very few signals were observed on the D sub genome (Dt). Sequencing and annotation of BAC clone 299N22, revealed a novel Ty3/gypsy transposon family, which was named as 'CICR'. This family is a significant contributor to size expansion in the A (sub) genome but not in the D (sub) genome. Further FISH analysis with the LTR of CICR as a probe revealed that CICR is lineage-specific, since massive repeats were found in A and B genomic groups, but not in C-G genomic groups within the Gossypium genus. Molecular evolutionary analysis of CICR suggested that tetraploid cottons evolved after silence of the transposon family 1-1.5 million years ago (Mya). Furthermore, A genomes are more homologous with B genomes, and the C, E, F, and G genomes likely diverged from a common ancestor prior to 3.5-4 Mya, the time when CICR appeared. The genomic variation caused by the insertion of CICR in the A (sub) genome may have played an important role in the speciation of organisms with A genomes.

CONCLUSIONS

The CICR family is highly repetitive in A and B genomes of Gossypium, but not amplified in the C-G genomes. The differential amount of CICR family in At and Dt will aid in partitioning sub genome sequences for chromosome assemblies during tetraploid genome sequencing and will act as a method for assessing the accuracy of tetraploid genomes by looking at the proportion of CICR elements in resulting pseudochromosome sequences. The timeline of the expansion of CICR family provides a new reference for cotton evolutionary analysis, while the impact on gene function caused by the insertion of CICR elements will be a target for further analysis of investigating phenotypic differences between A genome and D genome species.

The Evolution of Small-RNA-Mediated Silencing of an Invading Transposable Element

Transposable elements (TEs) are genomic parasites that impose fitness costs on their hosts by producing deleterious mutations and disrupting gametogenesis. Host genomes avoid these costs by regulating TE activity, particularly in germline cells where new insertions are heritable and TEs are exceptionally active. However, the capacity of different TE-associated fitness costs to select for repression in the host, and the role of selection in the evolution of TE regulation more generally remain controversial. In this study, we use forward, individual-based simulations to examine the evolution of small-RNA-mediated TE regulation, a conserved mechanism for TE repression that is employed by both prokaryotes and eukaryotes. To design and parameterize a biologically realistic model, we drew on an extensive survey of empirical studies of the transposition and regulation of P-element DNA transposons in Drosophila melanogaster. We observed that even under conservative assumptions, where small-RNA-mediated regulation reduces transposition only, repression evolves rapidly and adaptively after the genome is invaded by a new TE in simulated populations. We further show that the spread of repressor alleles through simulated populations is greatly enhanced by two additional TE-imposed fitness costs: dysgenic sterility and ectopic recombination. Finally, we demonstrate that the adaptive mutation rate to repression is a critical parameter that influences both the evolutionary trajectory of host repression and the associated proliferation of TEs after invasion in simulated populations. Our findings suggest that adaptive evolution of TE regulation may be stronger and more prevalent than previously appreciated, and provide a framework for interpreting empirical data.

MITE Tracker: an accurate approach to identify miniature inverted-repeat transposable elements in large genomes

BACKGROUND

Miniature inverted-repeat transposable elements (MITEs) are short, non-autonomous class II transposable elements present in a high number of conserved copies in eukaryote genomes. An accurate identification of these elements can help to shed light on the mechanisms controlling genome evolution and gene regulation. The structure and distribution of these elements are well-defined and therefore computational approaches can be used to identify MITEs sequences.

RESULTS

Here we describe MITE Tracker, a novel, open source software program that finds and classifies MITEs using an efficient alignment strategy to retrieve nearby inverted-repeat sequences from large genomes. This program groups them into high sequence homology families using a fast clustering algorithm and finally filters only those elements that were likely transposed from different genomic locations because of their low scoring flanking sequence alignment.

CONCLUSIONS

Many programs have been proposed to find MITEs hidden in genomes. However, none of them are able to process large-scale genomes such as that of bread wheat. Furthermore, in many cases the existing methods perform high false-positive rates (or miss rates). The rice genome was used as reference to compare MITE Tracker against known tools. Our method turned out to be the most reliable in our tests. Indeed, it revealed more known elements, presented the lowest false-positive number and was the only program able to run with the bread wheat genome as input. In wheat, MITE Tracker discovered 6013 MITE families and allowed the first structural exploration of MITEs in the complete bread wheat genome.

Transposable Elements in the Organization and Diversification of the Genome of Aegilops speltoides Tausch (Poaceae, Triticeae)

Repetitive DNA-specifically, transposable elements (TEs)-is a prevailing genomic fraction in cereals that underlies extensive genome reshuffling and intraspecific diversification in the wild. Although large amounts of data have been accumulated, the effect of TEs on the genome architecture and functioning is not fully understood. Here, plant genome organization was addressed by means of cloning and sequencing TE fragments of different types, which compose the largest portion of the Aegilops speltoides genome. Individual genotypes were analyzed cytogenetically using the cloned TE fragments as the DNA probes for fluorescence in situ hybridization (FISH). The obtained TE sequences of the Ty1-copia, Ty3-gypsy, LINE, and CACTA superfamilies showed the relatedness of the Ae. speltoides genome to the Triticeae tribe and similarities to evolutionarily distant species. A significant number of clones consisted of intercalated fragments of TEs of various types, in which Fatima (Ty3-gypsy) sequences predominated. At the chromosomal level, different TE clones demonstrated sequence-specific patterning, emphasizing the effect of the TE fraction on the Ae. speltoides genome architecture and intraspecific diversification. Altogether, the obtained data highlight the current species-specific organization and patterning of the mobile element fraction and point to ancient evolutionary events in the genome of Ae. speltoides.

Birth and Death of LTR-Retrotransposons in Aegilops tauschii

Long terminal repeat-retrotransposons (LTR-RTs) are a major component of all flowering plant genomes. To analyze the time dynamics of LTR-RTs, we modeled the insertion rates of the 35 most abundant LTR-RT families in the genome of Aegilops tauschii, one of the progenitors of wheat. Our model of insertion rate (birth) takes into account random variation in LTR divergence and the deletion rate (death) of LTR-RTs. Modeling the death rate is crucial because ignoring it would underestimate insertion rates in the distant past. We rejected the hypothesis of constancy of insertion rates for all 35 families and showed by simulations that our hypothesis test controlled the false-positive rate. LTR-RT insertions peaked from 0.064 to 2.39 MYA across the 35 families. Among other effects, the average age of elements within a family was negatively associated with recombination rate along a chromosome, with proximity to the closest gene, and weakly associated with the proximity to its 5' end. Elements within a family that were near genes colinear with genes in the genome of tetraploid emmer wheat tended to be younger than those near noncolinear genes. We discuss these associations in the context of genome evolution and stability of genome sizes in the tribe Triticeae. We demonstrate the general utility of our models by analyzing the two most abundant LTR-RT families in Arabidopsis lyrata, and show that these families differed in their insertion dynamics. Our estimation methods are available in the R package TE on CRAN.

Intraspecific and intraorganismal copy number dynamics of retrotransposons and tandem repeat in Aegilops speltoides Tausch (Poaceae, Triticeae)

Transposable elements (TE) and tandem repeats (TR) compose the largest fraction of the plant genome. The abundance and repatterning of repetitive DNA underlie intrapopulation polymorphisms and intraspecific diversification; however, the dynamics of repetitive elements in ontogenesis is not fully understood. Here, we addressed the genotype-specific and tissue-specific abundances and dynamics of the Ty1-copia, Ty3-gypsy, and LINE retrotransposons and species-specific Spelt1 tandem repeat in wild diploid goatgrass, Aegilops speltoides Tausch. Copy numbers of TEs and TR were estimated by real-time quantitative PCR in vegetative and generative tissues in original plants from contrasting allopatric populations and artificial intraspecific hybrids. The results showed that between leaves and somatic spike tissues as well as in progressive microsporogenesis of individual genotypes, the copy numbers of three TEs correlatively oscillated between 2- to 4-fold and the TR copy numbers fluctuated by 18- to 440-fold. Inter-individual and intraorganismal TEs and TR copy number dynamics demonstrate large-scale parallelism with extensive chromosomal rearrangements that were detected using fluorescent in situ hybridization in parental and hybrid genotypes. The data obtained indicate that tissue-specific differences in the abundance and pattern of repetitive sequences emerge during cell proliferation and differentiation in ontogenesis and reflect the reorganization of individual genomes in changing environments, especially in small peripheral population(s) under the influence of rapid climatic changes.

Epigenetic Variance, Performing Cooperative Structure with Genetics, Is Associated with Leaf Shape Traits in Widely Distributed Populations of Ornamental Tree Prunus mume

Increasing evidence shows that epigenetics plays an important role in phenotypic variance. However, little is known about epigenetic variation in the important ornamental tree Prunus mume. We used amplified fragment length polymorphism (AFLP) and methylation-sensitive amplified polymorphism (MSAP) techniques, and association analysis and sequencing to investigate epigenetic variation and its relationships with genetic variance, environment factors, and traits. By performing leaf sampling, the relative total methylation level (29.80%) was detected in 96 accessions of P. mume. And the relative hemi-methylation level (15.77%) was higher than the relative full methylation level (14.03%). The epigenetic diversity (I^∗ = 0.575, h^∗ = 0.393) was higher than the genetic diversity (I = 0.484, h = 0.319). The cultivated population displayed greater epigenetic diversity than the wild populations in both southwest and southeast China. We found that epigenetic variance and genetic variance, and environmental factors performed cooperative structures, respectively. In particular, leaf length, width and area were positively correlated with relative full methylation level and total methylation level, indicating that the DNA methylation level played a role in trait variation. In total, 203 AFLP and 423 MSAP associated markers were detected and 68 of them were sequenced. Homologous analysis and functional prediction suggested that the candidate marker-linked genes were essential for leaf morphology development and metabolism, implying that these markers play critical roles in the establishment of leaf length, width, area, and ratio of length to width.

Massive Loss of DNA Methylation in Nitrogen-, but Not in Phosphorus-Deficient Zea mays Roots Is Poorly Correlated With Gene Expression Differences

DNA methylation in plants plays a role in transposon silencing, genome stability and gene expression regulation. Environmental factors alter the methylation pattern of DNA and recently nutrient stresses, such as phosphate starvation, were shown to alter DNA methylation. Furthermore, DNA methylation had been frequently addressed in plants with notably small genomes that are poor in transposons. Here, we compare part of the DNA methylome of nitrogen- and phosphorus-deficient maize roots by reduced representation sequencing and analyze their relationship with gene expression under prolonged stresses. Tremendous DNA methylation loss was encountered in maize under nitrogen-deficiency, but much less with phosphorus-deficiency. This occurred only in the symmetrical cytosine context, predominantly in CG context, but also in the CHG context. In contrast to other plants, differential methylation in the more flexible CHH context was essentially absent. In both deficiency conditions a similar number of differentially expressed genes were found and differentially methylated regions (DMRs) were predominantly identified in transposable elements (TEs). A minor fraction of such DMRs was associated with altered gene expression of nearby genes. Interestingly, although these TEs were mostly hypomethylated, they were associated with both up- or down regulated gene expression. Our results suggest a different methylome regulation in maize compared to rice and Arabidopsis upon nutrient deficiencies and point to highly nutrient- and species-specific dynamics of genomic DNA methylation. Description of Significance: DNA methylation suppresses transposons in plant genomes, but was also associated with other genome protective functions and gene expression regulation. Recently it was shown that DNA methylation dynamically responds to several abiotic and biotic environmental factors, but to a large instance, DNA methylation is also heritable. DNA methylation changes have also been reported under phosphorus starvation in rice and Arabidopsis, but its relation with other nutrients and its importance for individual gene expression remains unclear. Here, DNA methylation changes upon the deficiency of two major essential nutrients, nitrogen and phosphorus, were studied in parallel with gene expression responses in maize roots. We show context, nutrient- and species-specific patterns in the methylome, as well as its relation with the nutrient-deficiency transcriptome. While cases of differentially methylated regions in the vicinity of differentially expressed genes were apparent, both positive and negative roles on the gene expression were identified, irrespective of the context.

tRNA-derived small RNAs target transposable element transcripts

tRNA-derived RNA fragments (tRFs) are 18–26 nucleotide small RNAs that are not random degradation products, but are rather specifically cleaved from mature tRNA transcripts. Abundant in stressed or viral-infected cells, the function and potential targets of tRFs are not known. We identified that in the unstressed wild-type male gamete containing pollen of flowering plants, and analogous reproductive structure in non-flowering plant species, tRFs accumulate to high levels. In the reference plant Arabidopsis thaliana, tRFs are processed by Dicer-like 1 and incorporated into Argonaute1 (AGO1), akin to a microRNA. We utilized the fact that many plant small RNAs direct cleavage of their target transcripts to demonstrate that the tRF–AGO1 complex acts to specifically target and cleave endogenous transposable element (TE) mRNAs produced from transcriptionally active TEs. The data presented here demonstrate that tRFs are bona-fide regulatory microRNA-like small RNAs involved in the regulation of genome stability through the targeting of TE transcripts.

Reexamining the P-Element Invasion of Drosophila melanogaster Through the Lens of piRNA Silencing

Transposable elements (TEs) are both important drivers of genome evolution and genetic parasites with potentially dramatic consequences for host fitness. The recent explosion of research on regulatory RNAs reveals that small RNA-mediated silencing is a conserved genetic mechanism through which hosts repress TE activity. The invasion of the Drosophila melanogaster genome by P elements, which happened on a historical timescale, represents an incomparable opportunity to understand how small RNA-mediated silencing of TEs evolves. Repression of P-element transposition emerged almost concurrently with its invasion. Recent studies suggest that this repression is implemented in part, and perhaps predominantly, by the Piwi-interacting RNA (piRNA) pathway, a small RNA-mediated silencing pathway that regulates TE activity in many metazoan germlines. In this review, I consider the P-element invasion from both a molecular and evolutionary genetic perspective, reconciling classic studies of P-element regulation with the new mechanistic framework provided by the piRNA pathway. I further explore the utility of the P-element invasion as an exemplar of the evolution of piRNA-mediated silencing. In light of the highly-conserved role for piRNAs in regulating TEs, discoveries from this system have taxonomically broad implications for the evolution of repression.

DNA methylation profiles of diverse Brachypodium distachyon align with underlying genetic diversity

DNA methylation, a common modification of genomic DNA, is known to influence the expression of transposable elements as well as some genes. Although commonly viewed as an epigenetic mark, evidence has shown that underlying genetic variation, such as transposable element polymorphisms, often associate with differential DNA methylation states. To investigate the role of DNA methylation variation, transposable element polymorphism, and genomic diversity, whole-genome bisulfite sequencing was performed on genetically diverse lines of the model cereal Brachypodium distachyon. Although DNA methylation profiles are broadly similar, thousands of differentially methylated regions are observed between lines. An analysis of novel transposable element indel variation highlighted hundreds of new polymorphisms not seen in the reference sequence. DNA methylation and transposable element variation is correlated with the genome-wide amount of genetic variation present between samples. However, there was minimal evidence that novel transposon insertions or deletions are associated with nearby differential methylation. This study highlights unique relationships between genetic variation and DNA methylation variation within Brachypodium and provides a valuable map of DNA methylation across diverse resequenced accessions of this model cereal species.

Chromosomal distribution and evolution of abundant retrotransposons in plants: gypsy elements in diploid and polyploid Brachiaria forage grasses

Like other eukaryotes, the nuclear genome of plants consists of DNA with a small proportion of low-copy DNA (genes and regulatory sequences) and very abundant DNA sequence motifs that are repeated thousands up to millions of times in the genomes including transposable elements (TEs) and satellite DNA. Retrotransposons, one class of TEs, are sequences that amplify via an RNA intermediate and reinsert into the genome, are often the major fraction of a genome. Here, we put research on retrotransposons into the larger context of plant repetitive DNA and genome behaviour, showing features of genome evolution in a grass genus, Brachiaria, in relation to other plant species. We show the contrasting amplification of different retroelement fractions across the genome with characteristics for various families and domains. The genus Brachiaria includes both diploid and polyploid species, with similar chromosome types and chromosome basic numbers x = 6, 7, 8 and 9. The polyploids reproduce asexually and are apomictic, but there are also sexual species. Cytogenetic studies and flow cytometry indicate a large variation in DNA content (C-value), chromosome sizes and genome organization. In order to evaluate the role of transposable elements in the genome and karyotype organization of species of Brachiaria, we searched for sequences similar to conserved regions of TEs in RNAseq reads library produced in Brachiaria decumbens. Of the 9649 TE-like contigs, 4454 corresponded to LTR-retrotransposons, and of these, 79.5 % were similar to members of the gypsy superfamily. Sequences of conserved protein domains of gypsy were used to design primers for producing the probes. The probes were used in FISH against chromosomes of accesses of B. decumbens, Brachiaria brizantha, Brachiaria ruziziensis and Brachiaria humidicola. Probes showed hybridization signals predominantly in proximal regions, especially those for retrotransposons of the clades CRM and Athila, while elements of Del and Tat exhibited dispersed signals, in addition to those proximal signals. These results show that the proximal region of Brachiaria chromosomes is a hotspot for retrotransposon insertion, particularly for the gypsy family. The combination of high-throughput sequencing and a chromosome-centric cytogenetic approach allows the abundance, organization and nature of transposable elements to be characterized in unprecedented detail. By their amplification and dispersal, retrotransposons can affect gene expression; they can lead to rapid diversification of chromosomes between species and, hence, are useful for studies of genome evolution and speciation in the Brachiaria genus. Centromeric regions can be identified and mapped, and retrotransposon markers can also assisting breeders in the developing and exploiting interspecific hybrids.

Repetitive sequences and epigenetic modification: inseparable partners play important roles in the evolution of plant sex chromosomes

The present review discusses the roles of repetitive sequences played in plant sex chromosome evolution, and highlights epigenetic modification as potential mechanism of repetitive sequences involved in sex chromosome evolution. Sex determination in plants is mostly based on sex chromosomes. Classic theory proposes that sex chromosomes evolve from a specific pair of autosomes with emergence of a sex-determining gene(s). Subsequently, the newly formed sex chromosomes stop recombination in a small region around the sex-determining locus, and over time, the non-recombining region expands to almost all parts of the sex chromosomes. Accumulation of repetitive sequences, mostly transposable elements and tandem repeats, is a conspicuous feature of the non-recombining region of the Y chromosome, even in primitive one. Repetitive sequences may play multiple roles in sex chromosome evolution, such as triggering heterochromatization and causing recombination suppression, leading to structural and morphological differentiation of sex chromosomes, and promoting Y chromosome degeneration and X chromosome dosage compensation. In this article, we review the current status of this field, and based on preliminary evidence, we posit that repetitive sequences are involved in sex chromosome evolution probably via epigenetic modification, such as DNA and histone methylation, with small interfering RNAs as the mediator.

Silencing of active transposable elements in plants

In plant genomes the vast majority of transposable elements (TEs) are found in a transcriptionally silenced state that is epigenetically propagated from generation to generation. Although the mechanism of this maintenance of silencing has been well studied, it is now clear that the pathways responsible for maintaining TEs in a silenced state differ from the pathways responsible for initially targeting the TE for silencing. Recently, attention in this field has focused on investigating the molecular mechanisms that initiate and establish TE silencing. Here we review the current models of how TEs are triggered for silencing, the data supporting each model, and the key future questions in this fast moving field.

A call for benchmarking transposable element annotation methods

DNA derived from transposable elements (TEs) constitutes large parts of the genomes of complex eukaryotes, with major impacts not only on genomic research but also on how organisms evolve and function. Although a variety of methods and tools have been developed to detect and annotate TEs, there are as yet no standard benchmarks-that is, no standard way to measure or compare their accuracy. This lack of accuracy assessment calls into question conclusions from a wide range of research that depends explicitly or implicitly on TE annotation. In the absence of standard benchmarks, toolmakers are impeded in improving their tools, annotators cannot properly assess which tools might best suit their needs, and downstream researchers cannot judge how accuracy limitations might impact their studies. We therefore propose that the TE research community create and adopt standard TE annotation benchmarks, and we call for other researchers to join the authors in making this long-overdue effort a success.

A Solution to the C-Value Paradox and the Function of Junk DNA: The Genome Balance Hypothesis

The Genome Balance Hypothesis originated from a recent study that provided a mechanism for the phenomenon of genome dominance in ancient polyploids: unique 24nt RNA coverage near genes is greater in genes on the recessive subgenome irrespective of differences in gene expression. 24nt RNAs target transposons. Transposon position effects are now hypothesized to balance the expression of networked genes and provide spring-like tension between pericentromeric heterochromatin and microtubules. The balance (coordination) of gene expression and centromere movement is under selection. Our hypothesis states that this balance can be maintained by many or few transposons about equally well. We explain known balanced distributions of junk DNA within genomes and between subgenomes in allopolyploids (and our hypothesis passes "the onion test" for any so-called solution to the C-value paradox). Importantly, when the allotetraploid maize chromosomes delete redundant genes, their nearby transposons are also lost; this result is explained if transposons near genes function. The Genome Balance Hypothesis is hypothetical because the position effect mechanisms implicated are not proved to apply to all junk DNA, and the continuous nature of the centromeric and gene position effects have not yet been studied as a single phenomenon.

The Ty1-copia LTR retroelement family PARTC is highly conserved in conifers over 200 MY of evolution

Long Terminal Repeat retroelements (LTR-RTs) are a major component of many plant genomes. Although well studied and described in angiosperms, their features and dynamics are poorly understood in gymnosperms. Representative complete copies of a Ty1-copia element isolate in Picea abies and named PARTC were identified in six other conifer species (Picea glauca, Pinus sylvestris, Pinus taeda, Abies sibirica, Taxus baccata and Juniperus communis) covering more than 200 million years of evolution. Here we characterized the structure of this element, assessed its abundance across conifers, studied the modes and timing of its amplification, and evaluated the degree of conservation of its extant copies at nucleotide level over distant species. We demonstrated that the element is ancient, abundant, widespread and its paralogous copies are present in the genera Picea, Pinus and Abies as an LTR-RT family. The amplification leading to the extant copies of PARTC occurred over long evolutionary times spanning 10s of MY and mostly took place after the speciation of the conifers analyzed. The level of conservation of PARTC is striking and may be explained by low substitution rates and limited removal mechanisms for LTR-RTs. These PARTC features and dynamics are representative of a more general scenario for LTR-RTs in gymnosperms quite different from that characterizing the vast majority of LTR-RT elements in angiosperms.