Identification of miniature inverted-repeat transposable elements (MITEs) and biogenesis of their siRNAs in the Solanaceae: new functional implications for MITEs

Genetic variability of aquaporin expression in maize: From eQTLs to a MITE insertion regulating PIP2;5 expression

Plant aquaporins are involved in numerous physiological processes, such as cellular homeostasis, tissue hydraulics, transpiration, and nutrient supply, and are key players of the response to environmental cues. While varying expression patterns of aquaporin genes have been described across organs, developmental stages, and stress conditions, the underlying regulation mechanisms remain elusive. Hence, this work aimed to shed light on the expression variability of 4 plasma membrane intrinsic protein (PIP) genes in maize (Zea mays) leaves, and its genetic causes, through expression quantitative trait locus (eQTL) mapping across a 252-hybrid diversity panel. Significant genetic variability in PIP transcript abundance was observed to different extents depending on the isoforms. The genome-wide association study mapped numerous eQTLs, both local and distant, thus emphasizing the existing natural diversity of PIP gene expression across the studied panel and the potential to reveal regulatory actors and mechanisms. One eQTL associated with PIP2;5 expression variation was characterized. Genomic sequence comparison and in vivo reporter assay attributed, at least partly, the local eQTL to a transposon-containing polymorphism in the PIP2;5 promoter. This work paves the way to the molecular understanding of PIP gene regulation and its possible integration into larger networks regulating physiological and stress adaptation processes.

A catalogue of recombination coldspots in interspecific tomato hybrids

Increasing natural resistance and resilience in plants is key for ensuring food security within a changing climate. Breeders improve these traits by crossing cultivars with their wild relatives and introgressing specific alleles through meiotic recombination. However, some genomic regions are devoid of recombination especially in crosses between divergent genomes, limiting the combinations of desirable alleles. Here, we used pooled-pollen sequencing to build a map of recombinant and non-recombinant regions between tomato and five wild relatives commonly used for introgressive tomato breeding. We detected hybrid-specific recombination coldspots that underscore the role of structural variations in modifying recombination patterns and maintaining genetic linkage in interspecific crosses. Crossover regions and coldspots show strong association with specific TE superfamilies exhibiting differentially accessible chromatin between somatic and meiotic cells. About two-thirds of the genome are conserved coldspots, located mostly in the pericentromeres and enriched with retrotransposons. The coldspots also harbor genes associated with agronomic traits and stress resistance, revealing undesired consequences of linkage drag and possible barriers to breeding. We presented examples of linkage drag that can potentially be resolved by pairing tomato with other wild species. Overall, this catalogue will help breeders better understand crossover localization and make informed decisions on generating new tomato varieties.

Conserved transcription factors NRZ1 and NRM1 regulate NLR receptor-mediated immunity

Plant innate immunity mediated by the nucleotide-binding leucine-rich repeat (NLR) class of immune receptors plays an important role in defense against various pathogens. Although key biochemical events involving NLR activation and signaling have been recently uncovered, we know very little about the transcriptional regulation of NLRs and their downstream signaling components. Here, we show that the Toll-Interleukin 1 receptor homology domain containing NLR (TNL) gene N (Necrosis), which confers resistance to Tobacco mosaic virus, is transcriptionally induced upon immune activation. We identified two conserved transcription factors, N required C3H zinc finger 1 (NRZ1) and N required MYB-like transcription factor 1 (NRM1), that activate N in an immune responsive manner. Genetic analyses indicated that NRZ1 and NRM1 positively regulate coiled-coil domain-containing NLR- and TNL-mediated immunity and function independently of the signaling component Enhanced Disease Susceptibility 1. Furthermore, NRZ1 functions upstream of NRM1 in cell death signaling, and their gene overexpression induces ectopic cell death and expression of NLR signaling components. Our findings uncovered a conserved transcriptional regulatory network that is central to NLR-mediated cell death and immune signaling in plants.

A transposable element-derived siRNAs involve DNA hypermethylation at the promoter of OsGSTZ4 for cadmium tolerance in rice

Environmental contaminants such as cadmium (Cd) pose high risks to crop production and human health. The genetic basis for regulation of Cd stress-responsive genes for plant adaptation to adverse environments remains poorly understood. In this study, we characterized a rice Zeta family glutathione-S-transferase (OsGSTZ4) gene for Cd detoxification. Heterologous expression of OsGSTZ4 in a yeast (Saccharomyces cerevisiae) conferred cellular Cd tolerance. Transgenic rice overexpressing OsGSTZ4 improved plant growth, attenuated Cd-induced toxicity, and accumulated more Cd in roots. OsGSTZ4 transcription was rapidly induced 3 h after Cd exposure and then declined to the basal level. This was followed by (days after Cd treatment) by CHH hypermethylation (by 41.2 %) at a MITE (Miniature Inverted-repeat Transposable Element) transposable element (TE) inserted in the 5'-untranscribed region (UTR) (-1,722 ∼ -1,392 bp) of OsGSTZ4. Meanwhile, three 24-nt siRNAs derived from the TE (-1,722 ∼ -1,471 bp) were detected and was also rapidly enriched under Cd stress. To validate the possibility that Cd-induced change in OsGSTZ4 expression correlates with the siRNAs-involved CHH methylation through an RdDM (RNA-directed DNA methylation) pathway, genetic analyses were performed. We found that the CHH methylation at the promoter and transcript level of OsGSTZ4 were compromised in the osdrm2 (loss of function for CHH methylation) and osrdr2i (defective in RNA-dependent RNA polymerase 2) but did not change in other types of methyltransferases such as osmet1, ossdg714 or osros1. Promoter deletion analyses confirmed that the siRNA target sequences were essential for the proper expression of OsGSTZ4. Our studies reveal an unusual feedback mechanism by which the Cd-induced rapid OsGSTZ4 expression for Cd tolerance would interplay with the late CHH hypermethylation to silence the TE through the 24-nt siRNAs- and Osdrm2-mediated RdDM pathway, and help understand the diversity of gene regulation via an epigenetic mechanism for rice adaptation to the environmental stress.

A multi-omic Nicotiana benthamiana resource for fundamental research and biotechnology

Nicotiana benthamiana is an invaluable model plant and biotechnology platform with a ~3 Gb allotetraploid genome. To further improve its usefulness and versatility, we have produced high-quality chromosome-level genome assemblies, coupled with transcriptome, epigenome, microRNA and transposable element datasets, for the ubiquitously used LAB strain and a related wild accession, QLD. In addition, single nucleotide polymorphism maps have been produced for a further two laboratory strains and four wild accessions. Despite the loss of five chromosomes from the ancestral tetraploid, expansion of intergenic regions, widespread segmental allopolyploidy, advanced diploidization and evidence of recent bursts of Copia pseudovirus (Copia) mobility not seen in other Nicotiana genomes, the two subgenomes of N. benthamiana show large regions of synteny across the Solanaceae. LAB and QLD have many genetic, metabolic and phenotypic differences, including disparate RNA interference responses, but are highly interfertile and amenable to genome editing and both transient and stable transformation. The LAB/QLD combination has the potential to be as useful as the Columbia-0/Landsberg errecta partnership, utilized from the early pioneering days of Arabidopsis genomics to today.

Miniature Inverted-Repeat Transposable Elements: Small DNA Transposons That Have Contributed to Plant MICRORNA Gene Evolution

Angiosperms form the largest phylum within the Plantae kingdom and show remarkable genetic variation due to the considerable difference in the nuclear genome size of each species. Transposable elements (TEs), mobile DNA sequences that can amplify and change their chromosome position, account for much of the difference in nuclear genome size between individual angiosperm species. Considering the dramatic consequences of TE movement, including the complete loss of gene function, it is unsurprising that the angiosperms have developed elegant molecular strategies to control TE amplification and movement. Specifically, the RNA-directed DNA methylation (RdDM) pathway, directed by the repeat-associated small-interfering RNA (rasiRNA) class of small regulatory RNA, forms the primary line of defense to control TE activity in the angiosperms. However, the miniature inverted-repeat transposable element (MITE) species of TE has at times avoided the repressive effects imposed by the rasiRNA-directed RdDM pathway. MITE proliferation in angiosperm nuclear genomes is due to their preference to transpose within gene-rich regions, a pattern of transposition that has enabled MITEs to gain further transcriptional activity. The sequence-based properties of a MITE results in the synthesis of a noncoding RNA (ncRNA), which, after transcription, folds to form a structure that closely resembles those of the precursor transcripts of the microRNA (miRNA) class of small regulatory RNA. This shared folding structure results in a MITE-derived miRNA being processed from the MITE-transcribed ncRNA, and post-maturation, the MITE-derived miRNA can be used by the core protein machinery of the miRNA pathway to regulate the expression of protein-coding genes that harbor homologous MITE insertions. Here, we outline the considerable contribution that the MITE species of TE have made to expanding the miRNA repertoire of the angiosperms.

Polymorphic inverted repeats near coding genes impact chromatin topology and phenotypic traits in Arabidopsis thaliana

Transposons are mobile elements that are commonly silenced to protect eukaryotic genome integrity. In plants, transposable element (TE)-derived inverted repeats (IRs) are commonly found near genes, where they affect host gene expression. However, the molecular mechanisms of such regulation are unclear in most cases. Expression of these IRs is associated with production of 24-nt small RNAs, methylation of the IRs, and drastic changes in local 3D chromatin organization. Notably, many of these IRs differ between Arabidopsis thaliana accessions, causing variation in short-range chromatin interactions and gene expression. CRISPR-Cas9-mediated disruption of two IRs leads to a switch in genome topology and gene expression with phenotypic consequences. Our data show that insertion of an IR near a gene provides an anchor point for chromatin interactions that profoundly impact the activity of neighboring loci. This turns IRs into powerful evolutionary agents that can contribute to rapid adaptation.

The role of transposon inverted repeats in balancing drought tolerance and yield-related traits in maize

The genomic basis underlying the selection for environmental adaptation and yield-related traits in maize remains poorly understood. Here we carried out genome-wide profiling of the small RNA (sRNA) transcriptome (sRNAome) and transcriptome landscapes of a global maize diversity panel under dry and wet conditions and uncover dozens of environment-specific regulatory hotspots. Transgenic and molecular studies of Drought-Related Environment-specific Super eQTL Hotspot on chromosome 8 (DRESH8) and ZmMYBR38, a target of DRESH8-derived small interfering RNAs, revealed a transposable element-mediated inverted repeats (TE-IR)-derived sRNA- and gene-regulatory network that balances plant drought tolerance with yield-related traits. A genome-wide scan revealed that TE-IRs associate with drought response and yield-related traits that were positively selected and expanded during maize domestication. These results indicate that TE-IR-mediated posttranscriptional regulation is a key molecular mechanism underlying the tradeoff between crop environmental adaptation and yield-related traits, providing potential genomic targets for the breeding of crops with greater stress tolerance but uncompromised yield.

A Conserved, Serine-Rich Protein Plays Opposite Roles in N-Mediated Immunity against TMV and N-Triggered Cell Death

Plant nucleotide-binding, leucine-rich, repeat-containing proteins (NLRs) play important roles in plant immunity. NLR expression and function are tightly regulated by multiple mechanisms. In this study, a conserved serine/arginine-rich protein (SR protein) was identified through the yeast one-hybrid screening of a tobacco cDNA library using DNA fragments from the N gene, an NLR that confers immunity to tobacco mosaic virus (TMV). This SR protein showed an interaction with a 3' genomic regulatory sequence (GRS) and has a potential role in regulating the alternative splicing of N. Thus, it was named SR regulator for N, abbreviated SR4N. Further study showed that SR4N plays a positive role in N-mediated cell death but a negative role in N protein accumulation. SR4N also promotes multiple virus replications in co-expression experiments, and this enhancement may not function through RNA silencing suppression, as it did not enhance 35S-GFP expression in co-infiltration experiments. Bioinformatic and molecular studies revealed that SR4N belongs to the SR2Z subtype of the SR protein family, which was conserved in both dicots and monocots, and its roles in repressing viral immunity and triggering cell death were also conserved. Our study revealed new roles for SR2Z family proteins in plant immunity against viruses.

Identification of Species-Specific MicroRNAs Provides Insights into Dynamic Evolution of MicroRNAs in Plants

MicroRNAs (miRNAs) are an important class of regulatory small RNAs that program gene expression, mainly at the post-transcriptional level. Although sporadic examples of species-specific miRNAs (termed SS-miRNAs) have been reported, a genome-scale study across a variety of distant species has not been assessed. Here, by comprehensively analyzing miRNAs in 81 plant species phylogenetically ranging from chlorophytes to angiosperms, we identified 8048 species-specific miRNAs from 5499 families, representing over 61.2% of the miRNA families in the examined species. An analysis of the conservation from different taxonomic levels supported the high turnover rate of SS-miRNAs, even over short evolutionary distances. A comparison of the intrinsic features between SS-miRNAs and NSS-miRNAs (non-species-specific miRNAs) indicated that the AU content of mature miRNAs was the most striking difference. Our data further illustrated a significant bias of the genomic coordinates towards SS-miRNAs lying close to or within genes. By analyzing the 125,267 putative target genes for the 7966 miRNAs, we found the preferentially regulated functions of SS-miRNAs related to diverse metabolic processes. Collectively, these findings underscore the dynamic evolution of miRNAs in the species-specific lineages.

The evolutionary history of small RNAs in Solanaceae

The Solanaceae or "nightshade" family is an economically important group with remarkable diversity. To gain a better understanding of how the unique biology of the Solanaceae relates to the family's small RNA (sRNA) genomic landscape, we downloaded over 255 publicly available sRNA data sets that comprise over 2.6 billion reads of sequence data. We applied a suite of computational tools to predict and annotate two major sRNA classes: (1) microRNAs (miRNAs), typically 20- to 22-nucleotide (nt) RNAs generated from a hairpin precursor and functioning in gene silencing and (2) short interfering RNAs (siRNAs), including 24-nt heterochromatic siRNAs typically functioning to repress repetitive regions of the genome via RNA-directed DNA methylation, as well as secondary phased siRNAs and trans-acting siRNAs generated via miRNA-directed cleavage of a polymerase II-derived RNA precursor. Our analyses described thousands of sRNA loci, including poorly understood clusters of 22-nt siRNAs that accumulate during viral infection. The birth, death, expansion, and contraction of these sRNA loci are dynamic evolutionary processes that characterize the Solanaceae family. These analyses indicate that individuals within the same genus share similar sRNA landscapes, whereas comparisons between distinct genera within the Solanaceae reveal relatively few commonalities.

Comprehensive analysis of both long and short read transcriptomes of a clonal and a seed-propagated model species reveal the prerequisites for transcriptional activation of autonomous and non-autonomous transposons in plants

Background

Transposable element (TE) transcription is a precursor to its mobilisation in host genomes. However, the characteristics of expressed TE loci, the identification of self-competent transposon loci contributing to new insertions, and the genomic conditions permitting their mobilisation remain largely unknown.

Results

Using Vitis vinifera embryogenic callus, we explored the impact of biotic stressors on transposon transcription through the exposure of the callus to live cultures of an endemic grapevine yeast, Hanseniaspora uvarum. We found that only 1.7–2.5% of total annotated TE loci were transcribed, of which 5–10% of these were full-length, and the expressed TE loci exhibited a strong location bias towards expressed genes. These trends in transposon transcription were also observed in RNA-seq data from Arabidopsis thaliana wild-type plants but not in epigenetically compromised Arabidopsis ddm1 mutants. Moreover, differentially expressed TE loci in the grapevine tended to share expression patterns with co-localised differentially expressed genes. Utilising nanopore cDNA sequencing, we found a strong correlation between the inclusion of intronic TEs in gene transcripts and the presence of premature termination codons in these transcripts. Finally, we identified low levels of full-length transcripts deriving from structurally intact TE loci in the grapevine model.

Conclusion

Our observations in two disparate plant models representing clonally and seed propagated plant species reveal a closely connected transcriptional relationship between TEs and co-localised genes, particularly when epigenetic silencing is not compromised. We found that the stress treatment alone was insufficient to induce large-scale full-length transcription from structurally intact TE loci, a necessity for non-autonomous and autonomous mobilisation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00271-5.

Genome size evolution in the beetle genus Diabrotica

Diabrocite corn rootworms are one of the most economically significant pests of maize in the United States and Europe and an emerging model for insect-plant interactions. Genome sizes of several species in the genus Diabrotica were estimated using flow cytometry along with that of Acalymma vittatum as an outgroup. Genome sizes ranged between 1.56 and 1.64 gigabase pairs and between 2.26 and 2.59 Gb, respectively, for the Diabrotica subgroups fucata and virgifera; the Acalymma vittatum genome size was around 1.65 Gb. This result indicated that a substantial increase in genome size occurred in the ancestor of the virgifera group. Further analysis of the fucata group and the virgifera group genome sequencing reads indicated that the genome size difference between the Diabrotica subgroups could be attributed to a higher content of transposable elements, mostly miniature inverted-transposable elements and gypsy-like long terminal repeat retroelements.

Recent amplification of microsatellite-associated miniature inverted-repeat transposable elements in the pineapple genome

BACKGROUND

Miniature inverted-repeat transposable elements (MITEs) are non-autonomous DNA transposable elements that play important roles in genome organization and evolution. Genome-wide identification and characterization of MITEs provide essential information for understanding genome structure and evolution.

RESULTS

We performed genome-wide identification and characterization of MITEs in the pineapple genome. The top two MITE families, accounting for 29.39% of the total MITEs and 3.86% of the pineapple genome, have insertion preference in (TA) n dinucleotide microsatellite regions. We therefore named these MITEs A. comosus microsatellite-associated MITEs (Ac-mMITEs). The two Ac-mMITE families, Ac-mMITE-1 and Ac-mMITE-2, shared sequence similarity in the terminal inverted repeat (TIR) regions, suggesting that these two Ac-mMITE families might be derived from a common or closely related autonomous elements. The Ac-mMITEs are frequently clustered via adjacent insertions. Among the 21,994 full-length Ac-mMITEs, 46.1% of them were present in clusters. By analyzing the Ac-mMITEs without (TA) n microsatellite flanking sequences, we found that Ac-mMITEs were likely derived from Mutator-like DNA transposon. Ac-MITEs showed highly polymorphic insertion sites between cultivated pineapples and their wild relatives. To better understand the evolutionary history of Ac-mMITEs, we filtered and performed comparative analysis on the two distinct groups of Ac-mMITEs, microsatellite-targeting MITEs (mt-MITEs) that are flanked by dinucleotide microsatellites on both sides and mutator-like MITEs (ml-MITEs) that contain 9/10 bp TSDs. Epigenetic analysis revealed a lower level of host-induced silencing on the mt-MITEs in comparison to the ml-MITEs, which partially explained the significantly higher abundance of mt-MITEs in pineapple genome. The mt-MITEs and ml-MITEs exhibited differential insertion preference to gene-related regions and RNA-seq analysis revealed their differential influences on expression regulation of nearby genes.

CONCLUSIONS

Ac-mMITEs are the most abundant MITEs in the pineapple genome and they were likely derived from Mutator-like DNA transposon. Preferential insertion in (TA) n microsatellite regions of Ac-mMITEs occurred recently and is likely the result of damage-limiting strategy adapted by Ac-mMITEs during co-evolution with their host. Insertion in (TA) n microsatellite regions might also have promoted the amplification of mt-MITEs. In addition, mt-MITEs showed no or negligible impact on nearby gene expression, which may help them escape genome control and lead to their amplification.

Distribution of MITE family Monkey King in rapeseed (Brassica napus L) and its influence on gene expression

Miniature inverted-repeat transposable elements (MITEs) are a group of class II transposable elements. The MITE Monkey King (MK) was first discovered upstream of BnFLC.A10. In this study, genome resequencing of four selected B. napus accessions, revealed more than 4000 distributed copies of MKs constituting ~2.4 Mb of the B. napus genomic sequence and caused 677 polymorphisms among the four accessions. MK -polymorphism-related markers across 128 natural and 58 synthetic accessions revealed more polymorphic MKs in natural than synthetic accessions. Ten MK -induced indels significantly affected the expression levels of the nearest gene based on RNAseq analysis, six of these effects were subsequently confirmed using qRT-PCR. Decreased expression pattern of MK -derived miRNA-bna-miR6031 was also observed under various stress treatments. Further research focused on the MITE families should promote not only our understanding of gene regulatory networks but also inform crop improvement efforts.

A Global Landscape of Miniature Inverted-Repeat Transposable Elements in the Carrot Genome

Miniature inverted-repeat transposable elements (MITEs) are the most abundant group of Class II mobile elements in plant genomes. Their presence in genic regions may alter gene structure and expression, providing a new source of functional diversity. Owing to their small size and lack of coding capacity, the identification of MITEs has been demanding. However, the increasing availability of reference genomes and bioinformatic tools provides better means for the genome-wide identification and analysis of MITEs and for the elucidation of their contribution to the evolution of plant genomes. We mined MITEs in the carrot reference genome DH1 using MITE-hunter and developed a curated carrot MITE repository comprising 428 families. Of the 31,025 MITE copies spanning 10.34 Mbp of the carrot genome, 54% were positioned in genic regions. Stowaways and Tourists were frequently present in the vicinity of genes, while Mutator-like MITEs were relatively more enriched in introns. hAT-like MITEs were relatively more frequently associated with transcribed regions, including untranslated regions (UTRs). Some carrot MITE families were shared with other Apiaceae species. We showed that hAT-like MITEs were involved in the formation of new splice variants of insertion-harboring genes. Thus, carrot MITEs contributed to the accretion of new diversity by altering transcripts and possibly affecting the regulation of many genes.

When junk DNA turns functional: transposon-derived non-coding RNAs in plants

Transposable elements (TEs) are major contributors to genome complexity in eukaryotes. TE mobilization may cause genome instability, although it can also drive genome diversity throughout evolution. TE transposition may influence the transcriptional activity of neighboring genes by modulating the epigenomic profile of the region or by altering the relative position of regulatory elements. Notably, TEs have emerged in the last few years as an important source of functional long and small non-coding RNAs. A plethora of small RNAs derived from TEs have been linked to the trans regulation of gene activity at the transcriptional and post-transcriptional levels. Furthermore, TE-derived long non-coding RNAs have been shown to modulate gene expression by interacting with protein partners, sequestering active small RNAs, and forming duplexes with DNA or other RNA molecules. In this review, we summarize our current knowledge of the functional and mechanistic paradigms of TE-derived long and small non-coding RNAs and discuss their role in plant development and evolution.

Secreted in Xylem Genes: Drivers of Host Adaptation in Fusarium oxysporum

Fusarium oxysporum (Fo) is a notorious pathogen that significantly contributes to yield losses in crops of high economic status. It is responsible for vascular wilt characterized by the browning of conductive tissue, wilting, and plant death. Individual strains of Fo are host specific (formae speciales), and approximately, 150 forms have been documented so far. The pathogen secretes small effector proteins in the xylem, termed as Secreted in Xylem (Six), that contribute to its virulence. Most of these proteins contain cysteine residues in even numbers. These proteins are encoded by SIX genes that reside on mobile pathogenicity chromosomes. So far, 14 proteins have been reported. However, formae speciales vary in SIX protein profile and their respective gene sequence. Thus, SIX genes have been employed as ideal markers for pathogen identification. Acquisition of SIX-encoding mobile pathogenicity chromosomes by non-pathogenic lines, through horizontal transfer, results in the evolution of new virulent lines. Recently, some SIX genes present on these pathogenicity chromosomes have been shown to be involved in defining variation in host specificity among formae speciales. Along these lines, the review entails the variability (formae speciales, races, and vegetative compatibility groups) and evolutionary relationships among members of F. oxysporum species complex (FOSC). It provides updated information on the diversity, structure, regulation, and (a)virulence functions of SIX genes. The improved understanding of roles of SIX in variability and virulence of Fo has significant implication in establishment of molecular framework and techniques for disease management. Finally, the review identifies the gaps in current knowledge and provides insights into potential research landscapes that can be explored to strengthen the understanding of functions of SIX genes.

Miniature inverted repeat transposable elements cis-regulate circular RNA expression and promote ethylene biosynthesis, reducing heat tolerance in Populus tomentosa

Transposable elements (TEs) and their reverse complementary sequence pairs (RCPs) are enriched around loci that produce circular RNAs (circRNAs) in plants. However, the function of these TE-RCP pairs in modulating circRNA expression remains elusive. Here, we identified 4609 circRNAs in poplar (Populus tomentosa) and showed that miniature inverted repeat transposable elements (MITEs)-RCPs were enriched in circRNA flanking regions. Moreover, we used expression quantitative trait nucleotide (eQTN) mapping to decipher the cis-regulatory role of MITEs. eQTN results showed that 14 single-nucleotide polymorphisms (SNPs) were significantly associated with Circ_0000408 and Circ_0003418 levels and the lead associated SNPs were located in MITE-RCP regions, indicating that MITE-RCP sequence variations affect exon circularization. Overexpression and knockdown analysis showed that Circ_0003418 positively modulated its parental gene, which encodes the RING-type E3 ligase XBAT32, and specifically increased the expression of the PtoXBAT32.5 transcript variant, which lacks the E3 ubiquitin ligase domain. Under heat stress, PtoXBAT32.5 expression was induced with up-regulation of Circ_0003418, resulting in increased production of ethylene and peroxidation of membrane lipids. Our findings thus reveal the cis-regulatory mechanism by which a MITE-RCP pair affects circRNA abundance in poplar and indicate that Circ_0003418 is a negative regulator of poplar heat tolerance via the ubiquitin-mediated protein modification pathway.

Screening of Helicoverpa armigera Mobilome Revealed Transposable Element Insertions in Insecticide Resistance Genes

Simple Summary

Transposable elements (TEs) are mobile DNA sequences that can copy themselves within a host genome. TE-mediated changes in regulation can lead to massive and rapid changes in expression, responses that are potentially highly adaptive when an organism is faced with a mortality agent in the environment, such as an insecticide. Helicoverpa armigera shows a hight number of reported cases of insecticide resistance worldwide, having evolved resistance against pyrethroids, organophosphates, carbamates, organochlorines, and recently to macrocyclic lactone spinosad and several Bacillus thuringiensis toxins. In the present study, we conducted a TE annotation using combined approaches, and the results revealed a total of 8521 TEs, representing 236,132 copies, covering 12.86% of the H. armigera genome. In addition, we underlined TE insertions in defensome genes and we successfully identified nine TE insertions belonging to the RTE, R2, CACTA, Mariner and hAT superfamilies.

Abstract

The cotton bollworm Helicoverpa armigera Hübner (Lepidoptera: Noctuidae) is an important pest of many crops that has developed resistance to almost all groups of insecticides used for its management. Insecticide resistance was often related to Transposable Element (TE) insertions near specific genes. In the present study, we deeply retrieve and annotate TEs in the H. armigera genome using the Pipeline to Retrieve and Annotate Transposable Elements, PiRATE. The results have shown that the TE library consists of 8521 sequences representing 236,132 TE copies, including 3133 Full-Length Copies (FLC), covering 12.86% of the H. armigera genome. These TEs were classified as 46.71% Class I and 53.29% Class II elements. Among Class I elements, Short and Long Interspersed Nuclear Elements (SINEs and LINEs) are the main families, representing 21.13% and 19.49% of the total TEs, respectively. Long Terminal Repeat (LTR) and Dictyostelium transposable element (DIRS) are less represented, with 5.55% and 0.53%, respectively. Class II elements are mainly Miniature Inverted Transposable Elements (MITEs) (49.11%), then Terminal Inverted Repeats (TIRs) (4.09%). Superfamilies of Class II elements, i.e., Transib, P elements, CACTA, Mutator, PIF-harbinger, Helitron, Maverick, Crypton and Merlin, were less represented, accounting for only 1.96% of total TEs. In addition, we highlighted TE insertions in insecticide resistance genes and we successfully identified nine TE insertions belonging to RTE, R2, CACTA, Mariner and hAT superfamilies. These insertions are hosted in genes encoding cytochrome P450 (CyP450), glutathione S-transferase (GST), and ATP-binding cassette (ABC) transporter belonging to the G and C1 family members. These insertions could therefore be involved in insecticide resistance observed in this pest.

GREEN STRIPE, encoding methylated TOMATO AGAMOUS-LIKE 1, regulates chloroplast development and Chl synthesis in fruit

Fruit development involves chloroplast development, carotenoid accumulation and fruit coloration. Although genetic regulation of fruit development has been extensively investigated, epigenetic regulation of fruit coloration remains largely unexplored. Here, we report a naturally occurring epigenetic regulation of TAGL1, and its impact on chloroplast development and fruit coloration. We used a genome-wide association study in combination with map-based cloning to identify the GREEN STRIPE (GS) locus, a methylated isoform of TAGL1 regulating diversified chloroplast development and carotenoid accumulation. Nonuniform pigmentation of fruit produced by GS was highly associated with methylation of the TAGL1 promoter, which is linked to a SNP at SL2.50ch07_63842838. High degrees of methylation of the TAGL1 promoter downregulated its expression, leading to green stripes. By contrast, low degrees of methylation led to light green stripes in gs. RNA-seq and ChIP collectively showed that the expression of genes involved with Chl synthesis and chloroplast development were significantly upregulated in green stripes relative to light green stripes. Quantitative PCR and dual luciferase assay confirmed that TAGL1 downregulates expression of SlMPEC, SlPsbQ, and SlCAB, and upregulates expression of PSY1 - genes which are associated with chloroplast development and carotenoid accumulation. Altogether, our findings regarding the GS locus demonstrate that naturally occurring methylation of TAGL1 has diverse effects on plastid development in fruit.

The Arabidopsis PHD-finger protein EDM2 has multiple roles in balancing NLR immune receptor gene expression

Plant NLR-type receptors serve as sensitive triggers of host immunity. Their expression has to be well-balanced, due to their interference with various cellular processes and dose-dependency of their defense-inducing activity. A genetic “arms race” with fast-evolving pathogenic microbes requires plants to constantly innovate their NLR repertoires. We previously showed that insertion of the COPIA-R7 retrotransposon into RPP7 co-opted the epigenetic transposon silencing signal H3K9me2 to a new function promoting expression of this Arabidopsis thaliana NLR gene. Recruitment of the histone binding protein EDM2 to COPIA-R7-associated H3K9me2 is required for optimal expression of RPP7. By profiling of genome-wide effects of EDM2, we now uncovered additional examples illustrating effects of transposons on NLR gene expression, strongly suggesting that these mobile elements can play critical roles in the rapid evolution of plant NLR genes by providing the “raw material” for gene expression mechanisms. We further found EDM2 to have a global role in NLR expression control. Besides serving as a positive regulator of RPP7 and a small number of other NLR genes, EDM2 acts as a suppressor of a multitude of additional NLR genes. We speculate that the dual functionality of EDM2 in NLR expression control arose from the need to compensate for fitness penalties caused by high expression of some NLR genes by suppression of others. Moreover, we are providing new insights into functional relationships of EDM2 with its interaction partner, the RNA binding protein EDM3/AIPP1, and its target gene IBM1, encoding an H3K9-demethylase.

We previously found the Arabidopsis thaliana PHD-finger protein EDM2 to serve as a chromatin-associated factor controlling expression of the NLR-type immune receptor gene RPP7. EDM2 binds to the transposon-silencing signal H3K9me2 and affects levels of this epigenetic mark at various loci. By genome-wide profiling of transcript- and H3K9me2-levels we now found EDM2 to have a broader role in controlling NLR gene expression. In order to mitigate fitness costs caused by its promoting effects on RPP7 expression and that of several other NLR genes, EDM2 seems to suppress expression of many additional members of this gene family. This observation is in line with multiple reports demonstrating the need for balanced expression of NLRs, which can substantially reduce overall plant fitness, but need to be present at certain minimal levels to confer sufficient immune protection. Our previous results demonstrated that the influence of EDM2 on RPP7 expression was co-opted to this immune receptor gene by the insertion of an EDM2-controlled transposon. Here, we are providing additional examples for transposon-associated effects on NLR gene expression, suggesting that these mobile elements play an important role for NLR genes by equipping members of this rapidly evolving gene family with regulatory mechanisms needed for balanced expression.

Genomic re-assessment of the transposable element landscape of the potato genome

We provide a comprehensive and reliable potato TE landscape, based on a wide variety of identification tools and integrative approaches, producing clear and ready-to-use outputs for the scientific community. Transposable elements (TEs) are DNA sequences with the ability to autoreplicate and move throughout the host genome. TEs are major drivers in stress response and genome evolution. Given their significance, the development of clear and efficient TE annotation pipelines has become essential for many species. The latest de novo TE discovery tools, along with available TEs from Repbase and sRNA-seq data, allowed us to perform a reliable potato TEs detection, classification and annotation through an open-source and freely available pipeline ( https://github.com/DiegoZavallo/TE_Discovery ). Using a variety of tools, approaches and rules, we were able to provide a clearly annotated of characterized TEs landscape. Additionally, we described the distribution of the different types of TEs across the genome, where LTRs and MITEs present a clear clustering pattern in pericentromeric and subtelomeric/telomeric regions respectively. Finally, we analyzed the insertion age and distribution of LTR retrotransposon families which display a distinct pattern between the two major superfamilies. While older Gypsy elements concentrated around heterochromatic regions, younger Copia elements located predominantly on euchromatic regions. Overall, we delivered not only a reliable, ready-to-use potato TE annotation files, but also all the necessary steps to perform de novo detection for other species.

Characterization of B-Genome Specific High Copy hAT MITE Families in Brassica nigra Genome

Miniature inverted-repeat transposable elements (MITEs) are non-autonomous class II transposons which have been shown to influence genome evolution. Brassica nigra L. (B-genome) is one of three Brassica diploids cultivated primarily as an oil crop, which harbors novel alleles important for breeding. Two new high copy hAT MITE families (BniHAT-1 and BniHAT-2) from the B-genome were characterized and their prevalence assessed in the genomes of the related diploids, rapa L. (A) and Brassica oleracea L. (C). Both novel MITE families were present at high copy numbers in the B-genome with 434 and 331 copies of BniHAT-1 and BniHAT-2, respectively. Yet less than 20 elements were identified in the genome assemblies of the A, and C -genomes, supporting B-genome specific proliferation of these MITE families. Although apparently randomly distributed across the genome, 68 and 70% of the B-genome MITEs were present within 2 kb flanking regions of annotated genes suggesting they might influence gene expression and/or function. In addition, MITE derived microRNAs and transcription factor binding sites suggested a putative role in gene regulation. Age of insertion analysis revealed that the major proliferation of these elements occurred during 2-3 million years ago. Additionally, site-specific polymorphism analyses showed that 44% MITEs were undergoing active amplification into the B-genome. Overall, this study provides a comprehensive analysis of two high copy MITE families, which were specifically amplified in the B-genome, suggesting a potential role in shaping the Brassica B-genome.

Evolution and domestication of Tc1/mariner transposons in the genome of African coelacanth (Latimeria chalumnae)

Here, we comprehensively analysed the abundance, diversity, and activity of Tc1/mariner transposons in African coelacanth (Latimeria chalumnae). Fifteen Tc1/mariner autonomous transposons were identified and grouped into six clades: DD34E/Tc1, DD34D/mariner, DD35D/Fot, DD31D/pogo, DD30-31D/pogo-like, and DD32-36D/Tigger, belonging to three known families: DD34E/Tc1, DD34D/mariner, and DD×D/pogo (DD35D/Fot, DD31D/pogo, DD30-31D/pogo-like, and DD32-36D/Tigger). Thirty-one miniature inverted-repeat transposable element (MITE) transposons of Tc1/mariner were also identified, and 20 of them display similarity to the identified autonomous transposons. The structural organization of these full Tc1/mariner elements includes a transposase gene flanked by terminal inverted repeats (TIRs) with TA dinucleotides. The transposases contain N-terminal DNA binding domain and a C-terminal catalytic domain characterized by the presence of a conservative D(Asp)DE(Glu)/D triad that is essential for transposase activity. The Tc1/mariner superfamily in coelacanth exhibited very low genome coverage (0.3%), but it experienced an extraordinary difference of proliferation dynamics among the six clades identified; moreover, most of them exhibited a very recent and current proliferation, suggesting that some copies of these transposons are putatively active. Additionally, at least four functional genes derived from Tc1/mariner transposons were found. We provide an up-to-date overview of Tc1/mariner in coelacanth, which may be helpful in determining genome and gene evolution in this living fossil.

Structure and diversity of 13S globulin zero-repeat subunit, the trypsin-resistant storage protein of common buckwheat (Fagopyrum esculentum M.) seeds

The zero-repeat subunit of 13S globulin, which lacks tandem repeat inserts, is trypsin-resistant and suggested to show higher allergenicity than the other subunits in common buckwheat (Fagopyrum esculentum Moench). To evaluate allelic variations and find novel alleles, the diversity of the zero-repeat genes was examined for two Japanese elite cultivars and 15 Pakistani landraces. The results demonstrated that two new alleles GlbNA1 and GlbNC1, plus three additional new alleles GlbNA2, GlbNA3, and GlbND, were identified besides the already-known GlbNA, GlbNB, and GlbNC alleles. In the Pakistani landraces, GlbNA was the most dominant allele (0.60-0.88 of allele frequency) in all except one landrace, where GlbNB was the most dominant allele (0.50 of allele frequency). Similar to GlbNC, the alleles GlbNA2 and GlbNA3 had extra ~200 bp MITE-like sequences around the stop codon. Secondary structure predictions of a sense strand demonstrated that the extra ~200 bp sequences of GlbNC, GlbNA2, and GlbNA3 can form rigid hairpin structures with free energies of -78.95, -67.06, and -29.90 kcal/mol, respectively. These structures may affect proper transcription and/or translation. In the GlbNC homozygous line, no transcript of a zero-repeat gene was detected, suggesting the material would be useful for developing hypoallergenic buckwheat.

Stowaway miniature inverted repeat transposable elements are important agents driving recent genomic diversity in wild and cultivated carrot

BACKGROUND

Miniature inverted repeat transposable elements (MITEs) are small non-autonomous DNA transposons that are ubiquitous in plant genomes, and are mobilised by their autonomous relatives. Stowaway MITEs are derived from and mobilised by elements from the mariner superfamily. Those elements constitute a significant portion of the carrot genome; however the variation caused by Daucus carota Stowaway MITEs (DcStos), their association with genes and their putative impact on genome evolution has not been comprehensively analysed.

RESULTS

Fourteen families of Stowaway elements DcStos occupy about 0.5% of the carrot genome. We systematically analysed 31 genomes of wild and cultivated Daucus carota, yielding 18.5 thousand copies of these elements, showing remarkable insertion site polymorphism. DcSto element demography differed based on the origin of the host populations, and corresponded with the four major groups of D. carota, wild European, wild Asian, eastern cultivated and western cultivated. The DcStos elements were associated with genes, and most frequently occurred in 5' and 3' untranslated regions (UTRs). Individual families differed in their propensity to reside in particular segments of genes. Most importantly, DcSto copies in the 2 kb regions up- and downstream of genes were more frequently associated with open reading frames encoding transcription factors, suggesting their possible functional impact. More than 1.5% of all DcSto insertion sites in different host genomes contained different copies in exactly the same position, indicating the existence of insertional hotspots. The DcSto7b family was much more polymorphic than the other families in cultivated carrot. A line of evidence pointed at its activity in the course of carrot domestication, and identified Dcmar1 as an active carrot mariner element and a possible source of the transposition machinery for DcSto7b.

CONCLUSION

Stowaway MITEs have made a substantial contribution to the structural and functional variability of the carrot genome.

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.

Genome-wide association analysis of the lipid and fatty acid metabolism regulatory network in the mesocarp of oil palm (Elaeis guineensis Jacq.) based on small noncoding RNA sequencing

Oil palm (Elaeis guineensis Jacq.) is the highest oil-yielding crop in the plant kingdom and accumulates 90% of palm oil in the mesocarp. However, the regulatory mechanisms of lipid and fatty acid (FA) metabolism in oil palm are just beginning to be understood, and more studies are needed, especially in the understanding of small noncoding RNA (ncRNA) and mRNA. Based on the deep sequencing of small noncoding RNAs and the degradome in five developmental mesocarp stages, 452 microRNAs (miRNAs), including 170 conserved known-miRNAs (kn-miRNAs) and 282 novel-miRNA (nov-miRNAs), were identified. After predicting the targets of those miRNAs to 37 FA synthesis-related genes, we found that 22 kn-miRNAs and 14 nov-miRNAs might be involved in FA metabolism pathways. Among them, eg-miR156c, eg-miR397, eg-miR444b and nov-miR129 regulated FA synthesis in plastids and the transport of FA-ACP from plastids to the endoplasmic reticulum by targeting acetyl-CoA carboxylase 1 (ACC1), long-chain acyl-CoA synthetase 9 (LACS9), LACS4 and enoyl-ACP reductase (ENR), respectively. Nov-miR138 and nov-miR59 targeted glycerol-3-phosphate acyltransferase (GPAT), and nov-miR274 targeted phosphatidate phosphatase 1 (PAP1). Both target genes are involved in triacylglycerol synthesis in the endoplasmic reticulum. Eg-miR156e and eg-miR156j played pivotal roles by targeting β-ketoacyl-CoA synthase 12 (KCS12), and nov-miR201 targets very-long-chain enoyl-CoA reductase (ECR). Several miRNAs were also predicted to indirectly regulate FA synthesis and lipid metabolism through the squamosa promoter-binding protein-like gene (SPL), NAC and MYB transcription factors. As a whole, indications of a complex and extensive miRNA-mRNA regulatory network associated with FA metabolism in the mesocarp of the oil palm is presented. The results help to broaden the knowledge of potential mechanisms that might be regulated by miRNAs through modulation of the expression of FA-related target gene metabolism in the oil palm.

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.

Parental experience modifies the Mimulus methylome

BACKGROUND

Transgenerational plasticity occurs when the environmental experience of an organism modifies the growth and development of its progeny. Leaf damage in Mimulus guttatus exhibits transgenerational plasticity mediated through differential expression of hundreds of genes. The epigenetic mechanisms that facilitate this response have yet to be described.

RESULTS

We performed whole genome bisulfite sequencing in the progeny of genetically identical damaged and control plants and developed a pipeline to compare differences in the mean and variance of methylation between treatment groups. We find that parental damage increases the variability of CG and CHG methylation among progeny, but does not alter the overall mean methylation. Instead it has positive effects in some regions and negative in others. We find 3,396 CHH, 203 CG, and 54 CHG Differentially Methylated Regions (DMRs) ranging from tens to thousands of base pairs scattered across the genome. CHG and CHH DMRs tended to overlap with transposable elements. CG DMRs tended to overlap with gene coding regions, many of which were previously found to be differentially expressed.

CONCLUSIONS

Genome-wide increases in methylome variation suggest that parental conditions can increase epigenetic diversity in response to stress. Additionally, the potential association between CG DMRs and differentially expressed genes supports the hypothesis that differential methylation is a mechanistic component of transgenerational plasticity in M. guttatus.

Characterization and chromosomal mapping of the DgmarMITE transposon in populations of Dichotomius (Luederwaldtinia) sericeus species complex (Coleoptera: Scarabaeidae)

Transposable elements are dispersed repetitive DNA sequences that can move within the genome and are related to genome and chromosome evolution, adaptation, and speciation. The aim of this study was to characterize and determine the chromosomal location and accumulation of a Mariner-like element in populations of four phylogenetically related species of the Dichotomius (Luederwaldtinia) sericeus complex. Mapping of the isolated element was performed by fluorescent in situ hybridization in different populations of analyzed species. Characterization of the isolated element revealed a degenerated transposon, named DgmarMITE. This transposon is 496-bp-long, AT rich (57%), and contains 24 bp terminal inverted repeats. In situ mapping revealed presence of this element only in two out of four species analyzed. DgmarMITE sites were located in heterochromatic and euchromatic regions and varied in location and number on the karyotypes of Dichotomius (L.) gilletti and D. (L.) guaribensis across different populations. These results demonstrate differential accumulation of the DgmarMITE in genomes of these species, which is probably due to the occurrence of ectopic recombination and cross-mobilization of the element mediated by the transposase of closely related or unrelated transposable elements.

Transcript-level expression control of plant NLR genes

Plant NLR genes encode sensitive immune receptors that can mediate the specific recognition of pathogen avirulence effectors and activate a strong defence response, termed effector-triggered immunity. The expression of NLRs requires strict regulation, as their ability to trigger immunity is dependent on their dose, and overexpression of NLRs results in autoimmunity and massive fitness costs. An elaborate interplay of different mechanisms controlling NLR transcript levels allows plants to maximize their defence capacity, whilst limiting negative impact on their fitness. Global suppression of NLR transcripts may be a prerequisite for the fast evolution of new NLR variants and the expansion of this gene family. Here, we summarize recent progress made towards a comprehensive understanding of NLR transcript-level expression control. Multiple mechanistic steps, including transcription as well as co-/post-transcriptional processing and transcript turn-over, contribute to balanced base levels of NLR transcripts and allow for dynamic adjustments to defence situations.

New insights into the genome repetitive fraction of the Antarctic bivalve Adamussium colbecki

Repetitive DNA represents the major component of the genome in both plant and animal species. It includes transposable elements (TEs), which are dispersed throughout the genome, and satellite DNAs (satDNAs), which are tandemly organized in long arrays. The study of the structure and organization of repetitive DNA contributes to our understanding of genome architecture and the mechanisms leading to its evolution. Molluscs represent one of the largest groups of invertebrates and include organisms with a wide variety of morphologies and lifestyles. To increase our knowledge of bivalves at the genome level, we analysed the Antarctic scallop Adamussium colbecki. The screening of the genomic library evidenced the presence of two novel satDNA elements and the CvA transposon. The interspecific investigation performed in this study demonstrated that one of the two satDNAs isolated in A. colbecki is widespread in polar molluscan species, indicating a possible link between repetitive DNA and abiotic factors. Moreover, the transcriptional activity of CvA and its presence in long-diverged bivalves suggests a possible role for this ancient element in shaping the genome architecture of this clade.

A role for small RNA in regulating innate immunity during plant growth

Plant genomes encode large numbers of nucleotide-binding (NB) leucine-rich repeat (LRR) immune receptors (NLR) that mediate effector triggered immunity (ETI) and play key roles in protecting crops from diseases caused by devastating pathogens. Fitness costs are associated with plant NLR genes and regulation of NLR genes by micro(mi)RNAs and phased small interfering RNAs (phasiRNA) is proposed as a mechanism for reducing these fitness costs. However, whether NLR expression and NLR-mediated immunity are regulated during plant growth is unclear. We conducted genome-wide transcriptome analysis and showed that NLR expression gradually increased while expression of their regulatory small RNAs (sRNA) gradually decreased as plants matured, indicating that sRNAs could play a role in regulating NLR expression during plant growth. We further tested the role of miRNA in the growth regulation of NLRs using the tobacco mosaic virus (TMV) resistance gene N, which was targeted by miR6019 and miR6020. We showed that N-mediated resistance to TMV effectively restricted this virus to the infected leaves of 6-week old plants, whereas TMV infection was lethal in 1- and 3-week old seedlings due to virus-induced systemic necrosis. We further found that N transcript levels gradually increased while miR6019 levels gradually decreased during seedling maturation that occurs in the weeks after germination. Analyses of reporter genes in transgenic plants showed that growth regulation of N expression was post-transcriptionally mediated by MIR6019/6020 whereas MIR6019/6020 was regulated at the transcriptional level during plant growth. TMV infection of MIR6019/6020 transgenic plants indicated a key role for miR6019-triggered phasiRNA production for regulation of N-mediated immunity. Together our results demonstrate a mechanistic role for miRNAs in regulating innate immunity during plant growth.

Genome of wild olive and the evolution of oil biosynthesis

Here we present the genome sequence and annotation of the wild olive tree (Olea europaea var. sylvestris), called oleaster, which is considered an ancestor of cultivated olive trees. More than 50,000 protein-coding genes were predicted, a majority of which could be anchored to 23 pseudochromosomes obtained through a newly constructed genetic map. The oleaster genome contains signatures of two Oleaceae lineage-specific paleopolyploidy events, dated at ∼28 and ∼59 Mya. These events contributed to the expansion and neofunctionalization of genes and gene families that play important roles in oil biosynthesis. The functional divergence of oil biosynthesis pathway genes, such as FAD2, SACPD, EAR, and ACPTE, following duplication, has been responsible for the differential accumulation of oleic and linoleic acids produced in olive compared with sesame, a closely related oil crop. Duplicated oleaster FAD2 genes are regulated by an siRNA derived from a transposable element-rich region, leading to suppressed levels of FAD2 gene expression. Additionally, neofunctionalization of members of the SACPD gene family has led to increased expression of SACPD2, 3, 5, and 7, consequently resulting in an increased desaturation of steric acid. Taken together, decreased FAD2 expression and increased SACPD expression likely explain the accumulation of exceptionally high levels of oleic acid in olive. The oleaster genome thus provides important insights into the evolution of oil biosynthesis and will be a valuable resource for oil crop genomics.

Wild tobacco genomes reveal the evolution of nicotine biosynthesis

Nicotine, the signature alkaloid of Nicotiana species responsible for the addictive properties of human tobacco smoking, functions as a defensive neurotoxin against attacking herbivores. However, the evolution of the genetic features that contributed to the assembly of the nicotine biosynthetic pathway remains unknown. We sequenced and assembled genomes of two wild tobaccos, Nicotiana attenuata (2.5 Gb) and Nicotiana obtusifolia (1.5 Gb), two ecological models for investigating adaptive traits in nature. We show that after the Solanaceae whole-genome triplication event, a repertoire of rapidly expanding transposable elements (TEs) bloated these Nicotiana genomes, promoted expression divergences among duplicated genes, and contributed to the evolution of herbivory-induced signaling and defenses, including nicotine biosynthesis. The biosynthetic machinery that allows for nicotine synthesis in the roots evolved from the stepwise duplications of two ancient primary metabolic pathways: the polyamine and nicotinamide adenine dinucleotide (NAD) pathways. In contrast to the duplication of the polyamine pathway that is shared among several solanaceous genera producing polyamine-derived tropane alkaloids, we found that lineage-specific duplications within the NAD pathway and the evolution of root-specific expression of the duplicated Solanaceae-specific ethylene response factor that activates the expression of all nicotine biosynthetic genes resulted in the innovative and efficient production of nicotine in the genus Nicotiana Transcription factor binding motifs derived from TEs may have contributed to the coexpression of nicotine biosynthetic pathway genes and coordinated the metabolic flux. Together, these results provide evidence that TEs and gene duplications facilitated the emergence of a key metabolic innovation relevant to plant fitness.

Genome-Wide Comparative Analysis of Miniature Inverted Repeat Transposable Elements in 19 Arabidopsis thaliana Ecotype Accessions

Miniature inverted repeat transposable elements (MITEs) are prevalent in eukaryotic genomes. They are known to critically influence the process of genome evolution and play a role in gene regulation. As the first study concentrated in the transposition activities of MITEs among different ecotype accessions within a species, we conducted a genome-wide comparative analysis by characterizing and comparing MITEs in 19 Arabidopsis thaliana accessions. A total of 343485 MITE putative sequences, including canonical, diverse and partial ones, were delineated from all 19 accessions. Within the entire population of MITEs sequences, 80.7% of them were previously unclassified MITEs, demonstrating a different genomic distribution and functionality compared to the classified MITEs. The interactions between MITEs and homologous genes across 19 accessions provided a fine source for analyzing MITE transposition activities and their impacts on genome evolution. Moreover, a significant proportion of MITEs were found located in the last exon of genes besides the ordinary intron locality, thus potentially modifying the end of genes. Finally, analysis of the impact of MITEs on gene expression suggests that migrations of MITEs have no detectable effect on the expression level for host genes across accessions.

iMITEdb: the genome-wide landscape of miniature inverted-repeat transposable elements in insects

Miniature inverted-repeat transposable elements (MITEs) have attracted much attention due to their widespread occurrence and high copy numbers in eukaryotic genomes. However, the systematic knowledge about MITEs in insects and other animals is still lacking. In this study, we identified 6012 MITE families from 98 insect species genomes. Comparison of these MITEs with known MITEs in the NCBI non-redundant database and Repbase showed that 5701(∼95%) of 6012 MITE families are novel. The abundance of MITEs varies drastically among different insect species, and significantly correlates with genome size. In general, larger genomes contain more MITEs than small genomes. Furthermore, all identified MITEs were included in a newly constructed database (iMITEdb) (http://gene.cqu.edu.cn/iMITEdb/), which has functions such as browse, search, BLAST and download. Overall, our results not only provide insight on insect MITEs but will also improve assembly and annotation of insect genomes. More importantly, the results presented in this study will promote studies of MITEs function, evolution and application in insects.

Database URL: http://gene.cqu.edu.cn/iMITEdb/

Clustering and evolutionary analysis of small RNAs identify regulatory siRNA clusters induced under drought stress in rice

Motivation

Drought tolerance is an important trait related to growth and yield in crop. Until now, drought related research has focused on coding genes. However, non-coding RNAs also respond significantly to environmental stimuli such as drought stress. Unfortunately, characterizing the role of siRNAs under drought stress is difficult since a large number of heterogenous siRNA species are expressed under drought stress and non-coding RNAs have very weak evolutionary conservation. Thus, to characterize the role of siRNAs, we need a well designed biological and bioinformatics strategy. In this paper, to characterize the function of siRNAs we developed and used a bioinformatics pipeline that includes a genomic-location based clustering technique and an evolutionary conservation tool.

Results

By comparing the wild type Nipponbare and two drought resistant rice varities, we found that 21 nt and 24 nt siRNAs are significantly expressed in the three rice plants but at different time points under a short-term (0, 1, and 6 hrs) drought treatment. siRNAs were up-regulated in the wild type at an early stage while the up-regulation was delayed in the two drought tolerant plants. Genes targeted by up-regulated siRNAs were related to oxidation reduction and proteolysis, which are well known to be associated with water deficit phenotypes. More interestingly, we found that siRNAs were located in intronic regions as clusters and were of high evolutionary conservation among monocot grass plants. In summary, we show that siRNAs are important respondents to drought stress and regulate genes related to the drought tolerance in water deficit conditions.

Electronic supplementary material

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

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.

Beta-amylase gene variability in introgressive wheat lines

Variability of the beta-amylase gene in bread wheat, artificial amphidiploids, and derived introgression wheat lines was analyzed. Variation in homeologous beta-amylase sequences caused by the presence of MITE (Miniature Inverted-Repeat Transposable Element) and its footprint has been identified in bread wheat. The previously unknown location of MITE in Triticum urartu and T. aestivum L. beta-amylase gene has been found. These species have a MITE sequence in the third intron of beta-amylase, as opposed to Aegilops comosa and a number of other Triticeae species, which have it in the fourth intron. These two MITEs from Ae. comosa and T. aestivum were shown to have low identity scores. Miosa, an artificial amphidiploid, which has the M genome from Ae. comosa was shown to lose the MITE sequences. This loss might be caused by genomic shock due to allopolyploidization.

The characteristics and functions of a miniature inverted-repeat transposable element TaMITE81 in the 5' UTR of TaCHS7BL from Triticum aestivum

Miniature inverted-repeat transposable elements (MITEs) are truncated derivatives of autonomous DNA transposons, and are dispersed abundantly in eukaryotic and prokaryotic genomes. In this article, a MITE, TaMITE81, was isolated from the 5' untranslated region (UTR) of TaCHS7BL, chalcone synthase (CHS) catalyzing the first committed step of anthocyanin biosynthesis, in the wheat cultivar 'Opata' with white grain. TaMITE81 was only 81 nucleotides, including a terminal inverted repeat with 39 nucleotides and was flanked by two nucleotides, "TA", target site duplications that were typical features of stowaway-like MITEs. Compared with the wheat cultivar 'Gy115' with purple grain, which is without the insertion, the expression of TaCHS7BL was lower in several organs of 'Opata'. The insertion of TaMITE81 into the 5' UTR of the GUS gene also reduced the transient expression of GUS on the coleoptiles of 'Opata', which means the insertion of TaMITE81 was the reason for the low expression of TaCHS7BL in 'Opata'. But the genotype of TaCHS7BL was not linked to phenotype of grain color in the RILs derived from a cross 'Gy115' and 'Opata'. The TaMITE81 density of the hexaploid variety of T. aestivum was more than 10 times that of diploid relatives, which implies that polyploidization caused the amplification of TaMITE81 homologous sequences. Further research should be conducted on decoding the relationship between TaCHS7BL and other traits relative to anthocyanin biosynthesis in wheat, and discovering the mechanism of TaMITE81 transposon action.

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.

A Signaling Cascade from miR444 to RDR1 in Rice Antiviral RNA Silencing Pathway

Plant RNA-DEPENDENT RNA POLYMERASE1 (RDR1) is a key component of the antiviral RNA-silencing pathway, contributing to the biogenesis of virus-derived small interfering RNAs. This enzyme also is responsible for producing virus-activated endogenous small interfering RNAs to stimulate the broad-spectrum antiviral activity through silencing host genes. The expression of RDR1 orthologs in various plants is usually induced by virus infection. However, the molecular mechanisms of activation of RDR1 expression in response to virus infection remain unknown. Here, we show that a monocot-specific microRNA, miR444, is a key factor in relaying the antiviral signaling from virus infection to OsRDR1 expression. The expression of miR444 is enhanced by infection with Rice stripe virus (RSV), and overexpression of miR444 improves rice (Oryza sativa) resistance against RSV infection accompanied by the up-regulation of OsRDR1 expression. We further show that three miR444 targets, the MIKC(C)-type MADS box proteins OsMADS23, OsMADS27a, and OsMADS57, form homodimers and heterodimers between them to repress the expression of OsRDR1 by directly binding to the CArG motifs of its promoter. Consequently, an increased level of miR444 diminishes the repressive roles of OsMADS23, OsMADS27a, and OsMADS57 on OsRDR1 transcription, thus activating the OsRDR1-dependent antiviral RNA-silencing pathway. We also show that overexpression of miR444-resistant OsMADS57 reduced OsRDR1 expression and rice resistance against RSV infection, and knockout of OsRDR1 reduced rice resistance against RSV infection. In conclusion, our results reveal a molecular cascade in the rice antiviral pathway in which miR444 and its MADS box targets directly control OsRDR1 transcription.

Evidence for the role of transposons in the recruitment of cis-regulatory motifs during the evolution of C4 photosynthesis

Background

C₄ photosynthesis evolved from C₃ photosynthesis and has higher light, water, and nitrogen use efficiencies. Several C₄ photosynthesis genes show cell-specific expression patterns, which are required for these high resource-use efficiencies. However, the mechanisms underlying the evolution of cis-regulatory elements that control these cell-specific expression patterns remain elusive.

Results

In the present study, we tested the hypothesis that the cis-regulatory motifs related to C₄ photosynthesis genes were recruited from non-photosynthetic genes and further examined potential mechanisms facilitating this recruitment. We examined 65 predicted bundle sheath cell-specific motifs, 17 experimentally validated cell-specific cis-regulatory elements, and 1,034 motifs derived from gene regulatory networks. Approximately 7, 5, and 1,000 of these three categories of motifs, respectively, were apparently recruited during the evolution of C₄ photosynthesis. In addition, we checked 1) the distance between the acceptors and the donors of potentially recruited motifs in a chromosome, and 2) whether the potentially recruited motifs reside within the overlapping region of transposable elements and the promoter of donor genes. The results showed that 7, 4, and 658 of the potentially recruited motifs might have moved via the transposable elements. Furthermore, the potentially recruited motifs showed higher binding affinity to transcription factors compared to randomly generated sequences of the same length as the motifs.

Conclusions

This study provides molecular evidence supporting the hypothesis that transposon-driven recruitment of pre-existing cis-regulatory elements from non-photosynthetic genes into photosynthetic genes plays an important role during C₄ evolution. The findings of the present study coincide with the observed repetitive emergence of C₄ during evolution.

Electronic supplementary material

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

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.