Comparison of class 2 transposable elements at superfamily resolution reveals conserved and distinct features in cereal grass genomes

Comparative Analysis of Transposable Elements in Strawberry Genomes of Different Ploidy Levels

Transposable elements (TEs) make up a large portion of plant genomes and play a vital role in genome structure, function, and evolution. Cultivated strawberry (Fragaria x ananassa) is one of the most important fruit crops, and its octoploid genome was formed through several rounds of genome duplications from diploid ancestors. Here, we built a pan-genome TE library for the Fragaria genus using ten published strawberry genomes at different ploidy levels, including seven diploids, one tetraploid, and two octoploids, and performed comparative analysis of TE content in these genomes. The TEs comprise 51.83% (F. viridis) to 60.07% (F. nilgerrensis) of the genomes. Long terminal repeat retrotransposons (LTR-RTs) are the predominant TE type in the Fragaria genomes (20.16% to 34.94%), particularly in F. iinumae (34.94%). Estimating TE content and LTR-RT insertion times revealed that species-specific TEs have shaped each strawberry genome. Additionally, the copy number of different LTR-RT families inserted in the last one million years reflects the genetic distance between Fragaria species. Comparing cultivated strawberry subgenomes to extant diploid ancestors showed that F. vesca and F. iinumae are likely the diploid ancestors of the cultivated strawberry, but not F. viridis. These findings provide new insights into the TE variations in the strawberry genomes and their roles in strawberry genome evolution.

Rates and spectra of de novo structural mutations in Chlamydomonas reinhardtii

Genetic variation originates from several types of spontaneous mutation, including single-nucleotide substitutions, short insertions and deletions (indels), and larger structural changes. Structural mutations (SMs) drive genome evolution and are thought to play major roles in evolutionary adaptation, speciation, and genetic disease, including cancers. Sequencing of mutation accumulation (MA) lines has provided estimates of rates and spectra of single-nucleotide and indel mutations in many species, yet the rate of new SMs is largely unknown. Here, we use long-read sequencing to determine the full mutation spectrum in MA lines derived from two strains (CC-1952 and CC-2931) of the green alga Chlamydomonas reinhardtii The SM rate is highly variable between strains and between MA lines, and SMs represent a substantial proportion of all mutations in both strains (CC-1952 6%; CC-2931 12%). The SM spectra differ considerably between the two strains, with almost all inversions and translocations occurring in CC-2931 MA lines. This variation is associated with heterogeneity in the number and type of active transposable elements (TEs), which comprise major proportions of SMs in both strains (CC-1952 22%; CC-2931 38%). In CC-2931, a Crypton and a previously undescribed type of DNA element have caused 71% of chromosomal rearrangements, whereas in CC-1952, a Dualen LINE is associated with 87% of duplications. Other SMs, notably large duplications in CC-2931, are likely products of various double-strand break repair pathways. Our results show that diverse types of SMs occur at substantial rates, and support prominent roles for SMs and TEs in evolution.

Resetting of the 24-nt siRNA landscape in rice zygotes

The zygote, a totipotent stem cell, is crucial to the life cycle of sexually reproducing organisms. It is produced by the fusion of two differentiated cells-the egg and sperm, which in plants have radically different siRNA transcriptomes from each other and from multicellular embryos. Owing to technical challenges, the epigenetic changes that accompany the transition from differentiated gametes to totipotent zygote are poorly understood. Because siRNAs serve as both regulators and outputs of the epigenome, we characterized small RNA transcriptomes of zygotes from rice. Zygote small RNAs exhibit extensive maternal carryover and an apparent lack of paternal contribution, indicated by absence of sperm signature siRNAs. Zygote formation is accompanied by widespread redistribution of 24-nt siRNAs relative to gametes, such that ∼70% of the zygote siRNA loci do not overlap any egg cell siRNA loci. Newly detected siRNA loci in zygote are gene-proximal and not associated with centromeric heterochromatin, similar to canonical siRNAs, in sharp contrast to gametic siRNA loci that are gene-distal and heterochromatic. In addition, zygote but not egg siRNA loci are associated with high DNA methylation in the mature embryo. Thus, the zygote begins transitioning before the first embryonic division to an siRNA profile that is associated with future RdDM in embryogenesis. These findings indicate that, in addition to changes in gene expression, the transition to totipotency in the plant zygote is accompanied by resetting of the epigenetic reprogramming that occurred during gamete formation.

The genomic ecosystem of transposable elements in maize

Transposable elements (TEs) constitute the majority of flowering plant DNA, reflecting their tremendous success in subverting, avoiding, and surviving the defenses of their host genomes to ensure their selfish replication. More than 85% of the sequence of the maize genome can be ascribed to past transposition, providing a major contribution to the structure of the genome. Evidence from individual loci has informed our understanding of how transposition has shaped the genome, and a number of individual TE insertions have been causally linked to dramatic phenotypic changes. Genome-wide analyses in maize and other taxa have frequently represented TEs as a relatively homogeneous class of fragmentary relics of past transposition, obscuring their evolutionary history and interaction with their host genome. Using an updated annotation of structurally intact TEs in the maize reference genome, we investigate the family-level dynamics of TEs in maize. Integrating a variety of data, from descriptors of individual TEs like coding capacity, expression, and methylation, as well as similar features of the sequence they inserted into, we model the relationship between attributes of the genomic environment and the survival of TE copies and families. In contrast to the wholesale relegation of all TEs to a single category of junk DNA, these differences reveal a diversity of survival strategies of TE families. Together these generate a rich ecology of the genome, with each TE family representing the evolution of a distinct ecological niche. We conclude that while the impact of transposition is highly family- and context-dependent, a family-level understanding of the ecology of TEs in the genome can refine our ability to predict the role of TEs in generating genetic and phenotypic diversity.

Assessing the regulatory potential of transposable elements using chromatin accessibility profiles of maize transposons

Transposable elements (TEs) have the potential to create regulatory variation both through the disruption of existing DNA regulatory elements and through the creation of novel DNA regulatory elements. In a species with a large genome, such as maize, many TEs interspersed with genes create opportunities for significant allelic variation due to TE presence/absence polymorphisms among individuals. We used information on putative regulatory elements in combination with knowledge about TE polymorphisms in maize to identify TE insertions that interrupt existing accessible chromatin regions (ACRs) in B73 as well as examples of polymorphic TEs that contain ACRs among four inbred lines of maize including B73, Mo17, W22, and PH207. The TE insertions in three other assembled maize genomes (Mo17, W22, or PH207) that interrupt ACRs that are present in the B73 genome can trigger changes to the chromatin, suggesting the potential for both genetic and epigenetic influences of these insertions. Nearly 20% of the ACRs located over 2 kb from the nearest gene are located within an annotated TE. These are regions of unmethylated DNA that show evidence for functional importance similar to ACRs that are not present within TEs. Using a large panel of maize genotypes, we tested if there is an association between the presence of TE insertions that interrupt, or carry, an ACR and the expression of nearby genes. While most TE polymorphisms are not associated with expression for nearby genes, the TEs that carry ACRs exhibit enrichment for being associated with higher expression of nearby genes, suggesting that these TEs may contribute novel regulatory elements. These analyses highlight the potential for a subset of TEs to rewire transcriptional responses in eukaryotic genomes.

Identification and characterization of inheritable structural variations induced by ion beam radiations in rice

High energy ion beams are effective physical mutagens for mutation induction in plants. Due to their high linear energy transfer (LET) property, they are known to generate single nucleotide variations (SNVs) and insertion/deletions (InDels, <50 bp) as well as structural variations (SVs). However, due to the technical difficulties to identify SVs, studies on ion beam induced SVs by genome sequencing have so far been limited in numbers and inadequate in nature, and knowledge of SVs is scarce with regards to their characteristics. In the present study, we identified and validated SVs in six M₄ plants (designated as Ar_50, Ar_100, C_150, C_200, Ne_50 and Ne_100 according to ion beam types and irradiation doses), two each induced by argon (⁴⁰Ar¹⁸⁺), carbon (¹²C⁶⁺) and neon (²⁰Ne¹⁰⁺) ion beams and performed in depth analyses of their characteristics. In total, 22 SVs were identified and validated, consisting of 11 deletions, 1 duplication, and 4 intra-chromosomal and 6 inter-chromosomal translocations. There were several SVs larger than 1 kbp. The SVs were distributed across the whole genome with an aggregation with SNVs and InDels only in the Ne_50 mutants. An enrichment of a 11-bp wide G-rich DNA motif 'GAAGGWGGRGG' was identified around the SV breakpoints. Three mechanisms might be involved in the SV formation, i.e., the expansion of tandem repeats, transposable element insertion, and non-allelic homologous recombination. Put together, the present study provides a preliminary view of SVs induced by Ar, C and Ne ion beam radiations, and as a pilot study, it contributes to our understanding of how SVs might form after ion beam irradiation in rice.

Perspective: 50 years of plant chromosome biology

The past 50 years has been the greatest era of plant science discovery, and most of the discoveries have emerged from or been facilitated by our knowledge of plant chromosomes. At last we have descriptive and mechanistic outlines of the information in chromosomes that programs plant life. We had almost no such information 50 years ago when few had isolated DNA from any plant species. The important features of genes have been revealed through whole genome comparative genomics and testing of variants using transgenesis. Progress has been enabled by the development of technologies that had to be invented and then become widely available. Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) have played extraordinary roles as model species. Unexpected evolutionary dramas were uncovered when learning that chromosomes have to manage constantly the vast numbers of potentially mutagenic families of transposons and other repeated sequences. The chromatin-based transcriptional and epigenetic mechanisms that co-evolved to manage the evolutionary drama as well as gene expression and 3-D nuclear architecture have been elucidated these past 20 years. This perspective traces some of the major developments with which I have become particularly familiar while seeking ways to improve crop plants. I draw some conclusions from this look-back over 50 years during which the scientific community has (i) exposed how chromosomes guard, readout, control, recombine, and transmit information that programs plant species, large and small, weed and crop, and (ii) modified the information in chromosomes for the purposes of genetic, physiological, and developmental analyses and plant improvement.

Genome-wide redistribution of 24-nt siRNAs in rice gametes

Gametes constitute a critical stage of the plant life cycle during which the genome undergoes reprogramming in preparation for embryogenesis. Here, we examined genome-wide distributions of small RNAs in the sperm and egg cells of rice. We found that 24-nt siRNAs, which are a hallmark of RNA-directed DNA methylation (RdDM) in plants, were depleted from heterochromatin boundaries in both gametes relative to vegetative tissues, reminiscent of siRNA patterns in DDM1-type nucleosome remodeler mutants. In sperm cells, 24-nt siRNAs were spread across heterochromatic regions, while in egg cells, 24-nt siRNAs were concentrated at a smaller number of heterochromatic loci throughout the genome, especially at loci which also produced siRNAs in other tissues. In both gametes, patterns of CHH methylation, typically a strong indicator of RdDM, were similar to vegetative tissues, although lower in magnitude. These findings indicate that the small RNA transcriptome undergoes large-scale redistribution in both male and female gametes, which is not correlated with recruitment of DNA methyltransferases in gametes and suggestive of unexplored regulatory activities of gamete small RNAs.

Transposition of a non-autonomous DNA transposon in the gene coding for a bHLH transcription factor results in a white bulb color of onions (Allium cepa L.)

A DNA transposon was found in the gene encoding a bHLH transcription factor. Genotypes of the marker tagging this DNA transposon perfectly co-segregated with color phenotypes in large F_2:3 populations A combined approach of bulked segregant analysis and RNA-Seq was used to isolate causal gene for C locus controlling white bulb color in onions (Allium cepa L.). A total of 114 contigs containing homozygous single nucleotide polymorphisms (SNPs) between white and yellow bulked RNAs were identified. Four of them showed high homologies with loci clustered in the middle of chromosome 5. SNPs in 34 contigs were confirmed by sequencing of PCR products. One of these contigs showed perfect linkage to the C locus in F_2:3 populations consisting of 2491 individuals. However, genotypes of molecular marker tagging this contig were inconsistent with color phenotypes of diverse breeding lines. A total of 146 contigs showed differential expression between yellow and white bulks. Among them, transcription levels of B2 gene encoding a bHLH transcription factor were significantly reduced in white RNA bulk and F_2:3 individuals, although there was no SNP in the coding region. Phylogenetic analysis showed that onion B2 was orthologous to bHLH-coding genes regulating anthocyanin biosynthesis pathway in other plant species. Promoter regions of B2 gene were obtained by genome walking and a 577-bp non-autonomous DNA transposon designated as AcWHITE was found in the white allele. Molecular marker tagging AcWHITE showed perfect linkage with the C locus. Marker genotypes of the white allele were detected in some white accessions. However, none of tested red or yellow onions contained AcWHITE insertion, implying that B2 gene was likely to be a casual gene for the C locus.

Benchmarking transposable element annotation methods for creation of a streamlined, comprehensive pipeline

BACKGROUND

Sequencing technology and assembly algorithms have matured to the point that high-quality de novo assembly is possible for large, repetitive genomes. Current assemblies traverse transposable elements (TEs) and provide an opportunity for comprehensive annotation of TEs. Numerous methods exist for annotation of each class of TEs, but their relative performances have not been systematically compared. Moreover, a comprehensive pipeline is needed to produce a non-redundant library of TEs for species lacking this resource to generate whole-genome TE annotations.

RESULTS

We benchmark existing programs based on a carefully curated library of rice TEs. We evaluate the performance of methods annotating long terminal repeat (LTR) retrotransposons, terminal inverted repeat (TIR) transposons, short TIR transposons known as miniature inverted transposable elements (MITEs), and Helitrons. Performance metrics include sensitivity, specificity, accuracy, precision, FDR, and F₁. Using the most robust programs, we create a comprehensive pipeline called Extensive de-novo TE Annotator (EDTA) that produces a filtered non-redundant TE library for annotation of structurally intact and fragmented elements. EDTA also deconvolutes nested TE insertions frequently found in highly repetitive genomic regions. Using other model species with curated TE libraries (maize and Drosophila), EDTA is shown to be robust across both plant and animal species.

CONCLUSIONS

The benchmarking results and pipeline developed here will greatly facilitate TE annotation in eukaryotic genomes. These annotations will promote a much more in-depth understanding of the diversity and evolution of TEs at both intra- and inter-species levels. EDTA is open-source and freely available: https://github.com/oushujun/EDTA.

Benchmarking transposable element annotation methods for creation of a streamlined, comprehensive pipeline

BACKGROUND

Sequencing technology and assembly algorithms have matured to the point that high-quality de novo assembly is possible for large, repetitive genomes. Current assemblies traverse transposable elements (TEs) and provide an opportunity for comprehensive annotation of TEs. Numerous methods exist for annotation of each class of TEs, but their relative performances have not been systematically compared. Moreover, a comprehensive pipeline is needed to produce a non-redundant library of TEs for species lacking this resource to generate whole-genome TE annotations.

RESULTS

We benchmark existing programs based on a carefully curated library of rice TEs. We evaluate the performance of methods annotating long terminal repeat (LTR) retrotransposons, terminal inverted repeat (TIR) transposons, short TIR transposons known as miniature inverted transposable elements (MITEs), and Helitrons. Performance metrics include sensitivity, specificity, accuracy, precision, FDR, and F1. Using the most robust programs, we create a comprehensive pipeline called Extensive de-novo TE Annotator (EDTA) that produces a filtered non-redundant TE library for annotation of structurally intact and fragmented elements. EDTA also deconvolutes nested TE insertions frequently found in highly repetitive genomic regions. Using other model species with curated TE libraries (maize and Drosophila), EDTA is shown to be robust across both plant and animal species.

CONCLUSIONS

The benchmarking results and pipeline developed here will greatly facilitate TE annotation in eukaryotic genomes. These annotations will promote a much more in-depth understanding of the diversity and evolution of TEs at both intra- and inter-species levels. EDTA is open-source and freely available: https://github.com/oushujun/EDTA.

The bracteatus pineapple genome and domestication of clonally propagated crops

Domestication of clonally propagated crops such as pineapple from South America was hypothesized to be a 'one-step operation'. We sequenced the genome of Ananas comosus var. bracteatus CB5 and assembled 513 Mb into 25 chromosomes with 29,412 genes. Comparison of the genomes of CB5, F153 and MD2 elucidated the genomic basis of fiber production, color formation, sugar accumulation and fruit maturation. We also resequenced 89 Ananas genomes. Cultivars 'Smooth Cayenne' and 'Queen' exhibited ancient and recent admixture, while 'Singapore Spanish' supported a one-step operation of domestication. We identified 25 selective sweeps, including a strong sweep containing a pair of tandemly duplicated bromelain inhibitors. Four candidate genes for self-incompatibility were linked in F153, but were not functional in self-compatible CB5. Our findings support the coexistence of sexual recombination and a one-step operation in the domestication of clonally propagated crops. This work guides the exploration of sexual and asexual domestication trajectories in other clonally propagated crops.

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.

Comparative analysis of miniature inverted-repeat transposable elements (MITEs) and long terminal repeat (LTR) retrotransposons in six Citrus species

BACKGROUND

Miniature inverted-repeat transposable elements (MITEs) and long terminal repeat (LTR) retrotransposons are ubiquitous in plants genomes, and highly important in their evolution and diversity. However, their mechanisms of insertion/amplification and roles in Citrus genome's evolution/diversity are still poorly understood.

RESULTS

To address this knowledge gap, we developed different computational pipelines to analyze, annotate and classify MITEs and LTR retrotransposons in six different sequenced Citrus species. We identified 62,010 full-length MITEs from 110 distinguished families. We observed MITEs tend to insert in gene related regions and enriched in promoters. We found that DTM63 is possibly an active Mutator-like MITE family in the traceable past and may still be active in Citrus. The insertion of MITEs resulted in massive polymorphisms and played an important role in Citrus genome diversity and gene structure variations. In addition, 6630 complete LTR retrotransposons and 13,371 solo-LTRs were identified. Among them, 12 LTR lineages separated before the differentiation of mono- and dicotyledonous plants. We observed insertion and deletion of LTR retrotransposons was accomplished with a dynamic balance, and their half-life in Citrus was ~ 1.8 million years.

CONCLUSIONS

These findings provide insights into MITEs and LTR retrotransposons and their roles in genome diversity in different Citrus genomes.

A miniature inverted-repeat transposable element, AddIn-MITE, located inside a WD40 gene is conserved in Andropogoneae grasses

Miniature inverted-repeat transposable elements (MITEs) have been associated with genic regions in plant genomes and may play important roles in the regulation of nearby genes via recruitment of small RNAs (sRNA) to the MITEs loci. We identified eight families of MITEs in the sugarcane genome assembly with MITE-Hunter pipeline. These sequences were found to be upstream, downstream or inserted into 67 genic regions in the genome. The position of the most abundant MITE (Stowaway-like) in genic regions, which we call AddIn-MITE, was confirmed in a WD40 gene. The analysis of four monocot species showed conservation of the AddIn-MITE sequence, with a large number of copies in their genomes. We also investigated the conservation of the AddIn-MITE’ position in the WD40 genes from sorghum, maize and, in sugarcane cultivars and wild Saccharum species. In all analyzed plants, AddIn-MITE has located in WD40 intronic region. Furthermore, the role of AddIn-MITE-related sRNA in WD40 genic region was investigated. We found sRNAs preferentially mapped to the AddIn-MITE than to other regions in the WD40 gene in sugarcane. In addition, the analysis of the small RNA distribution patterns in the WD40 gene and the structure of AddIn-MITE, suggests that the MITE region is a proto-miRNA locus in sugarcane. Together, these data provide insights into the AddIn-MITE role in Andropogoneae grasses.

QTL Mapping Combined With Comparative Analyses Identified Candidate Genes for Reduced Shattering in Setaria italica

Setaria (L.) P. Beauv is a genus of grasses that belongs to the Poaceae (grass) family, subfamily Panicoideae. Two members of the Setaria genus, Setaria italica (foxtail millet) and S. viridis (green foxtail), have been studied extensively over the past few years as model species for C4-photosynthesis and to facilitate genome studies in complex Panicoid bioenergy grasses. We exploited the available genetic and genomic resources for S. italica and its wild progenitor, S. viridis, to study the genetic basis of seed shattering. Reduced shattering is a key trait that underwent positive selection during domestication. Phenotyping of F2:3 and recombinant inbred line (RIL) populations generated from a cross between S. italica accession B100 and S. viridis accession A10 identified the presence of additive main effect quantitative trait loci (QTL) on chromosomes V and IX. As expected, enhanced seed shattering was contributed by the wild S. viridis. Comparative analyses pinpointed Sh1 and qSH1, two shattering genes previously identified in sorghum and rice, as potentially underlying the QTL on Setaria chromosomes IX and V, respectively. The Sh1 allele in S. italica was shown to carry a PIF/Harbinger MITE in exon 2, which gave rise to an alternatively spliced transcript that lacked exon 2. This MITE was universally present in S. italica accessions around the world and absent from the S. viridis germplasm tested, strongly suggesting a single origin of foxtail millet domestication. The qSH1 gene carried two MITEs in the 5'UTR. Presence of one or both MITEs was strongly associated with cultivated germplasm. If the MITE insertion(s) in qSH1 played a role in reducing shattering in S. italica accessions, selection for the variants likely occurred after the domestication of foxtail millet.

Distribution, Diversity, and Long-Term Retention of Grass Short Interspersed Nuclear Elements (SINEs)

Instances of highly conserved plant short interspersed nuclear element (SINE) families and their enrichment near genes have been well documented, but little is known about the general patterns of such conservation and enrichment and underlying mechanisms. Here, we perform a comprehensive investigation of the structure, distribution, and evolution of SINEs in the grass family by analyzing 14 grass and 5 other flowering plant genomes using comparative genomics methods. We identify 61 SINE families composed of 29,572 copies, in which 46 families are first described. We find that comparing with other grass TEs, grass SINEs show much higher level of conservation in terms of genomic retention: The origin of at least 26% families can be traced to early grass diversification and these families are among most abundant SINE families in 86% species. We find that these families show much higher level of enrichment near protein coding genes than families of relatively recent origin (51%:28%), and that 40% of all grass SINEs are near gene and the percentage is higher than other types of grass TEs. The pattern of enrichment suggests that differential removal of SINE copies in gene-poor regions plays an important role in shaping the genomic distribution of these elements. We also identify a sequence motif located at 3' SINE end which is shared in 17 families. In short, this study provides insights into structure and evolution of SINEs in the grass family.

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.

Duplication of host genes by transposable elements

The availability of large amounts of genomic and transcriptome sequences have allowed systematic surveys about the host gene sequences that have been duplicated by transposable elements. It is now clear that all super-families of transposons are capable of duplicating genes or gene fragments, and such incidents have been detected in a wide spectrum of organisms. Emerging evidence suggests that a considerable portion of them function as coding or non-coding sequences, driving innovations at molecular and phenotypic levels. Interestingly, the duplication events not only have to occur in the reproductive tissues to become heritable, but the duplicated copies are also preferentially expressed in those tissues. As a result, reproductive tissues may serve as the 'incubator' for genes generated by transposable elements.

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.

The repetitive landscape of the 5100 Mbp barley genome

Background

While transposable elements (TEs) comprise the bulk of plant genomic DNA, how they contribute to genome structure and organization is still poorly understood. Especially in large genomes where TEs make the majority of genomic DNA, it is still unclear whether TEs target specific chromosomal regions or whether they simply accumulate where they are best tolerated.

Results

Here, we present an analysis of the repetitive fraction of the 5100 Mb barley genome, the largest angiosperm genome to have a near-complete sequence assembly. Genes make only about 2% of the genome, while over 80% is derived from TEs. The TE fraction is composed of at least 350 different families. However, 50% of the genome is comprised of only 15 high-copy TE families, while all other TE families are present in moderate or low copy numbers. We found that the barley genome is highly compartmentalized with different types of TEs occupying different chromosomal “niches”, such as distal, interstitial, or proximal regions of chromosome arms. Furthermore, gene space represents its own distinct genomic compartment that is enriched in small non-autonomous DNA transposons, suggesting that these TEs specifically target promoters and downstream regions. Furthermore, their presence in gene promoters is associated with decreased methylation levels.

Conclusions

Our data show that TEs are major determinants of overall chromosome structure. We hypothesize that many of the the various chromosomal distribution patterns are the result of TE families targeting specific niches, rather than them accumulating where they have the least deleterious effects.

Electronic supplementary material

The online version of this article (10.1186/s13100-017-0102-3) contains supplementary material, which is available to authorized users.

Tracking the genome-wide outcomes of a transposable element burst over decades of amplification

Rice (Oryza sativa) has a unique combination of attributes that made it an ideal host to track the natural behavior of very active transposable elements (TEs) over generations. In this study, we have exploited its small genome and propagation by self or sibling pollination to identify and characterize two strain pairs, EG4/HEG4 and A119/A123, undergoing bursts of the nonautonomous miniature inverted repeat transposable element mPing. Comparative sequence analyses of these strains have advanced our understanding of (i) factors that contribute to sustaining a TE burst for decades, (ii) features that distinguish a natural TE burst from bursts in cell culture or mutant backgrounds, and (iii) the extent to which TEs can rapidly diversify the genome of an inbred organism.

To understand the success strategies of transposable elements (TEs) that attain high copy numbers, we analyzed two pairs of rice (Oryza sativa) strains, EG4/HEG4 and A119/A123, undergoing decades of rapid amplification (bursts) of the class 2 autonomous Ping element and the nonautonomous miniature inverted repeat transposable element (MITE) mPing. Comparative analyses of whole-genome sequences of the two strain pairs validated that each pair has been maintained for decades as inbreds since divergence from their respective last common ancestor. Strains EG4 and HEG4 differ by fewer than 160 SNPs and a total of 264 new mPing insertions. Similarly, strains A119 and A123 exhibited about half as many SNPs (277) as new mPing insertions (518). Examination of all other potentially active TEs in these genomes revealed only a single new insertion out of ∼40,000 loci surveyed. The virtual absence of any new TE insertions in these strains outside the mPing bursts demonstrates that the Ping/mPing family gradually attains high copy numbers by maintaining activity and evading host detection for dozens of generations. Evasion is possible because host recognition of mPing sequences appears to have no impact on initiation or maintenance of the burst. Ping is actively transcribed, and both Ping and mPing can transpose despite methylation of terminal sequences. This finding suggests that an important feature of MITE success is that host recognition does not lead to the silencing of the source of transposase.

Transposon-mediated epigenetic regulation contributes to phenotypic diversity and environmental adaptation in rice

Transposable elements (TEs) have long been regarded as 'selfish DNA', and are generally silenced by epigenetic mechanisms. However, work in the past decade has identified positive roles for TEs in generating genomic novelty and diversity in plants. In particular, recent studies suggested that TE-induced epigenetic alterations and modification of gene expression contribute to phenotypic variation and adaptation to geography or stress. These findings have led many to regard TEs, not as junk DNA, but as sources of control elements and genomic diversity. As a staple food crop and model system for genomic research on monocot plants, rice (Oryza sativa) has a modest-sized genome that harbors massive numbers of DNA transposons (class II transposable elements) scattered across the genome, which may make TE regulation of genes more prevalent. In this review, we summarize recent progress in research on the functions of rice TEs in modulating gene expression and creating new genes. We also examine the contributions of TEs to phenotypic diversity and adaptation to environmental conditions.

Epigenetic regulation and epigenomic landscape in rice

Epigenetic regulation has been implicated in the control of complex agronomic traits in rice (Oryza sativa), a staple food crop and model monocot plant. Recent advances in high-throughput sequencing and the moderately complex genome of rice have made it possible to study epigenetic regulation in rice on a genome-wide scale. This review discusses recent advances in our understanding of epigenetic regulation in rice, with an emphasis on the roles of key epigenetic regulators, the epigenomic landscape, epigenetic variation, transposon repression, and plant development.

Transposable element influences on gene expression in plants

Transposable elements (TEs) comprise a major portion of many plant genomes and bursts of TE movements cause novel genomic variation within species. In order to maintain proper gene function, plant genomes have evolved a variety of mechanisms to tolerate the presence of TEs within or near genes. Here, we review our understanding of the interactions between TEs and gene expression in plants by assessing three ways that transposons can influence gene expression. First, there is growing evidence that TE insertions within introns or untranslated regions of genes are often tolerated and have minimal impact on expression level or splicing. However, there are examples in which TE insertions within genes can result in aberrant or novel transcripts. Second, TEs can provide novel alternative promoters, which can lead to new expression patterns or original coding potential of an alternate transcript. Third, TE insertions near genes can influence regulation of gene expression through a variety of mechanisms. For example, TEs may provide novel cis-acting regulatory sites behaving as enhancers or insert within existing enhancers to influence transcript production. Alternatively, TEs may change chromatin modifications in regions near genes, which in turn can influence gene expression levels. Together, the interactions of genes and TEs provide abundant evidence for the role of TEs in changing basic functions within plant genomes beyond acting as latent genomic elements or as simple insertional mutagens. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

The First Rule of Plant Transposable Element Silencing: Location, Location, Location

Transposable elements (TEs) are mobile units of DNA that comprise large portions of plant genomes. Besides creating mutations via transposition and contributing to genome size, TEs play key roles in chromosome architecture and gene regulation. TE activity is repressed by overlapping mechanisms of chromatin condensation, epigenetic transcriptional silencing, and targeting by small interfering RNAs. The specific regulation of different TEs, as well as their different roles in chromosome architecture and gene regulation, is specified by where on the chromosome the TE is located: near a gene, within a gene, in a pericentromere/TE island, or at the centromere core. In this Review, we investigate the silencing mechanisms responsible for inhibiting TE activity for each of these chromosomal contexts, emphasizing that chromosomal location is the first rule dictating the specific regulation of each TE.

What makes up plant genomes: The vanishing line between transposable elements and genes

The ultimate source of evolution is mutation. As the largest component in plant genomes, transposable elements (TEs) create numerous types of mutations that cannot be mimicked by other genetic mechanisms. When TEs insert into genomic sequences, they influence the expression of nearby genes as well as genes unlinked to the insertion. TEs can duplicate, mobilize, and recombine normal genes or gene fragments, with the potential to generate new genes or modify the structure of existing genes. TEs also donate their transposase coding regions for cellular functions in a process called TE domestication. Despite the host defense against TE activity, a subset of TEs survived and thrived through discreet selection of transposition activity, target site, element size, and the internal sequence. Finally, TEs have established strategies to reduce the efficacy of host defense system by increasing the cost of silencing TEs. This review discusses the recent progress in the area of plant TEs with a focus on the interaction between TEs and genes.

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.

Contrasting allelic distribution of CO/Hd1 homologues in Miscanthus sinensis from the East Asian mainland and the Japanese archipelago

The genus Miscanthus is a perennial C4 grass native to eastern Asia and is a promising candidate bioenergy crop for cool temperate areas. Flowering time is a crucial factor governing regional and seasonal adaptation; in addition, it is also a key target trait for extending the vegetative phase to improve biomass potential. Homologues of CONSTANS (CO)/Heading date 1(Hd1) were cloned from Miscanthus sinensis and named MsiHd1. Sequences of MsiHd1 homologues were compared among 24 wild M. sinensis accessions from Japan, 14 from China, and three from South Korea. Two to five MsiHd1 alleles in each accession were identified, suggesting that MsiHd1 consists of at least three loci in the Miscanthus genome. Verifying the open reading frame in MsiHd1, they were classified as putative functional alleles without mutations or non-functional alleles caused by indels. The Neighbor-Joining tree indicated that one of the multiple MsiHd1 loci is a pseudogene locus without any functional alleles. The pseudogene locus was named MsiHd1b, and the other loci were considered to be part of the MsiHd1a multi-locus family. Interestingly, in most Japanese accessions 50% or more of the MsiHd1a alleles were non-functional, whereas accessions from the East Asian mainland harboured only functional alleles. Five novel miniature inverted transposable elements (MITEs) (MsiMITE1-MsiMITE5) were observed in MsiHd1a/b. MsiMITE1, detected in exon 1 of MsiHd1a, was only observed in Japanese accessions and its revertant alleles derived from retransposition were predominantly in Chinese accessions. These differences in MsiHd1a show that the dependency on functional MsiHd1a alleles is different between accessions from the East Asian mainland and Japan.

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

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

An RNA sequencing transcriptome analysis of the high-temperature stressed tall fescue reveals novel insights into plant thermotolerance

BACKGROUND

Tall fescue (Festuca arundinacea Schreb.) is major cool-season forage and turf grass species worldwide, but high-temperature is a major environmental stress that dramatically threaten forage production and turf management of tall fescue. However, very little is known about the whole-genome molecular mechanisms contributing to thermotolerance. The objectives of this study were to analyzed genome-wide gene expression profiles in the leaves of two tall fescue genotypes, heat tolerant 'PI578718' and heat sensitive 'PI234881' using high-throughput RNA sequencing.

RESULTS

A total of 262 million high-quality paired-end reads were generated and assembled into 31,803 unigenes with an average length of 1,840 bp. Of these, 12,974 unigenes showed different expression patterns in response to heat stress and were categorized into 49 Gene Ontology functional subcategories. In addition, the variance of enrichment degree in each functional subcategory between PI578718 and PI234881 increased with increasing treatment time. Cell division and cell cycle genes showed a massive increase in transcript abundance in heat-stressed plants and more activated genes were detected in PI 578718 by Kyoto Encyclopedia of Genes and Genomes pathways analysis. Low molecular weight heat shock protein (LMW-HSP, HSP20) showed activated in two stressed genotypes and high molecular weight HSP (HMW-HSP, HSP90) just in PI578718. Assimilation such as photosynthesis, carbon fixation, CH4, N, S metabolism decreased along with increased dissimilation such as oxidative phosphorylation.

CONCLUSIONS

The assembled transcriptome of tall fescue could serve as a global description of expressed genes and provide more molecular resources for future functional characterization analysis of genomics in cool-season turfgrass in response to high-temperature. Increased cell division, LMW/HMW-HSP, dissimilation and antioxidant transcript amounts in tall fescue were correlated with successful resistance to high temperature stress.

Accessible DNA and relative depletion of H3K9me2 at maize loci undergoing RNA-directed DNA methylation

RNA-directed DNA methylation (RdDM) in plants is a well-characterized example of RNA interference-related transcriptional gene silencing. To determine the relationships between RdDM and heterochromatin in the repeat-rich maize (Zea mays) genome, we performed whole-genome analyses of several heterochromatic features: dimethylation of lysine 9 and lysine 27 (H3K9me2 and H3K27me2), chromatin accessibility, DNA methylation, and small RNAs; we also analyzed two mutants that affect these processes, mediator of paramutation1 and zea methyltransferase2. The data revealed that the majority of the genome exists in a heterochromatic state defined by inaccessible chromatin that is marked by H3K9me2 and H3K27me2 but that lacks RdDM. The minority of the genome marked by RdDM was predominantly near genes, and its overall chromatin structure appeared more similar to euchromatin than to heterochromatin. These and other data indicate that the densely staining chromatin defined as heterochromatin differs fundamentally from RdDM-targeted chromatin. We propose that small interfering RNAs perform a specialized role in repressing transposons in accessible chromatin environments and that the bulk of heterochromatin is incompatible with small RNA production.

Building the sugarcane genome for biotechnology and identifying evolutionary trends

Background

Sugarcane is the source of sugar in all tropical and subtropical countries and is becoming increasingly important for bio-based fuels. However, its large (10 Gb), polyploid, complex genome has hindered genome based breeding efforts. Here we release the largest and most diverse set of sugarcane genome sequences to date, as part of an on-going initiative to provide a sugarcane genomic information resource, with the ultimate goal of producing a gold standard genome.

Results

Three hundred and seventeen chiefly euchromatic BACs were sequenced. A reference set of one thousand four hundred manually-annotated protein-coding genes was generated. A small RNA collection and a RNA-seq library were used to explore expression patterns and the sRNA landscape. In the sucrose and starch metabolism pathway, 16 non-redundant enzyme-encoding genes were identified. One of the sucrose pathway genes, sucrose-6-phosphate phosphohydrolase, is duplicated in sugarcane and sorghum, but not in rice and maize. A diversity analysis of the s6pp duplication region revealed haplotype-structured sequence composition. Examination of hom(e)ologous loci indicate both sequence structural and sRNA landscape variation. A synteny analysis shows that the sugarcane genome has expanded relative to the sorghum genome, largely due to the presence of transposable elements and uncharacterized intergenic and intronic sequences.

Conclusion

This release of sugarcane genomic sequences will advance our understanding of sugarcane genetics and contribute to the development of molecular tools for breeding purposes and gene discovery.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-540) contains supplementary material, which is available to authorized users.

Early embryogenesis-specific expression of the rice transposon Ping enhances amplification of the MITE mPing

Miniature inverted-repeat transposable elements (MITEs) are numerically predominant transposable elements in the rice genome, and their activities have influenced the evolution of genes. Very little is known about how MITEs can rapidly amplify to thousands in the genome. The rice MITE mPing is quiescent in most cultivars under natural growth conditions, although it is activated by various stresses, such as tissue culture, gamma-ray irradiation, and high hydrostatic pressure. Exceptionally in the temperate japonica rice strain EG4 (cultivar Gimbozu), mPing has reached over 1000 copies in the genome, and is amplifying owing to its active transposition even under natural growth conditions. Being the only active MITE, mPing in EG4 is an appropriate material to study how MITEs amplify in the genome. Here, we provide important findings regarding the transposition and amplification of mPing in EG4. Transposon display of mPing using various tissues of a single EG4 plant revealed that most de novo mPing insertions arise in embryogenesis during the period from 3 to 5 days after pollination (DAP), and a large majority of these insertions are transmissible to the next generation. Locus-specific PCR showed that mPing excisions and insertions arose at the same time (3 to 5 DAP). Moreover, expression analysis and in situ hybridization analysis revealed that Ping, an autonomous partner for mPing, was markedly up-regulated in the 3 DAP embryo of EG4, whereas such up-regulation of Ping was not observed in the mPing-inactive cultivar Nipponbare. These results demonstrate that the early embryogenesis-specific expression of Ping is responsible for the successful amplification of mPing in EG4. This study helps not only to elucidate the whole mechanism of mPing amplification but also to further understand the contribution of MITEs to genome evolution.

Transposable elements are major components of eukaryotic genomes, comprising a large portion of the genome in some species. Miniature inverted-repeat transposable elements (MITEs), which belong to the class II DNA transposable elements, are abundant in gene-rich regions, and their copy numbers are very high; therefore, they have been considered to contribute to genome evolution. Because MITEs are short and have no coding capacity, they cannot transpose their positions without the aid of transposase, provided in trans by their autonomous element(s). It has been unknown how MITEs amplify themselves to high copy numbers in the genome. Our results demonstrate that the rice active MITE mPing is mobilized in the embryo by the developmental stage-specific up-regulation of an autonomous element, Ping, and thereby successfully amplifies itself to a high copy number in the genome. The short-term expression of Ping is thought to be a strategy of the mPing family for amplifying mPing by escaping the silencing mechanism of the host genome.

mRNA and Small RNA Transcriptomes Reveal Insights into Dynamic Homoeolog Regulation of Allopolyploid Heterosis in Nascent Hexaploid Wheat

Nascent allohexaploid wheat may represent the initial genetic state of common wheat (Triticum aestivum), which arose as a hybrid between Triticum turgidum (AABB) and Aegilops tauschii (DD) and by chromosome doubling and outcompeted its parents in growth vigor and adaptability. To better understand the molecular basis for this success, we performed mRNA and small RNA transcriptome analyses in nascent allohexaploid wheat and its following generations, their progenitors, and the natural allohexaploid cultivar Chinese Spring, with the assistance of recently published A and D genome sequences. We found that nonadditively expressed protein-coding genes were rare but relevant to growth vigor. Moreover, a high proportion of protein-coding genes exhibited parental expression level dominance, with genes for which the total homoeolog expression level in the progeny was similar to that in T. turgidum potentially participating in development and those with similar expression to that in Ae. tauschii involved in adaptation. In addition, a high proportion of microRNAs showed nonadditive expression upon polyploidization, potentially leading to differential expression of important target genes. Furthermore, increased small interfering RNA density was observed for transposable element-associated D homoeologs in the allohexaploid progeny, which may account for biased repression of D homoeologs. Together, our data provide insights into small RNA-mediated dynamic homoeolog regulation mechanisms that may contribute to heterosis in nascent hexaploid wheat.

Advances in genome studies in plants and animals

The area of plant and animal genomics covers the entire suite of issues in biology because it aims to determine the structure and function of genetic material. Although specific issues define research advances at an organism level, it is evident that many of the fundamental features of genome structure and the translation of encoded information to function share common ground. The Plant and Animal Genome (PAG) conference held in San Diego (California), in January each year provides an overview across all organisms at the genome level, and often it is evident that investments in the human area provide leadership, applications, and discoveries for researchers studying other organisms. This mini-review utilizes the plenary lectures as a basis for summarizing the trends in the genome-level studies of organisms, and the lectures include presentations by Ewan Birney (EBI, UK), Eric Green (NIH, USA), John Butler (NIST, USA), Elaine Mardis (Washington, USA), Caroline Dean (John Innes Centre, UK), Trudy Mackay (NC State University, USA), Sue Wessler (UC Riverside, USA), and Patrick Wincker (Genoscope, France). The work reviewed is based on published papers. Where unpublished information is cited, permission to include the information in this manuscript was obtained from the presenters.

Dicer-like 3 produces transposable element-associated 24-nt siRNAs that control agricultural traits in rice

Transposable elements (TEs) and repetitive sequences make up over 35% of the rice (Oryza sativa) genome. The host regulates the activity of different TEs by different epigenetic mechanisms, including DNA methylation, histone H3K9 methylation, and histone H3K4 demethylation. TEs can also affect the expression of host genes. For example, miniature inverted repeat TEs (MITEs), dispersed high copy-number DNA TEs, can influence the expression of nearby genes. In plants, 24-nt small interfering RNAs (siRNAs) are mainly derived from repeats and TEs. However, the extent to which TEs, particularly MITEs associated with 24-nt siRNAs, affect gene expression remains elusive. Here, we show that the rice Dicer-like 3 homolog OsDCL3a is primarily responsible for 24-nt siRNA processing. Impairing OsDCL3a expression by RNA interference caused phenotypes affecting important agricultural traits; these phenotypes include dwarfism, larger flag leaf angle, and fewer secondary branches. We used small RNA deep sequencing to identify 535,054 24-nt siRNA clusters. Of these clusters, ∼82% were OsDCL3a-dependent and showed significant enrichment of MITEs. Reduction of OsDCL3a function reduced the 24-nt siRNAs predominantly from MITEs and elevated expression of nearby genes. OsDCL3a directly targets genes involved in gibberellin and brassinosteroid homeostasis; OsDCL3a deficiency may affect these genes, thus causing the phenotypes of dwarfism and enlarged flag leaf angle. Our work identifies OsDCL3a-dependent 24-nt siRNAs derived from MITEs as broadly functioning regulators for fine-tuning gene expression, which may reflect a conserved epigenetic mechanism in higher plants with genomes rich in dispersed repeats or TEs.

Organization and evolution of transposable elements along the bread wheat chromosome 3B

BACKGROUND

The 17 Gb bread wheat genome has massively expanded through the proliferation of transposable elements (TEs) and two recent rounds of polyploidization. The assembly of a 774 Mb reference sequence of wheat chromosome 3B provided us with the opportunity to explore the impact of TEs on the complex wheat genome structure and evolution at a resolution and scale not reached so far.

RESULTS

We develop an automated workflow, CLARI-TE, for TE modeling in complex genomes. We delineate precisely 56,488 intact and 196,391 fragmented TEs along the 3B pseudomolecule, accounting for 85% of the sequence, and reconstruct 30,199 nested insertions. TEs have been mostly silent for the last one million years, and the 3B chromosome has been shaped by a succession of bursts that occurred between 1 to 3 million years ago. Accelerated TE elimination in the high-recombination distal regions is a driving force towards chromosome partitioning. CACTAs overrepresented in the high-recombination distal regions are significantly associated with recently duplicated genes. In addition, we identify 140 CACTA-mediated gene capture events with 17 genes potentially created by exon shuffling and show that 19 captured genes are transcribed and under selection pressure, suggesting the important role of CACTAs in the recent wheat adaptation.

CONCLUSION

Accurate TE modeling uncovers the dynamics of TEs in a highly complex and polyploid genome. It provides novel insights into chromosome partitioning and highlights the role of CACTA transposons in the high level of gene duplication in wheat.

P-MITE: a database for plant miniature inverted-repeat transposable elements

Miniature inverted-repeat transposable elements (MITEs) are prevalent in eukaryotic species including plants. MITE families vary dramatically and usually cannot be identified based on homology. In this study, we de novo identified MITEs from 41 plant species, using computer programs MITE Digger, MITE-Hunter and/or Repetitive Sequence with Precise Boundaries (RSPB). MITEs were found in all, but one (Cyanidioschyzon merolae), species. Combined with the MITEs identified previously from the rice genome, >2.3 million sequences from 3527 MITE families were obtained from 41 plant species. In general, higher plants contain more MITEs than lower plants, with a few exceptions such as papaya, with only 538 elements. The largest number of MITEs is found in apple, with 237 302 MITE sequences. The number of MITE sequences in a genome is significantly correlated with genome size. A series of databases (plant MITE databases, P-MITE), available online at http://pmite.hzau.edu.cn/django/mite/, was constructed to host all MITE sequences from the 41 plant genomes. The databases are available for sequence similarity searches (BLASTN), and MITE sequences can be downloaded by family or by genome. The databases can be used to study the origin and amplification of MITEs, MITE-derived small RNAs and roles of MITEs on gene and genome evolution.