Transposable elements contribute to activation of maize genes in response to abiotic stress

Retrotransposon-based genetic variation of Poa annua populations from contrasting climate conditions

Background

Poa annua L. is an example of a plant characterized by abundant, worldwide distribution from polar to equatorial regions. Due to its high plasticity and extraordinary expansiveness, P. annua is considered an invasive species capable of occupying and surviving in a wide range of habitats including pioneer zones, areas intensively transformed by human activities, remote subarctic meadows and even the Antarctic Peninsula region.

Methods

In the present study, we evaluated the utility of inter-primer binding site (iPBS) markers for assessing the genetic variation of P. annua populations representing contrasting environments from the worldwide range of this species. The electrophoretic patterns of polymerase chain reaction products obtained for each individual were used to estimate the genetic diversity and differentiation between populations.

Results

iPBS genotyping revealed a pattern of genetic variation differentiating the six studied P. annua populations characterized by their different climatic conditions. According to the analysis of molecular variance, the greatest genetic variation was recorded among populations, whereas 41.75% was observed between individuals within populations. The results of principal coordinates analysis (PCoA) and model-based clustering analysis showed a clear subdivision of analyzed populations. According to PCoA, populations from Siberia and the Kola Peninsula were the most different from each other and showed the lowest genetic variability. The application of STRUCTURE software confirmed the unique character of the population from the Kola Peninsula.

Discussion

The lowest variability of the Siberia population suggested that it was subjected to genetic drift. However, although demographic expansion was indicated by negative values of Fu’s F_S statistic and analysis of mismatch distribution, it was not followed by significant traces of a bottleneck or a founder effect. For the Antarctic population, the observed level of genetic variation was surprisingly high, despite the observed significant traces of bottleneck/founder effect following demographic expansion, and was similar to that observed in populations from Poland and the Balkans. For the Antarctic population, the multiple introduction events from different sources are considered to be responsible for such an observation. Moreover, the results of STRUCTURE and PCoA showed that the P. annua from Antarctica has the highest genetic similarity to populations from Europe.

Conclusions

The observed polymorphism should be considered as a consequence of the joint influence of external abiotic stress and the selection process. Environmental changes, due to their ability to induce transposon activation, lead to the acceleration of evolutionary processes through the production of genetic variability.

Structural variants in 3000 rice genomes

Investigation of large structural variants (SVs) is a challenging yet important task in understanding trait differences in highly repetitive genomes. Combining different bioinformatic approaches for SV detection, we analyzed whole-genome sequencing data from 3000 rice genomes and identified 63 million individual SV calls that grouped into 1.5 million allelic variants. We found enrichment of long SVs in promoters and an excess of shorter variants in 5′ UTRs. Across the rice genomes, we identified regions of high SV frequency enriched in stress response genes. We demonstrated how SVs may help in finding causative variants in genome-wide association analysis. These new insights into rice genome biology are valuable for understanding the effects SVs have on gene function, with the prospect of identifying novel agronomically important alleles that can be utilized to improve cultivated rice.

Altered Promoter and G-Box Binding Factor for 1-Deoxy-d-Xylulose-5-Phosphate Synthase Gene Grown from Poa pratensis Seeds after Spaceflight

In plant cells, the nucleus DNA is considered the primary site of injury by the space environment, which could generate genetic alteration. As the part of genomic mutation, genetic variation in the promoter region could regulate gene expression. In the study, it is observed that there is a deletion in the upstream regulatory region of the 1-deoxy-d-xylulose-5-phosphate synthase 1 gene (PpDXS1) of Poa pratensis dwarf mutant and the PpDXS1 transcript abundance is lower in the dwarf mutant. It is indicated that the deletion in the promoter region between wild type and dwarf mutant could be responsible for the regulation of PpDXS1 gene expression. The PpDXS1 promoter of dwarf mutant shows a lower activity as determined by dual luciferase assay in Poa pratensis protoplast, as well as the GUS activity is lower in transgenic Poa pratensis plant. To further investigate the effect of the deletion in the promoter region on PpDXS1 transcript accumulation, the transient assay and yeast one-hybrid experiment demonstrate that the deletion comprises a motif which is a target of G-box binding factor (GBF1), and the motif correlates with an increase in transactivation by GBF1 protein. Taken together, these results indicate that the deletion in the promoter of PpDXS1 isolated from dwarf mutant is sufficient to account for the decrease in PpDXS1 transcript level and GBF1 can regulate the PpDXS1 gene expression, and subsequently affect accumulation of various isoprenoids throughout the plant.

Novel Insights into Plant Genome Evolution and Adaptation as Revealed through Transposable Elements and Non-Coding RNAs in Conifers

Plant genomes are punctuated by repeated bouts of proliferation of transposable elements (TEs), and these mobile bursts are followed by silencing and decay of most of the newly inserted elements. As such, plant genomes reflect TE-related genome expansion and shrinkage. In general, these genome activities involve two mechanisms: small RNA-mediated epigenetic repression and long-term mutational decay and deletion, that is, genome-purging. Furthermore, the spatial relationships between TE insertions and genes are an important force in shaping gene regulatory networks, their downstream metabolic and physiological outputs, and thus their phenotypes. Such cascading regulations finally set up a fitness differential among individuals. This brief review demonstrates factual evidence that unifies most updated conceptual frameworks covering genome size, architecture, epigenetic reprogramming, and gene expression. It aims to give an overview of the impact that TEs may have on genome and adaptive evolution and to provide novel insights into addressing possible causes and consequences of intimidating genome sizes (20⁻30 Gb) in a taxonomic group, conifers.

Centromere Repeats: Hidden Gems of the Genome

Satellite DNAs are now regarded as powerful and active contributors to genomic and chromosomal evolution. Paired with mobile transposable elements, these repetitive sequences provide a dynamic mechanism through which novel karyotypic modifications and chromosomal rearrangements may occur. In this review, we discuss the regulatory activity of satellite DNA and their neighboring transposable elements in a chromosomal context with a particular emphasis on the integral role of both in centromere function. In addition, we discuss the varied mechanisms by which centromeric repeats have endured evolutionary processes, producing a novel, species-specific centromeric landscape despite sharing a ubiquitously conserved function. Finally, we highlight the role these repetitive elements play in the establishment and functionality of de novo centromeres and chromosomal breakpoints that underpin karyotypic variation. By emphasizing these unique activities of satellite DNAs and transposable elements, we hope to disparage the conventional exemplification of repetitive DNA in the historically-associated context of ‘junk’.

Technological advances in maize breeding: past, present and future

Maize has for many decades been both one of the most important crops worldwide and one of the primary genetic model organisms. More recently, maize breeding has been impacted by rapid technological advances in sequencing and genotyping technology, transformation including genome editing, doubled haploid technology, parallelled by progress in data sciences and the development of novel breeding approaches utilizing genomic information. Herein, we report on past, current and future developments relevant for maize breeding with regard to (1) genome analysis, (2) germplasm diversity characterization and utilization, (3) manipulation of genetic diversity by transformation and genome editing, (4) inbred line development and hybrid seed production, (5) understanding and prediction of hybrid performance, (6) breeding methodology and (7) synthesis of opportunities and challenges for future maize breeding.

An updated gene atlas for maize reveals organ-specific and stress-induced genes

Maize (Zea mays L.), a model species for genetic studies, is one of the two most important crop species worldwide. The genome sequence of the reference genotype, B73, representative of the stiff stalk heterotic group was recently updated (AGPv4) using long-read sequencing and optical mapping technology. To facilitate the use of AGPv4 and to enable functional genomic studies and association of genotype with phenotype, we determined expression abundances for replicated mRNA-sequencing datasets from 79 tissues and five abiotic/biotic stress treatments revealing 36 207 expressed genes. Characterization of the B73 transcriptome across six organs revealed 4154 organ-specific and 7704 differentially expressed (DE) genes following stress treatment. Gene co-expression network analyses revealed 12 modules associated with distinct biological processes containing 13 590 genes providing a resource for further association of gene function based on co-expression patterns. Presence-absence variants (PAVs) previously identified using whole genome resequencing data from 61 additional inbred lines were enriched in organ-specific and stress-induced DE genes suggesting that PAVs may function in phenological variation and adaptation to environment. Relative to core genes conserved across the 62 profiled inbreds, PAVs have lower expression abundances which are correlated with their frequency of dispersion across inbreds and on average have significantly fewer co-expression network connections suggesting that a subset of PAVs may be on an evolutionary path to pseudogenization. To facilitate use by the community, we developed the Maize Genomics Resource website (maize.plantbiology.msu.edu) for viewing and data-mining these resources and deployed two new views on the maize electronic Fluorescent Pictograph Browser (bar.utoronto.ca/efp_maize).

Retrotransposon methylation and activity in wild fish (A. anguilla): A matter of size

Understanding how organisms cope with global change is a major question in many fields of biology. Mainly, understanding the molecular mechanisms supporting rapid phenotypic changes of organisms in response to stress and linking stress-induced molecular events to adaptive or adverse outcomes at the individual or population levels remain a major challenge in evolutionary biology, ecology or ecotoxicology. In this view, the present study aimed to test (i) whether environmental factors, especially pollutants, can trigger changes in the activity of retrotransposons (RTs) in wild fish and (ii) if changes in RT DNA methylation or transcription levels can be linked to modifications at the individual level. RTs are genetic elements that have the ability to replicate and integrate elsewhere in the genome. Although RTs are mainly quiescent during normal development, they can be experimentally activated under life-threatening conditions, affecting the fitness of their host. Wild eels were collected in four sampling sites presenting differing levels of contamination. The methylation level and the transcriptional activity of two RTs and two genes involved in development and cell differentiation were analyzed in fish liver in addition to the determination of fish contaminants levels and diverse growth and morphometric indices. An up-regulation of RTs associated to lower methylation levels and lower growth indices were observed in highly contaminated fish. Our results suggest that RT activation in fish experiencing stress conditions could have both detrimental and beneficial implications, affecting fish growth but promoting resistance to environmental stressors such as pollutants.

Computational aspects underlying genome to phenome analysis in plants

Recent advances in genomics technologies have greatly accelerated the progress in both fundamental plant science and applied breeding research. Concurrently, high-throughput plant phenotyping is becoming widely adopted in the plant community, promising to alleviate the phenotypic bottleneck. While these technological breakthroughs are significantly accelerating quantitative trait locus (QTL) and causal gene identification, challenges to enable even more sophisticated analyses remain. In particular, care needs to be taken to standardize, describe and conduct experiments robustly while relying on plant physiology expertise. In this article, we review the state of the art regarding genome assembly and the future potential of pangenomics in plant research. We also describe the necessity of standardizing and describing phenotypic studies using the Minimum Information About a Plant Phenotyping Experiment (MIAPPE) standard to enable the reuse and integration of phenotypic data. In addition, we show how deep phenotypic data might yield novel trait-trait correlations and review how to link phenotypic data to genomic data. Finally, we provide perspectives on the golden future of machine learning and their potential in linking phenotypes to genomic features.

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.

Epigenetic responses to abiotic stresses during reproductive development in cereals

KEY MESSAGE

Overview of current understanding of epigenetic alterations after abiotic stresses during reproductive development in cereals. Abiotic stresses, including heat, drought, cold, flooding, and salinity, negatively impact crop productivity. Various stages during reproductive development are especially sensitive to environmental stresses, which may lead to complete sterility and severe yield losses. Plants exhibit diverse responses to ameliorate stress damage. Changes in DNA methylation, histone modification as well as regulation of small RNA and long noncoding RNA pathways have been shown to represent key modulators in plant stress responses. During reproductive development in cereals, various protein complexes controlling histone and DNA methylation have been identified, revealing conserved and novel mechanisms regulating abiotic stress responses in cereals and other plant species. New findings highlight the role of transposable elements during stress periods. Here, we review our current understanding of epigenetic stress responses during male and female gametophyte formation (germline development), fertilization, early seed devolvement, and seed maturation in cereals. An integrative model of epigenetic responses during reproductive development in cereals is proposed, emphasizing the role of DNA methylation and histone modifications during abiotic stresses.

Ten things you should know about transposable elements

Transposable elements (TEs) are major components of eukaryotic genomes. However, the extent of their impact on genome evolution, function, and disease remain a matter of intense interrogation. The rise of genomics and large-scale functional assays has shed new light on the multi-faceted activities of TEs and implies that they should no longer be marginalized. Here, we introduce the fundamental properties of TEs and their complex interactions with their cellular environment, which are crucial to understanding their impact and manifold consequences for organismal biology. While we draw examples primarily from mammalian systems, the core concepts outlined here are relevant to a broad range of organisms.

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.

Low-temperature tolerance in land plants: Are transcript and membrane responses conserved?

Plants' tolerance of low temperatures is an economically and ecologically important limitation on geographic distributions and growing seasons. Tolerance for low temperatures varies significantly across different plant species, and different mechanisms likely act in different species. In order to survive low-temperature stress, plant membranes must maintain their fluidity in increasingly cold and oxidative cellular environments. The responses of different species to low-temperature stress include changes to the types and desaturation levels of membrane lipids, though the precise lipids affected tend to vary by species. Regulation of membrane dynamics and other low-temperature tolerance factors are controlled by both transcriptional and post-transcriptional mechanisms. Here, we review low-temperature induced changes in both membrane lipid composition and gene transcription across multiple related plant species with differing degrees of low-temperature tolerance. We attempt to define a core set of changes for transcripts and lipids across species and treatment variations. Some responses appear to be consistent across all species for which data are available, while many others appear likely to be species or family-specific. Potential rationales are presented, including variance in testing, reporting and the importance of considering the level of stress perceived by the plant.

Recent Activity in Expanding Populations and Purifying Selection Have Shaped Transposable Element Landscapes across Natural Accessions of the Mediterranean Grass Brachypodium distachyon

Transposable element (TE) activity has emerged as a major cause of variation in genome size and structure among species. To what extent TEs contribute to genetic variation and divergence within species, however, is much less clear, mainly because population genomic data have so far only been available for the classical model organisms. In this study, we use the annual Mediterranean grass Brachypodium distachyon to investigate TE dynamics in natural populations. Using whole-genome sequencing data for 53 natural accessions, we identified more than 5,400 TE polymorphisms across the studied genomes. We found, first, that while population bottlenecks and expansions have shaped genetic diversity in B. distachyon, these events did not lead to lineage-specific activations of TE families, as observed in other species. Instead, the same families have been active across the species range and TE activity is homogeneous across populations, indicating the presence of conserved regulatory mechanisms. Second, almost half of the TE insertion polymorphisms are accession-specific, most likely because of recent activity in expanding populations and the action of purifying selection. And finally, although TE insertion polymorphisms are underrepresented in and around genes, more than 1,000 of them occur in genic regions and could thus contribute to functional divergence. Our study shows that while TEs in B. distachyon are “well-behaved” compared with TEs in other species with larger genomes, they are an abundant source of lineage-specific genetic variation and may play an important role in population divergence and adaptation.

Shared and genetically distinct Zea mays transcriptome responses to ongoing and past low temperature exposure

BACKGROUND

Cold temperatures and their alleviation affect many plant traits including the abundance of protein coding gene transcripts. Transcript level changes that occur in response to cold temperatures and their alleviation are shared or vary across genotypes. In this study we identify individual transcripts and groups of functionally related transcripts that consistently respond to cold and its alleviation. Genes that respond differently to temperature changes across genotypes may have limited functional importance. We investigate if these genes share functions, and if their genotype-specific gene expression levels change in magnitude or rank across temperatures.

RESULTS

We estimate transcript abundances from over 22,000 genes in two unrelated Zea mays inbred lines during and after cold temperature exposure. Genotype and temperature contribute to many genes' abundances. Past cold exposure affects many fewer genes. Genes up-regulated in cold encode many cytokinin glucoside biosynthesis enzymes, transcription factors, signalling molecules, and proteins involved in diverse environmental responses. After cold exposure, protease inhibitors and cuticular wax genes are newly up-regulated, and environmentally responsive genes continue to be up-regulated. Genes down-regulated in response to cold include many photosynthesis, translation, and DNA replication associated genes. After cold exposure, DNA replication and translation genes are still preferentially downregulated. Lignin and suberin biosynthesis are newly down-regulated. DNA replication, reactive oxygen species response, and anthocyanin biosynthesis genes have strong, genotype-specific temperature responses. The ranks of genotypes' transcript abundances often change across temperatures.

CONCLUSIONS

We report a large, core transcriptome response to cold and the alleviation of cold. In cold, many of the core suite of genes are up or downregulated to control plant growth and photosynthesis and limit cellular damage. In recovery, core responses are in part to prepare for future stress. Functionally related genes are consistently and greatly up-regulated in a single genotype in response to cold or its alleviation, suggesting positive selection has driven genotype-specific temperature responses in maize.

Impact of transposable elements on genome structure and evolution in bread wheat

BACKGROUND

Transposable elements (TEs) are major components of large plant genomes and main drivers of genome evolution. The most recent assembly of hexaploid bread wheat recovered the highly repetitive TE space in an almost complete chromosomal context and enabled a detailed view into the dynamics of TEs in the A, B, and D subgenomes.

RESULTS

The overall TE content is very similar between the A, B, and D subgenomes, although we find no evidence for bursts of TE amplification after the polyploidization events. Despite the near-complete turnover of TEs since the subgenome lineages diverged from a common ancestor, 76% of TE families are still present in similar proportions in each subgenome. Moreover, spacing between syntenic genes is also conserved, even though syntenic TEs have been replaced by new insertions over time, suggesting that distances between genes, but not sequences, are under evolutionary constraints. The TE composition of the immediate gene vicinity differs from the core intergenic regions. We find the same TE families to be enriched or depleted near genes in all three subgenomes. Evaluations at the subfamily level of timed long terminal repeat-retrotransposon insertions highlight the independent evolution of the diploid A, B, and D lineages before polyploidization and cases of concerted proliferation in the AB tetraploid.

CONCLUSIONS

Even though the intergenic space is changed by the TE turnover, an unexpected preservation is observed between the A, B, and D subgenomes for features like TE family proportions, gene spacing, and TE enrichment near genes.

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.

Arabidopsis Serrate Coordinates Histone Methyltransferases ATXR5/6 and RNA Processing Factor RDR6 to Regulate Transposon Expression

Serrate (SE) is a key component in RNA metabolism. Little is known about whether and how it can regulate epigenetic silencing. Here, we report histone methyltransferases ATXR5 and ATXR6 (ATXR5/6) as novel partners of SE. ATXR5/6 deposit histone 3 lysine 27 monomethylation (H3K27me1) to promote heterochromatin formation, repress transposable elements (TEs), and control genome stability in Arabidopsis. SE binds to ATXR5/6-regulated TE loci and promotes H3K27me1 accumulation in these regions. Furthermore, SE directly enhances ATXR5 enzymatic activity in vitro. Unexpectedly, se mutation suppresses the TE reactivation and DNA re-replication phenotypes in the atxr5 atxr6 mutant. The suppression of TE expression results from triggering RNA-dependent RNA polymerase 6 (RDR6)-dependent RNA silencing in the se atxr5 atxr6 mutant. We propose that SE facilitates ATXR5/6-mediated deposition of the H3K27me1 mark while inhibiting RDR6-mediated RNA silencing to protect TE transcripts. Hence, SE coordinates epigenetic silencing and RNA processing machineries to fine-tune the TE expression.

Analysis of retrotransposon abundance, diversity and distribution in holocentric Eleocharis (Cyperaceae) genomes

BACKGROUND AND AIMS

Long terminal repeat-retrotransposons (LTR-RTs) comprise a large portion of plant genomes, with massive repeat blocks distributed across the chromosomes. Eleocharis species have holocentric chromosomes, and show a positive correlation between chromosome numbers and the amount of nuclear DNA. To evaluate the role of LTR-RTs in karyotype diversity in members of Eleocharis (subgenus Eleocharis), the occurrence and location of different members of the Copia and Gypsy superfamilies were compared, covering interspecific variations in ploidy levels (considering chromosome numbers), DNA C-values and chromosomal arrangements.

METHODS

The DNA C-value was estimated by flow cytometry. Genomes of Eleocharis elegans and E. geniculata were partially sequenced using Illumina MiSeq assemblies, which were a source for searching for conserved proteins of LTR-RTs. POL domains were used for recognition, comparing families and for probe production, considering different families of Copia and Gypsy superfamilies. Probes were obtained by PCR and used in fluorescence in situ hybridization (FISH) against chromosomes of seven Eleocharis species.

KEY RESULTS

A positive correlation between ploidy levels and the amount of nuclear DNA was observed, but with significant variations between samples with the same ploidy levels, associated with repetitive DNA fractions. LTR-RTs were abundant in E. elegans and E. geniculata genomes, with a predominance of Copia Sirevirus and Gypsy Athila/Tat clades. FISH using LTR-RT probes exhibited scattered and clustered signals, but with differences in the chromosomal locations of Copia and Gypsy. The diversity in LTR-RT locations suggests that there is no typical chromosomal distribution pattern for retrotransposons in holocentric chromosomes, except the CRM family with signals distributed along chromatids.

CONCLUSIONS

These data indicate independent fates for each LTR-RT family, including accumulation between and within chromosomes and genomes. Differential activity and small changes in LTR-RTs suggest a secondary role in nuclear DNA variation, when compared with ploidy changes.

Transposable element insertions shape gene regulation and melanin production in a fungal pathogen of wheat

BACKGROUND

Fungal plant pathogens pose major threats to crop yield and sustainable food production if they are highly adapted to their host and the local environment. Variation in gene expression contributes to phenotypic diversity within fungal species and affects adaptation. However, very few cases of adaptive regulatory changes have been reported in fungi and the underlying mechanisms remain largely unexplored. Fungal pathogen genomes are highly plastic and harbor numerous insertions of transposable elements, which can potentially contribute to gene expression regulation. In this work, we elucidated how transposable elements contribute to variation in melanin accumulation, a quantitative trait in fungi that affects survival under stressful conditions.

RESULTS

We demonstrated that differential transcriptional regulation of the gene encoding the transcription factor Zmr1, which controls expression of the genes in the melanin biosynthetic gene cluster, is responsible for variation in melanin accumulation in the fungal plant pathogen Zymoseptoria tritici. We show that differences in melanin levels between two strains of Z. tritici are due to two levels of transcriptional regulation: (1) variation in the promoter sequence of Zmr1 and (2) an insertion of transposable elements upstream of the Zmr1 promoter. Remarkably, independent insertions of transposable elements upstream of Zmr1 occurred in 9% of Z. tritici strains from around the world and negatively regulated Zmr1 expression, contributing to variation in melanin accumulation.

CONCLUSIONS

Our studies identified two levels of transcriptional control that regulate the synthesis of melanin. We propose that these regulatory mechanisms evolved to balance the fitness costs associated with melanin production against its positive contribution to survival in stressful environments.

Single-base resolution methylome analysis shows epigenetic changes in Arabidopsis seedlings exposed to microgravity spaceflight conditions on board the SJ-10 recoverable satellite

DNA methylation is a very important epigenetic modification that participates in many biological functions. Although many studies of DNA methylation have been reported in various plant species, few studies have assessed the global DNA methylation pattern in plants challenged by exposure to microgravity conditions. In this report, we mapped the Arabidopsis genome methylation pattern changes associated with microgravity conditions on board the Chinese recoverable scientific satellite SJ-10 at single-base resolution. Interestingly, we found epigenetic differences in Arabidopsis seedlings exposed to microgravity in that the Arabidopsis genome exhibits lower methylation levels in the CHG, CHH, and CpG contexts under microgravity conditions. Microgravity stimulation was related to altered methylation of a number of genes, including DNA methylation-associated genes, hormone signaling related genes, cell-wall modification genes and transposable elements (TEs). Relatively unstable DNA methylation of TEs was responsible for the induction of active transposons. These observations suggest that DNA demethylation within TEs may affect the transcription of transposons in response to microgravity conditions. In summary, the results of this investigation are beneficial for understanding the mechanism of plant adaptation to microgravity and improve strategies to allow plants to adapt to space.

Spaceflight alters the pattern of chemical tags that adorn DNA in plant seedlings. Weiming Cai and colleagues from the Shanghai Institute of Plant Physiology and Ecology, China, profiled the genome-wide epigenetic patterns of Arabidopsis thaliana seedlings that spent 60 h in the microgravity of low Earth orbit aboard the Shijian-10 recoverable satellite. They analyzed the distribution of methyl tags across the genome — an epigenetic mark that affects gene expression levels — and found that seedlings exposed to microgravity had lower methylation on average than control plants grown on the ground, although certain genes related to methylation, transcription factors and hormones tended to be more methylated. Epigenetic differences were also observed among genes involved in cell-wall modification and in transposable elements. The findings could help inform the design of plants optimized for growth in space.

Inpactor, Integrated and Parallel Analyzer and Classifier of LTR Retrotransposons and Its Application for Pineapple LTR Retrotransposons Diversity and Dynamics

One particular class of Transposable Elements (TEs), called Long Terminal Repeats (LTRs), retrotransposons, comprises the most abundant mobile elements in plant genomes. Their copy number can vary from several hundreds to up to a few million copies per genome, deeply affecting genome organization and function. The detailed classification of LTR retrotransposons is an essential step to precisely understand their effect at the genome level, but remains challenging in large-sized genomes, requiring the use of optimized bioinformatics tools that can take advantage of supercomputers. Here, we propose a new tool: Inpactor, a parallel and scalable pipeline designed to classify LTR retrotransposons, to identify autonomous and non-autonomous elements, to perform RT-based phylogenetic trees and to analyze their insertion times using High Performance Computing (HPC) techniques. Inpactor was tested on the classification and annotation of LTR retrotransposons in pineapple, a recently-sequenced genome. The pineapple genome assembly comprises 44% of transposable elements, of which 23% were classified as LTR retrotransposons. Exceptionally, 16.4% of the pineapple genome assembly corresponded to only one lineage of the Gypsy superfamily: Del, suggesting that this particular lineage has undergone a significant increase in its copy numbers. As demonstrated for the pineapple genome, Inpactor provides comprehensive data of LTR retrotransposons’ classification and dynamics, allowing a fine understanding of their contribution to genome structure and evolution. Inpactor is available at https://github.com/simonorozcoarias/Inpactor.

Transposable elements: all mobile, all different, some stress responsive, some adaptive?

Transposable elements (TEs) were first identified through the polymorphisms they induced in plants and animals. Genomic studies have later revealed that TEs were highly abundant in eukaryotic genomes. Recently, more precise single individual genomic analyses have unravelled the huge diversity of TE insertions in many plant and animal species. In most cases the stress conditions behind this diversity are not known and neither is the adaptive capacity of these natural TE-induced variants. Here, we review some of the most recent examples of TE-related impacts on gene expression at the locus or the genome level and discuss the rich diversity of the TE repertoire and its potential role in adaptive evolution.

Contribution of transposable elements in the plant's genome

Plants maintain extensive growth flexibility under different environmental conditions, allowing them to continuously and rapidly adapt to alterations in their environment. A large portion of many plant genomes consists of transposable elements (TEs) that create new genetic variations within plant species. Different types of mutations may be created by TEs in plants. Many TEs can avoid the host's defense mechanisms and survive alterations in transposition activity, internal sequence and target site. Thus, plant genomes are expected to utilize a variety of mechanisms to tolerate TEs that are near or within genes. TEs affect the expression of not only nearby genes but also unlinked inserted genes. TEs can create new promoters, leading to novel expression patterns or alternative coding regions to generate alternate transcripts in plant species. TEs can also provide novel cis-acting regulatory elements that act as enhancers or inserts within original enhancers that are required for transcription. Thus, the regulation of plant gene expression is strongly managed by the insertion of TEs into nearby genes. TEs can also lead to chromatin modifications and thereby affect gene expression in plants. TEs are able to generate new genes and modify existing gene structures by duplicating, mobilizing and recombining gene fragments. They can also facilitate cellular functions by sharing their transposase-coding regions. Hence, TE insertions can not only act as simple mutagens but can also alter the elementary functions of the plant genome. Here, we review recent discoveries concerning the contribution of TEs to gene expression in plant genomes and discuss the different mechanisms by which TEs can affect plant gene expression and reduce host defense mechanisms.

Proliferation of Regulatory DNA Elements Derived from Transposable Elements in the Maize Genome

Genomic regions free of nucleosomes, which are hypersensitive to DNase I digestion, are known as DNase I hypersensitive sites (DHSs) and frequently contain cis-regulatory DNA elements. To investigate their prevalence and characteristics in maize (Zea mays), we developed high-resolution genome-wide DHS maps using a modified DNase-seq technique. Maize DHSs exhibit depletion of nucleosomes and low levels of DNA methylation and are enriched with conserved noncoding sequences (CNSs). We developed a protoplast-based transient transformation assay to assess the potential gene expression enhancer and/or promoter functions associated with DHSs, which showed that more than 80% of DHSs overlapping with CNSs showed an enhancer function. Strikingly, nearly 25% of maize DHSs were derived from transposable elements (TEs), including both class I and class II transposons. Interestingly, TE-derived DHSs (teDHSs) homologous to retrotransposons were enriched with sequences related to the intrinsic cis-regulatory elements within the long terminal repeats of retrotransposons. We demonstrate that more than 80% of teDHSs can drive transcription of a reporter gene in protoplast assays. These results reveal the widespread occurrence of TE-derived cis-regulatory sequences and suggest that teDHSs play a major role in transcriptional regulation in maize.

Transposons: a blessing curse

The genomes of most plant species are dominated by transposable elements (TEs). Once considered as 'junk DNA', TEs are now known to have a major role in driving genome evolution. Over the last decade, it has become apparent that some stress conditions and other environmental stimuli can drive bursts of activity of certain TE families and consequently new TE insertions. These can give rise to altered gene expression patterns and phenotypes, with new TE insertions sometimes causing flanking genes to become transcriptionally responsive to the same stress conditions that activated the TE in the first place. Such connections between TE-mediated increases in diversity and an accelerated rate of genome evolution provide powerful mechanisms for plants to adapt more rapidly to new environmental conditions. This review will focus on environmentally induced transposition, the mechanisms by which it alters gene expression, and the consequences for plant genome evolution and breeding.

Exaptation of transposable element coding sequences

Transposable elements (TEs) are mobile genetic elements that were once perceived as merely selfish, but are now recognized as potent agents of adaptation. One way TEs contribute to genome evolution is through TE exaptation, a process whereby TEs, which usually persist by replicating in the genome, transform into novel host genes, which thereafter persist by conferring phenotypic benefits. Exapted TEs are known to contribute diverse and vital functions, and may facilitate punctuated equilibrium, yet we have little understanding about the process of TE exaptation. In order to facilitate our understanding of how TE coding sequences may become exapted, here we incorporate the findings of recent publications into a framework and six-step model.

Functional Divergence between Subgenomes and Gene Pairs after Whole Genome Duplications

Gene loss following whole genome duplication (WGD) is often biased, with one subgenome retaining more ancestral genes and the other sustaining more gene deletions. While bias toward the greater expression of gene copies on one subgenome can explain bias in gene loss, this raises the question to what drives differences in gene expression levels between subgenomes. Differences in chromatin modifications and epigenetic markers between subgenomes in several model species are now being identified, providing an explanation for bias in gene expression between subgenomes. WGDs can be classified into duplications with higher, biased gene loss and bias in gene expression between subgenomes versus those with lower, unbiased rates of gene loss and an absence of detectable bias between subgenomes; however, the originally proposed link between these two classes and whether WGD results from an allo- or autopolyploid event is inconsistent with recent data from the allopolyploid Capsella bursa-pastoris. The gene balance hypothesis can explain bias in the functional categories of genes retained following WGD, the difference in gene loss rates between unbiased and biased WGDs, and how plant genomes have avoided being overrun with genes encoding dose-sensitive subunits of multiprotein complexes. Comparisons of gene expression patterns between retained transcription factor pairs in maize suggest the high degree of retention for WGD-derived pairs of transcription factors may instead be explained by the older duplication-degeneration-complementation model.

Horizontal acquisition of transposable elements and viral sequences: patterns and consequences

It is becoming clear that most eukaryotic transposable elements (TEs) owe their evolutionary success in part to horizontal transfer events, which enable them to invade new species. Recent large-scale studies are beginning to unravel the mechanisms and ecological factors underlying this mode of transmission. Viruses are increasingly recognized as vectors in the process but also as a direct source of genetic material horizontally acquired by eukaryotic organisms. Because TEs and endogenous viruses are major catalysts of variation and innovation in genomes, we argue that horizontal inheritance has had a more profound impact in eukaryotic evolution than is commonly appreciated. To support this proposal, we compile a list of examples, including some previously unrecognized, whereby new host functions and phenotypes can be directly attributed to horizontally acquired TE or viral sequences. We predict that the number of examples will rapidly grow in the future as the prevalence of horizontal transfer in the life cycle of TEs becomes even more apparent, firmly establishing this form of non-Mendelian inheritance as a consequential facet of eukaryotic evolution.

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.

Transposable element activity, genome regulation and human health

A convergence of novel genome analysis technologies is enabling population genomic studies of human transposable elements (TEs). Population surveys of human genome sequences have uncovered thousands of individual TE insertions that segregate as common genetic variants, i.e. TE polymorphisms. These recent TE insertions provide an important source of naturally occurring human genetic variation. Investigators are beginning to leverage population genomic data sets to execute genome-scale association studies for assessing the phenotypic impact of human TE polymorphisms. For example, the expression quantitative trait loci (eQTL) analytical paradigm has recently been used to uncover hundreds of associations between human TE insertion variants and gene expression levels. These include population-specific gene regulatory effects as well as coordinated changes to gene regulatory networks. In addition, analyses of linkage disequilibrium patterns with previously characterized genome-wide association study (GWAS) trait variants have uncovered TE insertion polymorphisms that are likely causal variants for a variety of common complex diseases. Gene regulatory mechanisms that underlie specific disease phenotypes have been proposed for a number of these trait associated TE polymorphisms. These new population genomic approaches hold great promise for understanding how ongoing TE activity contributes to functionally relevant genetic variation within and between human populations.

Single-base methylome analysis reveals dynamic epigenomic differences associated with water deficit in apple

Cytosine methylation is an essential feature of epigenetic regulation and is involved in various biological processes. Although cytosine methylation has been analysed at the genomic scale for several plant species, there is a general lack of understanding of the dynamics of global and genic DNA methylation in plants growing in environments challenged with biotic and abiotic stresses. In this study, we mapped cytosine methylation at single-base resolution in the genome of commercial apple (Malus x domestica), and analysed changes in methylation patterns associated with water deficit in representative drought-sensitive and drought-tolerant cultivars. We found that the apple genome exhibits ~54%, ~38% and ~8.5% methylation at CG, CHG and CHH sequence contexts, respectively. We additionally documented changes in gene expression associated with water deficit in an attempt to link methylation and gene expression changes. Global methylation and transcription analysis revealed that promoter-unmethylated genes showed higher expression levels than promoter-methylated genes. Gene body methylation appears to be positively correlated with gene expression. Water deficit stress was associated with changes in methylation at a multitude of genes, including those encoding transcription factors (TFs) and transposable elements (TEs). These results present a methylome map of the apple genome and reveal widespread DNA methylation alterations in response to water deficit stress. These data will be helpful for understanding potential linkages between DNA methylation and gene expression in plants growing in natural environments and challenged with abiotic and biotic stresses.

Deep RNAseq indicates protective mechanisms of cold-tolerant indica rice plants during early vegetative stage

Cold-tolerance in rice may be related to increased cellulose deposition in the cell wall, membrane fatty acids unsaturation and differential expression of several newly identified genes. Low temperature exposure during early vegetative stages limits rice plant's growth and development. Most genes previously related to cold tolerance in rice are from the japonica subspecies. To help clarify the mechanisms that regulate cold tolerance in young indica rice plants, comparative transcriptome analysis of 6 h cold-treated (10 °C) leaves from two genotypes, cold-tolerant (CT) and cold-sensitive (CS), was performed. Differentially expressed genes were identified: 831 and 357 sequences more expressed in the tolerant and in the sensitive genotype, respectively. The genes with higher expression in the CT genotype were used in systems biology analyses to identify protein-protein interaction (PPI) networks and nodes (proteins) that are hubs and bottlenecks in the PPI. From the genes more expressed in the tolerant plants, 60% were reported as affected by cold in previous transcriptome experiments and 27% are located within QTLs related to cold tolerance during the vegetative stage. Novel cold-responsive genes were identified. Quantitative RT-PCR confirmed the high-quality of RNAseq libraries. Several genes related to cell wall assembly or reinforcement are cold-induced or constitutively highly expressed in the tolerant genotype. Cold-tolerant plants have increased cellulose deposition under cold. Genes related to lipid metabolism are more expressed in the tolerant genotype, which has higher membrane fatty acids unsaturation, with increasing levels of linoleic acid under cold. The CT genotype seems to have higher photosynthetic efficiency and antioxidant capacity, as well as more effective ethylene, Ca²⁺ and hormone signaling than the CS. These genes could be useful in future biotechnological approaches aiming to increase cold tolerance in rice.

LTR_retriever: A Highly Accurate and Sensitive Program for Identification of Long Terminal Repeat Retrotransposons

Long terminal repeat retrotransposons (LTR-RTs) are prevalent in plant genomes. The identification of LTR-RTs is critical for achieving high-quality gene annotation. Based on the well-conserved structure, multiple programs were developed for the de novo identification of LTR-RTs; however, these programs are associated with low specificity and high false discovery rates. Here, we report LTR_retriever, a multithreading-empowered Perl program that identifies LTR-RTs and generates high-quality LTR libraries from genomic sequences. LTR_retriever demonstrated significant improvements by achieving high levels of sensitivity (91%), specificity (97%), accuracy (96%), and precision (90%) in rice (Oryza sativa). LTR_retriever is also compatible with long sequencing reads. With 40k self-corrected PacBio reads equivalent to 4.5× genome coverage in Arabidopsis (Arabidopsis thaliana), the constructed LTR library showed excellent sensitivity and specificity. In addition to canonical LTR-RTs with 5'-TG…CA-3' termini, LTR_retriever also identifies noncanonical LTR-RTs (non-TGCA), which have been largely ignored in genome-wide studies. We identified seven types of noncanonical LTRs from 42 out of 50 plant genomes. The majority of noncanonical LTRs are Copia elements, with which the LTR is four times shorter than that of other Copia elements, which may be a result of their target specificity. Strikingly, non-TGCA Copia elements are often located in genic regions and preferentially insert nearby or within genes, indicating their impact on the evolution of genes and their potential as mutagenesis tools.

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.

De Novo Identification of sRNA Loci and Non-coding RNAs by High-Throughput Sequencing

Non-coding RNA transcripts, such as long non-coding RNAs, miRNAs, siRNAs, and transposon-originating transcripts, are involved in the regulation of RNA stability, protein translation, and/or the modulation of chromatin states. RNA-Seq can be used to catalog this diversity of novel transcripts and a joint analysis of these transcriptomic data can provide useful insights into epigenetic regulation of dynamic responses such as the stress response, which may not be deciphered from individual analysis of single transcript categories. Here, we present a protocol that allows the identification and analysis of small RNAs and long non-coding RNAs, together with the comparison of these species between different sample types.

Cis-Effects Condition the Induction of a Major Unfolded Protein Response Factor, ZmbZIP60, in Response to Heat Stress in Maize

Adverse environmental conditions such as heat and salt stress create endoplasmic reticulum (ER) stress in maize and set off the unfolded protein response (UPR). A key feature of the UPR is the upregulation of ZmbZIP60 and the splicing of its messenger RNA. We conducted an association analysis of a recombinant inbred line (RIL) derived from a cross of a tropical founder line, CML52 with a standard temperate line, B73. We found a major QTL conditioning heat-induced ZmbZIP60 expression located cis to the gene. Based on the premise that the QTL might be associated with the ZmbZIP60 promoter, we evaluated various maize inbred lines for their ability to upregulate the expression of ZmbZIP60 in response to heat stress. In general, tropical lines with promoter regions similar to CML52 were more robust in upregulating ZmbZIP60 in response to heat stress. This finding was confirmed by comparing the strength of the B73 and CML52 ZmbZIP60 promoters in transient maize protoplast assays. We concluded that the upstream region of ZmbZIP60 is important in conditioning the response to heat stress and was under selection in maize when adapted to different environments. Summary: Heat stress has large negative effects on maize grain yield. Heat stress creates ER stress in maize and sets off the UPR. We searched for factors conditioning heat induction of the UPR in maize seedlings by conducting an association analysis based on the upregulation of unspliced and spliced forms of ZmbZIP60 mRNA (ZmbZIP60u and ZmbZIP60s, respectively). ZmbZIP60u was upregulated more robustly by heat stress in the tropical maize line, CML52, than in B73, and a major QTL derived from the analysis of RILs from a cross of these two lines mapped in the vicinity of ZmbZIP60. We conducted a cis/trans test to determine whether the QTL was acting as a cis regulatory element or in trans, as might be expected for a transcription factor. We found that the QTL was acting in cis, likely involving the ZmbZIP60 promoter. ZmbZIP60 promoters in other temperate and tropical lines similar to CML52 showed enhanced expression of ZmbZIP60u by heat. The contribution of the CML52 promoter to heat induction of ZmbZIP60 was confirmed by analyzing the CML52 and B73 promoters linked to a luciferase reporter and assayed in heat-treated maize protoplasts.

Molecular characterization of 5' UTR of the lycopene epsilon cyclase (lcyE) gene among exotic and indigenous inbreds for its utilization in maize biofortification

Maize grains are the important source of food and energy, but possess very low proA (< 2.5 µg/g) compared to target level of 15 µg/g set by HarvestPlus to alleviate VAD. Favorable allele having variation in 5' untranslated region (UTR) of lycopene epsilon cyclase (lcyE) gene enhances concentration of proA in maize. To identify the sequence variation in 5' UTR of lcyE, a set of diverse 13 inbreds of indigenous and exotic origin was characterized for allelic constitution of lcyE. Inbreds possessed wide variation in proA (1.62-23.12 µg/g) with a mean of 9.64 µg/g. The proA in CIMMYT-HarvestPlus genotypes having favorable allele of lcyE was very high (22.28 µg/g), whereas the Indian inbreds with the same allele possessed very low proA (2.48 µg/g). Eight genotypes viz., HKI161, HKI163, HKI161-PV, HKI163-PV, HKI193-1-PV, HKI193-2-PV, HP704-22 and HP704-23 revealed the presence of favorable allele, while VQL1, DMRIL47, MGU-PV-123/C6, HKI193-1 and HKI193-2 showed the presence of unfavorable allele of lcyE gene. Sequence comparison of favorable allele of Indian (HKI161 and HKI163) and exotic genotypes (HP704-22 and HP704-23) revealed seven SNPs having three transitions (SNP1 and SNP3: G to A, SNP2: C to T) and four transversions (SNP4: C to G, SNP5: T to G, SNP6: G to C and SNP7: G to T). Four SNPs (SNP1: position 446, SNP2: position 458, SNP3: position 459 and SNP4: position 483) discriminated the low- and high- proA lines having favorable allele of lcyE 5'TE. These SNPs hold significance in enrichment of proA in maize for marker development and their use in marker-assisted selection.

Recombination rate plasticity: revealing mechanisms by design

For over a century, scientists have known that meiotic recombination rates can vary considerably among individuals, and that environmental conditions can modify recombination rates relative to the background. A variety of external and intrinsic factors such as temperature, age, sex and starvation can elicit 'plastic' responses in recombination rate. The influence of recombination rate plasticity on genetic diversity of the next generation has interesting and important implications for how populations evolve. Further, many questions remain regarding the mechanisms and molecular processes that contribute to recombination rate plasticity. Here, we review 100 years of experimental work on recombination rate plasticity conducted in Drosophila melanogaster We categorize this work into four major classes of experimental designs, which we describe via classic studies in D. melanogaster Based on these studies, we highlight molecular mechanisms that are supported by experimental results and relate these findings to studies in other systems. We synthesize lessons learned from this model system into experimental guidelines for using recent advances in genotyping technologies, to study recombination rate plasticity in non-model organisms. Specifically, we recommend (1) using fine-scale genome-wide markers, (2) collecting time-course data, (3) including crossover distribution measurements, and (4) using mixed effects models to analyse results. To illustrate this approach, we present an application adhering to these guidelines from empirical work we conducted in Drosophila pseudoobscuraThis article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Patterns and Consequences of Subgenome Differentiation Provide Insights into the Nature of Paleopolyploidy in Plants

Polyploidy is an important feature of plant genomes, but the nature of many polyploidization events remains to be elucidated. Here, we demonstrate that the evolutionary fates of the subgenomes in maize (Zea mays) and soybean (Glycine max) have followed different trajectories. One subgenome has been subject to relaxed selection, lower levels of gene expression, higher rates of transposable element accumulation, more small interfering RNAs and DNA methylation around genes, and higher rates of gene loss in maize, whereas none of these features were observed in soybean. Nevertheless, individual gene pairs exhibit differentiation with respect to these features in both species. In addition, we observed a higher number of chromosomal rearrangements and higher frequency of retention of duplicated genes in soybean than in maize. Furthermore, soybean "singletons" were found to be more frequently tandemly duplicated than "duplicates" in soybean, which may, to some extent, counteract the genome imbalance caused by gene loss. We propose that unlike in maize, in which two subgenomes were distinct prior to the allotetraploidization event and thus experienced global differences in selective constraints, in soybean, the two subgenomes were far less distinct prior to polyploidization, such that individual gene pairs, rather than subgenomes, experienced stochastic differences over longer periods of time, resulting in retention of the majority of duplicates.

Cut-and-Paste Transposons in Fungi with Diverse Lifestyles

Transposable elements (TEs) shape genomes via recombination and transposition, lead to chromosomal rearrangements, create new gene neighborhoods, and alter gene expression. They play key roles in adaptation either to symbiosis in Amanita genus or to pathogenicity in Pyrenophora tritici-repentis. Despite growing evidence of their importance, the abundance and distribution of mobile elements replicating in a "cut-and-paste" fashion is barely described so far. In order to improve our knowledge on this old and ubiquitous class of transposable elements, 1,730 fungal genomes were scanned using both de novo and homology-based approaches. DNA TEs have been identified across the whole data set and display uneven distribution from both DNA TE classification and fungal taxonomy perspectives. DNA TE content correlates with genome size, which confirms that many transposon families proliferate simultaneously. In contrast, it is independent from intron density, average gene distance and GC content. TE count is associated with species' lifestyle and tends to be elevated in plant symbionts and decreased in animal parasites. Lastly, we found that fungi with both RIP and RNAi systems have more total DNA TE sequences but less elements retaining a functional transposase, what reflects stringent control over transposition.

Abiotic Stress Phenotypes Are Associated with Conserved Genes Derived from Transposable Elements

Plant phenomics offers unique opportunities to accelerate our understanding of gene function and plant response to different environments, and may be particularly useful for studying previously uncharacterized genes. One important type of poorly characterized genes is those derived from transposable elements (TEs), which have departed from a mobility-driven lifestyle to attain new adaptive roles for the host (exapted TEs). We used phenomics approaches, coupled with reverse genetics, to analyze T-DNA insertion mutants of both previously reported and novel protein-coding exapted TEs in the model plant Arabidopsis thaliana. We show that mutations in most of these exapted TEs result in phenotypes, particularly when challenged by abiotic stress. We built statistical multi-dimensional phenotypic profiles and compared them to wild-type and known stress responsive mutant lines for each particular stress condition. We found that these exapted TEs may play roles in responses to phosphate limitation, tolerance to high salt concentration, freezing temperatures, and arsenic toxicity. These results not only experimentally validate a large set of putative functional exapted TEs recently discovered through computational analysis, but also uncover additional novel phenotypes for previously well-characterized exapted TEs in A. thaliana.

Analyses of methylomes of upland and lowland switchgrass (Panicum virgatum) ecotypes using MeDIP-seq and BS-seq

Background

Switchgrass is a crop with many desirable traits for bioenergy production. Plant genomes have high DNA methylation levels throughout genes and transposable elements and DNA methylation is known to play a role in silencing transposable elements. Here we analyzed methylomes in two switchgrass genotypes AP13 and VS16. AP13 is derived from a lowland ecotype and VS16, typically considered drought-tolerant, is derived from an upland ecotype, both genotypes are tetraploid (2n = 4× = 36).

Results

Methylated DNA immunoprecipitation-sequencing (MeDIP-seq) and bisulfite-sequencing (BS-seq) were used to profile DNA methylation in genomic features of AP13 and VS16. The methylation patterns in genes and transposable elements were similar to other plants, however, overall CHH methylation levels were comparatively low. Differentially methylated regions (DMRs) were assessed and a total of 1777 CG-DMRs, 573 CHG-DMRs, and 3 CHH-DMRs were detected between the two genotypes. TEs and their flanking regions were higher than that of genic regions. Different types of TEs had different methylation patterns, but the two LTRs (Copia and Gypsy) were similarly methylated, while LINEs and DNA transposons typically had different methylation patterns. MeDIP-seq data was compared to BS-seq data and most of the peaks generated by MeDIP-seq were confirmed to be highly methylated by BS-seq.

Conclusions

DNA methylation in switchgrass genotypes obtained from the two ecotypes were found similar. Collinear gene pairs in two subgenomes (A and B) were not significantly differentially methylated. Both BS-seq and MeDIP-seq methodologies were found effective. Methylation levels were highest at CG and least in CHH. Increased DNA methylation was seen in TEs compared to genic regions. Exploitation of TE methylations can be a viable option in future crop improvement.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4218-0) contains supplementary material, which is available to authorized users.

Pearl Millet Genome: Lessons from a Tough Crop

Pearl millet is an important cereal for food security in the arid regions of Africa and India. The recently published genome of this tough cereal crop has shed new light on its history and adaptation to dry, hot climates and paves the way for much-needed genomic-based breeding efforts.

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.