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

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.

Epigenetic Landscape Is Largely Shaped by Diversiform Transposons in Aegilops tauschii

Transposons (TEs) account for more than 80% of the wheat genome, the highest among all known crop species. They play an important role in shaping the elaborate genomic landscape, which is the key to the speciation of wheat. In this study, we analyzed the association between TEs, chromatin states, and chromatin accessibility in Aegilops tauschii, the D genome donor of bread wheat. We found that TEs contributed to the complex but orderly epigenetic landscape as chromatin states showed diverse distributions on TEs of different orders or superfamilies. TEs also contributed to the chromatin state and openness of potential regulatory elements, affecting the expression of TE-related genes. Some TE superfamilies, such as hAT-Ac, carry active/open chromatin regions. In addition, the histone mark H3K9ac was found to be associated with the accessibility shaped by TEs. These results suggest the role of diversiform TEs in shaping the epigenetic landscape and in gene expression regulation in Aegilops tauschii. This has positive implications for understanding the transposon roles in Aegilops tauschii or the wheat D genome.

THERMOSENSITIVE BARREN PANICLE (TAP) is required for rice panicle and spikelet development at high ambient temperature

In cereal plants, the size of the panicle (inflorescence) is a critical factor for yield. Panicle size is determined by a complex interplay of genetic and environmental factors, but the mechanisms underlying adaptations to temperature stress during panicle development remain largely unknown. We identify the rice THERMOSENSITIVE BARREN PANICLE (TAP) gene, which encodes a transposase-derived FAR1-RELATED SEQUENCE (FRS) protein and is responsible for regulating panicle and spikelet development at high ambient temperature. The tap mutants display high temperature-dependent reproductive abnormalities, including compromised secondary branch and spikelet initiation and pleiotropic floral organ defects. Consistent with its thermosensitive phenotype, TAP expression is induced by high temperature. TAP directly promotes the expression of OsYABBY3 (OsYAB3), OsYAB4, and OsYAB5, which encode key transcriptional regulators in panicle and spikelet development. In addition, TAP physically interacts with OsYAB4 and OsYAB5 proteins; phenotypic analysis of osyab4 tap-1 and osyab5 tap-1 double mutants indicates that TAP-OsYAB4/OsYAB5 complexes act to maintain normal panicle and spikelet development. Taken together, our study reveals the novel role of a TE-derived transcription factor in controlling rice panicle development under high ambient temperatures, shedding light on the molecular mechanism underlying the adaptation of cereal crops to increasing environmental temperatures.

Single-base resolution methylomes of somatic embryogenesis in Theobroma cacao L. reveal epigenome modifications associated with somatic embryo abnormalities

Propagation by somatic embryogenesis in Theobroma cacao has some issues to be solved, as many morphologically abnormal somatic embryos that do not germinate into plants are frequently observed, thus hampering plant production on a commercial scale. For the first time the methylome landscape of T. cacao somatic embryogenesis was examined, using whole-genome bisulfite sequencing technique, with the aim to understand the epigenetic basis of somatic embryo abnormalities. We identified 873 differentially methylated genes (DMGs) in the CpG context between zygotic embryos, normal and abnormal somatic embryos, with important roles in development, programmed cell death, oxidative stress, and hypoxia induction, which can help to explain the morphological abnormalities of somatic embryos. We also identified the role of ethylene and its precursor 1-aminocyclopropane-1-carboxylate in several biological processes, such as hypoxia induction, cell differentiation and cell polarity, that could be associated to the development of abnormal somatic embryos. The biological processes and the hypothesis of ethylene and its precursor involvement in the somatic embryo abnormalities in cacao are discussed.

Whole-genome DNA methylation patterns of Oryza sativa (L.) and Oryza glumaepatula (Steud) genotypes associated with aluminum response

Epigenetic mechanisms in crops have emerged as a fundamental factor in plant adaptation and acclimation to biotic and abiotic stresses. Among described epigenetic mechanisms, DNA methylation has been defined as the most studied epigenetic modification involved in several developmental processes. It has been shown that contrasting methylation marks are associated with gene expression variations between cultivated and wild crop species. In this study, we analyzed single-base resolution methylome maps for Oryza sativa (a cultivated species) and Oryza glumaepatula (a wild species) genotypes grown under control conditions. Our results showed that overall, genome-wide methylation profiles are mainly conserved between both species, nevertheless, there are several differentially methylated regions with species-specific methylation patterns. In addition, we analyzed the association of identified DNA methylation marks in relation with Aluminum-tolerance levels of studied genotypes. We found several differentially methylated regions (DMRs) and DMR-associated genes (DAGs) that are linked with Al tolerance. Some of these DAGs have been previously reported as differentially expressed under Al exposure in O. sativa. Complementarily a Transposable Elements (TE) analysis revealed that specific aluminum related genes have associated-TEs potentially regulated by DNA methylation. Interestingly, the DMRs and DAGs between Al-tolerant and susceptible genotypes were different between O. sativa and O. glumaepatula, suggesting that methylation patterns related to Al responses are unique for each rice species. Our findings provide novel insights into DNA methylation patterns in wild and cultivated rice genotypes and their possible role in the regulation of plant stress responses.

A null allele of the pol IV second subunit impacts stature and reproductive development in Oryza sativa

All eukaryotes possess three DNA-dependent RNA polymerases, Pols I-III, while land plants possess two additional polymerases, Pol IV and Pol V. Derived through duplication of Pol II subunits, Pol IV produces 24-nt short interfering RNAs that interact with Pol V transcripts to target de novo DNA methylation and silence transcription of transposons. Members of the grass family encode additional duplicated subunits of Pol IV and V, raising questions regarding the function of each paralog. In this study, we identify a null allele of the putative Pol IV second subunit, NRPD2, and demonstrate that NRPD2 is the sole subunit functioning with NRPD1 in small RNA production and CHH methylation in leaves. Homozygous nrpd2 mutants have neither gametophytic defects nor embryo lethality, although adult plants are dwarf and sterile.

Identification and Fine Mapping of the Candidate Gene Controlling Multi-Inflorescence in Brassica napus

As a desirable agricultural trait, multi-inflorescence (MI) fulfills the requirement of mechanized harvesting and yield increase in rapeseed (Brassica napus L.). However, the genetic mechanism underlying the multi-inflorescence trait remain poorly understood. We previously identified a difference of one pair of dominant genes between the two mapping parental materials. In this study, phenotype and expression analysis indicated that the imbalance of the CLAVATA (CLV)-WUSCHEL (WUS) feedback loop may contribute to the abnormal development of the shoot apical meristem (SAM). BnaMI was fine-mapped to a 55 kb genomic region combining with genotype and phenotype of 5768 BCF₁ individuals using a traditional mapping approach. Through comparative and expression analyses, combined with the annotation in Arabidopsis, five genes in this interval were identified as candidate genes. The present findings may provide assistance in functional analysis of the mechanism associated with multi-inflorescence and yield increase in rapeseed.

Temporal changes in transcripts of miniature inverted-repeat transposable elements during rice endosperm development

The repression of transcription from transposable elements (TEs) by DNA methylation is necessary to maintain genome integrity and prevent harmful mutations. However, under certain circumstances, TEs may escape from the host defense system and reactivate their transcription. In Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), DNA demethylases target the sequences derived from TEs in the central cell, the progenitor cell for the endosperm in the female gametophyte. Genome-wide DNA demethylation is also observed in the endosperm after fertilization. In the present study, we used a custom microarray to survey the transcripts generated from TEs during rice endosperm development and at selected time points in the embryo as a control. The expression patterns of TE transcripts are dynamically up- and downregulated during endosperm development, especially those of miniature inverted-repeat TEs (MITEs). Some TE transcripts were directionally controlled, whereas the other DNA transposons and retrotransposons were not. We also discovered the NUCLEAR FACTOR Y binding motif, CCAAT, in the region near the 5' terminal inverted repeat of Youren, one of the transcribed MITEs in the endosperm. Our results uncover dynamic changes in TE activity during endosperm development in rice.

Transposons Hidden in Arabidopsis thaliana Genome Assembly Gaps and Mobilization of Non-Autonomous LTR Retrotransposons Unravelled by Nanotei Pipeline

Long-read data is a great tool to discover new active transposable elements (TEs). However, no ready-to-use tools were available to gather this information from low coverage ONT datasets. Here, we developed a novel pipeline, nanotei, that allows detection of TE-contained structural variants, including individual TE transpositions. We exploited this pipeline to identify TE insertion in the Arabidopsis thaliana genome. Using nanotei, we identified tens of TE copies, including ones for the well-characterized ONSEN retrotransposon family that were hidden in genome assembly gaps. The results demonstrate that some TEs are inaccessible for analysis with the current A. thaliana (TAIR10.1) genome assembly. We further explored the mobilome of the ddm1 mutant with elevated TE activity. Nanotei captured all TEs previously known to be active in ddm1 and also identified transposition of non-autonomous TEs. Of them, one non-autonomous TE derived from (AT5TE33540) belongs to TR-GAG retrotransposons with a single open reading frame (ORF) encoding the GAG protein. These results provide the first direct evidence that TR-GAGs and other non-autonomous LTR retrotransposons can transpose in the plant genome, albeit in the absence of most of the encoded proteins. In summary, nanotei is a useful tool to detect active TEs and their insertions in plant genomes using low-coverage data from Nanopore genome sequencing.

Extensively Current Activity of Transposable Elements in Natural Rice Accessions Revealed by Singleton Insertions

Active transposable elements (TEs) have drawn more attention as they continue to create new insertions and contribute to genetic diversity of the genome. However, only a few have been discovered in rice up to now, and their activities are mostly induced by artificial treatments (e.g., tissue culture, hybridization etc.) rather than under normal growth conditions. To systematically survey the current activity of TEs in natural rice accessions and identify rice accessions carrying highly active TEs, the transposon insertion polymorphisms (TIPs) profile was used to identify singleton insertions, which were unique to a single accession and represented the new insertion of TEs in the genome. As a result, 10,924 high-confidence singletons from 251 TE families were obtained, covering all investigated TE types. The number of singletons varied substantially among different superfamilies/families, perhaps reflecting distinct current activity. Particularly, eight TE families maintained potentially higher activity in 3,000 natural rice accessions. Sixty percent of rice accessions were detected to contain singletons, indicating the extensive activity of TEs in natural rice accessions. Thirty-five TE families exhibited potentially high activity in at least one rice accession, and the majority of them showed variable activity among different rice groups/subgroups. These naturally active TEs would be ideal candidates for elucidating the molecular mechanisms underlying the transposition and activation of TEs, as well as investigating the interactions between TEs and the host genome.

Plant transfer RNA-derived fragments: Biogenesis and functions

Processing of mature transfer RNAs (tRNAs) produces complex populations of tRNA-derived fragments (tRFs). Emerging evidence shows that tRFs have important functions in bacteria, animals, and plants. Here, we review recent advances in understanding plant tRFs, focusing on their biological and cellular functions, such as regulating stress responses, mediating plant-pathogen interactions, and modulating post-transcriptional gene silencing and translation. We also review sequencing strategies and bioinformatics resources for studying tRFs in plants. Finally, we discuss future directions for plant tRF research, which will expand our knowledge of plant non-coding RNAs.

Deciphering the Epigenetic Alphabet Involved in Transgenerational Stress Memory in Crops

Although epigenetic modifications have been intensely investigated over the last decade due to their role in crop adaptation to rapid climate change, it is unclear which epigenetic changes are heritable and therefore transmitted to their progeny. The identification of epigenetic marks that are transmitted to the next generations is of primary importance for their use in breeding and for the development of new cultivars with a broad-spectrum of tolerance/resistance to abiotic and biotic stresses. In this review, we discuss general aspects of plant responses to environmental stresses and provide an overview of recent findings on the role of transgenerational epigenetic modifications in crops. In addition, we take the opportunity to describe the aims of EPI-CATCH, an international COST action consortium composed by researchers from 28 countries. The aim of this COST action launched in 2020 is: (1) to define standardized pipelines and methods used in the study of epigenetic mechanisms in plants, (2) update, share, and exchange findings in epigenetic responses to environmental stresses in plants, (3) develop new concepts and frontiers in plant epigenetics and epigenomics, (4) enhance dissemination, communication, and transfer of knowledge in plant epigenetics and epigenomics.

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.

Clustered and dispersed chromosomal distribution of the two classes of Revolver transposon family in rye (Secale cereale)

The chromosomal locations of a new class of Revolver transposon-like elements were analyzed by using FISH method on the metaphase chromosome in somatic cell division of the rye cultivar Petkus. First, the Revolver standard element probe λ2 was weakly hybridized throughout the rye chromosome, and comparatively large interstitial signals spotted with a dot shape were detected together with several telomeric regions. The dot shape interstitial signal was stably detected at one site on Chromosome (Chr) 1R (middle part of the interstitial region of the short arm), three sites on Chr 2R (distal part of the interstitial region and adjacent to the centromere on the short arm, middle part of the interstitial region of the long arm), and two sites on Chr 5R (middle part of the interstitial region and adjacent to the centromere on the long arm). The Revolver λ2 probe was effective for identification of 1R, 2R, and 5R chromosomes. On the other hand, Revolver nonautonomous element-specific L626-BARE-100 probe was strongly distributed throughout the rye chromosomes, and considerable numbers and diverse lengths of transcripts were detected by RT-PCR. Although the standard elements were found in localized clusters, the nonautonomous elements tended to be dispersed throughout the genome. Clustered nature of Revolver is a significantly rare case in genomics.

A copia-like retrotransposon insertion in the upstream region of the SHATTERPROOF1 gene, BnSHP1.A9, is associated with quantitative variation in pod shattering resistance in oilseed rape

Seed loss resulting from pod shattering is a major constraint in production of oilseed rape (Brassica napus L.). However, the molecular mechanisms underlying pod shatter resistance are not well understood. Here, we show that the pod shatter resistance at quantitative trait locus qSRI.A9.1 is controlled by one of the B. napus SHATTERPROOF1 homologs, BnSHP1.A9, in a doubled haploid population generated from parents designated R1 and R2 as well as in a diverse panel of oilseed rape. The R1 maternal parental line of the doubled haploid population carried the allele for shattering at qSRI.A9.1, while the R2 parental line carried the allele for shattering resistance. Quantitative RT-PCR showed that BnSHP1.A9 was expressed specifically in flower buds, flowers, and developing siliques in R1, while it was not expressed in any tissue of R2. Transgenic plants constitutively expressing either of the BnSHP1.A9 alleles from the R1 and R2 parental lines showed that both alleles are responsible for pod shattering, via a mechanism that promotes lignification of the enb layer. These findings indicated that the allelic differences in the BnSHP1.A9 gene per se are not the causal factor for quantitative variation in shattering resistance at qSRI.A9.1. Instead, a highly methylated copia-like long terminal repeat retrotransposon insertion (4803 bp) in the promotor region of the R2 allele of BnSHP1.A9 repressed the expression of BnSHP1.A9, and thus contributed to pod shatter resistance. Finally, we showed a copia-like retrotransposon-based marker, BnSHP1.A9R2, can be used for marker-assisted breeding targeting the pod shatter resistance trait in oilseed rape.

Epigenetic Landmarks of Leaf Senescence and Crop Improvement

This review synthesizes knowledge on epigenetic regulation of leaf senescence and discusses the possibility of using this knowledge to improve crop quality. This control level is implemented by different but interacting epigenetic mechanisms, including DNA methylation, covalent histone modifications, and non-covalent chromatin remodeling. The genetic and epigenetic changes may act alone or together and regulate the gene expression, which may result in heritable (stress memory) changes and may lead to crop survival. In the review, the question also arises whether the mitotically stable epigenetic information can be used for crop improvement. The barley crop model for early and late events of dark-induced leaf senescence (DILS), where the point of no return was defined, revealed differences in DNA and RNA modifications active in DILS compared to developmental leaf senescence. This suggests the possibility of a yet-to-be-discovered epigenetic-based switch between cell survival and cell death. Conclusions from the analyzed research contributed to the hypothesis that chromatin-remodeling mechanisms play a role in the control of induced leaf senescence. Understanding this mechanism in crops might provide a tool for further exploitation toward sustainable agriculture: so-called epibreeding.

Transcriptional regulation of genes bearing intronic heterochromatin in the rice genome

Intronic regions of eukaryotic genomes accumulate many Transposable Elements (TEs). Intronic TEs often trigger the formation of transcriptionally repressive heterochromatin, even within transcription-permissive chromatin environments. Although TE-bearing introns are widely observed in eukaryotic genomes, their epigenetic states, impacts on gene regulation and function, and their contributions to genetic diversity and evolution, remain poorly understood. In this study, we investigated the genome-wide distribution of intronic TEs and their epigenetic states in the Oryza sativa genome, where TEs comprise 35% of the genome. We found that over 10% of rice genes contain intronic heterochromatin, most of which are associated with TEs and repetitive sequences. These heterochromatic introns are longer and highly enriched in promoter-proximal positions. On the other hand, introns also accumulate hypomethylated short TEs. Genes with heterochromatic introns are implicated in various biological functions. Transcription of genes bearing intronic heterochromatin is regulated by an epigenetic mechanism involving the conserved factor OsIBM2, mutation of which results in severe developmental and reproductive defects. Furthermore, we found that heterochromatic introns evolve rapidly compared to non-heterochromatic introns. Our study demonstrates that heterochromatin is a common epigenetic feature associated with actively transcribed genes in the rice genome.

Gene duplication and transposition of mobile elements drive evolution of the Rpv3 resistance locus in grapevine

A wild grape haplotype (Rpv3-1) confers resistance to Plasmopara viticola. We mapped the causal factor for resistance to an interval containing a TIR-NB-LRR (TNL) gene pair that originated 1.6-2.6 million years ago by a tandem segmental duplication. Transient coexpression of the TNL pair in Vitis vinifera leaves activated pathogen-induced necrosis and reduced sporulation compared with control leaves. Even though transcripts of the TNL pair from the wild haplotype appear to be partially subject to nonsense-mediated mRNA decay, mature mRNA levels in a homozygous resistant genotype were individually higher than the mRNA trace levels observed for the orthologous single-copy TNL in sensitive genotypes. Allelic expression imbalance in a resistant heterozygote confirmed that cis-acting regulatory variation promotes expression in the wild haplotype. The movement of transposable elements had a major impact on the generation of haplotype diversity, altering the DNA context around similar TNL coding sequences and the GC-content in their proximal 5'-intergenic regions. The wild and domesticated haplotypes also diverged in conserved single-copy intergenic DNA, but the highest divergence was observed in intraspecific and not in interspecific comparisons. In this case, introgression breeding did not transgress the genetic boundaries of the domesticated species, because haplotypes present in modern varieties sometimes predate speciation events between wild and cultivated species.

CG hypomethylation leads to complex changes in DNA methylation and transpositional burst of diverse transposable elements in callus cultures of rice

CG methylation (^m CG) is essential for preserving genome stability in mammals, but this link remains obscure in plants. OsMET1-2, a major rice DNA methyltransferase, plays critical roles in maintaining ^m CG in rice. Null mutation of OsMET1-2 causes massive CG hypomethylation, rendering the mutant suitable to address the role of ^m CG in maintaining genome integrity in plants. Here, we analyzed ^m CG dynamics and genome stability in tissue cultures of OsMET1-2 homozygous (-/-) and heterozygous (+/-) mutants, and isogenic wild-type (WT). We found ^m CG levels in cultures of -/- were substantially lower than in those of WT and +/-, as expected. Unexpectedly, ^m CG levels in 1- and 3-year cultures of -/- were 77.6% and 48.7% higher, respectively, than in shoot, from which the cultures were initiated, suggesting substantial regain of ^m CG in -/- cultures, which contrasts to the general trend of ^m CG loss in all WT plant tissue cultures hitherto studied. Transpositional burst of diverse transposable elements (TEs) occurred only in -/- cultures, although no elevation of genome-wide mutation rate in the form of single nucleotide polymorphisms was detected. Altogether, our results establish an essential role of ^m CG in retaining TE immobility and hence genome stability in rice and likely in plants in general.

Interspecific variation across angiosperms in global DNA methylation: phylogeny, ecology and plant features in tropical and Mediterranean communities

The interspecific range of epigenetic variation and the degree to which differences between angiosperm species are related to geography, evolutionary history, ecological settings or species-specific traits, remain essentially unexplored. Genome-wide global DNA cytosine methylation is a tractable 'epiphenotypic' feature suitable for exploring these relationships. Global cytosine methylation was estimated in 279 species from two distant, ecologically disparate geographical regions: Mediterranean Spain and tropical México. At each region, four distinct plant communities were analyzed. Global methylation spanned a 10-fold range among species (4.8-42.2%). Interspecific differences were related to evolutionary trajectories, as denoted by a strong phylogenetic signal. Genomes of tropical species were on average less methylated than those of Mediterranean ones. Woody plants have genomes with lower methylation than perennial herbs, and genomes of widespread species were less methylated than those of species with restricted geographical distribution. The eight communities studied exhibited broad and overlapping interspecific variances in global cytosine methylation and only two of them differed in average methylation. Altogether, our broad taxonomic survey supported global methylation as a plant 'epiphenotypic' trait largely associated with species evolutionary history, genome size, range size and woodiness. Additional studies are required for better understanding the environmental components underlying local and geographical variation.

DNA methylation repatterning accompanying hybridization, whole genome doubling and homoeolog exchange in nascent segmental rice allotetraploids

Allopolyploidization, which entails interspecific hybridization and whole genome duplication (WGD), is associated with emergent genetic and epigenetic instabilities that are thought to contribute to adaptation and evolution. One frequent genomic consequence of nascent allopolyploidization is homoeologous exchange (HE), which arises from compromised meiotic fidelity and generates genetically and phenotypically variable progenies. Here, we used a genetically tractable synthetic rice segmental allotetraploid system to interrogate genome-wide DNA methylation and gene expression responses and outcomes to the separate and combined effects of hybridization, WGD and HEs. Progenies of the tetraploid rice were genomically diverse due to genome-wide HEs that affected all chromosomes, yet they exhibited overall methylome stability. Nonetheless, regional variation of cytosine methylation states was widespread in the tetraploids. Transcriptome profiling revealed genome-wide alteration of gene expression, which at least in part associates with changes in DNA methylation. Intriguingly, changes of DNA methylation and gene expression could be decoupled from hybridity and sustained and amplified by HEs. Our results suggest that HEs, a prominent genetic consequence of nascent allopolyploidy, can exacerbate, diversify and perpetuate the effects of allopolyploidization on epigenetic and gene expression variation, and hence may contribute to allopolyploid evolution.

Identification of an active miniature inverted-repeat transposable element mJing in rice

Miniature inverted-repeat transposable elements (MITEs) are structurally homogeneous non-autonomous DNA transposons with high copy numbers that play important roles in genome evolution and diversification. Here, we analyzed the rice high-tillering dwarf (htd) mutant in an advanced backcross population between cultivated and wild rice, and identified an active MITE named miniature Jing (mJing). The mJing element belongs to the PIF/Harbinger superfamily. japonica rice var. Nipponbare and indica var. 93-11 harbor 72 and 79 mJing family members, respectively, have undergone multiple rounds of amplification bursts during the evolution of Asian cultivated rice (Oryza sativa L.). A heterologous transposition experiment in Arabidopsis thaliana indicated that the autonomous element Jing is likely to have provides the transposase needed for mJing mobilization. We identified 297 mJing insertion sites and their presence/absence polymorphism among 71 rice samples through targeted high-throughput sequencing. The results showed that the copy number of mJing varies dramatically among Asian cultivated rice (O. sativa), its wild ancestor (O. rufipogon), and African cultivated rice (O. glaberrima) and that some mJing insertions are subject to directional selection. These findings suggest that the amplification and removal of mJing elements have played an important role in rice genome evolution and species diversification.

Prunus genetics and applications after de novo genome sequencing: achievements and prospects

Prior to the availability of whole-genome sequences, our understanding of the structural and functional aspects of Prunus tree genomes was limited mostly to molecular genetic mapping of important traits and development of EST resources. With public release of the peach genome and others that followed, significant advances in our knowledge of Prunus genomes and the genetic underpinnings of important traits ensued. In this review, we highlight key achievements in Prunus genetics and breeding driven by the availability of these whole-genome sequences. Within the structural and evolutionary contexts, we summarize: (1) the current status of Prunus whole-genome sequences; (2) preliminary and ongoing work on the sequence structure and diversity of the genomes; (3) the analyses of Prunus genome evolution driven by natural and man-made selection; and (4) provide insight into haploblocking genomes as a means to define genome-scale patterns of evolution that can be leveraged for trait selection in pedigree-based Prunus tree breeding programs worldwide. Functionally, we summarize recent and ongoing work that leverages whole-genome sequences to identify and characterize genes controlling 22 agronomically important Prunus traits. These include phenology, fruit quality, allergens, disease resistance, tree architecture, and self-incompatibility. Translationally, we explore the application of sequence-based marker-assisted breeding technologies and other sequence-guided biotechnological approaches for Prunus crop improvement. Finally, we present the current status of publically available Prunus genomics and genetics data housed mainly in the Genome Database for Rosaceae (GDR) and its updated functionalities for future bioinformatics-based Prunus genetics and genomics inquiry.

A high-quality apple genome assembly reveals the association of a retrotransposon and red fruit colour

A complete and accurate genome sequence provides a fundamental tool for functional genomics and DNA-informed breeding. Here, we assemble a high-quality genome (contig N50 of 6.99 Mb) of the apple anther-derived homozygous line HFTH1, including 22 telomere sequences, using a combination of PacBio single-molecule real-time (SMRT) sequencing, chromosome conformation capture (Hi-C) sequencing, and optical mapping. In comparison to the Golden Delicious reference genome, we identify 18,047 deletions, 12,101 insertions and 14 large inversions. We reveal that these extensive genomic variations are largely attributable to activity of transposable elements. Interestingly, we find that a long terminal repeat (LTR) retrotransposon insertion upstream of MdMYB1, a core transcriptional activator of anthocyanin biosynthesis, is associated with red-skinned phenotype. This finding provides insights into the molecular mechanisms underlying red fruit coloration, and highlights the utility of this high-quality genome assembly in deciphering agriculturally important trait in apple.

On the Importance to Acknowledge Transposable Elements in Epigenomic Analyses

Eukaryotic genomes comprise a large proportion of repeated sequences, an important fraction of which are transposable elements (TEs). TEs are mobile elements that have a significant impact on genome evolution and on gene functioning. Although some TE insertions could provide adaptive advantages to species, transposition is a highly mutagenic event that has to be tightly controlled to ensure its viability. Genomes have evolved sophisticated mechanisms to control TE activity, the most important being epigenetic silencing. However, the epigenetic control of TEs can also affect genes located nearby that can become epigenetically regulated. It has been proposed that the combination of TE mobilization and the induced changes in the epigenetic landscape could allow a rapid phenotypic adaptation to global environmental changes. In this review, we argue the crucial need to take into account the repeated part of genomes when studying the global impact of epigenetic modifications on an organism. We emphasize more particularly why it is important to carefully consider TEs and what bioinformatic tools can be used to do so.

PePIF1, a P-lineage of PIF-like transposable element identified in protocorm-like bodies of Phalaenopsis orchids

BACKGROUND

Orchids produce a colorless protocorm by symbiosis with fungi upon seed germination. For mass production of orchids, the prevailing approaches are both generation of protocorm-like bodies (PLBs) from callus and multiplication of adventitious buds on inflorescence. However, somaclonal variations occur during micropropagation.

RESULTS

We isolated the two most expressed transposable elements belonging to P Instability Factor (PIF)-like transposons. Among them, a potential autonomous element was identified by similarity analysis against the whole-genome sequence of Phalaenopsis equestris and named PePIF1. It contains a 19-bp terminal inverted repeat flanked by a 3-bp target site duplication and two coding regions encoding ORF1- and transposase-like proteins. Phylogenetic analysis revealed that PePIF1 belongs to a new P-lineage of PIF. Furthermore, two distinct families, PePIF1a and PePIF1b, with 29 and 37 putative autonomous elements, respectively, were isolated, along with more than 3000 non-autonomous and miniature inverted-repeat transposable element (MITE)-like elements. Among them, 828 PePIF1-related elements were inserted in 771 predicted genes. Intriguingly, PePIF1 was transposed in the somaclonal variants of Phalaenopsis cultivars, as revealed by transposon display, and the newly inserted genes were identified and sequenced.

CONCLUSION

A PIF-like element, PePIF1, was identified in the Phalaenopsis genome and actively transposed during micropropagation. With the identification of PePIF1, we have more understanding of the Phalaenopsis genome structure and somaclonal variations during micropropagation for use in orchid breeding and production.

CsMYB60 is a key regulator of flavonols and proanthocyanidans that determine the colour of fruit spines in cucumber

Spine colour is an important fruit quality trait that influences the commercial value of cucumber (Cucumis sativus). However, little is known about the metabolites and the regulatory mechanisms of their biosynthesis in black spine varieties. In this study, we determined that the pigments of black spines are flavonoids, including flavonols and proanthocyanidins (PAs). We identified CsMYB60 as the best candidate for the previously identified B (Black spine) locus. Expression levels of CsMYB60 and the key genes involved in flavonoid biosynthesis were higher in black-spine inbred lines than that in white-spine lines at different developmental stages. The insertion of a Mutator-like element (CsMULE) in the second intron of CsMYB60 decreased its expression in a white-spine line. Transient overexpression assays indicated that CsMYB60 is a key regulatory gene and Cs4CL is a key structural gene in the pigmentation of black spines. In addition, the DNA methylation level in the CsMYB60 promoter was much lower in the black-spine line compared with white-spine line. The CsMULE insert may decrease the expression level of CsMYB60, causing hindered synthesis of flavonols and PAs in cucumber fruit spines.

Exaptation at the molecular genetic level

The realization that body parts of animals and plants can be recruited or coopted for novel functions dates back to, or even predates the observations of Darwin. S.J. Gould and E.S. Vrba recognized a mode of evolution of characters that differs from adaptation. The umbrella term aptation was supplemented with the concept of exaptation. Unlike adaptations, which are restricted to features built by selection for their current role, exaptations are features that currently enhance fitness, even though their present role was not a result of natural selection. Exaptations can also arise from nonaptations; these are characters which had previously been evolving neutrally. All nonaptations are potential exaptations. The concept of exaptation was expanded to the molecular genetic level which aided greatly in understanding the enormous potential of neutrally evolving repetitive DNA-including transposed elements, formerly considered junk DNA-for the evolution of genes and genomes. The distinction between adaptations and exaptations is outlined in this review and examples are given. Also elaborated on is the fact that such distinctions are sometimes more difficult to determine; this is a widespread phenomenon in biology, where continua abound and clear borders between states and definitions are rare.

MITE Tracker: an accurate approach to identify miniature inverted-repeat transposable elements in large genomes

BACKGROUND

Miniature inverted-repeat transposable elements (MITEs) are short, non-autonomous class II transposable elements present in a high number of conserved copies in eukaryote genomes. An accurate identification of these elements can help to shed light on the mechanisms controlling genome evolution and gene regulation. The structure and distribution of these elements are well-defined and therefore computational approaches can be used to identify MITEs sequences.

RESULTS

Here we describe MITE Tracker, a novel, open source software program that finds and classifies MITEs using an efficient alignment strategy to retrieve nearby inverted-repeat sequences from large genomes. This program groups them into high sequence homology families using a fast clustering algorithm and finally filters only those elements that were likely transposed from different genomic locations because of their low scoring flanking sequence alignment.

CONCLUSIONS

Many programs have been proposed to find MITEs hidden in genomes. However, none of them are able to process large-scale genomes such as that of bread wheat. Furthermore, in many cases the existing methods perform high false-positive rates (or miss rates). The rice genome was used as reference to compare MITE Tracker against known tools. Our method turned out to be the most reliable in our tests. Indeed, it revealed more known elements, presented the lowest false-positive number and was the only program able to run with the bread wheat genome as input. In wheat, MITE Tracker discovered 6013 MITE families and allowed the first structural exploration of MITEs in the complete bread wheat genome.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Zebularine treatment is associated with deletion of FT-B1 leading to an increase in spikelet number in bread wheat

The number of rachis nodes (spikelets) on a wheat spike is a component of grain yield that correlates with flowering time. The genetic basis regulating flowering in cereals is well understood, but there are reports that flowering time can be modified at a high frequency by selective breeding, suggesting that it may be regulated by both epigenetic and genetic mechanisms. We investigated the role of DNA methylation in regulating spikelet number and flowering time by treating a semi-spring wheat with the demethylating agent, Zebularine. Three lines with a heritable increase in spikelet number were identified. The molecular basis for increased spikelet number was not determined in 2 lines, but the phenotype showed non-Mendelian inheritance, suggesting that it could have an epigenetic basis. In the remaining line, the increased spikelet phenotype behaved as a Mendelian recessive trait and late flowering was associated with a deletion encompassing the floral promoter, FT-B1. Deletion of FT-B1 delayed the transition to reproductive growth, extended the duration of spike development, and increased spikelet number under different temperature regimes and photoperiod. Transiently disrupting DNA methylation can generate novel flowering behaviour in wheat, but these changes may not be sufficiently stable for use in breeding programs.

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

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

Transposons: a blessing curse

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

Small RNA pathways responsible for non-cell-autonomous regulation of plant reproduction

In angiosperms, germline precursors and germ cells are always attached to or engulfed within somatic companion cells until just before fertilization. This is because sperm and egg cells develop as part of the multicellular gametophyte. Thus, the non-cell-autonomous regulation by somatic companions plays important roles in efficient reproduction, in addition to the cell-autonomous regulation. Epigenetic silencing of transposable elements is one of the central events by which the germline transmits the error-free genome to the next generation. This review focuses on small RNA-mediated epigenetic regulation of meiosis, spore formation and pollen development. Besides microRNA (miRNA) and small interfering RNA (siRNA), animals express PIWI-interacting RNA (piRNA), a germline-specific class of small RNAs. Plants lack piRNA-like RNAs and, instead, express unique classes of small RNAs: trans-acting siRNA (tasiRNA) and phased secondary siRNA (phasiRNA). Especially in grass species, 21- and 24-nucleotide phasiRNAs are abundant in anthers during premeiosis and meiosis. This review also describes recent progress in reproductive phasiRNA research.

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

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

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

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

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

Assessing and Exploiting Functional Diversity in Germplasm Pools to Enhance Abiotic Stress Adaptation and Yield in Cereals and Food Legumes

There is a need to accelerate crop improvement by introducing alleles conferring host plant resistance, abiotic stress adaptation, and high yield potential. Elite cultivars, landraces and wild relatives harbor useful genetic variation that needs to be more easily utilized in plant breeding. We review genome-wide approaches for assessing and identifying alleles associated with desirable agronomic traits in diverse germplasm pools of cereals and legumes. Major quantitative trait loci and single nucleotide polymorphisms (SNPs) associated with desirable agronomic traits have been deployed to enhance crop productivity and resilience. These include alleles associated with variation conferring enhanced photoperiod and flowering traits. Genetic variants in the florigen pathway can provide both environmental flexibility and improved yields. SNPs associated with length of growing season and tolerance to abiotic stresses (precipitation, high temperature) are valuable resources for accelerating breeding for drought-prone environments. Both genomic selection and genome editing can also harness allelic diversity and increase productivity by improving multiple traits, including phenology, plant architecture, yield potential and adaptation to abiotic stresses. Discovering rare alleles and useful haplotypes also provides opportunities to enhance abiotic stress adaptation, while epigenetic variation has potential to enhance abiotic stress adaptation and productivity in crops. By reviewing current knowledge on specific traits and their genetic basis, we highlight recent developments in the understanding of crop functional diversity and identify potential candidate genes for future use. The storage and integration of genetic, genomic and phenotypic information will play an important role in ensuring broad and rapid application of novel genetic discoveries by the plant breeding community. Exploiting alleles for yield-related traits would allow improvement of selection efficiency and overall genetic gain of multigenic traits. An integrated approach involving multiple stakeholders specializing in management and utilization of genetic resources, crop breeding, molecular biology and genomics, agronomy, stress tolerance, and reproductive/seed biology will help to address the global challenge of ensuring food security in the face of growing resource demands and climate change induced stresses.