Plant genome organisation and diversity: the year of the junk!

Identification of a genome-specific repetitive element in the Gossypium D genome

The activity of genome-specific repetitive sequences is the main cause of genome variation between Gossypium A and D genomes. Through comparative analysis of the two genomes, we retrieved a repetitive element termed ICRd motif, which appears frequently in the diploid Gossypium raimondii (D₅) genome but rarely in the diploid Gossypium arboreum (A₂) genome. We further explored the existence of the ICRd motif in chromosomes of G. raimondii, G. arboreum, and two tetraploid (AADD) cotton species, Gossypium hirsutum and Gossypium barbadense, by fluorescence in situ hybridization (FISH), and observed that the ICRd motif exists in the D₅ and D-subgenomes but not in the A₂ and A-subgenomes. The ICRd motif comprises two components, a variable tandem repeat (TR) region and a conservative sequence (CS). The two constituents each have hundreds of repeats that evenly distribute across 13 chromosomes of the D₅genome. The ICRd motif (and its repeats) was revealed as the common conservative region harbored by ancient Long Terminal Repeat Retrotransposons. Identification and investigation of the ICRd motif promotes the study of A and D genome differences, facilitates research on Gossypium genome evolution, and provides assistance to subgenome identification and genome assembling.

Identification of a genome-specific repetitive element in the Gossypium D genome

The activity of genome-specific repetitive sequences is the main cause of genome variation between Gossypium A and D genomes. Through comparative analysis of the two genomes, we retrieved a repetitive element termed ICRd motif, which appears frequently in the diploid Gossypium raimondii (D₅) genome but rarely in the diploid Gossypium arboreum (A₂) genome. We further explored the existence of the ICRd motif in chromosomes of G. raimondii, G. arboreum, and two tetraploid (AADD) cotton species, Gossypium hirsutum and Gossypium barbadense, by fluorescence in situ hybridization (FISH), and observed that the ICRd motif exists in the D₅ and D-subgenomes but not in the A₂ and A-subgenomes. The ICRd motif comprises two components, a variable tandem repeat (TR) region and a conservative sequence (CS). The two constituents each have hundreds of repeats that evenly distribute across 13 chromosomes of the D₅genome. The ICRd motif (and its repeats) was revealed as the common conservative region harbored by ancient Long Terminal Repeat Retrotransposons. Identification and investigation of the ICRd motif promotes the study of A and D genome differences, facilitates research on Gossypium genome evolution, and provides assistance to subgenome identification and genome assembling.

Rolling-circle transposons catalyze genomic innovation in a mammalian lineage

Rolling-circle transposons (Helitrons) are a newly discovered group of mobile DNA widespread in plant and invertebrate genomes but limited to the bat family Vespertilionidae among mammals. Little is known about the long-term impact of Helitron activity because the genomes where Helitron activity has been extensively studied are predominated by young families. Here, we report a comprehensive catalog of vetted Helitrons from the 7× Myotis lucifugus genome assembly. To estimate the timing of transposition, we scored presence/absence across related vespertilionid genome sequences with estimated divergence times. This analysis revealed that the Helibat family has been a persistent source of genomic innovation throughout the vespertilionid diversification from approximately 30–36 Ma to as recently as approximately 1.8–6 Ma. This is the first report of persistent Helitron transposition over an extended evolutionary timeframe. These findings illustrate that the pattern of Helitron activity is akin to the vertical persistence of LINE retrotransposons in primates and other mammalian lineages. Like retrotransposition in primates, rolling-circle transposition has generated lineage-specific variation and accounts for approximately 110 Mb, approximately 6% of the genome of M. lucifugus. The Helitrons carry a heterogeneous assortment of host sequence including retroposed messenger RNAs, retrotransposons, DNA transposons, as well as introns, exons and regulatory regions (promoters, 5′-untranslated regions [UTRs], and 3′-UTRs) of which some are evolving in a pattern suggestive of purifying selection. Evidence that Helitrons have contributed putative promoters, exons, splice sites, polyadenylation sites, and microRNA-binding sites to transcripts otherwise conserved across mammals is presented, and the implication of Helitron activity to innovation in these unique mammals is discussed.

Scanning of transposable elements and analyzing expression of transposase genes of sweet potato [Ipomoea batatas]

BACKGROUND

Transposable elements (TEs) are the most abundant genomic components in eukaryotes and affect the genome by their replications and movements to generate genetic plasticity. Sweet potato performs asexual reproduction generally and the TEs may be an important genetic factor for genome reorganization. Complete identification of TEs is essential for the study of genome evolution. However, the TEs of sweet potato are still poorly understood because of its complex hexaploid genome and difficulty in genome sequencing. The recent availability of the sweet potato transcriptome databases provides an opportunity for discovering and characterizing the expressed TEs.

METHODOLOGY/PRINCIPAL FINDINGS

We first established the integrated-transcriptome database by de novo assembling four published sweet potato transcriptome databases from three cultivars in China. Using sequence-similarity search and analysis, a total of 1,405 TEs including 883 retrotransposons and 522 DNA transposons were predicted and categorized. Depending on mapping sets of RNA-Seq raw short reads to the predicted TEs, we compared the quantities, classifications and expression activities of TEs inter- and intra-cultivars. Moreover, the differential expressions of TEs in seven tissues of Xushu 18 cultivar were analyzed by using Illumina digital gene expression (DGE) tag profiling. It was found that 417 TEs were expressed in one or more tissues and 107 in all seven tissues. Furthermore, the copy number of 11 transposase genes was determined to be 1-3 copies in the genome of sweet potato by Real-time PCR-based absolute quantification.

CONCLUSIONS/SIGNIFICANCE

Our result provides a new method for TE searching on species with transcriptome sequences while lacking genome information. The searching, identification and expression analysis of TEs will provide useful TE information in sweet potato, which are valuable for the further studies of TE-mediated gene mutation and optimization in asexual reproduction. It contributes to elucidating the roles of TEs in genome evolution.

Identification of multiple binding sites for the THAP domain of the Galileo transposase in the long terminal inverted-repeats

Galileo is a DNA transposon responsible for the generation of several chromosomal inversions in Drosophila. In contrast to other members of the P-element superfamily, it has unusually long terminal inverted-repeats (TIRs) that resemble those of Foldback elements. To investigate the function of the long TIRs we derived consensus and ancestral sequences for the Galileo transposase in three species of Drosophilids. Following gene synthesis, we expressed and purified their constituent THAP domains and tested their binding activity towards the respective Galileo TIRs. DNase I footprinting located the most proximal DNA binding site about 70 bp from the transposon end. Using this sequence we identified further binding sites in the tandem repeats that are found within the long TIRs. This suggests that the synaptic complex between Galileo ends may be a complicated structure containing higher-order multimers of the transposase. We also attempted to reconstitute Galileo transposition in Drosophila embryos but no events were detected. Thus, although the limited numbers of Galileo copies in each genome were sufficient to provide functional consensus sequences for the THAP domains, they do not specify a fully active transposase. Since the THAP recognition sequence is short, and will occur many times in a large genome, it seems likely that the multiple binding sites within the long, internally repetitive, TIRs of Galileo and other Foldback-like elements may provide the transposase with its binding specificity.

Highly diverse chromoviruses of Beta vulgaris are classified by chromodomains and chromosomal integration

Background

Chromoviruses are one of the three genera of Ty3-gypsy long terminal repeat (LTR) retrotransposons, and are present in high copy numbers in plant genomes. They are widely distributed within the plant kingdom, with representatives even in lower plants such as green and red algae. Their hallmark is the presence of a chromodomain at the C-terminus of the integrase. The chromodomain exhibits structural characteristics similar to proteins of the heterochromatin protein 1 (HP1) family, which mediate the binding of each chromovirus type to specific histone variants. A specific integration via the chromodomain has been shown for only a few chromoviruses. However, a detailed study of different chromoviral clades populating a single plant genome has not yet been carried out.

Results

We conducted a comprehensive survey of chromoviruses within the Beta vulgaris (sugar beet) genome, and found a highly diverse chromovirus population, with significant differences in element size, primarily caused by their flanking LTRs. In total, we identified and annotated full-length members of 16 families belonging to the four plant chromoviral clades: CRM, Tekay, Reina, and Galadriel. The families within each clade are structurally highly conserved; in particular, the position of the chromodomain coding region relative to the polypurine tract is clade-specific. Two distinct groups of chromodomains were identified. The group II chromodomain was present in three chromoviral clades, whereas families of the CRM clade contained a more divergent motif. Physical mapping using representatives of all four clades identified a clade-specific integration pattern. For some chromoviral families, we detected the presence of expressed sequence tags, indicating transcriptional activity.

Conclusions

We present a detailed study of chromoviruses, belonging to the four major clades, which populate a single plant genome. Our results illustrate the diversity and family structure of B. vulgaris chromoviruses, and emphasize the role of chromodomains in the targeted integration of these viruses. We suggest that the diverse sets of plant chromoviruses with their different localization patterns might help to facilitate plant-genome organization in a structural and functional manner.

Characterization of transposable elements in the ectomycorrhizal fungus Laccaria bicolor

Background

The publicly available Laccaria bicolor genome sequence has provided a considerable genomic resource allowing systematic identification of transposable elements (TEs) in this symbiotic ectomycorrhizal fungus. Using a TE-specific annotation pipeline we have characterized and analyzed TEs in the L. bicolor S238N-H82 genome.

Methodology/Principal Findings

TEs occupy 24% of the 60 Mb L. bicolor genome and represent 25,787 full-length and partial copy elements distributed within 171 families. The most abundant elements were the Copia-like. TEs are not randomly distributed across the genome, but are tightly nested or clustered. The majority of TEs exhibits signs of ancient transposition except some intact copies of terminal inverted repeats (TIRS), long terminal repeats (LTRs) and a large retrotransposon derivative (LARD) element. There were three main periods of TE expansion in L. bicolor: the first from 57 to 10 Mya, the second from 5 to 1 Mya and the most recent from 0.5 Mya ago until now. LTR retrotransposons are closely related to retrotransposons found in another basidiomycete, Coprinopsis cinerea.

Conclusions

This analysis 1) represents an initial characterization of TEs in the L. bicolor genome, 2) contributes to improve genome annotation and a greater understanding of the role TEs played in genome organization and evolution and 3) provides a valuable resource for future research on the genome evolution within the Laccaria genus.

Involvement of disperse repetitive sequences in wheat/rye genome adjustment

The union of different genomes in the same nucleus frequently results in hybrid genotypes with improved genome plasticity related to both genome remodeling events and changes in gene expression. Most modern cereal crops are polyploid species. Triticale, synthesized by the cross between wheat and rye, constitutes an excellent model to study polyploidization functional implications. We intend to attain a deeper knowledge of dispersed repetitive sequence involvement in parental genome reshuffle in triticale and in wheat-rye addition lines that have the entire wheat genome plus each rye chromosome pair. Through Random Amplified Polymorphic DNA (RAPD) analysis with OPH20 10-mer primer we unraveled clear alterations corresponding to the loss of specific bands from both parental genomes. Moreover, the sequential nature of those events was revealed by the increased absence of rye-origin bands in wheat-rye addition lines in comparison with triticale. Remodeled band sequencing revealed that both repetitive and coding genome domains are affected in wheat-rye hybrid genotypes. Additionally, the amplification and sequencing of pSc20H internal segments showed that the disappearance of parental bands may result from restricted sequence alterations and unraveled the involvement of wheat/rye related repetitive sequences in genome adjustment needed for hybrid plant stabilization.

Emerging knowledge from genome sequencing of crop species

Extensive insights into the genome composition, organization, and evolution have been gained from the plant genome sequencing and annotation ongoing projects. The analysis of crop genomes provided surprising evidences with important implications in plant origin and evolution: genome duplication, ancestral re-arrangements and unexpected polyploidization events opened new doors to address fundamental questions related to species proliferation, adaptation, and functional modulations. Detailed paleogenomic analysis led to many speculation on how chromosomes have been shaped over time in terms of gene content and order. The completion of the genome sequences of several major crops, prompted to a detailed identification and annotation of transposable elements: new hypothesis related to their composition, chromosomal distribution, insertion models, amplification rate, and evolution patterns are coming up. Availability of full genome sequence of several crop species as well as from many accessions within species is providing new keys for biodiversity exploitation and interpretation. Re-sequencing is enabling high-throughput genotyping to identify a wealth of SNP and afterward to produce haplotype maps necessary to accurately associate molecular variation to phenotype. Conservation genomics is emerging as a powerful tool to explain adaptation, genetic drift, natural selection, hybridization and to estimate genetic variation, fitness and population's viability.

MASiVEdb: the Sirevirus Plant Retrotransposon Database

Background

Sireviruses are an ancient genus of the Copia superfamily of LTR retrotransposons, and the only one that has exclusively proliferated within plant genomes. Based on experimental data and phylogenetic analyses, Sireviruses have successfully infiltrated many branches of the plant kingdom, extensively colonizing the genomes of grass species. Notably, it was recently shown that they have been a major force in the make-up and evolution of the maize genome, where they currently occupy ~21% of the nuclear content and ~90% of the Copia population. It is highly likely, therefore, that their life dynamics have been fundamental in the genome composition and organization of a plethora of plant hosts. To assist studies into their impact on plant genome evolution and also facilitate accurate identification and annotation of transposable elements in sequencing projects, we developed MASiVEdb (Mapping and Analysis of SireVirus Elements Database), a collective and systematic resource of Sireviruses in plants.

Description

Taking advantage of the increasing availability of plant genomic sequences, and using an updated version of MASiVE, an algorithm specifically designed to identify Sireviruses based on their highly conserved genome structure, we populated MASiVEdb (http://bat.infspire.org/databases/masivedb/) with data on 16,243 intact Sireviruses (total length >158Mb) discovered in 11 fully-sequenced plant genomes. MASiVEdb is unlike any other transposable element database, providing a multitude of highly curated and detailed information on a specific genus across its hosts, such as complete set of coordinates, insertion age, and an analytical breakdown of the structure and gene complement of each element. All data are readily available through basic and advanced query interfaces, batch retrieval, and downloadable files. A purpose-built system is also offered for detecting and visualizing similarity between user sequences and Sireviruses, as well as for coding domain discovery and phylogenetic analysis.

Conclusion

MASiVEdb is currently the most comprehensive directory of Sireviruses, and as such complements other efforts in cataloguing plant transposable elements and elucidating their role in host genome evolution. Such insights will gradually deepen, as we plan to further improve MASiVEdb by phylogenetically mapping Sireviruses into families, by including data on fragments and solo LTRs, and by incorporating elements from newly-released genomes.

Gene space dynamics during the evolution of Aegilops tauschii, Brachypodium distachyon, Oryza sativa, and Sorghum bicolor genomes

Nine different regions totaling 9.7 Mb of the 4.02 Gb Aegilops tauschii genome were sequenced using the Sanger sequencing technology and compared with orthologous Brachypodium distachyon, Oryza sativa (rice), and Sorghum bicolor (sorghum) genomic sequences. The ancestral gene content in these regions was inferred and used to estimate gene deletion and gene duplication rates along each branch of the phylogenetic tree relating the four species. The total gene number in the extant Ae. tauschii genome was estimated to be 36,371. The gene deletion and gene duplication rates and total gene numbers in the four genomes were used to estimate the total gene number in each node of the phylogenetic tree. The common ancestor of the Brachypodieae and Triticeae lineages was estimated to have had 28,558 genes, and the common ancestor of the Panicoideae, Ehrhartoideae, and Pooideae subfamilies was estimated to have had 27,152 or 28,350 genes, depending on the ancestral gene scenario. Relative to the Brachypodieae and Triticeae common ancestor, the gene number was reduced in B. distachyon by 3,026 genes and increased in Ae. tauschii by 7,813 genes. The sum of gene deletion and gene duplication rates, which reflects the rate of gene synteny loss, was correlated with the rate of structural chromosome rearrangements and was highest in the Ae. tauschii lineage and lowest in the rice lineage. The high rate of gene space evolution in the Ae. tauschii lineage accounts for the fact that, contrary to the expectations, the level of synteny between the phylogenetically more related Ae. tauschii and B. distachyon genomes is similar to the level of synteny between the Ae. tauschii genome and the genomes of the less related rice and sorghum. The ratio of gene duplication to gene deletion rates in these four grass species closely parallels both the total number of genes in a species and the overall genome size. Because the overall genome size is to a large extent a function of the repeated sequence content in a genome, we suggest that the amount and activity of repeated sequences are important factors determining the number of genes in a genome.

Transpositional landscape of the rice genome revealed by paired-end mapping of high-throughput re-sequencing data

Transposable elements (TEs) are mobile entities that densely populate most eukaryotic genomes and contribute to both their structural and functional dynamics. However, most TE-related sequences in both plant and animal genomes correspond to inactive, degenerated elements, due to the combined effect of silencing pathways and elimination through deletions. One of the major difficulties in fully characterizing the molecular basis of genetic diversity of a given species lies in establishing its genome-wide transpositional activity. Here, we provide an extensive survey of the transpositional landscape of a plant genome using a deep sequencing strategy. This was achieved through paired-end mapping of a fourfold coverage of the genome of rice mutant line derived from an in vitro callus culture using Illumina technology. Our study shows that at least 13 TE families are active in this genotype, causing 34 new insertions. This next-generation sequencing-based strategy provides new opportunities to quantify the impact of TEs on the genome dynamics of the species.

Are we Genomic Mosaics? Variations of the Genome of Somatic Cells can Contribute to Diversify our Phenotypes

Theoretical and experimental evidences support the hypothesis that the genomes and the epigenomes may be different in the somatic cells of complex organisms. In the genome, the differences range from single base substitutions to chromosome number; in the epigenome, they entail multiple postsynthetic modifications of the chromatin. Somatic genome variations (SGV) may accumulate during development in response both to genetic programs, which may differ from tissue to tissue, and to environmental stimuli, which are often undetected and generally irreproducible. SGV may jeopardize physiological cellular functions, but also create novel coding and regulatory sequences, to be exposed to intraorganismal Darwinian selection. Genomes acknowledged as comparatively poor in genes, such as humans', could thus increase their pristine informational endowment. A better understanding of SGV will contribute to basic issues such as the "nature vs nurture" dualism and the inheritance of acquired characters. On the applied side, they may explain the low yield of cloning via somatic cell nuclear transfer, provide clues to some of the problems associated with transdifferentiation, and interfere with individual DNA analysis. SGV may be unique in the different cells types and in the different developmental stages, and thus explain the several hundred gaps persisting in the human genomes "completed" so far. They may compound the variations associated to our epigenomes and make of each of us an "(epi)genomic" mosaic. An ensuing paradigm is the possibility that a single genome (the ephemeral one assembled at fertilization) has the capacity to generate several different brains in response to different environments.

Codon usage biases of transposable elements and host nuclear genes in Arabidopsis thaliana and Oryza sativa

Transposable elements (TEs) are mobile genetic entities ubiquitously distributed in nearly all genomes. High frequency of codons ending in A/T in TEs has been previously observed in some species. In this study, the biases in nucleotide composition and codon usage of TE transposases and host nuclear genes were investigated in the AT-rich genome of Arabidopsis thaliana and the GC-rich genome of Oryza sativa. Codons ending in A/T are more frequently used by TEs compared with their host nuclear genes. A remarkable positive correlation between highly expressed nuclear genes and C/G-ending codons were detected in O. sativa (r=0.944 and 0.839, respectively, P<0.0001) but not in A. thaliana, indicating a close association between the GC content and gene expression level in monocot species. In both species, TE codon usage biases are similar to that of weakly expressed genes. The expression and activity of TEs may be strictly controlled in plant genomes. Mutation bias and selection pressure have simultaneously acted on the TE evolution in A. thaliana and O. sativa. The consistently observed biases of nucleotide composition and codon usage of TEs may also provide a useful clue to accurately detect TE sequences in different species.

Absence of close-facing retrotransposons: a comparison of molecular data and theory

Retrotransposons occur in extremely large numbers in many eukaryotic genomes. However, little is known of the factors which affect the distribution of close proximity elements. In this work we investigate the frequency of close facing retrotransposons in a plant species with extremely high numbers of retrotransposons. Molecular observations are compared with predictions of a mathematical model that assumes a uniform probability of retrotransposon insertion into the genome. The mathematical model plays the role of a null hypothesis. We find that compared with the predictions of the model, there is a statistically significant deficit of identical copies of facing retroelements that are close to one another. This suggests that an efficient mechanism exists that removes or limits close facing retroelements.

Dehydrin variants associated with superior freezing tolerance in alfalfa (Medicago sativa L.)

A cDNA (msaCIG) encoding a cold-inducible Y(2)K(4) dehydrin in alfalfa (Medicago sativa spp. sativa) was shown to share extensive homology with sequences from other species and subspecies of Medicago. Differences were mainly the result of the occurrence of large indels, amino acids substitutions/deletions and sequence duplications. Using a combination of a bulk segregant analysis and RFLP hybridization, we uncovered an msaCIG polymorphism that increases in frequency in response to recurrent selection for superior freezing tolerance. Progenies from crosses between genotypes with (D+) or without (D-) the polymorphic dehydrin significantly differed in their tolerance to subfreezing temperatures. Based on the msaCIG sequence, we looked for intragenic variations that could be associated to the polymorphism detected on Southern blots. Amplifications with primers targeting the 3' half side of msaCIG revealed fragment size variations between pools of genotypes with (+) or without (-) the polymorphism. Three major groups of amplicons of approximately 370 nt (G1), 330 nt (G2), and 290 nt (G3) were distinguished. The G2 group was more intensively amplified in pools of genotypes with the polymorphic dehydrin and was associated to a superior freezing tolerance phenotype. Sequences analysis revealed that size variation in the 3' half was attributable to the variable occurrence of large indels. Single amino acid substitutions and/or deletions caused major differences in the prediction of the secondary structure of the polypeptides. The identification of dehydrin variants associated to superior freezing tolerance paves the way to the development of functional markers and the fixation of favorable alleles in various genetic backgrounds.

Molecular characterization of the Sasanda LTR copia retrotransposon family uncovers their recent amplification in Triticum aestivum (L.) genome

Retrotransposons constitute a major proportion of the Triticeae genomes. Genome-scale studies have revealed their role in evolution affecting both genome structure and function and their potential for the development of novel markers. In this study, family members of an LTR copia retrotransposon which mediated the duplication of the gene encoding the high molecular weight glutenin subunit Bx7 in cultivar Glenlea were characterized. This novel element was named Sasanda_EU157184-1 (TREP3516). High density filters of the Glenlea hexaploid wheat BAC library were screened with a Sasanda long terminal repeat (LTR)-specific probe and approximately 1,075 positive clones representing an estimated copy number of 347 elements per haploid genome were identified. The 242 BAC clones with the strongest hybridization signal were selected. To maximize isolation of complete elements, this subset of clones was screened with a reverse transcriptase (RT) domain probe and DNA was isolated from the 133 clones that produced a strong hybridization signal. Left (5') and right (3') LTRs as well as the RT domains were PCR amplified and sequencing was carried out on the final subset of 121 clones. Evolutionary relationships were inferred from a data set consisting of 100 RT, 102 5' LTR and 100 3' LTR sequences representing 233, 451 and 495 informative sites for comparison, respectively. Neighbour-joining tree indicated that the element is at least 1.8 million years old and has evolved into a minimum of five sub-families. The insertion times of the 89 complete elements were estimated based on the divergence between their LTRs. Corroborating the inference from the RT domain, analysis of the LTR domains also indicated bursts of amplification from 2.6 million years ago (MYA) to now, except for one member dated to 4.6 +/- 0.7 MYA, which corresponds to the interval of divergence of Triticum and Aegilops (3 MYA) and divergence of Triticum and Rye (7 MYA). In 44 elements, the 5' and 3' LTRs were identical indicating recent transposition activity. The element can be used to develop retrotransposon-based markers such as sequence-specific amplified polymorphism, retrotransposon microsatellite amplified polymorphism and inter-retrotransposon amplified polymorphism, all of which are well suited for genotyping studies.

The contribution of recombination to heterozygosity differs among plant evolutionary lineages and life-forms

BACKGROUND

Despite its role as a generator of haplotypic variation, little is known about how the rates of recombination evolve across taxa. Recombination is a very labile force, susceptible to evolutionary and life trait related processes, which have also been correlated with general levels of genetic diversity. For example, in plants, it has been shown that long-lived outcrossing taxa, such as trees, have higher heterozygosity (He) at SSRs and allozymes than selfing or annual species. However, some of these tree taxa have surprisingly low levels of nucleotide diversity at the DNA sequence level, which points to recombination as a potential generator of genetic diversity in these organisms. In this study, we examine how genome-wide and within-gene rates of recombination evolve across plant taxa, determine whether such rates are influenced by the life-form adopted by species, and evaluate if higher genome-wide rates of recombination translate into higher He values, especially in trees.

RESULTS

Estimates of genome-wide (cM/Mb) recombination rates from 81 higher plants showed a significant phylogenetic signal. The use of different comparative phylogenetic models demonstrated that there is a positive correlation between recombination rate and He (0.83 +/- 0.29), and that trees have higher rates of genome-wide recombination than short-lived herbs and shrubs. A significant taxonomic component was further made evident by our models, as conifers exhibited lower recombination rates than angiosperms. This trend was also found at the within-gene level.

CONCLUSIONS

Altogether, our results illustrate how both common ancestry and life-history traits have to be taken into account for understanding the evolution of genetic diversity and genomic rates of recombination across plant species, and highlight the relevance of species life forms to explain general levels of diversity and recombination.

Gene expression profiles deciphering rice phenotypic variation between Nipponbare (Japonica) and 93-11 (Indica) during oxidative stress

Rice is a very important food staple that feeds more than half the world's population. Two major Asian cultivated rice (Oryza sativa L.) subspecies, japonica and indica, show significant phenotypic variation in their stress responses. However, the molecular mechanisms underlying this phenotypic variation are still largely unknown. A common link among different stresses is that they produce an oxidative burst and result in an increase of reactive oxygen species (ROS). In this study, methyl viologen (MV) as a ROS agent was applied to investigate the rice oxidative stress response. We observed that 93-11 (indica) seedlings exhibited leaf senescence with severe lesions under MV treatment compared to Nipponbare (japonica). Whole-genome microarray experiments were conducted, and 1,062 probe sets were identified with gene expression level polymorphisms between the two rice cultivars in addition to differential expression under MV treatment, which were assigned as Core Intersectional Probesets (CIPs). These CIPs were analyzed by gene ontology (GO) and highlighted with enrichment GO terms related to toxin and oxidative stress responses as well as other responses. These GO term-enriched genes of the CIPs include glutathine S-transferases (GSTs), P450, plant defense genes, and secondary metabolism related genes such as chalcone synthase (CHS). Further insertion/deletion (InDel) and regulatory element analyses for these identified CIPs suggested that there may be some eQTL hotspots related to oxidative stress in the rice genome, such as GST genes encoded on chromosome 10. In addition, we identified a group of marker genes individuating the japonica and indica subspecies. In summary, we developed a new strategy combining biological experiments and data mining to study the possible molecular mechanism of phenotypic variation during oxidative stress between Nipponbare and 93-11. This study will aid in the analysis of the molecular basis of quantitative traits.

Comparative analysis between homoeologous genome segments of Brassica napus and its progenitor species reveals extensive sequence-level divergence

Homoeologous regions of Brassica genomes were analyzed at the sequence level. These represent segments of the Brassica A genome as found in Brassica rapa and Brassica napus and the corresponding segments of the Brassica C genome as found in Brassica oleracea and B. napus. Analysis of synonymous base substitution rates within modeled genes revealed a relatively broad range of times (0.12 to 1.37 million years ago) since the divergence of orthologous genome segments as represented in B. napus and the diploid species. Similar, and consistent, ranges were also identified for single nucleotide polymorphism and insertion-deletion variation. Genes conserved across the Brassica genomes and the homoeologous segments of the genome of Arabidopsis thaliana showed almost perfect collinearity. Numerous examples of apparent transduplication of gene fragments, as previously reported in B. oleracea, were observed in B. rapa and B. napus, indicating that this phenomenon is widespread in Brassica species. In the majority of the regions studied, the C genome segments were expanded in size relative to their A genome counterparts. The considerable variation that we observed, even between the different versions of the same Brassica genome, for gene fragments and annotated putative genes suggest that the concept of the pan-genome might be particularly appropriate when considering Brassica genomes.

The problem of the eukaryotic genome size

The current state of knowledge concerning the unsolved problem of the huge interspecific eukaryotic genome size variations not correlating with the species phenotypic complexity (C-value enigma also known as C-value paradox) is reviewed. Characteristic features of eukaryotic genome structure and molecular mechanisms that are the basis of genome size changes are examined in connection with the C-value enigma. It is emphasized that endogenous mutagens, including reactive oxygen species, create a constant nuclear environment where any genome evolves. An original quantitative model and general conception are proposed to explain the C-value enigma. In accordance with the theory, the noncoding sequences of the eukaryotic genome provide genes with global and differential protection against chemical mutagens and (in addition to the anti-mutagenesis and DNA repair systems) form a new, third system that protects eukaryotic genetic information. The joint action of these systems controls the spontaneous mutation rate in coding sequences of the eukaryotic genome. It is hypothesized that the genome size is inversely proportional to functional efficiency of the anti-mutagenesis and/or DNA repair systems in a particular biological species. In this connection, a model of eukaryotic genome evolution is proposed.

Coding DNA repeated throughout intergenic regions of the Arabidopsis thaliana genome: evolutionary footprints of RNA silencing

Pyknons are non-random sequence patterns significantly repeated throughout non-coding genomic DNA that also appear at least once among coding genes. They are interesting because they portend an unforeseen connection between coding and non-coding DNA. Pyknons have only been discovered in the human genome, so it is unknown whether pyknons have wider biological relevance or are simply a phenomenon of the human genome. To address this, DNA sequence patterns from the Arabidopsis thaliana genome were detected using a probability-based method. 24 654 statistically significant sequence patterns, 16 to 24 nucleotides long, repeating 10 or more times in non-coding DNA also appeared in 46% of A. thaliana protein-coding genes. A. thaliana pyknons exhibit features similar to human pyknons, including being distinct sequence patterns, having multiple instances in genes and having remarkable similarity to small RNA sequences with roles in gene silencing. Chromosomal position mapping revealed that genomic pyknon density has concordance with siRNA and transposable element positioning density. Because the A. thaliana and human genomes have approximately the same number of genes but drastically different amounts of non-coding DNA, these data reveal that pyknons represent a biologically important link between coding and non-coding DNA. Because of the association of pyknons with siRNAs and localization to silenced regions of heterochromatin, we postulate that RNA-mediated gene silencing leads to the accumulation of gene sequences in non-coding DNA regions.

Plant origin and ploidy influence gene expression and life cycle characteristics in an invasive weed

BACKGROUND

Ecological, evolutionary and physiological studies have thus far provided an incomplete picture of why some plants become invasive; therefore we used genomic resources to complement and advance this field. In order to gain insight into the invasive mechanism of Centaurea stoebe we compared plants of three geo-cytotypes, native Eurasian diploids, native Eurasian tetraploids and introduced North American tetraploids, grown in a common greenhouse environment. We monitored plant performance characteristics and life cycle habits and characterized the expression of genes related to constitutive defense and genome stability using quantitative PCR.

RESULTS

Plant origin and ploidy were found to have a significant effect on both life cycle characteristics and gene expression, highlighting the importance of comparing appropriate taxonomic groups in studies of native and introduced plant species. We found that introduced populations of C. stoebe exhibit reduced expression of transcripts related to constitutive defense relative to their native tetraploid counterparts, as might be expected based on ideas of enemy release and rapid evolution. Measurements of several vegetative traits were similar for all geo-cytotypes; however, fecundity of tetraploids was significantly greater than diploids, due in part to their polycarpic nature. A simulation of seed production over time predicts that introduced tetraploids have the highest fecundity of the three geo-cytotypes.

CONCLUSION

Our results suggest that characterizing gene expression in an invasive species using populations from both its native and introduced range can provide insight into the biology of plant invasion that can complement traditional measurements of plant performance. In addition, these results highlight the importance of using appropriate taxonomic units in ecological genomics investigations.

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.

Invasion of the Arabidopsis genome by the tobacco retrotransposon Tnt1 is controlled by reversible transcriptional gene silencing

Long terminal repeat (LTR) retrotransposons are generally silent in plant genomes. However, they often constitute a large proportion of repeated sequences in plants. This suggests that their silencing is set up after a certain copy number is reached and/or that it can be released in some circumstances. We introduced the tobacco (Nicotiana tabacum) LTR retrotransposon Tnt1 into Arabidopsis (Arabidopsis thaliana), thus mimicking the horizontal transfer of a retrotransposon into a new host species and allowing us to study the regulatory mechanisms controlling its amplification. Tnt1 is transcriptionally silenced in Arabidopsis in a copy number-dependent manner. This silencing is associated with 24-nucleotide short-interfering RNAs targeting the promoter localized in the LTR region and with the non-CG site methylation of these sequences. Consequently, the silencing of Tnt1 is not released in methyltransferase1 mutants, in contrast to decrease in DNA methylation1 or polymerase IVa mutants. Stable reversion of Tnt1 silencing is obtained when the number of Tnt1 elements is reduced to two by genetic segregation. Our results support a model in which Tnt1 silencing in Arabidopsis occurs via an RNA-directed DNA methylation process. We further show that silencing can be partially overcome by some stresses.

DNA transposons and the evolution of eukaryotic genomes

Transposable elements are mobile genetic units that exhibit broad diversity in their structure and transposition mechanisms. Transposable elements occupy a large fraction of many eukaryotic genomes and their movement and accumulation represent a major force shaping the genes and genomes of almost all organisms. This review focuses on DNA-mediated or class 2 transposons and emphasizes how this class of elements is distinguished from other types of mobile elements in terms of their structure, amplification dynamics, and genomic effect. We provide an up-to-date outlook on the diversity and taxonomic distribution of all major types of DNA transposons in eukaryotes, including Helitrons and Mavericks. We discuss some of the evolutionary forces that influence their maintenance and diversification in various genomic environments. Finally, we highlight how the distinctive biological features of DNA transposons have contributed to shape genome architecture and led to the emergence of genetic innovations in different eukaryotic lineages.

A unified classification system for eukaryotic transposable elements

Our knowledge of the structure and composition of genomes is rapidly progressing in pace with their sequencing. The emerging data show that a significant portion of eukaryotic genomes is composed of transposable elements (TEs). Given the abundance and diversity of TEs and the speed at which large quantities of sequence data are emerging, identification and annotation of TEs presents a significant challenge. Here we propose the first unified hierarchical classification system, designed on the basis of the transposition mechanism, sequence similarities and structural relationships, that can be easily applied by non-experts. The system and nomenclature is kept up to date at the WikiPoson web site.

Case study for identification of potentially indel-caused alternative expression isoforms in the rice subspecies japonica and indica by integrative genome analysis

Alternative splicing (AS) is one of the most significant components of the functional complexity of the eukaryote genome, increasing protein diversity, creating isoforms, and affecting mRNA stability. Recently, whole genome sequences and large microarray data sets have become available, making data integration feasible and allowing the study of the possible regulatory mechanism of AS in rice (Oryza sativa) by erecting and testing hypotheses before doing bench studies. We have developed a new strategy and have identified 215 rice genes with alternative expression isoforms related to insertion and deletion (indel) between subspecies indica and subspecies japonica. We did a case study for alternative expression isoforms of the rice peroxidase gene LOC_Os06g48030 to investigate possible mechanisms by which indels caused alternative splicing between the indica and the japonica varieties by mining of array data together with validation by RT-PCR and genome sequencing analysis. Multiple poly(A) signals were detected in the specific indel region for LOC_Os06g48030. We present a new methodology to promote more discoveries of potentially indel-caused AS genes in rice, which may serve as the foundation for research into the regulatory mechanism of alternative expression isoforms between subspecies.

Rapid evolution and complex structural organization in genomic regions harboring multiple prolamin genes in the polyploid wheat genome

Genes encoding wheat prolamins belong to complicated multi-gene families in the wheat genome. To understand the structural complexity of storage protein loci, we sequenced and analyzed orthologous regions containing both gliadin and LMW-glutenin genes from the A and B genomes of a tetraploid wheat species, Triticum turgidum ssp. durum. Despite their physical proximity to one another, the gliadin genes and LMW-glutenin genes are organized quite differently. The gliadin genes are found to be more clustered than the LMW-glutenin genes which are separated from each other by much larger distances. The separation of the LMW-glutenin genes is the result of both the insertion of large blocks of repetitive DNA owing to the rapid amplification of retrotransposons and the presence of genetic loci interspersed between them. Sequence comparisons of the orthologous regions reveal that gene movement could be one of the major factors contributing to the violation of microcolinearity between the homoeologous A and B genomes in wheat. The rapid sequence rearrangements and differential insertion of repetitive DNA has caused the gene islands to be not conserved in compared regions. In addition, we demonstrated that the i-type LMW-glutenin originated from a deletion of 33-bps in the 5' coding region of the m-type gene. Our results show that multiple rounds of segmental duplication of prolamin genes have driven the amplification of the omega-gliadin genes in the region; such segmental duplication could greatly increase the repetitive DNA content in the genome depending on the amount of repetitive DNA present in the original duplicate region.

Common sequence polymorphisms shaping genetic diversity in Arabidopsis thaliana

The genomes of individuals from the same species vary in sequence as a result of different evolutionary processes. To examine the patterns of, and the forces shaping, sequence variation in Arabidopsis thaliana, we performed high-density array resequencing of 20 diverse strains (accessions). More than 1 million nonredundant single-nucleotide polymorphisms (SNPs) were identified at moderate false discovery rates (FDRs), and approximately 4% of the genome was identified as being highly dissimilar or deleted relative to the reference genome sequence. Patterns of polymorphism are highly nonrandom among gene families, with genes mediating interaction with the biotic environment having exceptional polymorphism levels. At the chromosomal scale, regional variation in polymorphism was readily apparent. A scan for recent selective sweeps revealed several candidate regions, including a notable example in which almost all variation was removed in a 500-kilobase window. Analyzing the polymorphisms we describe in larger sets of accessions will enable a detailed understanding of forces shaping population-wide sequence variation in A. thaliana.

Genome plasticity a key factor in the success of polyploid wheat under domestication

Wheat was domesticated about 10,000 years ago and has since spread worldwide to become one of the major crops. Its adaptability to diverse environments and end uses is surprising given the diversity bottlenecks expected from recent domestication and polyploid speciation events. Wheat compensates for these bottlenecks by capturing part of the genetic diversity of its progenitors and by generating new diversity at a relatively fast pace. Frequent gene deletions and disruptions generated by a fast replacement rate of repetitive sequences are buffered by the polyploid nature of wheat, resulting in subtle dosage effects on which selection can operate.

A first step in understanding an invasive weed through its genes: an EST analysis of invasive Centaurea maculosa

BACKGROUND

The economic and biological implications of plant invasion are overwhelming; however, the processes by which plants become successful invaders are not well understood. Limited genetic resources are available for most invasive and weedy species, making it difficult to study molecular and genetic aspects that may be associated with invasion.

RESULTS

As an initial step towards understanding the molecular mechanisms by which plants become invasive, we have generated a normalized Expressed Sequence Tag (EST) library comprising seven invasive populations of Centaurea maculosa, an invasive aster in North America. Seventy-seven percent of the 4423 unique transcripts showed significant similarity to existing proteins in the NCBI database and could be grouped based on gene ontology assignments.

CONCLUSION

The C. maculosa EST library represents an initial step towards looking at gene-specific expression in this species, and will pave the way for creation of other resources such as microarray chips that can help provide a view of global gene expression in invasive C. maculosa and its native counterparts. To our knowledge, this is the first published set of ESTs derived from an invasive weed that will be targeted to study invasive behavior. Understanding the genetic basis of evolution for increased invasiveness in exotic plants is critical to understanding the mechanisms through which exotic invasions occur.

Mechanisms and rates of birth and death of dispersed duplicated genes during the evolution of a multigene family in diploid and tetraploid wheats

A family of 5 genes that evolved within the past 1.9 Myr in diploid wheat was characterized. The ancestral gene, ALP-A1, is on chromosome 1A and encodes an aci-reductone dioxygenase-like protein. The duplicated genes ALP-A2, ALP-A3, ALP-A4.1, and ALP-A4.2 acquired complete coding sequences but lost the original promoter. They are on chromosomes 4A, 2A, 6A and 6A, respectively, and evolved sequentially, the youngest duplicated gene always producing the next duplicate. It is shown that dispersed gene duplication rate consists of the primary rate (duplications of ancestral genes) and the secondary rate (duplications of genes that had been generated by recent duplications). The primary rate was 2.5 x 10(-3) gene(-1) Myr(-1) in diploid wheat. The secondary rate was 5.2 x 10(-2) gene(-1) Myr(-1) in the ALP family. The 20-fold acceleration of the secondary rate was caused by the insertion of the ALP-A2 gene into a novel type transposon. Only the ALP-A1 and ALP-A3 genes are transcribed. The transcription of ALP-A3 is directed by a promoter within a DNA fragment similar to a CACTA type of DNA transposons, making ALP-A3 a new gene. The ALP-A3 transcript is longer than that of the ALP-A1. The half-life of ALP duplicated genes was estimated to be 0.87 Myr. Strong purifying selection acting on the ancestral gene ALP-A1 was undiminished by the evolution of duplicated genes. The evolution of the ALP family shows that repeated elements facilitate both gene duplication and expression of duplicated genes and highlights their importance for the evolution of gene repertoire in large plant genomes.

Retrotransposons: Metaparasites and Agents of Genome Evolution

Transposable elements comprise the bulk of higher plant genomes. The majority of these elements are the Class I LTR retrotransposons, which transpose via an RNA intermediate in a "Copy-and-Paste" mechanism. Because retrotransposons use cellular resources and their own enzymes to replicate independently of the genome as a whole, and have thereby become in many cases more predominant than the cellular genes, they have been considered "selfish DNA" and nuclear parasites. They are thought to share many features of the internal life cycle of retroviruses such as HIV (lentiviruses). However, whereas at least some of the retroviruses arriving in an organism during an infection must be functional in order for the infection to proceed, some LTR retrotransposon families appear to completely lack active members even though they remain mobile. Furthermore, the process of retrotransposition is inherently error-prone and mutagenic, giving rise to "pseudospecies," or clusters of imperfect copies. The non-autonomous retrotransposons are able to cis- and trans-parasitize host retrotransposons to gain mobility, much as do defective interfering particles of RNA viruses. Hence, a complex dynamic is set up, whereby the impact of retrotransposons on genomes can be under selection on the organismal level; the impact of non-autonomous retrotransposons on autonomous ones can likewise be under selection if there is selection on the autonomous elements themselves. We are exploring the retrotransposon life cycle and the causes and possible consequences of non-autonomy at each stage regarding genome evolution.