Contrasting evolutionary patterns and target specificities among three Tourist-like MITE families in the maize genome

Insertion of a miniature inverted-repeat transposable element into the promoter of OsTCP4 results in more tillers and a lower grain size in rice

A class I PCF type protein, TCP4, was identified as a transcription factor associated with both grain size and tillering through a DNA pull-down-MS assay combined with a genome-wide association study. This transcription factor was found to have a significant role in the variations among the 533 rice accessions, dividing them into two main subspecies. A Tourist-like miniature inverted-repeat transposable element (MITE) was discovered in the promoter of TCP4 in japonica/geng accessions (TCP4M+), which was found to suppress the expression of TCP4 at the transcriptional level. The MITE-deleted haplotype (TCP4M-) was mainly found in indica/xian accessions. ChIP-qPCR and EMSA demonstrated the binding of TCP4 to promoters of grain reservoir genes such as SSIIa and Amy3D in vivo and in vitro, respectively. The introduction of the genomic sequence of TCP4M+ into different TCP4M- cultivars was found to affect the expression of TCP4 in the transgenic rice, resulting in decreased expression of its downstream target gene SSIIa, increased tiller number, and decreased seed length. This study revealed that a Tourist-like MITE contributes to subspecies divergence by regulating the expression of TCP4 in response to environmental pressure, thus influencing source-sink balance by regulating starch biosynthesis in rice.

Pan-genome: A promising resource for noncoding RNA discovery in plants

Plant genomes contain both protein-coding and noncoding sequences including transposable elements (TEs) and noncoding RNAs (ncRNAs). The ncRNAs are recognized as important elements that play fundamental roles in the structural organization and function of plant genomes. Despite various hypotheses, TEs are believed to be a major precursor of ncRNAs. Transposable elements are also prime factors that cause genomic variation among members of a species. Hence, TEs pose a major challenge in the discovery and analysis of ncRNAs. With the increase in the number of sequenced plant genomes, it is now accepted that a single reference genome is insufficient to represent the complete genomic diversity and contents of a species, and exploring the pan-genome of a species is critical. In this review, we summarize the recent progress in the field of plant pan-genomes. We also discuss TEs and their roles in ncRNA biogenesis and present our perspectives on the application of pan-genomes for the discovery of ncRNAs to fully explore and exploit their biological roles in plants.

Amplification of miniature inverted-repeat transposable elements and the associated impact on gene regulation and alternative splicing in mulberry (Morus notabilis)

Background

Miniature inverted-repeat transposable elements (MITEs) are common in eukaryotic genomes, and are important for genomic evolution.

Results

In the present study, the identification of MITEs in the mulberry genome revealed 286,122 MITE-related sequences, including 90,789 full-length elements. The amplification of mulberry MITEs and the influence of MITEs on the evolution of the mulberry genome were analyzed. The timing of MITE amplifications varied considerably among the various MITE families. Fifty-one MITE families have undergone a single round of amplification, while the other families developed from multiple amplifications. Most mulberry MITEs were inserted near genes and some could regulate gene expression through small RNAs. An analysis of transgenic plants indicated that MITE insertions can upregulate the expression of a target gene. Moreover, MITEs are frequently associated with alternative splicing events (exonizations).

Conclusion

The data presented herein provide insights into the generation of MITEs as well as their impact on gene regulation and evolution in mulberry.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0169-0) contains supplementary material, which is available to authorized users.

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.

Genome-wide identification and evolution of TC1/Mariner in the silkworm (Bombyx mori) genome

TC1/Mariner transposons belong to class II transposable elements (TEs) that use DNA-mediated "cut and paste" mechanism to transpose, and they have been identified in almost all organisms. Although silkworm (Bombyx mori) has a large amount of TC1/Mariner elements, the genome wide information of this superfamily in the silkworm is unknown. In this study, we have identified 2670 TC1/Mariner (Bmmar) elements in the silkworm genome. All the TEs were classified into 22 families by means of fgclust, a tool of repetitive sequence classification, seven of which was first reported in this study. Phylogenetic and structure analyses based on the catalytic domain (DDxD/E) of transposase sequences indicated that all members of TC1/Mariner were grouped into five subgroups: Mariner, Tc1, maT, DD40D and DD41D/E. Of these five subgroups, maT rather than Mariner possessed most members of TC1/Mariner (51.23%) in the silkworm genome. In particular, phylogenetic analysis and structure analysis revealed that Bmmar15 (DD40D) formed a new basal subgroup of TC1/Mariner element in insects, which was referred to as bmori. Furthermore, we concluded that DD40D appeared to intermediate between mariner and Tc1. Finally, we estimated the insertion time for each copy of TC1/Mariner in the silkworm and found that most of members were dramatically amplified during a period from 0 to 1 mya. Moreover, the detailed functional data analysis showed that Bmmar1, Bmmar6 and Bmmar9 had EST evidence and intact transposases. These implied that TC1/Mariner might have potential transpositional activity. In conclusion, this study provides some new insights into the landscape, origin and evolution of TC1/Mariner in the insect genomes.

Genome-wide characterization and evolution analysis of miniature inverted-repeat transposable elements (MITEs) in moso bamboo (Phyllostachys heterocycla)

Moso bamboo MITEs were genome-wide identified first time, and data shows that MITEs contribute to the genomic diversity and differentiation of bamboo. Miniature inverted-repeat transposable elements (MITEs) are widespread in animals and plants. There are a large number of transposable elements in moso bamboo (Phyllostachys heterocycla var. pubescens) genome, but the genome-wide information of moso bamboo MITEs is not known yet. Here we identified 362 MITE families with a total of 489,592 MITE-related sequences, accounting for 4.74 % of the moso bamboo genome. The 362 MITE families are clustered into six known and one unknown super-families. Our analysis indicated that moso bamboo MITEs preferred to reside in or near the genes that might be involved in regulation of host gene expression. Of the seven super-families, three might undergo major expansion event twice, respectively, during 8-11 million years ago (mya) ago and 22-28 mya ago; two might experience a long expansion period from 6 to 13 mya. Almost 1/3 small RNAs might be derived from the MITE sequences. Some MITE families generate small RNAs mainly from the terminals, while others predominantly from the central region. Given the high copy number of MITEs, many siRNAs and miRNAs derived from MITE sequences and the preferential insertion of MITE into gene regions, MITEs may contribute to the genomic diversity and differentiation of bamboo.

Widespread and evolutionary analysis of a MITE family Monkey King in Brassicaceae

BACKGROUND

Miniature inverted repeat transposable elements (MITEs) are important components of eukaryotic genomes, with hundreds of families and many copies, which may play important roles in gene regulation and genome evolution. However, few studies have investigated the molecular mechanisms involved. In our previous study, a Tourist-like MITE, Monkey King, was identified from the promoter region of a flowering time gene, BnFLC.A10, in Brassica napus. Based on this MITE, the characteristics and potential roles on gene regulation of the MITE family were analyzed in Brassicaceae.

RESULTS

The characteristics of the Tourist-like MITE family Monkey King in Brassicaceae, including its distribution, copies and insertion sites in the genomes of major Brassicaceae species were analyzed in this study. Monkey King was actively amplified in Brassica after divergence from Arabidopsis, which was indicated by the prompt increase in copy number and by phylogenetic analysis. The genomic variations caused by Monkey King insertions, both intra- and inter-species in Brassica, were traced by PCR amplification. Genomic sequence analysis showed that most complete Monkey King elements are located in gene-rich regions, less than 3kb from genes, in both the B. rapa and A. thaliana genomes. Sixty-seven Brassica expressed sequence tags carrying Monkey King fragments were also identified from the NCBI database. Bisulfite sequencing identified specific DNA methylation of cytosine residues in the Monkey King sequence. A fragment containing putative TATA-box motifs in the MITE sequence could bind with nuclear protein(s) extracted from leaves of B. napus plants. A Monkey King-related microRNA, bna-miR6031, was identified in the microRNA database. In transgenic A. thaliana, when the Monkey King element was inserted upstream of 35S promoter, the promoter activity was weakened.

CONCLUSION

Monkey King, a Brassicaceae Tourist-like MITE family, has amplified relatively recently and has induced intra- and inter-species genomic variations in Brassica. Monkey King elements are most abundant in the vicinity of genes and may have a substantial effect on genome-wide gene regulation in Brassicaceae. Monkey King insertions potentially regulate gene expression and genome evolution through epigenetic modification and new regulatory motif production.

Identification and evolution of the silkworm helitrons and their contribution to transcripts

In this study, we developed a structure-based approach to identify Helitrons in four lepidopterans and systematically analysed Helitrons in the silkworm genome. We found that the content of Helitrons varied greatly among genomes. The silkworm genome harboured 67,555 Helitron-related sequences that could be classified into 21 families and accounted for ≈ 4.23% of the genome. Thirteen of the families were new. Three families were putatively autonomous and included the replication initiator motif and helicase domain. The silkworm Helitrons were widely and randomly distributed in the genome. Most Helitron families radiated within the past 2 million years and experienced a single burst of expansion. These Helitron families captured 3724 gene fragments and contributed to at least 1.4% of the silkworm full-length cDNAs, suggesting important roles of Helitrons in the evolution of the silkworm genes. In addition, we found that some new Helitrons were generated by combinations of other Helitrons. Overall, the results presented in this study provided insights into the generation and evolution of Helitron transposons and their contribution to transcripts.

Miniature inverted-repeat transposable elements: discovery, distribution, and activity

Eukaryotic organisms have dynamic genomes, with transposable elements (TEs) as a major contributing factor. Although the large autonomous TEs can significantly shape genomic structures during evolution, genomes often harbor more miniature nonautonomous TEs that can infest genomic niches where large TEs are rare. In spite of their cut-and-paste transposition mechanisms that do not inherently favor copy number increase, miniature inverted-repeat transposable elements (MITEs) are abundant in eukaryotic genomes and exist in high copy numbers. Based on the large number of MITE families revealed in previous studies, accurate annotation of MITEs, particularly in newly sequenced genomes, will identify more genomes highly rich in these elements. Novel families identified from these analyses, together with the currently known families, will further deepen our understanding of the origins, transposase sources, and dramatic amplification of these elements.

Reconstructing the evolutionary history of transposable elements

The impact of transposable elements (TEs) on genome structure, plasticity, and evolution is still not well understood. The recent availability of complete genome sequences makes it possible to get new insights on the evolutionary dynamics of TEs from the phylogenetic analysis of their multiple copies in a wide range of species. However, this source of information is not always fully exploited. Here, we show how the history of transposition activity may be qualitatively and quantitatively reconstructed by considering the distribution of transposition events in the phylogenetic tree, along with the tree topology. Using statistical models developed to infer speciation and extinction rates in species phylogenies, we demonstrate that it is possible to estimate the past transposition rate of a TE family, as well as how this rate varies with time. This methodological framework may not only facilitate the interpretation of genomic data, but also serve as a basis to develop new theoretical and statistical models.

msechBari, a new MITE-like element in Drosophila sechellia related to the Bari transposon

A few occurrences of miniature inverted-repeat transposable elements (MITEs) have been reported in species of the genus Drosophila. Here, we describe msechBari, a MITE-like element in Drosophila sechellia. The element is short, approximately 90 bp in length, AT-rich and occurs in association with, or close to, genes, characteristics that are typical for MITEs. The identification was performed in silico using the sequenced genome of D. sechellia and confirmed in a laboratory strain. This short element is related to the Bari_DM transposon of Drosophila melanogaster, having terminal inverted repeats (TIRs) of a similar length and a high identity with the full-length Bari_DM element. The estimated recent origin of the element and the homogeneity observed between copies found in the genome suggests that msechBari could be active in D. sechellia.

Maize genetic diversity and association mapping using transposable element insertion polymorphisms

Transposable elements are the major component of the maize genome and presumably highly polymorphic yet they have not been used in population genetics and association analyses. Using the Transposon Display method, we isolated and converted into PCR-based markers 33 Miniature Inverted Repeat Transposable Elements (MITE) polymorphic insertions. These polymorphisms were genotyped on a population-based sample of 26 American landraces for a total of 322 plants. Genetic diversity was high and partitioned within and among landraces. The genetic groups identified using Bayesian clustering were in agreement with published data based on SNPs and SSRs, indicating that MITE polymorphisms reflect maize genetic history. To explore the contribution of MITEs to phenotypic variation, we undertook an association mapping approach in a panel of 367 maize lines phenotyped for 26 traits. We found a highly significant association between the marker ZmV1-9, on chromosome 1, and male flowering time. The variance explained by this association is consistent with a flowering delay of +123 degree-days. This MITE insertion is located at only 289 nucleotides from the 3' end of a Cytochrome P450-like gene, a region that was never identified in previous association mapping or QTL surveys. Interestingly, we found (i) a non-synonymous mutation located in the exon 2 of the gene in strong linkage disequilibrium with the MITE polymorphism, and (ii) a perfect sequence homology between the MITE sequence and a maize siRNA that could therefore potentially interfere with the expression of the Cytochrome P450-like gene. Those two observations among others offer exciting perspectives to validate functionally the role of this region on phenotypic variation.

Whole genome resequencing reveals natural target site preferences of transposable elements in Drosophila melanogaster

Transposable elements are mobile DNA sequences that integrate into host genomes using diverse mechanisms with varying degrees of target site specificity. While the target site preferences of some engineered transposable elements are well studied, the natural target preferences of most transposable elements are poorly characterized. Using population genomic resequencing data from 166 strains of Drosophila melanogaster, we identified over 8,000 new insertion sites not present in the reference genome sequence that we used to decode the natural target preferences of 22 families of transposable element in this species. We found that terminal inverted repeat transposon and long terminal repeat retrotransposon families present clade-specific target site duplications and target site sequence motifs. Additionally, we found that the sequence motifs at transposable element target sites are always palindromes that extend beyond the target site duplication. Our results demonstrate the utility of population genomics data for high-throughput inference of transposable element targeting preferences in the wild and establish general rules for terminal inverted repeat transposon and long terminal repeat retrotransposon target site selection in eukaryotic genomes.

The turbulent life of Sirevirus retrotransposons and the evolution of the maize genome: more than ten thousand elements tell the story

Sireviruses are one of the three genera of Copia long terminal repeat (LTR) retrotransposons, exclusive to and highly abundant in plants, and with a unique, among retrotransposons, genome structure. Yet, perhaps due to the few references to the Sirevirus origin of some families, compounded by the difficulty in correctly assigning retrotransposon families into genera, Sireviruses have hardly featured in recent research. As a result, analysis at this key level of classification and details of their colonization and impact on plant genomes are currently lacking. Recently, however, it became possible to accurately assign elements from diverse families to this genus in one step, based on highly conserved sequence motifs. Hence, Sirevirus dynamics in the relatively obese maize genome can now be comprehensively studied. Overall, we identified >10 600 intact and approximately 28 000 degenerate Sirevirus elements from a plethora of families, some brought into the genus for the first time. Sireviruses make up approximately 90% of the Copia population and it is the only genus that has successfully infiltrated the genome, possibly by experiencing intense amplification during the last 600 000 years, while being constantly recycled by host mechanisms. They accumulate in chromosome-distal gene-rich areas, where they insert in between gene islands, mainly in preferred zones within their own genomes. Sirevirus LTRs are heavily methylated, while there is evidence for a palindromic consensus target sequence. This work brings Sireviruses in the spotlight, elucidating their lifestyle and history, and suggesting their crucial role in the current genomic make-up of maize, and possibly other plant hosts.

Miniature inverted-repeat transposable elements (MITEs) have been accumulated through amplification bursts and play important roles in gene expression and species diversity in Oryza sativa

Miniature inverted–repeat transposable elements (MITEs) are predicted to play important roles on genome evolution. We developed a BLASTN-based approach for de novo identification of MITEs and systematically analyzed MITEs in rice genome. The genome of rice cultivar Nipponbare (Oryza sativa ssp. japonica) harbors 178,533 MITE-related sequences classified into 338 families. Pairwise nucleotide diversity and phylogenetic tree analysis indicated that individual MITE families were resulted from one or multiple rounds of amplification bursts. The timing of amplification burst varied considerably between different MITE families or subfamilies. MITEs are associated with 23,623 (58.2%) genes in rice genome. At least 7,887 MITEs are transcribed and more than 3,463 were transcribed with rice genes. The MITE sequences transcribed with rice coding genes form 1,130 pairs of potential natural sense/antisense transcripts. MITEs generate 23.5% (183,837 of 781,885) of all small RNAs identified from rice. Some MITE families generated small RNAs mainly from the terminals, while other families generated small RNAs predominantly from the central region. More than half (51.8%) of the MITE-derived small RNAs were generated exclusively by MITEs located away from genes. Genome-wide analysis showed that genes associated with MITEs have significantly lower expression than genes away from MITEs. Approximately 14.8% of loci with full-length MITEs have presence/absence polymorphism between rice cultivars 93-11 (O. sativa ssp. indica) and Nipponbare. Considering that different sets of genes may be regulated by MITE-derived small RNAs in different genotypes, MITEs provide considerable diversity for O. sativa.

BraSto, a Stowaway MITE from Brassica: recently active copies preferentially accumulate in the gene space

We characterized a Brassica miniature inverted repeat transposable element (MITE) from the Stowaway superfamily, designated BraSto (Bra ssica Sto waway). BraSto copy number was assessed using real-time quantitative PCR in the two diploid species B. rapa (genome A) and B. oleracea (genome C) and the corresponding allotetraploid species B. napus (genome AC). Phylogenetic relationships among a set of 131 BraSto copies were then analyzed. BraSto appears to have been only moderately amplified in the Brassica genome and was still active recently with marks of proliferation in both diploid Brassica species, which diverged 3.75 million years ago, but also in the allotetraploid species after reuniting of the two differentiated genomes. We characterized insertion sites for low-divergence BraSto copies among the gene space of the B. rapa genome using bioinformatics approaches. For BraSto copies localized nearby or within genes, we observed frequent associations of BraSto with putative promoters and regulatory regions of genes, but exclusion from coding regions. In addition, BraSto was significantly similar to several Brassica expressed sequence tags (ESTs), including stress-induced ESTs. We also demonstrated the enrichment of BraSto sequences in binding sites for transcription factors and other regulatory elements. Our results lead to the question of a role for BraSto in the regulation of gene expression: this putative role, if further confirmed experimentally, would help to obtain a new insight into the significance of MITEs in the functional plant genome.

Genome size and transposable element content as determined by high-throughput sequencing in maize and Zea luxurians

The genome of maize (Zea mays ssp. mays) consists mostly of transposable elements (TEs) and varies in size among lines. This variation extends to other species in the genus Zea: although maize and Zea luxurians diverged only ∼140,000 years ago, their genomes differ in size by ∼50%. We used paired-end Illumina sequencing to evaluate the potential contribution of TEs to the genome size difference between these two species. We aligned the reads both to a filtered gene set and to an exemplar database of unique repeats representing 1,514 TE families; ∼85% of reads mapped against TE repeats in both species. The relative contribution of TE families to the B73 genome was highly correlated with previous estimates, suggesting that reliable estimates of TE content can be obtained from short high-throughput sequencing reads, even at low coverage. Because we used paired-end reads, we could assess whether a TE was near a gene by determining if one paired read mapped to a TE and the second read mapped to a gene. Using this method, Class 2 DNA elements were found significantly more often in genic regions than Class 1 RNA elements, but Class 1 elements were found more often near other TEs. Overall, we found that both Class 1 and 2 TE families account for ∼70% of the genome size difference between B73 and luxurians. Interestingly, the relative abundance of TE families was conserved between species (r = 0.97), suggesting genome-wide control of TE content rather than family-specific effects.

Burst expansion, distribution and diversification of MITEs in the silkworm genome

Background

Miniature inverted-repeat transposable elements (MITEs) are widespread in plants and animals. Although silkworm (Bombyx mori) has a large amount of and a variety of transposable elements, the genome-wide information of the silkworm MITEs is unknown.

Results

We used structure-based and homology approaches to search for MITEs in the silkworm genome. We identified 17 MITE families with a total of 5785 members, accounting for ~0.4% of the genome. 7 of 17 MITE families are completely novel based on the nucleotide composition of target site duplication (TSD) and/or terminal inverted repeats (TIR). Silkworm MITEs were widely and nonrandom distributed in the genome. One family named BmMITE-2 might experience a recent burst expansion. Network and diversity analyses for each family revealed different diversification patterns of the silkworm MITEs, reflecting the signatures of genome-shocks that silkworm experienced. Most silkworm MITEs preferentially inserted into or near genes and BmMITE-11 that encodes a germline-restricted small RNA might silence its the closest genes in silkworm ovary through a small RNA pathway.

Conclusions

Silkworm harbors 17 MITE families. The silkworm MITEs preferred to reside in or near genes and one MITE might be involved in gene silence. Our results emphasize the exceptional role of MITEs in transcriptional regulation of genes and have general implications to understand interaction between MITEs and their host genome.