Plant transposable elements: where genetics meets genomics

Novel Centromeric Loci of the Wine and Beer Yeast Dekkera bruxellensis CEN1 and CEN2

The wine and beer yeast Dekkera bruxellensis thrives in environments that are harsh and limiting, especially in concentrations with low oxygen and high ethanol. Its different strains’ chromosomes greatly vary in number (karyotype). This study isolates two novel centromeric loci (CEN1 and CEN2), which support both the yeast’s autonomous replication and the stable maintenance of plasmids. In the sequenced genome of the D. bruxellensis strain CBS 2499, CEN1 and CEN2 are each present in one copy. They differ from the known “point” CEN elements, and their biological activity is retained within ~900–1300 bp DNA segments. CEN1 and CEN2 have features of both “point” and “regional” centromeres: They contain conserved DNA elements, ARSs, short repeats, one tRNA gene, and transposon-like elements within less than 1 kb. Our discovery of a miniature inverted-repeat transposable element (MITE) next to CEN2 is the first report of such transposons in yeast. The transformants carrying circular plasmids with cloned CEN1 and CEN2 undergo a phenotypic switch: They form fluffy colonies and produce three times more biofilm. The introduction of extra copies of CEN1 and CEN2 promotes both genome rearrangements and ploidy shifts, with these effects mediated by homologous recombination (between circular plasmid and genome centromere copy) or by chromosome breakage when integrated. Also, the proximity of the MITE-like transposon to CEN2 could translocate CEN2 within the genome or cause chromosomal breaks, so promoting genome dynamics. With extra copies of CEN1 and CEN2, the yeast’s enhanced capacities to rearrange its genome and to change its gene expression could increase its abilities for exploiting new and demanding niches.

Long Terminal Repeat Retrotransposon Content in Eight Diploid Sunflower Species Inferred from Next-Generation Sequence Data

The most abundant transposable elements (TEs) in plant genomes are Class I long terminal repeat (LTR) retrotransposons represented by superfamilies gypsy and copia. Amplification of these superfamilies directly impacts genome structure and contributes to differential patterns of genome size evolution among plant lineages. Utilizing short-read Illumina data and sequence information from a panel of Helianthus annuus (sunflower) full-length gypsy and copia elements, we explore the contribution of these sequences to genome size variation among eight diploid Helianthus species and an outgroup taxon, Phoebanthus tenuifolius. We also explore transcriptional dynamics of these elements in both leaf and bud tissue via RT-PCR. We demonstrate that most LTR retrotransposon sublineages (i.e., families) display patterns of similar genomic abundance across species. A small number of LTR retrotransposon sublineages exhibit lineage-specific amplification, particularly in the genomes of species with larger estimated nuclear DNA content. RT-PCR assays reveal that some LTR retrotransposon sublineages are transcriptionally active across all species and tissue types, whereas others display species-specific and tissue-specific expression. The species with the largest estimated genome size, H. agrestis, has experienced amplification of LTR retrotransposon sublineages, some of which have proliferated independently in other lineages in the Helianthus phylogeny.

Characterization of the Poplar Pan-Genome by Genome-Wide Identification of Structural Variation

Many recent studies have emphasized the important role of structural variation (SV) in determining human genetic and phenotypic variation. In plants, studies aimed at elucidating the extent of SV are still in their infancy. Evidence has indicated a high presence and an active role of SV in driving plant genome evolution in different plant species.

With the aim of characterizing the size and the composition of the poplar pan-genome, we performed a genome-wide analysis of structural variation in three intercrossable poplar species: Populus nigra, Populus deltoides, and Populus trichocarpa. We detected a total of 7,889 deletions and 10,586 insertions relative to the P. trichocarpa reference genome, covering respectively 33.2 Mb and 62.9 Mb of genomic sequence, and 3,230 genes affected by copy number variation (CNV). The majority of the detected variants are inter-specific in agreement with a recent origin following separation of species.

Insertions and deletions (INDELs) were preferentially located in low-gene density regions of the poplar genome and were, for the majority, associated with the activity of transposable elements. Genes affected by SV showed lower-than-average expression levels and higher levels of dN/dS, suggesting that they are subject to relaxed selective pressure or correspond to pseudogenes.

Functional annotation of genes affected by INDELs showed over-representation of categories associated with transposable elements activity, while genes affected by genic CNVs showed enrichment in categories related to resistance to stress and pathogens. This study provides a genome-wide catalogue of SV and the first insight on functional and structural properties of the poplar pan-genome.

High-throughput sequencing and graph-based cluster analysis facilitate microsatellite development from a highly complex genome

Despite recent advances in high-throughput sequencing, difficulties are often encountered when developing microsatellites for species with large and complex genomes. This probably reflects the close association in many species of microsatellites with cryptic repetitive elements. We therefore developed a novel approach for isolating polymorphic microsatellites from the club-legged grasshopper (Gomphocerus sibiricus), an emerging quantitative genetic and behavioral model system. Whole genome shotgun Illumina MiSeq sequencing was used to generate over three million 300 bp paired-end reads, of which 67.75% were grouped into 40,548 clusters within RepeatExplorer. Annotations of the top 468 clusters, which represent 60.5% of the reads, revealed homology to satellite DNA and a variety of transposable elements. Evaluating 96 primer pairs in eight wild-caught individuals, we found that primers mined from singleton reads were six times more likely to amplify a single polymorphic microsatellite locus than primers mined from clusters. Our study provides experimental evidence in support of the notion that microsatellites associated with repetitive elements are less likely to successfully amplify. It also reveals how advances in high-throughput sequencing and graph-based repetitive DNA analysis can be leveraged to isolate polymorphic microsatellites from complex genomes.

The evolutionary ecology of molecular replicators

By reasonable criteria, life on the Earth consists mainly of molecular replicators. These include viruses, transposons, transpovirons, coviruses and many more, with continuous new discoveries like Sputnik Virophage. Their study is inherently multidisciplinary, spanning microbiology, genetics, immunology and evolutionary theory, and the current view is that taking a unified approach has great power and promise. We support this with a new, unified, model of their evolutionary ecology, using contemporary evolutionary theory coupling the Price equation with game theory, studying the consequences of the molecular replicators' promiscuous use of each others' gene products for their natural history and evolutionary ecology. Even at this simple expository level, we can make a firm prediction of a new class of replicators exploiting viruses such as lentiviruses like SIVs, a family which includes HIV: these have been explicitly stated in the primary literature to be non-existent. Closely connected to this departure is the view that multicellular organism immunology is more about the management of chronic infections rather than the elimination of acute ones and new understandings emerging are changing our view of the kind of theatre we ourselves provide for the evolutionary play of molecular replicators. This study adds molecular replicators to bacteria in the emerging field of sociomicrobiology.

Retrozymes are a unique family of non-autonomous retrotransposons with hammerhead ribozymes that propagate in plants through circular RNAs

Background

Catalytic RNAs, or ribozymes, are regarded as fossils of a prebiotic RNA world that have remained in the genomes of modern organisms. The simplest ribozymes are the small self-cleaving RNAs, like the hammerhead ribozyme, which have been historically considered biological oddities restricted to some RNA pathogens. Recent data, however, indicate that small self-cleaving ribozymes are widespread in genomes, although their functions are still unknown.

Results

We reveal that hammerhead ribozyme sequences in plant genomes form part of a new family of small non-autonomous retrotransposons with hammerhead ribozymes, referred to as retrozymes. These elements contain two long terminal repeats of approximately 350 bp, each harbouring a hammerhead ribozyme that delimitates a variable region of 600–1000 bp with no coding capacity. Retrozymes are actively transcribed, which gives rise to heterogeneous linear and circular RNAs that accumulate differentially depending on the tissue or developmental stage of the plant. Genomic and transcriptomic retrozyme sequences are highly heterogeneous and share almost no sequence homology among species except the hammerhead ribozyme motif and two small conserved domains typical of Ty3-gypsy long terminal repeat retrotransposons. Moreover, we detected the presence of RNAs of both retrozyme polarities, which suggests events of independent RNA-RNA rolling-circle replication and evolution, similarly to that of infectious circular RNAs like viroids and viral satellite RNAs.

Conclusions

Our work reveals that circular RNAs with hammerhead ribozymes are frequently occurring molecules in plant and, most likely, metazoan transcriptomes, which explains the ubiquity of these genomic ribozymes and suggests a feasible source for the emergence of circular RNA plant pathogens.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-1002-4) contains supplementary material, which is available to authorized users.

Identification of Bari Transposons in 23 Sequenced Drosophila Genomes Reveals Novel Structural Variants, MITEs and Horizontal Transfer

Bari elements are members of the Tc1-mariner superfamily of DNA transposons, originally discovered in Drosophila melanogaster, and subsequently identified in silico in 11 sequenced Drosophila genomes and as experimentally isolated in four non-sequenced Drosophila species. Bari-like elements have been also studied for their mobility both in vivo and in vitro. We analyzed 23 Drosophila genomes and carried out a detailed characterization of the Bari elements identified, including those from the heterochromatic Bari1 cluster in D. melanogaster. We have annotated 401 copies of Bari elements classified either as putatively autonomous or inactive according to the structure of the terminal sequences and the presence of a complete transposase-coding region. Analyses of the integration sites revealed that Bari transposase prefers AT-rich sequences in which the TA target is cleaved and duplicated. Furthermore evaluation of transposon’s co-occurrence near the integration sites of Bari elements showed a non-random distribution of other transposable elements. We also unveil the existence of a putatively autonomous Bari1 variant characterized by two identical long Terminal Inverted Repeats, in D. rhopaloa. In addition, we detected MITEs related to Bari transposons in 9 species. Phylogenetic analyses based on transposase gene and the terminal sequences confirmed that Bari-like elements are distributed into three subfamilies. A few inconsistencies in Bari phylogenetic tree with respect to the Drosophila species tree could be explained by the occurrence of horizontal transfer events as also suggested by the results of dS analyses. This study further clarifies the Bari transposon’s evolutionary dynamics and increases our understanding on the Tc1-mariner elements’ biology.

Structural organization of fatty acid desaturase loci in linseed lines with contrasting linolenic acid contents

Flax (Linum usitatissimum L.), the richest crop source of omega-3 fatty acids (FAs), is a diploid plant with an estimated genome size of ~370 Mb and is well suited for studying genomic organization of agronomically important traits. In this study, 12 bacterial artificial chromosome clones harbouring the six FA desaturase loci sad1, sad2, fad2a, fad2b, fad3a and fad3b from the conventional variety CDC Bethune and the high linolenic acid line M5791 were sequenced, analysed and compared to determine the structural organization of these loci and to gain insights into the genetic mechanisms underlying FA composition in flax. With one gene every 3.2-4.6 kb, the desaturase loci have a higher gene density than the genome's average of one gene per 7.8-8.2 kb. The gene order and orientation across the two genotypes were generally conserved with the exception of the sad1 locus that was predicted to have additional genes in CDC Bethune. High sequence conservation in both genic and intergenic regions of the sad and fad2b loci contrasted with the significant level of variation of the fad2a and fad3 loci, with SNPs being the most frequently observed mutation type. The fad2a locus had 297 SNPs and 36 indels over ~95 kb contrasting with the fad2b locus that had a mere seven SNPs and four indels in ~110 kb. Annotation of the gene-rich loci revealed other genes of known role in lipid or carbohydrate metabolic/catabolic pathways. The organization of the fad2b locus was particularly complex with seven copies of the fad2b gene in both genotypes. The presence of Gypsy, Copia, MITE, Mutator, hAT and other novel repeat elements at the desaturase loci was similar to that of the whole genome. This structural genomic analysis provided some insights into the genomic organization and composition of the main desaturase loci of linseed and of their complex evolution through both tandem and whole genome duplications.

Evolutionary origin of Rosaceae-specific active non-autonomous hAT elements and their contribution to gene regulation and genomic structural variation

Transposable elements account for approximately 30 % of the Prunus genome; however, their evolutionary origin and functionality remain largely unclear. In this study, we identified a hAT transposon family, termed Moshan, in Prunus. The Moshan elements consist of three types, aMoshan, tMoshan, and mMoshan. The aMoshan and tMoshan types contain intact or truncated transposase genes, respectively, while the mMoshan type is miniature inverted-repeat transposable element (MITE). The Moshan transposons are unique to Rosaceae, and the copy numbers of different Moshan types are significantly correlated. Sequence homology analysis reveals that the mMoshan MITEs are direct deletion derivatives of the tMoshan progenitors, and one kind of mMoshan containing a MuDR-derived fragment were amplified predominately in the peach genome. The mMoshan sequences contain cis-regulatory elements that can enhance gene expression up to 100-fold. The mMoshan MITEs can serve as potential sources of micro and long noncoding RNAs. Whole-genome re-sequencing analysis indicates that mMoshan elements are highly active, and an insertion into S-haplotype-specific F-box gene was reported to cause the breakdown of self-incompatibility in sour cherry. Taken together, all these results suggest that the mMoshan elements play important roles in regulating gene expression and driving genomic structural variation in Prunus.

The concerted impact of domestication and transposon insertions on methylation patterns between dogs and grey wolves

The process of domestication can exert intense trait-targeted selection on genes and regulatory regions. Specifically, rapid shifts in the structure and sequence of genomic regulatory elements could provide an explanation for the extensive, and sometimes extreme, variation in phenotypic traits observed in domesticated species. Here, we explored methylation differences from >24 000 cytosines distributed across the genomes of the domesticated dog (Canis familiaris) and the grey wolf (Canis lupus). PCA and model-based cluster analyses identified two primary groups, domestic vs. wild canids. A scan for significantly differentially methylated sites (DMSs) revealed species-specific patterns at 68 sites after correcting for cell heterogeneity, with weak yet significant hypermethylation typical of purebred dogs when compared to wolves (59% and 58%, P < 0.05, respectively). Additionally, methylation patterns at eight genes significantly deviated from neutrality, with similar trends of hypermethylation in purebred dogs. The majority (>66%) of differentially methylated regions contained or were associated with repetitive elements, indicative of a genotype-mediated trend. However, DMSs were also often linked to functionally relevant genes (e.g. neurotransmitters). Finally, we utilized known genealogical relationships among Yellowstone wolves to survey transmission stability of methylation marks, from which we found a substantial fraction that demonstrated high heritability (both H(2) and h(2 ) > 0.99). These analyses provide a unique epigenetic insight into the molecular consequences of recent selection and radiation of our most ancient domesticated companion, the dog. These findings suggest selection has acted on methylation patterns, providing a new genomic perspective on phenotypic diversification in domesticated species.

Principal components analysis--K-means transposon element based foxtail millet core collection selection method

Background

Core collections are important tools in genetic resources research and administration. At present, most core collection selection criteria are based on one of the following item characteristics: passport data, genetic markers, or morphological traits, which may lead to inadequate representations of variability in the complete collection. The development of a comprehensive methodology that includes as much element data as possible has been explored poorly. Using a collection of (Setaria italica sbsp. italica (L.) P. Beauv.) as a model, we developed a method for core collection construction based on genotype data and numerical representations of agromorphological traits, thereby improving the selection process.

Results

Principal component analysis allows the selection of the most informative discriminators among the various elements evaluated, regardless of whether they are genetic or morphological, thereby providing an adequate criterion for further K-mean clustering. Overall, the core collections of S. italica constructed using only genotype data demonstrated overall better validation scores than other core collections that we generated. However, core collection based on both genotype and agromorphological characteristics represented the overall diversity adequately.

Conclusions

The inclusion of both genotype and agromorphological characteristics as a comprehensive dataset in this methodology ensures that agricultural traits are considered in the core collection construction. This approach will be beneficial for genetic resources management and research activities for S. italica as well as other genetic resources.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0343-z) contains supplementary material, which is available to authorized users.

Identification and characterization of abundant repetitive sequences in Eragrostis tef cv. Enatite genome

BACKGROUND

Eragrostis tef is an allotetraploid (2n = 4 × = 40) annual, C4 grass with an estimated nuclear genome size of 730 Mbp. It is widely grown in Ethiopia, where it provides basic nutrition for more than half of the population. Although a draft assembly of the E. tef genome was made available in 2014, characterization of the repetitive portion of the E. tef genome has not been a subject of a detailed analysis. Repetitive sequences constitute most of the DNA in eukaryotic genomes. Transposable elements are usually the most abundant repetitive component in plant genomes. They contribute to genome size variation, cause mutations, can result in chromosomal rearrangements, and influence gene regulation. An extensive and in depth characterization of the repetitive component is essential in understanding the evolution and function of the genome.

RESULTS

Using new paired-end sequence data and a de novo repeat identification strategy, we identified the most repetitive elements in the E. tef genome. Putative repeat sequences were annotated based on similarity to known repeat groups in other grasses. Altogether we identified 1,389 medium/highly repetitive sequences that collectively represent about 27% of the teff genome. Phylogenetic analyses of the most important classes of TEs were carried out in a comparative framework including paralog elements from rice and maize. Finally, an abundant tandem repeat accounting for more than 4% of the whole genome was identified and partially characterized.

CONCLUSIONS

Analyzing a large sample of randomly sheared reads we obtained a library of the repetitive sequences of E. tef. The approach we used was designed to avoid underestimation of repeat contribution; such underestimation is characteristic of whole genome assembly projects. The data collected represent a valuable resource for further analysis of the genome of this important orphan crop.

Isolation and characterization of putative functional long terminal repeat retrotransposons in the Pyrus genome

Background

Long terminal repeat (LTR)-retrotransposons constitute 42.4 % of the genome of the ‘Suli’ pear (Pyrus pyrifolia white pear group), implying that retrotransposons have played important roles in Pyrus evolution. Therefore, further analysis of retrotransposons will enhance our understanding of the evolutionary history of Pyrus.

Results

We identified 1836 LTR-retrotransposons in the ‘Suli’ pear genome, of which 440 LTR-retrotransposons were predicted to contain at least two of three gene models (gag, integrase and reverse transcriptase). Because these were most likely to be functional transposons, we focused our analyses on this set of 440. Most of the LTR-retrotransposons were estimated to have inserted into the genome less than 2.5 million years ago. Sequence analysis showed that the reverse transcriptase component of the identified LTR-retrotransposons was highly heterogeneous. Analyses of transcripts assembled from RNA-Seq databases of two cultivars of Pyrus species showed that LTR-retrotransposons were expressed in the buds and fruit of Pyrus. A total of 734 coding sequences in the ‘Suli’ genome were disrupted by the identified LTR-retrotransposons. Five high-copy-number LTR-retrotransposon families were identified in Pyrus. These families were rarely found in the genomes of Malus and Prunus, but were distributed extensively in Pyrus and abundance varied between species.

Conclusions

We identified potentially functional, full-length LTR-retrotransposons with three gene models in the ‘Suli’ genome. The analysis of RNA-seq data demonstrated that these retrotransposons are expressed in the organs of pears. The differential copy number of LTR-retrotransposon families between Pyrus species suggests that the transposition of retrotransposons is an important evolutionary force driving the genetic divergence of species within the genus.

Electronic supplementary material

The online version of this article (doi:10.1186/s13100-016-0058-8) contains supplementary material, which is available to authorized users.

The effect of transposable elements on phenotypic variation: insights from plants to humans

Transposable elements (TEs), originally discovered in maize as controlling elements, are the main components of most eukaryotic genomes. TEs have been regarded as deleterious genomic parasites due to their ability to undergo massive amplification. However, TEs can regulate gene expression and alter phenotypes. Also, emerging findings demonstrate that TEs can establish and rewire gene regulatory networks by genetic and epigenetic mechanisms. In this review, we summarize the key roles of TEs in fine-tuning the regulation of gene expression leading to phenotypic plasticity in plants and humans, and the implications for adaption and natural selection.

Effects of genome structure variation, homeologous genes and repetitive DNA on polyploid crop research in the age of genomics

Compared to diploid species, allopolyploid crop species possess more complex genomes, higher productivity, and greater adaptability to changing environments. Next generation sequencing techniques have produced high-density genetic maps, whole genome sequences, transcriptomes and epigenomes for important polyploid crops. However, several problems interfere with the full application of next generation sequencing techniques to these crops. Firstly, different types of genomic variation affect sequence assembly and QTL mapping. Secondly, duplicated or homoeologous genes can diverge in function and then lead to emergence of many minor QTL, which increases difficulties in fine mapping, cloning and marker assisted selection. Thirdly, repetitive DNA sequences arising in polyploid crop genomes also impact sequence assembly, and are increasingly being shown to produce small RNAs to regulate gene expression and hence phenotypic traits. We propose that these three key features should be considered together when analyzing polyploid crop genomes. It is apparent that dissection of genomic structural variation, elucidation of the function and mechanism of interaction of homoeologous genes, and investigation of the de novo roles of repeat sequences in agronomic traits are necessary for genomics-based crop breeding in polyploids.

A Comprehensive Prescription for Plant miRNA Identification

microRNAs (miRNAs) are tiny ribo-regulatory molecules involved in various essential pathways for persistence of cellular life, such as development, environmental adaptation, and stress response. In recent years, miRNAs have become a major focus in molecular biology because of their functional and diagnostic importance. This interest in miRNA research has resulted in the development of many specific software and pipelines for the identification of miRNAs and their specific targets, which is the key for the elucidation of miRNA-modulated gene expression. While the well-recognized importance of miRNAs in clinical research pushed the emergence of many useful computational identification approaches in animals, available software and pipelines are fewer for plants. Additionally, existing approaches suffers from mis-identification and annotation of plant miRNAs since the miRNA mining process for plants is highly prone to false-positives, particularly in cereals which have a highly repetitive genome. Our group developed a homology-based in silico miRNA identification approach for plants, which utilizes two Perl scripts "SUmirFind" and "SUmirFold" and since then, this method helped identify many miRNAs particularly from crop species such as Triticum or Aegliops. Herein, we describe a comprehensive updated guideline by the implementation of two new scripts, "SUmirPredictor" and "SUmirLocator," and refinements to our previous method in order to identify genuine miRNAs with increased sensitivity in consideration of miRNA identification problems in plants. Recent updates enable our method to provide more reliable and precise results in an automated fashion in addition to solutions for elimination of most false-positive predictions, miRNA naming and miRNA mis-annotation. It also provides a comprehensive view to genome/transcriptome-wide location of miRNA precursors as well as their association with transposable elements. The "SUmirPredictor" and "SUmirLocator" scripts are freely available together with a reference high-confidence plant miRNA list.

The loss-of-function GLABROUS 3 mutation in cucumber is due to LTR-retrotransposon insertion in a class IV HD-ZIP transcription factor gene CsGL3 that is epistatic over CsGL1

BACKGROUND

Trichomes, developed from the protodermal cells (the outermost cell layer of the embryo), are hair-like structures covering the aerial parts of plants. The genetic network regulating trichome development has been extensively studied and well understood in the model species Arabidopsis thaliana, which bears unicellular, non-glandular and branched trichomes. However, little is known about the genetic and molecular basis of organogenesis of multi-cellular trichomes in plant species like cucumber (Cucumis sativus L.), which are likely different from Arabidopsis.

RESULTS

We identified a new trichome mutant in cucumber which exhibited a completely glabrous phenotype on all aerial organs. Genetic analysis indicated that the glabrous phenotype was inherited as a single recessive gene, csgl3. Fine genetic mapping delimited the csgl3 locus into a 68.4 kb region with 12 predicted genes. Genetic analysis, sequence alignment and allelic variation survey in natural populations identified Csa6G514870 encoding a class IV homeodomain-associated leucine zipper (HD-ZIP) transcription factor as the only candidate for CsGL3, which was 5188 bp in length with 10 predicted exons. Gene expression analysis revealed the loss-of-function of CsGL3 in the mutant due to the insertion of a 5-kb long terminal repeat (LTR) retrotransposon in the 4th exon of CsGL3. Linkage analysis in a segregating population and gene expression analysis of the CsGL1 and CsGL3 genes in csgl1, csgl3, and csgl1 + 3 genetic backgrounds uncovered interactions between the two genes. Phylogenetic analysis among 28 class IV HD-ZIP protein sequences from five species placed cucumber CsGL3 into the same clade with 7 other members that play important roles in trichome initiation.

CONCLUSIONS

The new glabrous mutation in cucumber was controlled by a single recessive locus csgl3, which was phenotypically and genetically distinct from two previously reported glabrous mutants csgl1 and csgl2. The glabrous phenotype in csgl3 was due to insertion of an autonomous, active, class I transposable element in CsGL3, a class IV HD-ZIP transcription factor. CsGL3 was epistatic to CsGL1. CsGL3 seemed to play important roles in cucumber trichome initiation whereas CsGL1 may act downstream in the trichome development pathway(s). Findings from the present study provide new insights into genetic control of trichome development in cucumber.

Genetic and genomic toolbox of Zea mays

Maize has a long history of genetic and genomic tool development and is considered one of the most accessible higher plant systems. With a fully sequenced genome, a suite of cytogenetic tools, methods for both forward and reverse genetics, and characterized phenotype markers, maize is amenable to studying questions beyond plant biology. Major discoveries in the areas of transposons, imprinting, and chromosome biology came from work in maize. Moving forward in the post-genomic era, this classic model system will continue to be at the forefront of basic biological study. In this review, we outline the basics of working with maize and describe its rich genetic toolbox.

The role of short RNA loops in recognition of a single-hairpin exon derived from a mammalian-wide interspersed repeat

Splice-site selection is controlled by secondary structure through sequestration or approximation of splicing signals in primary transcripts but the exact role of even the simplest and most prevalent structural motifs in exon recognition remains poorly understood. Here we took advantage of a single-hairpin exon that was activated in a mammalian-wide interspersed repeat (MIR) by a mutation stabilizing a terminal triloop, with splice sites positioned close to each other in a lower stem of the hairpin. We first show that the MIR exon inclusion in mRNA correlated inversely with hairpin stabilities. Employing a systematic manipulation of unpaired regions without altering splice-site configuration, we demonstrate a high correlation between exon inclusion of terminal tri- and tetraloop mutants and matching tri-/tetramers in splicing silencers/enhancers. Loop-specific exon inclusion levels and enhancer/silencer associations were preserved across primate cell lines, in 4 hybrid transcripts and also in the context of a distinct stem, but only if its loop-closing base pairs were shared with the MIR hairpin. Unlike terminal loops, splicing activities of internal loop mutants were predicted by their intramolecular Watson-Crick interactions with the antiparallel strand of the MIR hairpin rather than by frequencies of corresponding trinucleotides in splicing silencers/enhancers. We also show that splicing outcome of oligonucleotides targeting the MIR exon depend on the identity of the triloop adjacent to their antisense target. Finally, we identify proteins regulating MIR exon recognition and reveal a distinct requirement of adjacent exons for C-terminal extensions of Tra2α and Tra2β RNA recognition motifs.

Identification and characterization of a minisatellite contained within a novel miniature inverted-repeat transposable element (MITE) of Porphyromonas gingivalis

Background

Repetitive regions of DNA and transposable elements have been found to constitute large percentages of eukaryotic and prokaryotic genomes. Such elements are known to be involved in transcriptional regulation, host-pathogen interactions and genome evolution.

Results

We identified a minisatellite contained within a miniature inverted-repeat transposable element (MITE) in Porphyromonas gingivalis. The P. gingivalis minisatellite and associated MITE, named ‘BrickBuilt’, comprises a tandemly repeating twenty-three nucleotide DNA sequence lacking spacer regions between repeats, and with flanking ‘leader’ and ‘tail’ subunits that include small inverted-repeat ends. Forms of the BrickBuilt MITE are found 19 times in the genome of P. gingivalis strain ATCC 33277, and also multiple times within the strains W83, TDC60, HG66 and JCVI SC001. BrickBuilt is always located intergenically ranging between 49 and 591 nucleotides from the nearest upstream and downstream coding sequences. Segments of BrickBuilt contain promoter elements with bidirectional transcription capabilities.

Conclusions

We performed a bioinformatic analysis of BrickBuilt utilizing existing whole genome sequencing, microarray and RNAseq data, as well as performing in vitro promoter probe assays to determine potential roles, mechanisms and regulation of the expression of these elements and their affect on surrounding loci. The multiplicity, localization and limited host range nature of MITEs and MITE-like elements in P. gingivalis suggest that these elements may play an important role in facilitating genome evolution as well as modulating the transcriptional regulatory system.

Electronic supplementary material

The online version of this article (doi:10.1186/s13100-015-0049-1) contains supplementary material, which is available to authorized users.

A transposable element in a NAC gene is associated with drought tolerance in maize seedlings

Drought represents a major constraint on maize production worldwide. Understanding the genetic basis for natural variation in drought tolerance of maize may facilitate efforts to improve this trait in cultivated germplasm. Here, using a genome-wide association study, we show that a miniature inverted-repeat transposable element (MITE) inserted in the promoter of a NAC gene (ZmNAC111) is significantly associated with natural variation in maize drought tolerance. The 82-bp MITE represses ZmNAC111 expression via RNA-directed DNA methylation and H3K9 dimethylation when heterologously expressed in Arabidopsis. Increasing ZmNAC111 expression in transgenic maize enhances drought tolerance at the seedling stage, improves water-use efficiency and induces upregulation of drought-responsive genes under water stress. The MITE insertion in the ZmNAC111 promoter appears to have occurred after maize domestication and spread among temperate germplasm. The identification of this MITE insertion provides insight into the genetic basis for natural variation in maize drought tolerance.

Drought is a major cause of yield loss in maize and understanding the genetic determinants of natural variation in drought tolerance may aid breeding programs produce more tolerant varieties. Here, Mao et al. identify a MITE transposon insertion in a NAC transcription factor, which is associated with natural variation in drought tolerance.

The Role of Vertical and Horizontal Transfer in the Evolutionary Dynamics of PIF-Like Transposable Elements in Triticeae

PIF-like transposable elements are members of the PIF/Harbinger superfamily of DNA transposons found in the genomes of many plants, animals, and fungi. The evolution of the gene that encodes the transposase responsible for mobilizing PIF-like elements has been studied in both plants and animals, but the elements' history in flowering plants remains poorly known. In this work, we describe the phylogenetic distribution and evolution of PIF-like elements in the genomes of 21 diploid species from the wheat tribe, Triticeae, and we present the first convincing evidence of horizontal transfer of PIF elements in plant genomes. A phylogenetic analysis of 240 PIF sequences based on the conserved region of the transposase domain revealed at least four main transposase lineages. Their complex evolutionary history can be best explained by a combination of vertical transmission with differential evolutionary success among lineages, and occasional horizontal transfer between phylogenetically distant Triticeae genera. In addition, we identified 127 potentially functional transposase sequences indicating possible recent activity of PIF.

Precise repair of mPing excision sites is facilitated by target site duplication derived microhomology

BACKGROUND

A key difference between the Tourist and Stowaway families of miniature inverted repeat transposable elements (MITEs) is the manner in which their excision alters the genome. Upon excision, Stowaway-like MITEs and the associated Mariner elements usually leave behind a small duplication and short sequences from the end of the element. These small insertions or deletions known as "footprints" can potentially disrupt coding or regulatory sequences. In contrast, Tourist-like MITEs and the associated PIF/Pong/Harbinger elements generally excise precisely, returning the genome to its original state. The purpose of this study was to determine the mechanisms underlying these excision differences, including the role of the host DNA repair mechanisms.

RESULTS

The transposition of the Tourist-like element, mPing, and the Stowaway-like element, 14T32, were evaluated using yeast transposition assays. Assays performed in yeast strains lacking non-homologous end joining (NHEJ) enzymes indicated that the excision sites of both elements were primarily repaired by NHEJ. Altering the target site duplication (TSD) sequences that flank these elements reduced the transposition frequency. Using yeast strains with the ability to repair the excision site by homologous repair showed that some TSD changes disrupt excision of the element. Changing the ends of mPing to produce non-matching TSDs drastically reduced repair of the excision site and resulted in increased generation of footprints.

CONCLUSIONS

Together these results indicate that the difference in Tourist and Stowaway excision sites results from transposition mechanism characteristics. The TSDs of both elements play a role in element excision, but only the mPing TSDs actively participate in excision site repair. Our data suggests that Tourist-like elements excise with staggered cleavage of the TSDs, which provides microhomology that facilitates precise repair. This slight modification in the transposition mechanism results in more efficient repair of the double stranded break, and thus, may be less harmful to host genomes by disrupting fewer genes.

Genome-Wide Dosage-Dependent and -Independent Regulation Contributes to Gene Expression and Evolutionary Novelty in Plant Polyploids

Polyploidy provides evolutionary and morphological novelties in many plants and some animals. However, the role of genome dosage and composition in gene expression changes remains poorly understood. Here, we generated a series of resynthesized Arabidopsis tetraploids that contain 0-4 copies of Arabidopsis thaliana and Arabidopsis arenosa genomes and investigated ploidy and hybridity effects on gene expression. Allelic expression can be defined as dosage dependent (expression levels correlate with genome dosages) or otherwise as dosage independent. Here, we show that many dosage-dependent genes contribute to cell cycle, photosynthesis, and metabolism, whereas dosage-independent genes are enriched in biotic and abiotic stress responses. Interestingly, dosage-dependent genes tend to be preserved in ancient biochemical pathways present in both plant and nonplant species, whereas many dosage-independent genes belong to plant-specific pathways. This is confirmed by an independent analysis using Arabidopsis phylostratigraphic map. For A. thaliana loci, the dosage-dependent alleles are devoid of TEs and tend to correlate with H3K9ac, H3K4me3, and CG methylation, whereas the majority of dosage-independent alleles are enriched with TEs and correspond to H3K27me1, H3K27me3, and CHG (H = A, T, or C) methylation. Furthermore, there is a parent-of-origin effect on nonadditively expressed genes in the reciprocal allotetraploids especially when A. arenosa is used as the pollen donor, leading to metabolic and morphological changes. Thus, ploidy, epigenetic modifications, and cytoplasmic-nuclear interactions shape gene expression diversity in polyploids. Dosage-dependent expression can maintain growth and developmental stability, whereas dosage-independent expression can facilitate functional divergence between homeologs (subfunctionalization and/or neofunctionalization) during polyploid evolution.

Encapsidation of Host RNAs by Cucumber Necrosis Virus Coat Protein during both Agroinfiltration and Infection

Next-generation sequence analysis of virus-like particles (VLPs) produced during agroinfiltration of cucumber necrosis virus (CNV) coat protein (CP) and of authentic CNV virions was conducted to assess if host RNAs can be encapsidated by CNV CP. VLPs containing host RNAs were found to be produced during agroinfiltration, accumulating to approximately 1/60 the level that CNV virions accumulated during infection. VLPs contained a variety of host RNA species, including the major rRNAs as well as cytoplasmic, chloroplast, and mitochondrial mRNAs. The most predominant host RNA species encapsidated in VLPs were chloroplast encoded, consistent with the efficient targeting of CNV CP to chloroplasts during agroinfiltration. Interestingly, droplet digital PCR analysis showed that the CNV CP mRNA expressed during agroinfiltration was the most efficiently encapsidated mRNA, suggesting that the CNV CP open reading frame may contain a high-affinity site or sites for CP binding and thus contribute to the specificity of CNV RNA encapsidation. Approximately 0.09% to 0.7% of the RNA derived from authentic CNV virions contained host RNA, with chloroplast RNA again being the most prominent species. This is consistent with our previous finding that a small proportion of CNV CP enters chloroplasts during the infection process and highlights the possibility that chloroplast targeting is a significant aspect of CNV infection. Remarkably, 6 to 8 of the top 10 most efficiently encapsidated nucleus-encoded RNAs in CNV virions correspond to retrotransposon or retrotransposon-like RNA sequences. Thus, CNV could potentially serve as a vehicle for horizontal transmission of retrotransposons to new hosts and thereby significantly influence genome evolution.

IMPORTANCE

Viruses predominantly encapsidate their own virus-related RNA species due to the possession of specific sequences and/or structures on viral RNA which serve as high-affinity binding sites for the coat protein. In this study, we show, using next-generation sequence analysis, that CNV also encapsidates host RNA species, which account for ∼0.1% of the RNA packaged in CNV particles. The encapsidated host RNAs predominantly include chloroplast RNAs, reinforcing previous observations that CNV CP enters chloroplasts during infection. Remarkably, the most abundantly encapsidated cytoplasmic mRNAs consisted of retrotransposon-like RNA sequences, similar to findings recently reported for flock house virus (A. Routh, T. Domitrovic, and J. E. Johnson, Proc Natl Acad Sci U S A 109:1907-1912, 2012). Encapsidation of retrotransposon sequences may contribute to their horizontal transmission should CNV virions carrying retrotransposons infect a new host. Such an event could lead to large-scale genomic changes in a naive plant host, thus facilitating host evolutionary novelty.

Active transposable elements recover species boundaries and geographic structure in Madagascan coffee species

The completion of the genome assembly for the economically important coffee plant Coffea canephora (Rubiaceae) has allowed the use of bioinformatic tools to identify and characterize a diverse array of transposable elements (TEs), which can be used in evolutionary studies of the genus. An overview of the copy number and location within the C. canephora genome of four TEs is presented. These are tested for their use as molecular markers to unravel the evolutionary history of the Millotii Complex, a group of six wild coffee (Coffea) species native to Madagascar. Two TEs from the Gypsy superfamily successfully recovered some species boundaries and geographic structure among samples, whereas a TE from the Copia superfamily did not. Notably, species occurring in evergreen moist forests of eastern and southeastern Madagascar were divergent with respect to species in other habitats and regions. Our results suggest that the peak of transpositional activity of the Gypsy and Copia TEs occurred, respectively, before and after the speciation events of the tested Madagascan species. We conclude that the utilization of active TEs has considerable potential to unravel the evolutionary history and delimitation of closely related Coffea species. However, the selection of TE needs to be experimentally tested, since each element has its own evolutionary history. Different TEs with similar copy number in a given species can render different dendrograms; thus copy number is not a good selection criterion to attain phylogenetic resolution.

A new approach for annotation of transposable elements using small RNA mapping

Transposable elements (TEs) are mobile genomic DNA sequences found in most organisms. They so densely populate the genomes of many eukaryotic species that they are often the major constituents. With the rapid generation of many plant genome sequencing projects over the past few decades, there is an urgent need for improved TE annotation as a prerequisite for genome-wide studies. Analogous to the use of RNA-seq for gene annotation, we propose a new method for de novo TE annotation that uses as a guide 24 nt-siRNAs that are a part of TE silencing pathways. We use this new approach, called TASR (for Transposon Annotation using Small RNAs), for de novo annotation of TEs in Arabidopsis, rice and soybean and demonstrate that this strategy can be successfully applied for de novo TE annotation in plants.Executable PERL is available for download from: http://tasr-pipeline.sourceforge.net/.

Diversity, abundance, and evolutionary dynamics of Pong-like transposable elements in Triticeae

Pong-like elements are members of the PIF/Harbinger superfamily of DNA transposons that has been described in many plants, animals, and fungi. Most Pong elements contain two open reading frames (ORFs). One encodes a transposase (ORF2) that catalyzes transposition of Pong and related non-autonomous elements, while the function of the second is unknown. Little is known about the evolutionary history of Pong elements in flowering plants. In this work, we present the first comprehensive analysis of the diversity, abundance, and evolution of the Pong-like transposase gene in the genomes of 21 diploid species from the wheat tribe, Triticeae, and we present the first convincing evidence of horizontal transfer of nuclear-encoded Pong elements in any organism. A phylogenetic analysis of nearly 300 Pong sequences based on a conserved region of the transposase domain revealed a complex evolutionary history of Pong elements that can be best explained by ancestral polymorphism, followed by differential evolutionary success of some transposase lineages, and by occasional horizontal transfer between phylogenetically distant genera. In addition, we used transposon display to estimate the abundance of the transposase gene within Triticeae genomes, and our results revealed varying levels of Pong proliferation, with numbers of transposase copies ranging from 22 to 92. Comparisons of Pong transposase abundance to flow cytometry estimates of genome size revealed that larger Triticeae genome size was not correlated with transposase abundance.

Short-term hybridisation activates Tnt1 and Tto1 Copia retrotransposons in wild tuber-bearing Solanum species

Interspecific hybridisation in tuber-bearing species of Solanum is a common phenomenon and represents an important source of variability, crucial for adaptation and speciation of potato species. In this regard, the effects of interspecific hybridisation on retrotransposon families present in the genomes, and their consequent effects on generation of genetic variability in wild tuber-bearing Solanum species, are poorly characterised. The aim of this study was to analyse the activity of retrotransposons in inter- and intraspecific hybrids between S. kurtzianum and S. microdontum, obtained by controlled crosses, and the effects on morphological, genetic and epigenetic variability. For genetic and epigenetic analysis, S-SAP (sequence-specific amplification polymorphism) and TMD (transposon methylation display) techniques were used, respectively, with specific primers for Tnt1 and Tto1 retrotransposon families (Order LTR, Superfamily Copia). The results indicate that at morphological level, interspecific hybrid genotypes differ from their parental species, whereas derived intraspecific hybrids do not. In both cases, we observed significant reductions in pollen grain viability, and a negative correlation with Tnt1 mobility. Both retrotransposons, Tto1 and Tnt1, were mobilised in the genotypes analysed, with mobility ranging from 0 to 7.8%. Furthermore, at the epigenetic level, demethylation was detected in the vicinity of Tnt1 and Tto1 in the hybrids compared with the parental genotypes. These patterns were positively correlated with the activity of the retrotransposons. The results suggest a possible mechanism through which hybridisation events generate genetic variability in tuber-bearing species of Solanum through retrotranposon activation.

The diversification of PHIS transposon superfamily in eukaryotes

BACKGROUND

PHIS transposon superfamily belongs to DNA transposons and includes PIF/Harbinger, ISL2EU, and Spy transposon groups. These three groups have similar DDE domain-containing transposases; however, their coding capacity, species distribution, and target site duplications (TSDs) are significantly different.

RESULTS

In this study, we systematically identified and analyzed PHIS transposons in 836 sequenced eukaryotic genomes using transposase homology search and structure approach. In total, 380 PHIS families were identified in 112 genomes and 168 of 380 families were firstly reported in this study. Besides previous identified PIF/Harbinger, ISL2EU, and Spy groups, three new types (called Pangu, NuwaI, and NuwaII) of PHIS superfamily were identified; each has its own distinctive characteristics, especially in TSDs. Pangu and NuwaII transposons are characterized by 5'-ANT-3' and 5'-C|TNA|G-3' TSDs, respectively. Both transposons are widely distributed in plants, fungi, and animals; the NuwaI transposons are characterized by 5'-CWG-3' TSDs and mainly distributed in animals.

CONCLUSIONS

Here, in total, 380 PHIS families were identified in eukaryotes. Among these 380 families, 168 were firstly reported in this study. Furthermore, three new types of PHIS superfamily were identified. Our results not only enrich the transposon diversity but also have extensive significance for improving genome sequence assembly and annotation of higher organisms.

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.

Identification, Diversity and Evolution of MITEs in the Genomes of Microsporidian Nosema Parasites

Miniature inverted-repeat transposable elements (MITEs) are short, non-autonomous DNA transposons, which are widespread in most eukaryotic genomes. However, genome-wide identification, origin and evolution of MITEs remain largely obscure in microsporidia. In this study, we investigated structural features for de novo identification of MITEs in genomes of silkworm microsporidia Nosema bombycis and Nosema antheraeae, as well as a honeybee microsporidia Nosema ceranae. A total of 1490, 149 and 83 MITE-related sequences from 89, 17 and five families, respectively, were found in the genomes of the above-mentioned species. Species-specific MITEs are predominant in each genome of microsporidian Nosema, with the exception of three MITE families that were shared by N. bombycis and N. antheraeae. One or multiple rounds of amplification occurred for MITEs in N. bombycis after divergence between N. bombycis and the other two species, suggesting that the more abundant families in N. bombycis could be attributed to the recent amplification of new MITEs. Significantly, some MITEs that inserted into the homologous protein-coding region of N. bombycis were recruited as introns, indicating that gene expansion occurred during the evolution of microsporidia. NbS31 and NbS24 had polymorphisms in different geographical strains of N. bombycis, indicating that they could still be active. In addition, several small RNAs in the MITEs in N. bombycis are mainly produced from both ends of the MITEs sequence.

Genetic and epigenetic changes in somatic hybrid introgression lines between wheat and tall wheatgrass

Broad phenotypic variations were induced in derivatives of an asymmetric somatic hybridization of bread wheat (Triticum aestivum) and tall wheatgrass (Thinopyrum ponticum Podp); however, how these variations occurred was unknown. We explored the nature of these variations by cytogenetic assays and DNA profiling techniques to characterize six genetically stable somatic introgression lines. Karyotyping results show the six lines similar to their wheat parent, but GISH analysis identified the presence of a number of short introgressed tall wheatgrass chromatin segments. DNA profiling revealed many genetic and epigenetic differences, including sequences deletions, altered regulation of gene expression, changed patterns of cytosine methylation, and the reactivation of retrotransposons. Phenotypic variations appear to result from altered repetitive sequences combined with the epigenetic regulation of gene expression and/or retrotransposon transposition. The extent of genetic and epigenetic variation due to the maintenance of parent wheat cells in tissue culture was assessed and shown to be considerably lower than had been induced in the introgression lines. Asymmetric somatic hybridization provides appropriate material to explore the nature of the genetic and epigenetic variations induced by genomic shock.

Construction of a linkage map based on retrotransposon insertion polymorphisms in sweetpotato via high-throughput sequencing

Sweetpotato (Ipomoea batatas L.) is an outcrossing hexaploid species with a large number of chromosomes (2n = 6x = 90). Although sweetpotato is one of the world's most important crops, genetic analysis of the species has been hindered by its genetic complexity combined with the lack of a whole genome sequence. In the present study, we constructed a genetic linkage map based on retrotransposon insertion polymorphisms using a mapping population derived from a cross between 'Purple Sweet Lord' (PSL) and '90IDN-47' cultivars. High-throughput sequencing and subsequent data analyses identified many Rtsp-1 retrotransposon insertion sites, and their allele dosages (simplex, duplex, triplex, or double-simplex) were determined based on segregation ratios in the mapping population. Using a pseudo-testcross strategy, 43 and 47 linkage groups were generated for PSL and 90IDN-47, respectively. Interestingly, most of these insertions (~90%) were present in a simplex manner, indicating their utility for linkage map construction in polyploid species. Additionally, our approach led to savings of time and labor for genotyping. Although the number of markers herein was insufficient for map-based cloning, our trial analysis exhibited the utility of retrotransposon-based markers for linkage map construction in sweetpotato.

Characterization of expressed sequence tags from Lilium longiflorum in vernalized and non-vernalized bulbs

In Lilium longiflorum, vernalization is both an obligatory requirement and the major factor affecting flowering time, however, little is known about the molecular regulation of this mechanism in Lilium and other flowering bulbs. Exposure of L. longiflorum bulbs to 9 weeks at 4°C greatly promoted stem elongation within the bulb, floral transition and flowering. Subtraction libraries of vernalized (V) and non-vernalized (NV) bulb meristems were constructed. 671 and 479 genes were sequenced, from which 72 and 82 proteins were inferred for the NV-V and the V-NV libraries, respectively. Much lower transcription levels and putative gene functions were recorded in the NV-V libraries compared the V-NV libraries. However, a large number of genes annotated to transposable elements (TEs), represented more than 20% of the sequenced cDNA were expressed in the NV-V libraries, as opposed to less than 2% in the V-NV libraries. The expression profile of several genes potentially involved in the vernalization pathway was assessed. Expression of LlSOC1, the lily homologue of SUPPRESSOR OF OVER-EXPRESSION OF CO1 (SOC1), an important flowering gene in several plant species, found in the V-NV library, was highly up-regulated during bulb meristem cold exposure. The subtraction libraries provided a fast tool for relevant gene isolation.

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

Transposable elements (TEs) account for a large portion of the genome in many eukaryotic species. Despite their reputation as “junk” DNA or genomic parasites deleterious for the host, TEs have complex interactions with host genes and the potential to contribute to regulatory variation in gene expression. It has been hypothesized that TEs and genes they insert near may be transcriptionally activated in response to stress conditions. The maize genome, with many different types of TEs interspersed with genes, provides an ideal system to study the genome-wide influence of TEs on gene regulation. To analyze the magnitude of the TE effect on gene expression response to environmental changes, we profiled gene and TE transcript levels in maize seedlings exposed to a number of abiotic stresses. Many genes exhibit up- or down-regulation in response to these stress conditions. The analysis of TE families inserted within upstream regions of up-regulated genes revealed that between four and nine different TE families are associated with up-regulated gene expression in each of these stress conditions, affecting up to 20% of the genes up-regulated in response to abiotic stress, and as many as 33% of genes that are only expressed in response to stress. Expression of many of these same TE families also responds to the same stress conditions. The analysis of the stress-induced transcripts and proximity of the transposon to the gene suggests that these TEs may provide local enhancer activities that stimulate stress-responsive gene expression. Our data on allelic variation for insertions of several of these TEs show strong correlation between the presence of TE insertions and stress-responsive up-regulation of gene expression. Our findings suggest that TEs provide an important source of allelic regulatory variation in gene response to abiotic stress in maize.

Transposable elements are mobile DNA elements that are a prevalent component of many eukaryotic genomes. While transposable elements can often have deleterious effects through insertions into protein-coding genes they may also contribute to regulatory variation of gene expression. There are a handful of examples in which specific transposon insertions contribute to regulatory variation of nearby genes, particularly in response to environmental stress. We sought to understand the genome-wide influence of transposable elements on gene expression responses to abiotic stress in maize, a plant with many families of transposable elements located in between genes. Our analysis suggests that a small number of maize transposable element families may contribute to the response of nearby genes to abiotic stress by providing stress-responsive enhancer-like functions. The specific insertions of transposable elements are often polymorphic within a species. Our data demonstrate that allelic variation for insertions of the transposable elements associated with stress-responsive expression can contribute to variation in the regulation of nearby genes. Thus novel insertions of transposable elements provide a potential mechanism for genes to acquire cis-regulatory influences that could contribute to heritable variation for stress response.

From genes to markers: exploiting gene sequence information to develop tools for plant breeding

Once the sequence is known for a gene of interest, it is usually possible to design markers to detect polymorphisms within the gene. Such markers can be particularly useful in plant breeding, especially if they detect the causal polymorphism within the gene and are diagnostic of the phenotype. In this chapter, we (1) discuss how gene sequences are obtained and aligned and how polymorphic sites can be identified or predicted; (2) explain the principles of PCR primer design and PCR amplification and provide guidelines for their application in the design and testing of markers; (3) discuss detection methods for presence/absence (dominant) polymorphisms, length polymorphisms and single nucleotide polymorphisms (SNPs); and (4) outline some of the factors that affect the utility of markers in plant breeding and explain how markers can be evaluated (validated) for use in plant breeding.

A small RNA mediated regulation of a stress-activated retrotransposon and the tissue specific transposition during the reproductive period in Arabidopsis

Transposable elements (TEs) are key elements that facilitate genome evolution of the host organism. A number of studies have assessed the functions of TEs, which change gene expression in the host genome. Activation of TEs is controlled by epigenetic modifications such as DNA methylation and histone modifications. Several recent studies have reported that TEs can also be activated by biotic or abiotic stress in some plants. We focused on a Ty1/copia retrotransposon, ONSEN, that is activated by heat stress (HS) in Arabidopsis. We found that transcriptional activation of ONSEN was regulated by a small interfering RNA (siRNA)-related pathway, and the activation could also be induced by oxidative stress. Mutants deficient in siRNA biogenesis that were exposed to HS at the initial stages of vegetative growth showed transgenerational transposition. The transposition was also detected in the progeny, which originated from tissue that had differentiated after exposure to the HS. The results indicated that in some undifferentiated cells, transpositional activity could be maintained quite long after exposure to the HS.

Evolutionary dynamics of hAT DNA transposon families in Saccharomycetaceae

Transposable elements (TEs) are widespread in eukaryotes but uncommon in yeasts of the Saccharomycotina subphylum, in terms of both host species and genome fraction. The class II elements are especially scarce, but the hAT element Rover is a noteworthy exception that deserves further investigation. Here, we conducted a genome-wide analysis of hAT elements in 40 ascomycota. A novel family, Roamer, was found in three species, whereas Rover was detected in 15 preduplicated species from Kluyveromyces, Eremothecium, and Lachancea genera, with up to 41 copies per genome. Rover acquisition seems to have occurred by horizontal transfer in a common ancestor of these genera. The detection of remote Rover copies in Naumovozyma dairenensis and in the sole Saccharomyces cerevisiae strain AWRI1631, without synteny, suggests that two additional independent horizontal transfers took place toward these genomes. Such patchy distribution of elements prevents any anticipation of TE presence in incoming sequenced genomes, even closely related ones. The presence of both putative autonomous and defective Rover copies, as well as their diversification into five families, indicate particular dynamics of Rover elements in the Lachancea genus. Especially, we discovered the first miniature inverted-repeat transposable elements (MITEs) to be described in yeasts, together with their parental autonomous copies. Evidence of MITE insertion polymorphism among Lachancea waltii strains suggests their recent activity. Moreover, 40% of Rover copies appeared to be involved in chromosome rearrangements, showing the large structural impact of TEs on yeast genome and opening the door to further investigations to understand their functional and evolutionary consequences.

Genome-wide development of transposable elements-based markers in foxtail millet and construction of an integrated database

Transposable elements (TEs) are major components of plant genome and are reported to play significant roles in functional genome diversity and phenotypic variations. Several TEs are highly polymorphic for insert location in the genome and this facilitates development of TE-based markers for various genotyping purposes. Considering this, a genome-wide analysis was performed in the model plant foxtail millet. A total of 30,706 TEs were identified and classified as DNA transposons (24,386), full-length Copia type (1,038), partial or solo Copia type (10,118), full-length Gypsy type (1,570), partial or solo Gypsy type (23,293) and Long- and Short-Interspersed Nuclear Elements (3,659 and 53, respectively). Further, 20,278 TE-based markers were developed, namely Retrotransposon-Based Insertion Polymorphisms (4,801, ∼24%), Inter-Retrotransposon Amplified Polymorphisms (3,239, ∼16%), Repeat Junction Markers (4,451, ∼22%), Repeat Junction-Junction Markers (329, ∼2%), Insertion-Site-Based Polymorphisms (7,401, ∼36%) and Retrotransposon-Microsatellite Amplified Polymorphisms (57, 0.2%). A total of 134 Repeat Junction Markers were screened in 96 accessions of Setaria italica and 3 wild Setaria accessions of which 30 showed polymorphism. Moreover, an open access database for these developed resources was constructed (Foxtail millet Transposable Elements-based Marker Database; http://59.163.192.83/ltrdb/index.html). Taken together, this study would serve as a valuable resource for large-scale genotyping applications in foxtail millet and related grass species.

TE-Tracker: systematic identification of transposition events through whole-genome resequencing

Background

Transposable elements (TEs) are DNA sequences that are able to move from their location in the genome by cutting or copying themselves to another locus. As such, they are increasingly recognized as impacting all aspects of genome function. With the dramatic reduction in cost of DNA sequencing, it is now possible to resequence whole genomes in order to systematically characterize novel TE mobilization in a particular individual. However, this task is made difficult by the inherently repetitive nature of TE sequences, which in some eukaryotes compose over half of the genome sequence. Currently, only a few software tools dedicated to the detection of TE mobilization using next-generation-sequencing are described in the literature. They often target specific TEs for which annotation is available, and are only able to identify families of closely related TEs, rather than individual elements.

Results

We present TE-Tracker, a general and accurate computational method for the de-novo detection of germ line TE mobilization from re-sequenced genomes, as well as the identification of both their source and destination sequences. We compare our method with the two classes of existing software: specialized TE-detection tools and generic structural variant (SV) detection tools. We show that TE-Tracker, while working independently of any prior annotation, bridges the gap between these two approaches in terms of detection power. Indeed, its positive predictive value (PPV) is comparable to that of dedicated TE software while its sensitivity is typical of a generic SV detection tool. TE-Tracker demonstrates the benefit of adopting an annotation-independent, de novo approach for the detection of TE mobilization events. We use TE-Tracker to provide a comprehensive view of transposition events induced by loss of DNA methylation in Arabidopsis. TE-Tracker is freely available at http://www.genoscope.cns.fr/TE-Tracker.

Conclusions

We show that TE-Tracker accurately detects both the source and destination of novel transposition events in re-sequenced genomes. Moreover, TE-Tracker is able to detect all potential donor sequences for a given insertion, and can identify the correct one among them. Furthermore, TE-Tracker produces significantly fewer false positives than common SV detection programs, thus greatly facilitating the detection and analysis of TE mobilization events.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0377-z) contains supplementary material, which is available to authorized users.

Miniature inverted-repeat transposable elements (MITEs) in rice were originated and amplified predominantly after the divergence of Oryza and Brachypodium and contributed considerable diversity to the species

Miniature inverted-repeat transposable elements (MITEs) are ubiquitous in high eukaryotic genomes. More than 178,000 MITE sequences of 338 families are present in the genome of rice (Oryza sativa) cultivar Nipponbare. Interestingly, only two of the 338 MITE families have homologous sequences in the genome of Brachypodium distachyon, a relative in the grass family. Therefore, the vast majority of MITEs in the rice genome were originated and amplified after the divergence of Oryza and Brachypodium. Comparison between rice cultivar Nipponbare and another rice cultivar 93–11 showed 14.8% of MITEs exhibit presence/absence (P/A) polymorphism. The P/A polymorphism was mainly attributed to recent MITE transpositions, while less than 10% of the P/A polymorphism was caused by MITE excisions. Therefore, the high P/A polymorphisms of MITEs may generate considerable gene expression and phenotypic diversity for O. sativa.

The origin and evolution of six miniature inverted-repeat transposable elements in Bombyx mori and Rhodnius prolixus

Miniature inverted-repeat transposable elements (MITEs) are a specific group of nonautonomous DNA transposons, and they are distributed in a wide range of hosts. However, the origin and evolutionary history of MITEs in eukaryotic genomes remain unclear. In this study, six MITEs were identified in the silkworm (Bombyx mori). Five elements are grouped into four known superfamilies of DNA transposons, and one represents a novel class of MITEs. Unexpectedly, six similar MITEs are also present in the triatomine bug (Rhodnius prolixus) that diverged from the common ancestor with the silkworm about 370 Ma. However, they show different lengths in two species, suggesting that they are different derivatives of progenitor transposons. Three direct progenitor transposons (Sola1, hobo/Ac/Tam [hAT], and Ginger2) are also identified in some other organisms, and several lines of evidence suggested that these autonomous elements might have been independently and horizontally transferred into their hosts. Furthermore, it is speculated that the twisted-wing parasites may be the candidate vectors for these horizontal transfers. The data presented in this study provide some new insights into the origin and evolutionary history of MITEs in the silkworm and triatomine bug.

Epigenetic regulation and functional exaptation of transposable elements in higher plants

Transposable elements (TEs) are mobile genetic elements that can proliferate in their host genomes. Because of their robust amplification, TEs have long been considered 'selfish DNA', harmful insertions that can threaten host genome integrity. The idea of TEs as junk DNA comes from analysis of epigenetic silencing of their mobility in plants and animals. This idea contrasts with McClintock's characterization of TEs as 'controlling elements'. Emerging studies on the regulatory functions of TEs in plant genomes have updated McClintock's characterization, indicating exaptation of TEs for genetic regulation. In this review, we summarize recent progress in TE silencing, particularly in Arabidopsis and rice, and show that TEs provide an abundant, natural source of regulation for the host genome.

A non-autonomous insect piggyBac transposable element is mobile in tobacco

The piggyBac transposable element, originally isolated from a virus in an insect cell line, is a valuable molecular tool for transgenesis and mutagenesis of invertebrates. For heterologous transgenesis in a variety of mammals, transfer of the piggyBac transposable element from an ectopic plasmid only requires expression of piggyBac transposase. To determine if piggyBac could function in dicotyledonous plants, a two-element system was developed in tobacco (Nicotiana tabacum) to test for transposable element excision and insertion. The first transgenic line constitutively expressed piggyBac transposase, while the second transgenic line contained at least two non-autonomous piggyBac transposable elements. Progeny from crosses of the two transgenic lines was analyzed for piggyBac excision and transposition. Several progeny displayed excision events, and all the sequenced excision sites exhibited evidence of the precise excision mechanism characteristic of piggyBac transposase. Two unique transposition insertion events were identified that each included diagnostic duplication of the target site. These data indicate that piggyBac transposase is active in a dicotyledonous plant, although at a low frequency.