Transposons and genome evolution in plants

Genome-wide study of glutathione transferases and their regulation in flufenacet susceptible and resistant black-grass (Alopecurus myosuroides Huds.)

BACKGROUND

Glutathione transferases (GSTs) are enzymes with a wide range of functions, including herbicide detoxification. Up-regulation of GSTs and their detoxification activity enables the grass weed black-grass (Alopecurus myosuroides Huds.) to metabolize the very-long-chain fatty acid synthesis inhibitor flufenacet and other herbicides leading to multiple herbicide resistance. However, the genomic organization and regulation of GSTs genes is still poorly understood.

RESULTS

In this genome-wide study the location and expression of 115 GSTs were investigated using a recently published black-grass genome. Particularly, the most abundant GSTs of class tau and phi were typically clustered and often followed similar expression patterns but possessed divergent upstream regulatory regions. Similarities were found in the promoters of the most up-regulated GSTs, which are located next to each other in a cluster. The binding motif of the E2F/DP transcription factor complex in the promoter of an up-regulated GST was identical in susceptible and resistant plants, however, adjacent sequences differed. This led to a stronger binding of proteins to the motif of the susceptible plant, indicating repressor activity.

CONCLUSIONS

This study constitutes the first analysis dealing with the genomic investigation of GST genes found in black-grass and their transcriptional regulation. It highlights the complexity of the evolution of GSTs in black-grass, their duplication and divergence over time. The large number of GSTs allows weeds to detoxify a broad spectrum of herbicides. Ultimately, more research is needed to fully elucidate the regulatory mechanisms of GST expression. © 2024 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Molecular mechanisms and regulation of recombination frequency and distribution in plants

KEY MESSAGE

Recent developments in understanding the distribution and distinctive features of recombination hotspots are reviewed and approaches are proposed to increase recombination frequency in coldspot regions. Recombination events during meiosis provide the foundation and premise for creating new varieties of crops. The frequency of recombination in different genomic regions differs across eukaryote species, with recombination generally occurring more frequently at the ends of chromosomes. In most crop species, recombination is rare in centromeric regions. If a desired gene variant is linked in repulsion with an undesired variant of a second gene in a region with a low recombination rate, obtaining a recombinant plant combining two favorable alleles will be challenging. Traditional crop breeding involves combining desirable genes from parental plants into offspring. Therefore, understanding the mechanisms of recombination and factors affecting the occurrence of meiotic recombination is important for crop breeding. Here, we review chromosome recombination types, recombination mechanisms, genes and proteins involved in the meiotic recombination process, recombination hotspots and their regulation systems and discuss how to increase recombination frequency in recombination coldspot regions.

A super-pangenome of the North American wild grape species

BACKGROUND

Capturing the genetic diversity of wild relatives is crucial for improving crops because wild species are valuable sources of agronomic traits that are essential to enhance the sustainability and adaptability of domesticated cultivars. Genetic diversity across a genus can be captured in super-pangenomes, which provide a framework for interpreting genomic variations.

RESULTS

Here we report the sequencing, assembly, and annotation of nine wild North American grape genomes, which are phased and scaffolded at chromosome scale. We generate a reference-unbiased super-pangenome using pairwise whole-genome alignment methods, revealing the extent of the genomic diversity among wild grape species from sequence to gene level. The pangenome graph captures genomic variation between haplotypes within a species and across the different species, and it accurately assesses the similarity of hybrids to their parents. The species selected to build the pangenome are a great representation of the genus, as illustrated by capturing known allelic variants in the sex-determining region and for Pierce's disease resistance loci. Using pangenome-wide association analysis, we demonstrate the utility of the super-pangenome by effectively mapping short reads from genus-wide samples and identifying loci associated with salt tolerance in natural populations of grapes.

CONCLUSIONS

This study highlights how a reference-unbiased super-pangenome can reveal the genetic basis of adaptive traits from wild relatives and accelerate crop breeding research.

Gene mining and genomics-assisted breeding empowered by the pangenome of tea plant Camellia sinensis

Tea is one of the world's oldest crops and is cultivated to produce beverages with various flavours. Despite advances in sequencing technologies, the genetic mechanisms underlying key agronomic traits of tea remain unclear. In this study, we present a high-quality pangenome of 22 elite cultivars, representing broad genetic diversity in the species. Our analysis reveals that a recent long terminal repeat burst contributed nearly 20% of gene copies, introducing functional genetic variants that affect phenotypes such as leaf colour. Our graphical pangenome improves the efficiency of genome-wide association studies and allows the identification of key genes controlling bud flush timing. We also identified strong correlations between allelic variants and flavour-related chemistries. These findings deepen our understanding of the genetic basis of tea quality and provide valuable genomic resources to facilitate its genomics-assisted breeding.

The blackgrass genome reveals patterns of non-parallel evolution of polygenic herbicide resistance

Globally, weedy plants are a major constraint to sustainable crop production. Much of the success of weeds rests with their ability to rapidly adapt in the face of human-mediated management of agroecosystems. Alopecurus myosuroides (blackgrass) is a widespread and impactful weed affecting agriculture in Europe. Here we report a chromosome-scale genome assembly of blackgrass and use this reference genome to explore the genomic/genetic basis of non-target site herbicide resistance (NTSR). Based on our analysis of F2 seed families derived from two distinct blackgrass populations with the same NTSR phenotype, we demonstrate that the trait is polygenic and evolves from standing genetic variation. We present evidence that selection for NTSR has signatures of both parallel and non-parallel evolution. There are parallel and non-parallel changes at the transcriptional level of several stress- and defence-responsive gene families. At the genomic level, however, the genetic loci underpinning NTSR are different (non-parallel) between seed families. We speculate that variation in the number, regulation and function of stress- and defence-related gene families enable weedy species to rapidly evolve NTSR via exaptation of genes within large multi-functional gene families. These results provide novel insights into the potential for, and nature of plant adaptation in rapidly changing environments.

Predicting Cloned Disease Resistance Gene Homologs (CDRHs) in Radish, Underutilised Oilseeds, and Wild Brassicaceae Species

Brassicaceae crops, including Brassica, Camelina and Raphanus species, are among the most economically important crops globally; however, their production is affected by several diseases. To predict cloned disease resistance (R) gene homologs (CDRHs), we used the protein sequences of 49 cloned R genes against fungal and bacterial diseases in Brassicaceae species. In this study, using 20 Brassicaceae genomes (17 wild and 3 domesticated species), 3172 resistance gene analogs (RGAs) (2062 nucleotide binding-site leucine-rich repeats (NLRs), 497 receptor-like protein kinases (RLKs) and 613 receptor-like proteins (RLPs)) were identified. CDRH clusters were also observed in Arabis alpina, Camelina sativa and Cardamine hirsuta with assigned chromosomes, consisting of 62 homogeneous (38 NLR, 17 RLK and 7 RLP clusters) and 10 heterogeneous RGA clusters. This study highlights the prevalence of CDRHs in the wild relatives of the Brassicaceae family, which may lay the foundation for rapid identification of functional genes and genomics-assisted breeding to develop improved disease-resistant Brassicaceae crop cultivars.

Genome sequencing and analysis uncover the regulatory elements involved in the development and oil biosynthesis of Pongamia pinnata (L.) - A potential biodiesel feedstock

Due to rapid industrialization, the consumption of petro-products has increased, while fossil fuel resources have been gradually depleted. There has been a resurgence of interest in plant-derived biofuels as a sustainable alternative to fossil fuels for the purpose of reducing greenhouse gas emissions. Pongamia pinnata L., which is also known as Millettia pinnata is an oil-yielding, leguminous tree with a large and complex genome. Despite its multiple industrial applications, this orphan tree species has inconsistent yields and a limited understanding of its functional genomics. We assessed physiological and morphological characteristics of five high-yielding pongamia accessions and deduced important yield descriptors. Furthermore, we sequenced the genome of this potential biofuel feedstock using Illumina HiSeq, NextSeq, and MiSeq platforms to generate paired-end reads. Around 173 million processed reads amounting to 65.2 Gb were assembled into a 685 Mb genome, with a gap rate of 0.02%. The sequenced scaffolds were used to identify 30,000 gene models, 406,385 Simple-Sequence-Repeat (SSR) markers, and 43.6% of repetitive sequences. We further analyzed the structural information of genes belonging to certain key metabolic pathways, including lipid metabolism, photosynthesis, circadian rhythms, plant-pathogen interactions, and karanjin biosynthesis, all of which are commercially significant for pongamia. A total of 2,219 scaffolds corresponding to 29 transcription factor families provided valuable information about gene regulation in pongamia. Similarity studies and phylogenetic analysis revealed a monophyletic group of Fabaceae members wherein pongamia out-grouped from Glycine max and Cajanus cajan, revealing its unique ability to synthesize oil for biodiesel. This study is the first step toward completing the genome sequence of this imminent biofuel tree species. Further attempts at re-sequencing with different read chemistry will certainly improve the genetic resources at the chromosome level and accelerate the molecular breeding programs.

Mining of Cloned Disease Resistance Gene Homologs (CDRHs) in Brassica Species and Arabidopsis thaliana

Simple Summary

Developing cultivars with resistance genes (R genes) is an effective strategy to support high yield and quality in Brassica crops. The availability of clone R gene and genomic sequences in Brassica species and Arabidopsis thaliana provide the opportunity to compare genomic regions and survey R genes across genomic databases. In this paper, we aim to identify genes related to cloned genes through sequence identity, providing a repertoire of species-wide related R genes in Brassica crops. The comprehensive list of candidate R genes can be used as a reference for functional analysis.

Abstract

Various diseases severely affect Brassica crops, leading to significant global yield losses and a reduction in crop quality. In this study, we used the complete protein sequences of 49 cloned resistance genes (R genes) that confer resistance to fungal and bacterial diseases known to impact species in the Brassicaceae family. Homology searches were carried out across Brassica napus, B. rapa, B. oleracea, B. nigra, B. juncea, B. carinata and Arabidopsis thaliana genomes. In total, 660 cloned disease R gene homologs (CDRHs) were identified across the seven species, including 431 resistance gene analogs (RGAs) (248 nucleotide binding site-leucine rich repeats (NLRs), 150 receptor-like protein kinases (RLKs) and 33 receptor-like proteins (RLPs)) and 229 non-RGAs. Based on the position and distribution of specific homologs in each of the species, we observed a total of 87 CDRH clusters composed of 36 NLR, 16 RLK and 3 RLP homogeneous clusters and 32 heterogeneous clusters. The CDRHs detected consistently across the seven species are candidates that can be investigated for broad-spectrum resistance, potentially providing resistance to multiple pathogens. The R genes identified in this study provide a novel resource for the future functional analysis and gene cloning of Brassicaceae R genes towards crop improvement.

The Genome of the CTG(Ser1) Yeast Scheffersomyces stipitis Is Plastic

Microorganisms need to adapt to environmental changes, and genome plasticity can lead to rapid adaptation to hostile environments by increasing genetic diversity. Here, we investigate genome plasticity in the CTG(Ser1) yeast Scheffersomyces stipitis, an organism with an enormous potential for second-generation biofuel production. We demonstrate that S. stipitis has an intrinsically plastic genome and that different S. stipitis isolates have genomes with distinct chromosome organizations. Real-time evolution experiments show that S. stipitis genome plasticity is common and rapid since extensive genomic changes with fitness benefits are detected following in vitro evolution experiments. Hybrid MinION Nanopore and Illumina genome sequencing identify retrotransposons as major drivers of genome diversity. Indeed, the number and position of retrotransposons are different in different S. stipitis isolates, and retrotransposon-rich regions of the genome are sites of chromosome rearrangements. Our findings provide important insights into the adaptation strategies of the CTG(Ser1) yeast clade and have critical implications in the development of second-generation biofuels. These data highlight that genome plasticity is an essential factor for developing sustainable S. stipitis platforms for second-generation biofuels production. IMPORTANCE Genomes contain genes encoding the information needed to build the organism and allow it to grow and develop. Genomes are described as stable structures where genes have specific positions within a chromosome. Changes in gene dosage and position are viewed as harmful. However, it is becoming increasingly clear that genome plasticity can benefit microbial organisms that need to adapt rapidly to environmental changes. Mechanisms of genome plasticity are still poorly understood. This study focuses on Scheffersomyces stipitis, a yeast that holds great potential for second-generation biofuel production generated from forestry and agriculture waste. We demonstrate that S. stipitis chromosomes are easily reshuffled and that chromosome reshuffling is linked to adaptation to hostile environments. Genome sequencing demonstrates that mobile genetic elements, called transposons, mediate S. stipitis genome reshuffling. These data highlight that understanding genome plasticity is important for developing sustainable S. stipitis platforms for second-generation biofuels production.

The Dynamism of Transposon Methylation for Plant Development and Stress Adaptation

Plant development processes are regulated by epigenetic alterations that shape nuclear structure, gene expression, and phenotypic plasticity; these alterations can provide the plant with protection from environmental stresses. During plant growth and development, these processes play a significant role in regulating gene expression to remodel chromatin structure. These epigenetic alterations are mainly regulated by transposable elements (TEs) whose abundance in plant genomes results in their interaction with genomes. Thus, TEs are the main source of epigenetic changes and form a substantial part of the plant genome. Furthermore, TEs can be activated under stress conditions, and activated elements cause mutagenic effects and substantial genetic variability. This introduces novel gene functions and structural variation in the insertion sites and primarily contributes to epigenetic modifications. Altogether, these modifications indirectly or directly provide the ability to withstand environmental stresses. In recent years, many studies have shown that TE methylation plays a major role in the evolution of the plant genome through epigenetic process that regulate gene imprinting, thereby upholding genome stability. The induced genetic rearrangements and insertions of mobile genetic elements in regions of active euchromatin contribute to genome alteration, leading to genomic stress. These TE-mediated epigenetic modifications lead to phenotypic diversity, genetic variation, and environmental stress tolerance. Thus, TE methylation is essential for plant evolution and stress adaptation, and TEs hold a relevant military position in the plant genome. High-throughput techniques have greatly advanced the understanding of TE-mediated gene expression and its associations with genome methylation and suggest that controlled mobilization of TEs could be used for crop breeding. However, development application in this area has been limited, and an integrated view of TE function and subsequent processes is lacking. In this review, we explore the enormous diversity and likely functions of the TE repertoire in adaptive evolution and discuss some recent examples of how TEs impact gene expression in plant development and stress adaptation.

Potential evidence for transgenerational epigenetic memory in Arabidopsis thaliana following spaceflight

Plants grown in spaceflight exhibited differential methylation responses and this is important because plants are sessile, they are constantly exposed to a variety of environmental pressures and respond to them in many ways. We previously showed that the Arabidopsis genome exhibited lower methylation level after spaceflight for 60 h in orbit. Here, using the offspring of the seedlings grown in microgravity environment in the SJ-10 satellite for 11 days and returned to Earth, we systematically studied the potential effects of spaceflight on DNA methylation, transcriptome, and phenotype in the offspring. Whole-genome methylation analysis in the first generation of offspring (F₁) showed that, although there was no significant difference in methylation level as had previously been observed in the parent plants, some residual imprints of DNA methylation differences were detected. Combined DNA methylation and RNA-sequencing analysis indicated that expression of many pathways, such as the abscisic acid-activated pathway, protein phosphorylation, and nitrate signaling pathway, etc. were enriched in the F₁ population. As some phenotypic differences still existed in the F₂ generation, it was suggested that these epigenetic DNA methylation modifications were partially retained, resulting in phenotypic differences in the offspring. Furthermore, some of the spaceflight-induced heritable differentially methylated regions (DMRs) were retained. Changes in epigenetic modifications caused by spaceflight affected the growth of two future seed generations. Altogether, our research is helpful in better understanding the adaptation mechanism of plants to the spaceflight environment.

In order to investigate whether the effects of spaceflight on plants persist in future generations, Xu et al studied the offspring of Arabidopsis thaliana seedlings that had been grown in a microgravity environment for 11 days. They found that epigenetic modifications caused by spaceflight potentially affected the growth of two future seed generations, shedding light on the longevity of the effects of spaceflight on plants.

Whole-genome microsynteny-based phylogeny of angiosperms

Plant genomes vary greatly in size, organization, and architecture. Such structural differences may be highly relevant for inference of genome evolution dynamics and phylogeny. Indeed, microsynteny-the conservation of local gene content and order-is recognized as a valuable source of phylogenetic information, but its use for the inference of large phylogenies has been limited. Here, by combining synteny network analysis, matrix representation, and maximum likelihood phylogenetic inference, we provide a way to reconstruct phylogenies based on microsynteny information. Both simulations and use of empirical data sets show our method to be accurate, consistent, and widely applicable. As an example, we focus on the analysis of a large-scale whole-genome data set for angiosperms, including more than 120 available high-quality genomes, representing more than 50 different plant families and 30 orders. Our 'microsynteny-based' tree is largely congruent with phylogenies proposed based on more traditional sequence alignment-based methods and current phylogenetic classifications but differs for some long-contested and controversial relationships. For instance, our synteny-based tree finds Vitales as early diverging eudicots, Saxifragales within superasterids, and magnoliids as sister to monocots. We discuss how synteny-based phylogenetic inference can complement traditional methods and could provide additional insights into some long-standing controversial phylogenetic relationships.

The Terrestrial Carnivorous Plant Utricularia reniformis Sheds Light on Environmental and Life-Form Genome Plasticity

Utricularia belongs to Lentibulariaceae, a widespread family of carnivorous plants that possess ultra-small and highly dynamic nuclear genomes. It has been shown that the Lentibulariaceae genomes have been shaped by transposable elements expansion and loss, and multiple rounds of whole-genome duplications (WGD), making the family a platform for evolutionary and comparative genomics studies. To explore the evolution of Utricularia, we estimated the chromosome number and genome size, as well as sequenced the terrestrial bladderwort Utricularia reniformis (2n = 40, 1C = 317.1-Mpb). Here, we report a high quality 304 Mb draft genome, with a scaffold NG50 of 466-Kb, a BUSCO completeness of 87.8%, and 42,582 predicted genes. Compared to the smaller and aquatic U. gibba genome (101 Mb) that has a 32% repetitive sequence, the U. reniformis genome is highly repetitive (56%). The structural differences between the two genomes are the result of distinct fractionation and rearrangements after WGD, and massive proliferation of LTR-retrotransposons. Moreover, GO enrichment analyses suggest an ongoing gene birth-death-innovation process occurring among the tandem duplicated genes, shaping the evolution of carnivory-associated functions. We also identified unique patterns of developmentally related genes that support the terrestrial life-form and body plan of U. reniformis. Collectively, our results provided additional insights into the evolution of the plastic and specialized Lentibulariaceae genomes.

Autopolyploidization in switchgrass alters phenotype and flowering time via epigenetic and transcription regulation

Polyploidization is a significant source of genomic and organism diversification during plant evolution, and leads to substantial alterations in plant phenotypes and natural fitness. To help understand the phenotypic and molecular impacts of autopolyploidization, we conducted epigenetic and full-transcriptomic analyses of a synthesized autopolyploid accession of switchgrass (Panicum virgatum) in order to interpret the molecular and phenotypic changes. We found that mCHH levels were decreased in both genic and transposable element (TE) regions, and that TE methylation near genes was decreased as well. Among 142 differentially expressed genes involved in cell division, cellulose biosynthesis, auxin response, growth, and reproduction processes, 75 of them were modified by 122 differentially methylated regions, 10 miRNAs, and 15 siRNAs. In addition, up-regulated PvTOE1 and suppressed PvFT probably contribute to later flowering time of the autopolyploid. The expression changes were probably associated with modification of nearby methylation sites and siRNAs. We also experimentally demonstrated that expression levels of PvFT and PvTOE1 were regulated by DNA methylation, supporting the link between alterations in methylation induced by polyploidization and the phenotypic changes that were observed. Collectively, our results show epigenetic modifications in synthetic autopolyploid switchgrass for the first time, and support the hypothesis that polyploidization-induced methylation is an important cause of phenotypic alterations and is potentially important for plant evolution and improved fitness.

In silico Phylogenetic Analysis of hAT Transposable Elements in Plants

Transposable elements of the hAT family exhibit a cross-kingdom distribution. The plant hAT transposons are proposed to play a critical role in plant adaptive evolution and DNA damage repair. The sequencing of an increasing number of plant genomes has facilitated the discovery of a plethora of hAT elements. This enabled us to perform an in-depth phylogenetic analysis of consensus hAT sequences in the fully-sequenced genomes of 11 plant species that represent diverse taxonomic divisions. Four putative nucleotide sequences were detected in cottonwood that were similar to the corresponding animal hAT elements, which are possibly sequence artifacts. Phylogenetic trees were constructed based both on the known and putative hAT sequences, by employing two different methods of phylogenetic inference. On the basis of the reconstructed phylogeny, plant hAT elements have rather evolved through kingdom-specific vertical gene transfer and gene amplifications within eudicotyledons, monocotyledons, and chlorophytes. Furthermore, the plant hAT sequences were searched for conserved DNA and amino acid sequence features. In this way, diagnostic sequence patterns were detected which allowed us to assign functional annotations to the plant hAT sequences.

On the Search for Retrotransposons: Alternative Protocols to Obtain Sequences to Learn Profile Hidden Markov Models

Profile hidden Markov models (pHMMs) have been used to search for transposable elements (TEs) in genomes. For the learning of pHMMs aimed to search for TEs of the retrotransposon class, the conventional protocol is to use the whole internal nucleotide portions of these elements as representative sequences. To further explore the potential of pHMMs in such a search, we propose five alternative ways to obtain the sets of representative sequences of TEs other than the conventional protocol. In this study, we are interested in Bel-PAO, Copia, Gypsy, and DIRS superfamilies from the retrotransposon class. We compared the pHMMs of all six protocols. The test results show that, for each TE superfamily, the pHMMs of at least two of the proposed protocols performed better than the conventional one and that the number of correct predictions provided by the latter can be improved by considering together the results of one or more of the alternative protocols.

Transcriptional activity of PIF and Pong-like Class II transposable elements in Triticeae

Background

Transposable elements are major contributors to genome size and variability, accounting for approximately 70–80% of the maize, barley, and wheat genomes. PIF and Pong-like elements belong to two closely-related element families within the PIF/Harbinger superfamily of Class II (DNA) transposons. Both elements contain two open reading frames; one encodes a transposase (ORF2) that catalyzes transposition of the functional elements and their related non-autonomous elements, while the function of the second is still debated. In this work, we surveyed for PIF- and Pong-related transcriptional activity in 13 diploid Triticeae species, all of which have been previously shown to harbor extensive within-genome diversity of both groups of elements.

Results

The results revealed that PIF elements have considerable transcriptional activity in Triticeae, suggesting that they can escape the initial levels of plant cell control and are regulated at the post-transcriptional level. Phylogenetic analysis of 156 PIF cDNA transposase fragments along with 240 genomic partial transposase sequences showed that most, if not all, PIF clades are transcriptionally competent, and that multiple transposases coexisting within a single genome have the potential to act simultaneously. In contrast, we did not detect any transcriptional activity of Pong elements in any sample.

Conclusions

The lack of Pong element transcription shows that even closely related transposon families can exhibit wide variation in their transposase transcriptional activity within the same genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1028-6) contains supplementary material, which is available to authorized users.

Electromagnetic Field Seems to Not Influence Transcription via CTCT Motif in Three Plant Promoters

It was proposed that magnetic fields (MFs) can influence gene transcription via CTCT motif located in human HSP70 promoter. To check the universality of this mechanism, we estimated the potential role of this motif on plant gene transcription in response to MFs using both bioinformatics and experimental studies. We searched potential promoter sequences (1000 bp upstream) in the potato Solanum tuberosum and thale cress Arabidopsis thaliana genomes for the CTCT sequence. The motif was found, on average, 3.6 and 4.3 times per promoter (148,487 and 134,361 motifs in total) in these two species, respectively; however, the CTCT sequences were not randomly distributed in the promoter regions but were preferentially located near the transcription initiation site and were closely packed. The closer these CTCT sequences to the transcription initiation site, the smaller distance between them in both plants. One can assume that genes with many CTCT motifs in their promoter regions can be potentially regulated by MFs. To check this assumption, we tested the influence of MFs on gene expression in a transgenic potato with three promoters (16R, 20R, and 5UGT) containing from 3 to 12 CTCT sequences and starting expression of β-glucuronidase as a reported gene. The potatoes were exposed to a 50 Hz 60-70 A/m MF for 30 min and the reporter gene activity was measured for up to 24 h. Although other factors induced the reporter gene activity, the MF did not. It implies the CTCT motif does not mediate in response to MF in the tested plant promoters.

The DNA methylation level against the background of the genome size and t-heterochromatin content in some species of the genus Secale L

Methylation of cytosine in DNA is one of the most important epigenetic modifications in eukaryotes and plays a crucial role in the regulation of gene activity and the maintenance of genomic integrity. DNA methylation and other epigenetic mechanisms affect the development, differentiation or the response of plants to biotic and abiotic stress. This study compared the level of methylation of cytosines on a global (ELISA) and genomic scale (MSAP) between the species of the genus Secale. We analyzed whether the interspecific variation of cytosine methylation was associated with the size of the genome (C-value) and the content of telomeric heterochromatin. MSAP analysis showed that S. sylvestre was the most distinct species among the studied rye taxa; however, the results clearly indicated that these differences were not statistically significant. The total methylation level of the studied loci was very similar in all taxa and ranged from 60% in S. strictum ssp. africanum to 66% in S. cereale ssp. segetale, which confirmed the lack of significant differences in the sequence methylation pattern between the pairs of rye taxa. The level of global cytosine methylation in the DNA was not significantly associated with the content of t-heterochromatin and did not overlap with the existing taxonomic rye relationships. The highest content of 5-methylcytosine was found in S. cereale ssp. segetale (83%), while very low in S. strictum ssp. strictum (53%), which was significantly different from the methylation state of all taxa, except for S. sylvestre. The other studied taxa of rye had a similar level of methylated cytosine ranging from 66.42% (S. vavilovii) to 74.41% in (S. cereale ssp. afghanicum). The results obtained in this study are evidence that the percentage of methylated cytosine cannot be inferred solely based on the genome size or t-heterochromatin. This is a significantly more complex issue.

Activation of Tag1 transposable elements in Arabidopsis dedifferentiating cells and their regulation by CHROMOMETHYLASE 3-mediated CHG methylation

Dedifferentiation, that is, the acquisition of stem cell-like state, commonly induced by stress (e.g., protoplasting), is characterized by open chromatin conformation, a chromatin state that could lead to activation of transposable elements (TEs). Here, we studied the activation of the Arabidopsis class II TE Tag1, in which two copies, situated close to each other (near genes) on chromosome 1 are found in Landsberg erecta (Ler) but not in Columbia (Col). We first transformed protoplasts with a construct in which a truncated Tag1 (ΔTag1 non-autonomous) blocks the expression of a reporter gene AtMBD5-GFP and found a relatively high ectopic excision of ΔTag1 accompanied by expression of AtMBD5-GFP in protoplasts derived from Ler compared to Col; further increase was observed in ddm1 (decrease in DNA methylation1) protoplasts (Ler background). Ectopic excision was associated with transcription of the endogenous Tag1 and changes in histone H3 methylation at the promoter region. Focusing on the endogenous Tag1 elements we found low level of excision in Ler protoplasts, which was slightly and strongly enhanced in ddm1 and cmt3 (chromomethylase3) protoplasts, respectively, concomitantly with reduction in Tag1 gene body (GB) CHG methylation and increased Tag1 transcription; strong activation of Tag1 was also observed in cmt3 leaves. Notably, in cmt3, but not in ddm1, Tag1 elements were excised out from their original sites and transposed elsewhere in the genome. Our results suggest that dedifferentiation is associated with Tag1 activation and that CMT3 rather than DDM1 plays a central role in restraining Tag1 activation via inducing GB CHG methylation.

The plasma membrane NADPH oxidase OsRbohA plays a crucial role in developmental regulation and drought-stress response in rice

Plasma membrane NADPH oxidases are major producers of reactive oxygen species (ROS) in plant cells under normal growth and stress conditions. In the present study the total activity of rice NADPH oxidases and the transcription of OsRbohA, which encodes an Oryza sativa plasma membrane NADPH oxidase, were stimulated by drought. OsRbohA was expressed in all tissues examined throughout development. Its mRNA was upregulated by a number of factors, including heat, drought, salt, oxidative stress and methyl jasmonate treatment. Compared with wild-type (WT), the OsRbohA-knockout mutant osrbohA exhibited upregulated expression of other respiratory burst oxidase homolog genes and multiple abnormal agronomic traits, including reduced biomass, low germination rate and decreased pollen viability and seed fertility. However, OsRbohA-overexpressing transgenic plants showed no differences in these traits compared with WT. Although osrbohA leaves and roots produced more ROS than WT, the mutant had lesser intracellular ROS. In contrast, OsRbohA-overexpressing transgenic plants exhibited higher ROS production at the intracellular level and in tissues. Ablation of OsRbohA impaired the tolerance of plants to various water stresses, whereas its overexpression enhanced the tolerance. In addition, a number of genes related to energy supply, substrate transport, stress response and transcriptional regulation were differentially expressed in osrbohA plants even under normal growth conditions, suggesting that OsRbohA has fundamental and broad functions in rice. These results indicate that OsRbohA-mediated processes are governed by complex signaling pathways that function during the developmental regulation and drought-stress response in rice.

Global DNA cytosine methylation as an evolving trait: phylogenetic signal and correlated evolution with genome size in angiosperms

DNA cytosine methylation is a widespread epigenetic mechanism in eukaryotes, and plant genomes commonly are densely methylated. Genomic methylation can be associated with functional consequences such as mutational events, genomic instability or altered gene expression, but little is known on interspecific variation in global cytosine methylation in plants. In this paper, we compare global cytosine methylation estimates obtained by HPLC and use a phylogenetically-informed analytical approach to test for significance of evolutionary signatures of this trait across 54 angiosperm species in 25 families. We evaluate whether interspecific variation in global cytosine methylation is statistically related to phylogenetic distance and also whether it is evolutionarily correlated with genome size (C-value). Global cytosine methylation varied widely between species, ranging between 5.3% (Arabidopsis) and 39.2% (Narcissus). Differences between species were related to their evolutionary trajectories, as denoted by the strong phylogenetic signal underlying interspecific variation. Global cytosine methylation and genome size were evolutionarily correlated, as revealed by the significant relationship between the corresponding phylogenetically independent contrasts. On average, a ten-fold increase in genome size entailed an increase of about 10% in global cytosine methylation. Results show that global cytosine methylation is an evolving trait in angiosperms whose evolutionary trajectory is significantly linked to changes in genome size, and suggest that the evolutionary implications of epigenetic mechanisms are likely to vary between plant lineages.

Differential genome evolution and speciation of Coix lacryma-jobi L. and Coix aquatica Roxb. hybrid guangxi revealed by repetitive sequence analysis and fine karyotyping

Abstract

Background

Coix, Sorghum and Zea are closely related plant genera in the subtribe Maydeae. Coix comprises 9–11 species with different ploidy levels (2n = 10, 20, 30, and 40). The exclusively cultivated C. lacryma-jobi L. (2n = 20) is widely used in East and Southeast Asia for food and medicinal applications. Three fertile cytotypes (2n = 10, 20, and 40) have been reported for C. aquatica Roxb. One sterile cytotype (2n = 30) closely related to C. aquatica has been recently found in Guangxi of China. This putative hybrid has been named C. aquatica HG (Hybrid Guangxi). The genome composition and the evolutionary history of C. lacryma-jobi and C. aquatica HG are largely unclear.

Results

About 76% of the genome of C. lacryma-jobi and 73% of the genome of C. aquatica HG are repetitive DNA sequences as shown by low coverage genome sequencing followed by similarity-based cluster analysis. In addition, long terminal repeat (LTR) retrotransposable elements are dominant repetitive sequences in these two genomes, and the proportions of many repetitive sequences in whole genome varied greatly between the two species, indicating evolutionary divergence of them. We also found that a novel 102 bp variant of centromeric satellite repeat CentX and two other satellites only appeared in C. aquatica HG. The results from FISH analysis with repeat probe cocktails and the data from chromosomes pairing in meiosis metaphase showed that C. lacryma-jobi is likely a diploidized paleotetraploid species and C. aquatica HG is possibly a recently formed hybrid. Furthermore, C. lacryma-jobi and C. aquatica HG shared more co-existing repeat families and higher sequence similarity with Sorghum than with Zea.

Conclusions

The composition and abundance of repetitive sequences are divergent between the genomes of C. lacryma-jobi and C. aquatica HG. The results from fine karyotyping analysis and chromosome pairing suggested diploidization of C. lacryma-jobi during evolution and C. aquatica HG is a recently formed hybrid. The genome-wide comparison of repetitive sequences indicated that the repeats in Coix were more similar to those in Sorghum than to those in Zea, which is consistent with the phylogenetic relationship reported by previous work.

Electronic supplementary material

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

Transcriptional activity, chromosomal distribution and expression effects of transposable elements in Coffea genomes

Plant genomes are massively invaded by transposable elements (TEs), many of which are located near host genes and can thus impact gene expression. In flowering plants, TE expression can be activated (de-repressed) under certain stressful conditions, both biotic and abiotic, as well as by genome stress caused by hybridization. In this study, we examined the effects of these stress agents on TE expression in two diploid species of coffee, Coffea canephora and C. eugenioides, and their allotetraploid hybrid C. arabica. We also explored the relationship of TE repression mechanisms to host gene regulation via the effects of exonized TE sequences. Similar to what has been seen for other plants, overall TE expression levels are low in Coffea plant cultivars, consistent with the existence of effective TE repression mechanisms. TE expression patterns are highly dynamic across the species and conditions assayed here are unrelated to their classification at the level of TE class or family. In contrast to previous results, cell culture conditions per se do not lead to the de-repression of TE expression in C. arabica. Results obtained here indicate that differing plant drought stress levels relate strongly to TE repression mechanisms. TEs tend to be expressed at significantly higher levels in non-irrigated samples for the drought tolerant cultivars but in drought sensitive cultivars the opposite pattern was shown with irrigated samples showing significantly higher TE expression. Thus, TE genome repression mechanisms may be finely tuned to the ideal growth and/or regulatory conditions of the specific plant cultivars in which they are active. Analysis of TE expression levels in cell culture conditions underscored the importance of nonsense-mediated mRNA decay (NMD) pathways in the repression of Coffea TEs. These same NMD mechanisms can also regulate plant host gene expression via the repression of genes that bear exonized TE sequences.

The Reference Genome of the Halophytic Plant Eutrema salsugineum

Halophytes are plants that can naturally tolerate high concentrations of salt in the soil, and their tolerance to salt stress may occur through various evolutionary and molecular mechanisms. Eutrema salsugineum is a halophytic species in the Brassicaceae that can naturally tolerate multiple types of abiotic stresses that typically limit crop productivity, including extreme salinity and cold. It has been widely used as a laboratorial model for stress biology research in plants. Here, we present the reference genome sequence (241 Mb) of E. salsugineum at 8× coverage sequenced using the traditional Sanger sequencing-based approach with comparison to its close relative Arabidopsis thaliana. The E. salsugineum genome contains 26,531 protein-coding genes and 51.4% of its genome is composed of repetitive sequences that mostly reside in pericentromeric regions. Comparative analyses of the genome structures, protein-coding genes, microRNAs, stress-related pathways, and estimated translation efficiency of proteins between E. salsugineum and A. thaliana suggest that halophyte adaptation to environmental stresses may occur via a global network adjustment of multiple regulatory mechanisms. The E. salsugineum genome provides a resource to identify naturally occurring genetic alterations contributing to the adaptation of halophytic plants to salinity and that might be bioengineered in related crop species.

Promoter variation and transcript divergence in Brassicaceae lineages of FLOWERING LOCUS T

Brassica napus (AACC, 2n = 38), an oil crop of world-wide importance, originated from interspecific hybridization of B. rapa (AA, 2n = 20) and B. oleracea (CC, 2n = 18), and has six FLOWERING LOCUS T (FT) paralogues. Two located on the homeologous chromosomes A2 and C2 arose from a lineage distinct from four located on A7 and C6. A set of three conserved blocks A, B and C, which were found to be essential for FT activation by CONSTANS (CO) in Arabidopsis, was identified within the FT upstream region in B. napus and its progenitor diploids. However, on chromosome C2, insertion of a DNA transposable element (TE) and a retro-element in FT upstream blocks A and B contributed to significant structural divergence between the A and C genome orthologues. Phylogenetic analysis of upstream block A indicated the conserved evolutionary relationships of distinct FT genes within Brassicaceae. We conclude that the ancient At-α whole genome duplication contributed to distinct ancestral lineages for this key adaptive gene, which co-exist within the same genus. FT-A2 was found to be transcribed in all leaf samples from different developmental stages in both B. rapa and B. napus, whereas FT-C2 was not transcribed in either B. napus or B. oleracea. Silencing of FT-C2 appeared to result from TE insertion and consequent high levels of cytosine methylation in TE sequences within upstream block A. Interestingly, FT-A7/C6 paralogues were specifically silenced in winter type B. napus but abundantly expressed in spring type cultivars under vernalization-free conditions. Motif prediction indicated the presence of two CO protein binding sites within all Brassica block A and additional sites for FT activation in block C. We propose that the ancestral whole genome duplications have contributed to more complex mechanisms of floral regulation and niche adaptation in Brassica compared to Arabidopsis.

Propidium iodide for making heterochromatin more evident in the C-banding technique

The detection of regions of heterochromatin has been the subject of intense investigation. We investigated an adaptation of the commonly used technique by replacing the nonfluorescent dye, Giemsa, by a fluorescent one, propidium iodide. This adaptation produces greater contrast of the heterochromatic bands in metaphase chromosomes and can be especially valuable when the organisms studied possess heterochromatin that is pale and difficult to visualize. We discuss the interactions of these two dyes with DNA and the excitation of the fluorescent dye when irradiated with ultraviolet light.

Recent retrotransposon insertions are methylated and phylogenetically clustered in japonica rice (Oryza sativa spp. japonica)

In plants, the genome of the host responds to the amplification of transposable elements (TEs) with DNA methylation. However, neither the factors involved in TE methylation nor the dynamics of the host-TE interaction are well resolved. Here, we identify 5,522 long terminal repeat retrotransposons (LTR-RT) in the genome of Oryza sativa ssp. japonica and then assess methylation for individual elements. Our analyses uncover three strong trends: long LTR-RTs are more highly methylated, the insertion times of LTR-RTs are negatively correlated with methylation, and young LTR-RTs tend to be closer to genes than older insertions. Additionally, a phylogenetic examination of the gypsy-like LTR-RT superfamily revealed that methylation is phylogenetically correlated. Given these observations, we present a model suggesting that the phylogenetic correlation among related LTR-RTs is a primary mechanism driving methylation. In this model, bursts of transposition produce new elements with high sequence similarity. The host machinery identifies proliferating elements as well as closely related LTR-RTs through cross-homology. In addition, our data are consistent with previous hypotheses that methylated LTR-RT elements are removed preferentially from regions near genes, explaining some of the observed age distribution.

Rapid divergence of repetitive DNAs in Brassica relatives

Centromeric, subtelomeric, and telomeric repetitive DNAs were characterized in Brassica species and the related Raphanus sativus and Arabidopsis thaliana. In general, rapid divergence of the repeats was found. The centromeric tandem satellite repeats were differentially distributed in the species studied, suggesting that centromeric repeats have diverged during the evolution of the A/C and B genome lineages. Sequence analysis of centromeric repeats suggested rapid evolution. Pericentromere-associated retrotransposons were identified and showed divergence during the evolution of the lineages as centromeric repeats. A novel subtelomeric tandem repeat from B. nigra was found to be conserved across the diploid Brassica genomes; however, this sequence was not identified in the related species. In contrast to previous studies, interstitial telomere-like repeats were identified in the pericentromeres of Brassica chromosomes, and these repeats may be associated with genomic stability. These results provide insight into genome evolution during polyploidization in Brassica and divergence within the Brassicaceae.

Physical organization of mixed protease inhibitor gene clusters, coordinated expression and association with resistance to late blight at the StKI locus on potato chromosome III

Protease inhibitors (PIs) play a role in plant defence against pests and pathogens as well as in plant development. Potato (Solanum tuberosum) contains abundant levels of diverse PIs. Most potato Kunitz-type inhibitor (KTI) genes map to the StKI locus on potato chromosome III, which is linked to a quantitative trait locus (QTL) for resistance to Phytophthora infestans. To elucidate the physical organization of PIs at the StKI locus, we screened bacterial artificial chromosome (BAC) libraries with KTI probes. Ten different clones were selected, sequenced and annotated. Of 100 putative genes, 22 corresponded to five PI classes. Expression analysis by quantitative real-time PCR (qRT-PCR) using PI class-specific primers in different tissues of the tetraploid potato cultivars 'Nikita' and 'Baltica' revealed different transcript levels, depending on PI type and genotype. During the compatible interaction with a complex race of P. infestans, four PI classes showed coordinated expression over 3 d after infection, a strong decrease in infected leaves and a transient induction in systemic leaves. Basal transcript levels in non-infected leaves differed strongly between the two genotypes examined. Two microsatellite markers located within the PI gene cluster were associated with resistance to P. infestans in a population of potato varieties and breeding clones.

Insight into durum wheat Lpx-B1: a small gene family coding for the lipoxygenase responsible for carotenoid bleaching in mature grains

BACKGROUND

The yellow colour of pasta products is one of the main criteria used by consumers to assess pasta quality. This character is due to the presence of carotenoid pigments in semolina. During pasta processing, oxidative degradation of carotenoid pigments occurs mainly due to lipoxygenase (LOX). In durum wheat (Triticum durum Desf.), two Lpx-1 genes have been identified on chromosome 4B, Lpx-B1.1 and Lpx-B1.2, and evidences have been reported that the deletion of Lpx-B1.1 is associated with a strong reduction in LOX activity in semolina. In the present study, we characterised the Lpx-B1 gene family identified in a durum wheat germplasm collection and related the distribution and expression of the Lpx-B1 genes and alleles to variations in LOX activity in the mature grains.

RESULTS

In addition to the already known Lpx-B1.1 and Lpx-B1.2 genes, a new gene was identified, Lpx-B1.3, along with three different Lpx-B1.1 alleles, Lpx-B1.1a, Lpx-B1.1b and the partially deleted Lpx-B1.1c. Screening of the germplasm collection showed that all of the genotypes have one of the three Lpx-B1.1 alleles, associated with either Lpx-B1.2 or Lpx-B1.3, thus showing that in this collection the two genes are alternatives. Therefore, based on Lpx-B1 distribution, three different haplotypes were distinguished: haplotype I, carrying Lpx-B1.3 and the Lpx-B1.1b allele; haplotype II carrying Lpx-B1.2 and the Lpx-B1.1a allele; and haplotype III carrying Lpx-B1.2 and the Lpx-B1.1c allele. Determination of Lpx-B1 transcript abundance and total LOX activity in mature grains revealed differences among these three haplotypes: haplotypes I, II and III showed high, intermediate and low levels, respectively, of functional Lpx-B1 transcripts and enzymatic activity.

CONCLUSIONS

In this germplasm collection, the Lpx-B1 gene family accounts for most of the total LOX activity in the mature grains. Information on these Lpx-B1 haplotypes provides significant improvement for prediction of LOX-1 activity levels in mature grains, and will therefore help in breeding programmes aimed at selection of new durum wheat genotypes with higher carotenoid contents in their end products.

Genome Size in Diploids, Allopolyploids, and Autopolyploids of Mediterranean Triticeae

Nuclear DNA amount, determined by the flow cytometry method, in diploids, natural and synthetic allopolyploids, and natural and synthetic autopolyploids of the tribe Triticeae (Poaceae) is reviewed here and discussed. In contrast to the very small and nonsignificant variation in nuclear DNA amount that was found at the intraspecific level, the variation at the interspecific level is very large. Evidently changes in genome size are either the cause or the result of speciation. Typical autopolyploids had the expected additive DNA amount of their diploid parents, whereas natural and synthetic cytologically diploidized autopolyploids and natural and synthetic allopolyploids had significantly less DNA than the sum of their parents. Thus, genome downsizing, occurring during or immediately after the formation of these polyploids, provides the physical basis for their cytological diploidization, that is, diploid-like meiotic behavior. Possible mechanisms that are involved in genome downsizing and the biological significance of this phenomenon are discussed.

Tracing the evolution of the floral homeotic B- and C-function genes through genome synteny

The evolution of the floral homeotic genes has been characterized using phylogenetic and functional studies. It is possible to enhance these studies by comparing gene content and order between species to determine the evolutionary history of the regulatory genes. Here, we use a synteny-based approach to trace the evolution of the floral B- and C-function genes that are required for specification of the reproductive organs. Consistent with previous phylogenetic studies, we show that the euAP3–TM6 split occurred after the monocots and dicots diverged. The Arabidopsis TM6 and papaya euAP3 genes are absent from the respective genomes, and we have detected loci from which these genes were lost. These data indicate that either the TM6 or the euAP3 lineage genes can be lost without detriment to flower development. In contrast, PI is essential for male reproductive organ development; yet, contrary to predictions, complex genomic rearrangements have resulted in almost complete breakdown of synteny at the PI locus. In addition to showing the evolution of B-function genes through the prediction of ancestral loci, similar reconstructions reveal the origins of the C-function AG and PLE lineages in dicots, and show the shared ancestry with the monocot C-function genes. During our studies, we found that transposable elements (TEs) present in sequenced Antirrhinum genomic clones limited comparative studies. A pilot survey of the Antirrhinum data revealed that gene-rich regions contain an unusually high degree of TEs of very varied types, which will be an important consideration for future genome sequencing efforts.

Genomic microstructure and differential expression of the genes encoding UDP-glucose:sinapate glucosyltransferase (UGT84A9) in oilseed rape (Brassica napus)

In oilseed rape (Brassica napus), the glucosyltransferase UGT84A9 catalyzes the formation of 1-O-sinapoyl-beta-glucose, which feeds as acyl donor into a broad range of accumulating sinapate esters, including the major antinutritive seed component sinapoylcholine (sinapine). Since down-regulation of UGT84A9 was highly efficient in decreasing the sinapate ester content, the genes encoding this enzyme were considered as potential targets for molecular breeding of low sinapine oilseed rape. B. napus harbors two distinguishable sequence types of the UGT84A9 gene designated as UGT84A9-1 and UGT84A9-2. UGT84A9-1 is the predominantly expressed variant, which is significantly up-regulated during the seed filling phase, when sinapate ester biosynthesis exhibits strongest activity. In the allotetraploid genome of B. napus, UGT84A9-1 is represented by two loci, one derived from the Brassica C-genome (UGT84A9a) and one from the Brassica A-genome (UGT84A9b). Likewise, for UGT84A9-2 two loci were identified in B. napus originating from both diploid ancestor genomes (UGT84A9c, Brassica C-genome; UGT84A9d, Brassica A-genome). The distinct UGT84A9 loci were genetically mapped to linkage groups N15 (UGT84A9a), N05 (UGT84A9b), N11 (UGT84A9c) and N01 (UGT84A9d). All four UGT84A9 genomic loci from B. napus display a remarkably low micro-collinearity with the homologous genomic region of Arabidopsis thaliana chromosome III, but exhibit a high density of transposon-derived sequence elements. Expression patterns indicate that the orthologous genes UGT84A9a and UGT84A9b should be considered for mutagenesis inactivation to introduce the low sinapine trait into oilseed rape.

Sequence-level comparative analysis of the Brassica napus genome around two stearoyl-ACP desaturase loci

We conducted a sequence-level comparative analyses, at the scale of complete bacterial artificial chromosome (BAC) clones, between the genome of the most economically important Brassica species, Brassica napus (oilseed rape), and those of Brassica rapa, the genome of which is currently being sequenced, and Arabidopsis thaliana. We constructed a new B. napus BAC library and identified and sequenced clones that contain homoeologous regions of the genome including stearoyl-ACP desaturase-encoding genes. We sequenced the orthologous region of the genome of B. rapa and conducted comparative analyses between the Brassica sequences and those of the orthologous region of the genome of A. thaliana. The proportion of genes conserved (approximately 56%) is lower than has been reported previously between A. thaliana and Brassica (approximately 66%). The gene models for sets of conserved genes were used to determine the extent of nucleotide conservation of coding regions. This was found to be 84.2 +/- 3.9% and 85.8 +/- 3.7% between the B. napus A and C genomes, respectively, and that of A. thaliana, which is consistent with previous results for other Brassica species, and 97.5 +/- 3.1% between the B. napus A genome and B. rapa, and 93.1 +/- 4.9% between the B. napus C genome and B. rapa. The divergence of the B. napus genes from the A genome and the B. rapa genes was greater than anticipated and indicates that the A genome ancestor of the B. napus cultivar studied was relatively distantly related to the cultivar of B. rapa selected for genome sequencing.

The protist Trichomonas vaginalis harbors multiple lineages of transcriptionally active Mutator-like elements

Background

For three decades the Mutator system was thought to be exclusive of plants, until the first homolog representatives were characterized in fungi and in early-diverging amoebas earlier in this decade.

Results

Here, we describe and characterize four families of Mutator-like elements in a new eukaryotic group, the Parabasalids. These Trichomonas vaginalis Mutator- like elements, or TvMULEs, are active in T. vaginalis and patchily distributed among 12 trichomonad species and isolates. Despite their relatively distinctive amino acid composition, the inclusion of the repeats TvMULE1, TvMULE2, TvMULE3 and TvMULE4 into the Mutator superfamily is justified by sequence, structural and phylogenetic analyses. In addition, we identified three new TvMULE-related sequences in the genome sequence of Candida albicans. While TvMULE1 is a member of the MuDR clade, predominantly from plants, the other three TvMULEs, together with the C. albicans elements, represent a new and quite distinct Mutator lineage, which we named TvCaMULEs. The finding of TvMULE1 sequence inserted into other putative repeat suggests the occurrence a novel TE family not yet described.

Conclusion

These findings expand the taxonomic distribution and the range of functional motif of MULEs among eukaryotes. The characterization of the dynamics of TvMULEs and other transposons in this organism is of particular interest because it is atypical for an asexual species to have such an extreme level of TE activity; this genetic landscape makes an interesting case study for causes and consequences of such activity. Finally, the extreme repetitiveness of the T. vaginalis genome and the remarkable degree of sequence identity within its repeat families highlights this species as an ideal system to characterize new transposable elements.

Change of gene structure and function by non-homologous end-joining, homologous recombination, and transposition of DNA

An important objective in genome research is to relate genome structure to gene function. Sequence comparisons among orthologous and paralogous genes and their allelic variants can reveal sequences of functional significance. Here, we describe a 379-kb region on chromosome 1 of maize that enables us to reconstruct chromosome breakage, transposition, non-homologous end-joining, and homologous recombination events. Such a high-density composition of various mechanisms in a small chromosomal interval exemplifies the evolution of gene regulation and allelic diversity in general. It also illustrates the evolutionary pace of changes in plants, where many of the above mechanisms are of somatic origin. In contrast to animals, somatic alterations can easily be transmitted through meiosis because the germline in plants is contiguous to somatic tissue, permitting the recovery of such chromosomal rearrangements. The analyzed region contains the P1-wr allele, a variant of the genetically well-defined p1 gene, which encodes a Myb-like transcriptional activator in maize. The P1-wr allele consists of eleven nearly perfect P1-wr 12-kb repeats that are arranged in a tandem head-to-tail array. Although a technical challenge to sequence such a structure by shotgun sequencing, we overcame this problem by subcloning each repeat and ordering them based on nucleotide variations. These polymorphisms were also critical for recombination and expression analysis in presence and absence of the trans-acting epigenetic factor Ufo1. Interestingly, chimeras of the p1 and p2 genes, p2/p1 and p1/p2, are framing the P1-wr cluster. Reconstruction of sequence amplification steps at the p locus showed the evolution from a single Myb-homolog to the multi-gene P1-wr cluster. It also demonstrates how non-homologous end-joining can create novel gene fusions. Comparisons to orthologous regions in sorghum and rice also indicate a greater instability of the maize genome, probably due to diploidization following allotetraploidization.

Mutator transposon in maize and MULEs in the plant genome

Mutator (Mu) is by far the most mutagenic plant transposon. The high frequency of transposition and the tendency to insert into low copy sequences for such transposon have made it the primary means by which genes are mutagenized in maize (Zea mays L.). Mus like elements (MULEs) are widespread among angiosperms and multiple-diverged functional variants can be present in a single genome. MULEs often capture genetic sequences. These Pack-MuLEs can mobilize thousands of gene fragments, which may have had a significant impact on host genome evolution. There is also evidence that MULEs can move between reproductively isolated species. Here we present an overview of the discovery, features and utility of Mu transposon. Classification of Mu elements and future directions of related research are also discussed. Understanding Mu will help us elucidate the dynamic genome.

Progress in plant CACTA elements

Transposable elements are DNA fragments that can insert new chromosomal locations. On the basis of the mechanism of transposition, transposable elements were divided into two classes. Class 1 elements were retroelements that used reverse transposase to transpose by an RNA intermediate. Class 2 elements or DNA transposons transposed directly from DNA to DNA. Of the Class 2 elements, CACTA superfamily, so far identified exclusively in plants and previously regarded as low-copy-transposon for the conserved mechanism of propagation, recently received considerable interest because of their increasing evidence reiterating their high copies in some plant genomes. This article aimed at outlining CACTA elements with regard to their structure, transposition, and utilization.

Identification of differentially expressed proteins between hybrid and parents in wheat (Triticum aestivum L.) seedling leaves

In spite of commercial use of heterosis in agriculture, the molecular basis of heterosis is poorly understood. To gain a better understanding of the molecular basis of wheat heterosis, we carried out a comparative proteomic analysis in seedling leaves between wheat hybrid and parents. Common wheat (Triticum aestivum L., 2n = 6x = 42, AABBDD) Line 3338 and spelt wheat (Triticum spelta L., 2n = 6x = 42, AABBDD) Line 2463 were used to produce a heterotic F(1) hybrid. The expression patterns of the total proteins were compared in seedling leaves between hybrid and its parents by using two-dimensional gel electrophoresis with two pH ranges for the first dimension separation. Among ~900 protein spots reproducibly detected, 49 protein spots were identified as being differentially expressed between hybrid and its parental lines (P < 0.05) for more than 1.5-folds. Six possible modes of differential expression were observed, including high- and low-parent dominance, underdominance, and overdominance, uniparent silencing and uniparent dominance. Moreover, 30 of the 49 differentially expressed protein spots were identified, which were involved in metabolism, signal transduction, energy, cell growth and division, disease and defense, secondary metabolism. These results indicated that wheat hybridization can cause protein expression differences between hybrid and its parents; these proteins were involved in diverse physiological process pathways, which might be responsible for the observed heterosis.

A transcriptionally active copia-like retroelement in Citrus limon

The plant nuclear genome is largely composed of mobile DNA, which can rearrange genomes and other individual gene structure and also affect gene regulation through various promoted activities: transposition, insertion, excision, chromosome breakage, and ectopic recombination. Ty1-copia-like retrotransposon is a widespread class of transposable elements in the plant kingdom, representing a large part of the total DNA content. Here, a novel retrotransposon-like sequence was isolated and identified as the Ty1-copia-like reverse transcriptase domain (named here CLCoy1), based on the homology of known elements. Fluorescence in situ hybridization, revealed that CLCoy1 was mainly located in telomeric and sub-telomeric regions along the Citrus chromosomes. CLCoy1 composes 3.6% of the genome and, interestingly, while transposons are mostly specific to a species, this element was identified in other Citrus species such as Citrus aurantium, Fortunella margarita and Citrus paradisi, but undetected in Poncirus trifoliata. We also determined that wounding, salt and cell culture stress produced transcriptional activation of this novel retroelement in Citrus limon. The novel Ty1-copia-like element CLCoy1 may have played a major role in shaping genome structure and size during Citrus species evolution.

The CACTA transposon Bot1 played a major role in Brassica genome divergence and gene proliferation

We isolated and characterized a Brassica C genome-specific CACTA element, which was designated Bot1 (Brassica oleracea transposon 1). After analysing phylogenetic relationships, copy numbers and sequence similarity of Bot1 and Bot1 analogues in B. oleracea (C genome) versus Brassica rapa (A genome), we concluded that Bot1 has encountered several rounds of amplification in the oleracea genome only, and has played a major role in the recent rapa and oleracea genome divergence. We performed in silico analyses of the genomic organization and internal structure of Bot1, and established which segment of Bot1 is C-genome specific. Our work reports a fully characterized Brassica repetitive sequence that can distinguish the Brassica A and C chromosomes in the allotetraploid Brassica napus, by fluorescent in situ hybridization. We demonstrated that Bot1 carries a host S locus-associated SLL3 gene copy. We speculate that Bot1 was involved in the proliferation of SLL3 around the Brassica genome. The present study reinforces the assumption that transposons are a major driver of genome and gene evolution in higher plants.

Genome size and genome evolution in diploid Triticeae species

One of the intriguing issues concerning the dynamics of plant genomes is the occurrence of intraspecific variation in nuclear DNA amount. The aim of this work was to assess the ranges of intraspecific, interspecific, and intergeneric variation in nuclear DNA content of diploid species of the tribe Triticeae (Poaceae) and to examine the relation between life form or habitat and genome size. Altogether, 438 plants representing 272 lines that belong to 22 species were analyzed. Nuclear DNA content was estimated by flow cytometry. Very small intraspecific variation in DNA amount was found between lines of Triticeae diploid species collected from different habitats or between different morphs. In contrast to the constancy in nuclear DNA amount at the intraspecific level, there are significant differences in genome size between the various diploid species. Within the genus Aegilops, the 1C DNA amount ranged from 4.84 pg in A. caudata to 7.52 pg in A. sharonensis; among genera, the 1C DNA amount ranged from 4.18 pg in Heteranthelium piliferum to 9.45 pg in Secale montanum. No evidence was found for a smaller genome size in annual, self-pollinating species relative to perennial, cross-pollinating ones. Diploids that grow in the southern part of the group's distribution have larger genomes than those growing in other parts of the distribution. The contrast between the low variation at the intraspecific level and the high variation at the interspecific one suggests that changes in genome size originated in close temporal proximity to the speciation event, i.e., before, during, or immediately after it. The possible effects of sudden changes in genome size on speciation processes are discussed.