Evolutionary dynamics of an ancient retrotransposon family provides insights into evolution of genome size in the genus Oryza

Oryza genome evolution through a tetraploid lens

Oryza is a remarkable genus comprising 27 species and 11 genome types, with ~3.4-fold genome size variation, that possesses a virtually untapped reservoir of genes that can be used for crop improvement and neodomestication. Here we present 11 chromosome-level assemblies (nine tetraploid, two diploid) in the context of ~15 million years of evolution and show that the core Oryza (sub)genome is only ~200 Mb and largely syntenic, whereas the remaining nuclear fractions (~80-600 Mb) are intermingled, plastic and rapidly evolving. For the halophyte Oryza coarctata, we found that despite detection of gene fractionation in the subgenomes, homoeologous genes were expressed at higher levels in one subgenome over the other in a mosaic form, demonstrating subgenome equivalence. The integration of these 11 new reference genomes with previously published genome datasets provides a nearly complete view of the consequences of evolution for genome diversification across the genus.

Complete mitochondrial genome assembly and comparative analysis of Colocasia esculenta

Colocasia esculenta ranks as the fifth most important tuber crop and is known for its high nutritional and medicinal value. However, there is no research on its mitochondrial genome, hindering in-depth exploration of its genomic resources and genetic relationships. Using second- and third-generation sequencing technologies, we assembled and annotated the mitogenome of C. esculenta. Its mitogenome mainly consists of five circular DNA molecules, with a total length of 594,811 bp and a GC content of 46.25%. A total of 55 genes, 157 simple sequence repeats, 29 tandem repeat sequences, 202 dispersed repeat sequences, and 625 RNA editing sites were detected. Most protein-coding genes use ATG as the start codon, and the third position of the codon tends to be A or T (U). GAA, AUU, and UUU are the most common codons in C. esculenta mitochondria. Finally, based on 28 representative plant species, a phylogenetic tree was constructed, revealing a close relationship between C. esculenta and Araceae. This study provides comprehensive information on C. esculenta, laying a foundation for crop genetics and molecular breeding.

Comparative analysis of the complete mitogenomes of Camellia sinensis var. sinensis and C. sinensis var. assamica provide insights into evolution and phylogeny relationship

Introduction

Among cultivated tea plants (Camellia sinensis), only four mitogenomes for C. sinensis var. assamica (CSA) have been reported so far but none for C. sinensis var. sinensis (CSS). Here, two mitogenomes of CSS (CSSDHP and CSSRG) have been sequenced and assembled.

Methods

Using a combination of Illumina and Nanopore data for the first time. Comparison between CSS and CSA mitogenomes revealed a huge heterogeneity.

Results

The number of the repetitive sequences was proportional to the mitogenome size and the repetitive sequences dominated the intracellular gene transfer segments (accounting for 88.7%- 92.8% of the total length). Predictive RNA editing analysis revealed that there might be significant editing in NADH dehydrogenase subunit transcripts. Codon preference analysis showed a tendency to favor A/T bases and T was used more frequently at the third base of the codon. ENc plots analysis showed that the natural selection play an important role in shaping the codon usage bias, and Ka/Ks ratios analysis indicated Nad1 and Sdh3 genes may have undergone positive selection. Further, phylogenetic analysis shows that six C. sinensis clustered together, with the CSA and CSS forming two distinct branches, suggesting two different evolutionary pathway.

Discussion

Altogether, this investigation provided an insight into evolution and phylogeny relationship of C. sinensis mitogenome, thereby enhancing comprehension of the evolutionary patterns within C. sinensis species.

Comparative Analysis of Transposable Elements in Strawberry Genomes of Different Ploidy Levels

Transposable elements (TEs) make up a large portion of plant genomes and play a vital role in genome structure, function, and evolution. Cultivated strawberry (Fragaria x ananassa) is one of the most important fruit crops, and its octoploid genome was formed through several rounds of genome duplications from diploid ancestors. Here, we built a pan-genome TE library for the Fragaria genus using ten published strawberry genomes at different ploidy levels, including seven diploids, one tetraploid, and two octoploids, and performed comparative analysis of TE content in these genomes. The TEs comprise 51.83% (F. viridis) to 60.07% (F. nilgerrensis) of the genomes. Long terminal repeat retrotransposons (LTR-RTs) are the predominant TE type in the Fragaria genomes (20.16% to 34.94%), particularly in F. iinumae (34.94%). Estimating TE content and LTR-RT insertion times revealed that species-specific TEs have shaped each strawberry genome. Additionally, the copy number of different LTR-RT families inserted in the last one million years reflects the genetic distance between Fragaria species. Comparing cultivated strawberry subgenomes to extant diploid ancestors showed that F. vesca and F. iinumae are likely the diploid ancestors of the cultivated strawberry, but not F. viridis. These findings provide new insights into the TE variations in the strawberry genomes and their roles in strawberry genome evolution.

Genome assembly and analysis of Lactuca virosa: implications for lettuce breeding

Lettuce (Lactuca sativa L.) is a leafy vegetable crop with ongoing breeding efforts related to quality, resilience, and innovative production systems. To breed resilient and resistant lettuce in the future, valuable genetic variation found in close relatives could be further exploited. Lactuca virosa (2x = 2n = 18), a wild relative assigned to the tertiary lettuce gene pool, has a much larger genome (3.7 Gbp) than Lactuca sativa (2.5 Gbp). It has been used in interspecific crosses and is a donor to modern crisphead lettuce cultivars. Here, we present a de novo reference assembly of L. virosa with high continuity and complete gene space. This assembly facilitated comparisons to the genome of L. sativa and to that of the wild species L. saligna, a representative of the secondary lettuce gene pool. To assess the diversity in gene content, we classified the genes of the 3 Lactuca species as core, accessory, and unique. In addition, we identified 3 interspecific chromosomal inversions compared to L. sativa, which each may cause recombination suppression and thus hamper future introgression breeding. Using 3-way comparisons in both reference-based and reference-free manners, we show that the proliferation of long-terminal repeat elements has driven the genome expansion of L. virosa. Further, we performed a genome-wide comparison of immune genes, nucleotide-binding leucine-rich repeat, and receptor-like kinases among Lactuca spp. and indicated the evolutionary patterns and mechanisms behind their expansions. These genome analyses greatly facilitate the understanding of genetic variation in L. virosa, which is beneficial for the breeding of improved lettuce varieties.

Comparative Analysis of the Complete Mitochondrial Genomes of Apium graveolens and Apium leptophyllum Provide Insights into Evolution and Phylogeny Relationships

The genus Apium, belonging to the family Apiaceae, comprises roughly 20 species. Only two species, Apium graveolens and Apium leptophyllum, are available in China and are both rich in nutrients and have favorable medicinal properties. However, the lack of genomic data has severely constrained the study of genetics and evolution in Apium plants. In this study, Illumina NovaSeq 6000 and Nanopore sequencing platforms were employed to identify the mitochondrial genomes of A. graveolens and A. leptophyllum. The complete lengths of the mitochondrial genomes of A. graveolens and A. leptophyllum were 263,017 bp and 260,164 bp, respectively, and contained 39 and 36 protein-coding genes, five and six rRNA genes, and 19 and 20 tRNA genes. Consistent with most angiosperms, both A. graveolens and A. leptophyllum showed a preference for codons encoding leucine (Leu). In the mitochondrial genome of A. graveolens, 335 SSRs were detected, which is higher than the 196 SSRs found in the mitochondrial genome of A. leptophyllum. Studies have shown that the most common RNA editing type is C-to-U, but, in our study, both A. graveolens and A. leptophyllum exhibited the U-C editing type. Furthermore, the transfer of the mitochondrial genomes of A. graveolens and A. leptophyllum into the chloroplast genomes revealed homologous sequences, accounting for 8.14% and 4.89% of the mitochondrial genome, respectively. Lastly, in comparing the mitochondrial genomes of 29 species, it was found that A. graveolens, A. leptophyllum, and Daucus carota form a sister group with a support rate of 100%. Overall, this investigation furnishes extensive insights into the mitochondrial genomes of A. graveolens and A. leptophyllum, thereby enhancing comprehension of the traits and evolutionary patterns within the Apium genus. Additionally, it offers supplementary data for evolutionary and comparative genomic analyses of other species within the Apiaceae family.

The final piece of the Triangle of U: Evolution of the tetraploid Brassica carinata genome

Ethiopian mustard (Brassica carinata) is an ancient crop with remarkable stress resilience and a desirable seed fatty acid profile for biofuel uses. Brassica carinata is one of six Brassica species that share three major genomes from three diploid species (AA, BB, and CC) that spontaneously hybridized in a pairwise manner to form three allotetraploid species (AABB, AACC, and BBCC). Of the genomes of these species, that of B. carinata is the least understood. Here, we report a chromosome scale 1.31-Gbp genome assembly with 156.9-fold sequencing coverage for B. carinata, completing the reference genomes comprising the classic Triangle of U, a classical theory of the evolutionary relationships among these six species. Our assembly provides insights into the hybridization event that led to the current B. carinata genome and the genomic features that gave rise to the superior agronomic traits of B. carinata. Notably, we identified an expansion of transcription factor networks and agronomically important gene families. Completion of the Triangle of U comparative genomics platform has allowed us to examine the dynamics of polyploid evolution and the role of subgenome dominance in the domestication and continuing agronomic improvement of B. carinata and other Brassica species.

Comparative repeatome analysis reveals new evidence on genome evolution in wild diploid Arachis (Fabaceae) species

Opposing changes in the abundance of satellite DNA and long terminal repeat (LTR) retroelements are the main contributors to the variation in genome size and heterochromatin amount in Arachis diploids. The South American genus Arachis (Fabaceae) comprises 83 species organized in nine taxonomic sections. Among them, section Arachis is characterized by species with a wide genome and karyotype diversity. Such diversity is determined mainly by the amount and composition of repetitive DNA. Here we performed computational analysis on low coverage genome sequencing to infer the dynamics of changes in major repeat families that led to the differentiation of genomes in diploid species (x = 10) of genus Arachis, focusing on section Arachis. Estimated repeat content ranged from 62.50 to 71.68% of the genomes. Species with different genome composition tended to have different landscapes of repeated sequences. Athila family retrotransposons were the most abundant and variable lineage among Arachis repeatomes, with peaks of transpositional activity inferred at different times in the evolution of the species. Satellite DNAs (satDNAs) were less abundant, but differentially represented among species. High rates of evolution of an AT-rich superfamily of satDNAs led to the differential accumulation of heterochromatin in Arachis genomes. The relationship between genome size variation and the repetitive content is complex. However, largest genomes presented a higher accumulation of LTR elements and lower contents of satDNAs. In contrast, species with lowest genome sizes tended to accumulate satDNAs in detriment of LTR elements. Phylogenetic analysis based on repetitive DNA supported the genome arrangement of section Arachis. Altogether, our results provide the most comprehensive picture on the repeatome dynamics that led to the genome differentiation of Arachis species.

Characterisation of LTR-Retrotransposons of Stevia rebaudiana and Their Use for the Analysis of Genetic Variability

Stevia rebaudiana is one of the most important crops belonging to the Asteraceae family. Stevia is cultivated all over the world as it represents a valid natural alternative to artificial sweeteners thanks to its leaves, which produce steviol glycosides that have high sweetening power and reduced caloric value. In this work, the stevia genome sequence was used to isolate and characterise full-length long-terminal repeat retrotransposons (LTR-REs), which account for more than half of the genome. The Gypsy retrotransposons were twice as abundant as the Copia ones. A disproportionate abundance of elements belonging to the Chromovirus/Tekay lineage was observed among the Gypsy elements. Only the SIRE and Angela lineages represented significant portions of the genome among the Copia elements. The dynamics with which LTR-REs colonised the stevia genome were also estimated; all isolated full-length elements turned out to be relatively young, with a proliferation peak around 1–2 million years ago. However, a different analysis conducted by comparing sequences encoding retrotranscriptase showed the occurrence of an older period in which there was a lot of LTR-RE proliferation. Finally, a group of isolated full-length elements belonging to the lineage Angela was used to analyse the genetic variability in 25 accessions of S. rebaudiana using the Inter-Retrotransposon Amplified Polymorphism (IRAP) protocol. The obtained fingerprints highlighted a high degree of genetic variability and were used to study the genomic structures of the different accessions. It was hypothesised that there are four ancestral subpopulations at the root of the analysed accessions, which all turned out to be admixed. Overall, these data may be useful for genome sequence annotations and for evaluating genetic variability in this species, which may be useful in stevia breeding.

The Singular Evolution of Olea Genome Structure

The current view of plant genome evolution proposes that genome size has mainly been determined by polyploidisation and amplification/loss of transposons, with a minor role played by other repeated sequences, such as tandem repeats. In cultivated olive (Olea europaea subsp. europaea var. europaea), available data suggest a singular model of genome evolution, in which a massive expansion of tandem-repeated sequences accompanied changes in nuclear architecture. This peculiar scenario highlights the importance of focusing on Olea genus evolution, to shed light on mechanisms that led to its present genomic structure. Next-generation sequencing technologies, bioinformatics and in situ hybridisation were applied to study the genomic structure of five related Olea taxa, which originated at different times from their last common ancestor. On average, repetitive DNA in the Olea taxa ranged from ~59% to ~73% of the total genome, showing remarkable differences in terms of composition. Among repeats, we identified 11 major families of tandem repeats, with different abundances in the analysed taxa, five of which were novel discoveries. Interestingly, overall tandem repeat abundance was inversely correlated to that of retrotransposons. This trend might imply a competition in the proliferation of these repeat classes. Indeed, O. paniculata, the species closest to the Olea common ancestor, showed very few tandem-repeated sequences, while it was rich in long terminal repeat retrotransposons, suggesting that the amplification of tandem repeats occurred after its divergence from the Olea ancestor. Furthermore, some tandem repeats were physically localised in closely related O. europaea subspecies (i.e., cultivated olive and O. europaea subsp. cuspidata), which showed a significant difference in tandem repeats abundance. For 4 tandem repeats families, a similar number of hybridisation signals were observed in both subspecies, apparently indicating that, after their dissemination throughout the olive genome, these tandem repeats families differentially amplified maintaining the same positions in each genome. Overall, our research identified the temporal dynamics shaping genome structure during Olea speciation, which represented a singular model of genome evolution in higher plants.

The complete chloroplast genomes of three Hamamelidaceae species: Comparative and phylogenetic analyses

Hamamelidaceae is an important group that represents the origin and early evolution of angiosperms. Its plants have many uses, such as timber, medical, spice, and ornamental uses. In this study, the complete chloroplast genomes of Loropetalum chinense (R. Br.) Oliver, Corylopsis glandulifera Hemsl., and Corylopsis velutina Hand.-Mazz. were sequenced using the Illumina NovaSeq 6000 platform. The sizes of the three chloroplast genomes were 159,402 bp (C. glandulifera), 159,414 bp (C. velutina), and 159,444 bp (L. chinense), respectively. These chloroplast genomes contained typical quadripartite structures with a pair of inverted repeat (IR) regions (26,283, 26,283, and 26,257 bp), a large single-copy (LSC) region (88,134, 88,146, and 88,160 bp), and a small single-copy (SSC) region (18,702, 18,702, and 18,770 bp). The chloroplast genomes encoded 132-133 genes, including 85-87 protein-coding genes, 37-38 tRNA genes, and 8 rRNA genes. The coding regions were composed of 26,797, 26,574, and 26,415 codons, respectively, most of which ended in A/U. A total of 37-43 long repeats and 175-178 simple sequence repeats (SSRs) were identified, and the SSRs contained a higher number of A + T than G + C bases. The genome comparison showed that the IR regions were more conserved than the LSC or SSC regions, while the noncoding regions contained higher variability than the gene coding regions. Phylogenetic analyses revealed that species in the same genus tended to cluster together. Chunia Hung T. Chang, Mytilaria Lecomte, and Disanthus Maxim. may have diverged early and Corylopsis Siebold & Zucc. was closely related to Loropetalum R. Br. This study provides valuable information for further species identification, evolution, and phylogenetic studies of Hamamelidaceae plants.

InpactorDB: A Classified Lineage-Level Plant LTR Retrotransposon Reference Library for Free-Alignment Methods Based on Machine Learning

Long terminal repeat (LTR) retrotransposons are mobile elements that constitute the major fraction of most plant genomes. The identification and annotation of these elements via bioinformatics approaches represent a major challenge in the era of massive plant genome sequencing. In addition to their involvement in genome size variation, LTR retrotransposons are also associated with the function and structure of different chromosomal regions and can alter the function of coding regions, among others. Several sequence databases of plant LTR retrotransposons are available for public access, such as PGSB and RepetDB, or restricted access such as Repbase. Although these databases are useful to identify LTR-RTs in new genomes by similarity, the elements of these databases are not fully classified to the lineage (also called family) level. Here, we present InpactorDB, a semi-curated dataset composed of 130,439 elements from 195 plant genomes (belonging to 108 plant species) classified to the lineage level. This dataset has been used to train two deep neural networks (i.e., one fully connected and one convolutional) for the rapid classification of these elements. In lineage-level classification approaches, we obtain up to 98% performance, indicated by the F1-score, precision and recall scores.

Improved Brassica oleracea JZS assembly reveals significant changing of LTR-RT dynamics in different morphotypes

Brassica oleracea is an important vegetable crop that has provided ancestor genomes of the two most important Brassica oil crops, Brassica napus and Brassica carinata. The current B. oleracea reference genome (JZS, also named 02-12) displays problems of large mis-assemblies, low sequence continuity, and low assembly integrity, thus limiting genomic analysis. We reported an updated assembly of the B. oleracea reference genome (JZS v2) obtained through single-molecule sequencing and chromosome conformation capture technologies. We assembled an additional 83.16 Mb of genomic sequences, and the updated genome features a contig N50 size of 2.37 Mb, representing an ~ 88-fold improvement. We detected a new round of long terminal repeat retrotransposon (LTR-RT) burst in the new assembly. Comparative analysis with the reported genome sequences of two other genomes of B. oleracea (TO1000 and HDEM) identified extensive gene order and gene structural variation. In addition, we found that the genome-specific amplification of Gypsy-like LTR-RTs occurred around 0-1 million years ago (MYA). In particular, the athila, tat, and Del families were extensively amplified in JZS around 0-1 MYA. Moreover, we identified that the syntenic genes were modified due to the insertion of genome-specific LTR-RTs. These results indicated that the genome-specific LTR-RT dynamics were associated with genome diversification in B. oleracea.

The impact of transposable elements on the structure, evolution and function of the rice genome

Transposable elements (TEs) are ubiquitous in plants and are the primary genomic component of the majority of taxa. Knowledge of their impact on the structure, function and evolution of plant genomes is therefore a priority in the field of genomics. Rice, as one of the most prevalent crops for food security worldwide, has been subjected to intense research efforts over recent decades. Consequently, a considerable amount of genomic resources has been generated and made freely available to the scientific community. These can be exploited both to improve our understanding of some basic aspects of genome biology of this species and to develop new concepts for crop improvement. In this review, we describe the current knowledge on how TEs have shaped rice chromosomes and propose a new strategy based on a genome-wide association study (GWAS) to address the important question of their functional impact on this crop.

Potential of Oryza officinalis to augment the cold tolerance genetic mechanisms of Oryza sativa by network complementation

Oryza officinalis is an accessible alien donor for genetic improvement of rice. Comparison across a representative panel of Oryza species showed that the wild O. officinalis and cultivated O. sativa ssp. japonica have similar cold tolerance potentials. The possibility that either distinct or similar genetic mechanisms are involved in the low temperature responses of each species was addressed by comparing their transcriptional networks. General similarities were supported by shared transcriptomic signatures indicative of equivalent metabolic, hormonal, and defense status. However, O. officinalis has maintained an elaborate cold-responsive brassinosteroid-regulated BES1-network that appeared to have been fragmented in O. sativa. BES1-network is potentially important for integrating growth-related responses with physiological adjustments and defenses through the protection of photosynthetic machinery and maintenance of stomatal aperture, oxidative defenses, and osmotic adjustment. Equivalent physiological processes are functional in O. sativa but their genetic mechanisms are under the direct control of ABA-dependent, DREB-dependent and/or oxidative-mediated networks uncoupled to BES1. While O. officinalis and O. sativa represent long periods of speciation and domestication, their comparable cold tolerance potentials involve equivalent physiological processes but distinct genetic networks. BES1-network represents a novel attribute of O. officinalis with potential applications in diversifying or complementing other mechanisms in the cultivated germplasm.

Transcriptome Analysis and Functional Identification of Xa13 and Pi-ta Orthologs in Oryza granulata

Nees & Arn. ex Watt, a perennial wild rice species with a GG genome, preserves many important genes for cultivated rice ( L.) improvement. At present, however, no genetic resource is available for studying . Here, we report 91,562 high-quality transcripts of assembled de novo. Moreover, comparative transcriptome analysis revealed that 1311 single-copy orthologous pairs shared by and (Zoll. & Moritzi) Baill. that may have undergone adaptive evolution. We performed an analysis of the genes potentially involved in plant-pathogen interactions to explore the molecular basis of disease resistance, and isolated the full-length cDNAs of () and () orthologs from . The overexpression of in Nipponbare and functional characterization showed enhanced the resistance of transgenic Nipponbare to rice blast resulting from the presence of the gene. , an alternatively spliced transcript of the blast resistance gene in encodes a 1024-amino acid polypeptide with a C-terminal thioredoxin domain. This study provides an important resource for functional and evolutionary studies of the genus .

Different histories of two highly variable LTR retrotransposons in sunflower species

In the Helianthus genus, very large intra- and interspecific variability related to two specific retrotransposons of Helianthus annuus (Helicopia and SURE) exists. When comparing these two sequences to sunflower sequence databases recently produced by our lab, the Helicopia family was shown to belong to the Maximus/SIRE lineage of the Sirevirus genus of the Copia superfamily, whereas the SURE element (whose superfamily was not even previously identified) was classified as a Gypsy element of the Ogre/Tat lineage of the Metavirus genus. Bioinformatic analysis of the two retrotransposon families revealed their genomic abundance and relative proliferation timing. The genomic abundance of these families differed significantly among 12 Helianthus species. The ratio between the abundance of long terminal repeats and their reverse transcriptases suggested that the SURE family has relatively more solo long terminal repeats than does Helicopia. Pairwise comparisons of Illumina reads encoding the reverse transcriptase domain indicated that SURE amplification may have occurred more recently than that of Helicopia. Finally, the analysis of population structure based on the SURE and Helicopia polymorphisms of 32 Helianthus species evidenced two subpopulations, which roughly corresponded to species of the Helianthus and Divaricati/Ciliares sections. However, a number of species showed an admixed structure, confirming the importance of interspecific hybridisation in the evolution of this genus. In general, these two retrotransposon families differentially contributed to interspecific variability, emphasising the need to refer to specific families when studying genome evolution.

Genome-wide analysis of LTR-retrotransposon diversity and its impact on the evolution of the genus Helianthus (L.)

Background

Genome divergence by mobile elements activity and recombination is a continuous process that plays a key role in the evolution of species. Nevertheless, knowledge on retrotransposon-related variability among species belonging to the same genus is still limited. Considering the importance of the genus Helianthus, a model system for studying the ecological genetics of speciation and adaptation, we performed a comparative analysis of the repetitive genome fraction across ten species and one subspecies of sunflower, focusing on long terminal repeat retrotransposons at superfamily, lineage and sublineage levels.

Results

After determining the relative genome size of each species, genomic DNA was isolated and subjected to Illumina sequencing. Then, different assembling and clustering approaches allowed exploring the repetitive component of all genomes. On average, repetitive DNA in Helianthus species represented more than 75% of the genome, being composed mostly by long terminal repeat retrotransposons. Also, the prevalence of Gypsy over Copia superfamily was observed and, among lineages, Chromovirus was by far the most represented. Although nearly all the same sublineages are present in all species, we found considerable variability in the abundance of diverse retrotransposon lineages and sublineages, especially between annual and perennial species.

Conclusions

This large variability should indicate that different events of amplification or loss related to these elements occurred following species separation and should have been involved in species differentiation. Our data allowed us inferring on the extent of interspecific repetitive DNA variation related to LTR-RE abundance, investigating the relationship between changes of LTR-RE abundance and the evolution of the genus, and determining the degree of coevolution of different LTR-RE lineages or sublineages between and within species. Moreover, the data suggested that LTR-RE abundance in a species was affected by the annual or perennial habit of that species.

Electronic supplementary material

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

Evolution of plant genome architecture

We have witnessed an explosion in our understanding of the evolution and structure of plant genomes in recent years. Here, we highlight three important emergent realizations: (1) that the evolutionary history of all plant genomes contains multiple, cyclical episodes of whole-genome doubling that were followed by myriad fractionation processes; (2) that the vast majority of the variation in genome size reflects the dynamics of proliferation and loss of lineage-specific transposable elements; and (3) that various classes of small RNAs help shape genomic architecture and function. We illustrate ways in which understanding these organism-level and molecular genetic processes can be used for crop plant improvement.

Structural Diversity of a Novel LTR Retrotransposon, RTPOSON, in the Genus Oryza

Retrotransposons with long terminal repeats (LTRs) are the most abundant transposable elements in plant genomes. A novel LTR retrotransposon named RTPOSON primarily occurs in the genus Oryza and in several species of the Poaceae family. RTPOSON has been identified in the Ty1-copia group of retrotransposons because two of its open reading frames encode an uncharacterized protein and UBN2_2 and zinc knuckle, respectively. More than 700 RTPOSONs were identified in Oryza genomes; 127 RTPOSONs with LTRs and gag-pol elements were classified into three subgroups. The subgroup RTPOSON_sub3 had the smallest DNA size and 97% (32/33) of RTPOSON elements from Oryza punctata are classified in this group. The insertion time of these RTPOSONs varied and their proliferation occurred within the last 8 Mya, with two bursting periods within the last 1.5–5.0 Mya. A total of 37 different orthologous insertions of RTPOSONs, with different nested transposable elements and gene fragments, were identified by comparing the genomes of ssp. japonica cv. Nipponbare and ssp. indica cv. 93–11. A part of intact RTPOSON elements was evolved independently after the divergence of indica and japonica. In addition, intact RTPOSONs and homologous fragments were preferentially retained or integrated in genic regions. This novel LTR retrotransposon, RTPOSON, might have an impact on genome evolution, genic innovation, and genetic variation.

Landscape and evolutionary dynamics of terminal repeat retrotransposons in miniature in plant genomes

Background

Terminal repeat retrotransposons in miniature (TRIMs) are a unique group of small long terminal repeat retrotransposons that are difficult to identify. Thus far, only a few TRIMs have been characterized in the euphyllophytes, and their evolutionary and biological significance as well as their transposition mechanisms are poorly understood.

Results

Using a combination of de novo and homology-based methods, we annotate TRIMs in 48 plant genome sequences, spanning land plants to algae. The TRIMs are grouped into 156 families including 145 that were previously undefined. Notably, we identify the first TRIMs in a lycophyte and non-vascular plants. The majority of the TRIM families are highly conserved and shared within and between plant families. Unlike other long terminal repeat retrotransposons, TRIMs are enriched in or near genes; they are also targeted by sRNAs between 21 and 24 nucleotides in length, and are frequently found in CG body-methylated genes. Importantly, we also identify putative autonomous retrotransposons and very recent transpositions of a TRIM element in Oryza sativa.

Conclusions

We perform the most comprehensive analysis of TRIM transposons thus far and report that TRIMs are ubiquitous across plant genomes. Our results show that TRIMs are more frequently associated with large and CG body-methylated genes that have undergone strong purifying selection. Our findings also indicate that TRIMs are likely derived from internal deletions of large long terminal repeat retrotransposons. Finally, our data and methodology are important resources for the characterization and evolutionary and genomic studies of long terminal repeat retrotransposons in other genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0867-y) contains supplementary material, which is available to authorized users.

The Ty1-copia LTR retroelement family PARTC is highly conserved in conifers over 200 MY of evolution

Long Terminal Repeat retroelements (LTR-RTs) are a major component of many plant genomes. Although well studied and described in angiosperms, their features and dynamics are poorly understood in gymnosperms. Representative complete copies of a Ty1-copia element isolate in Picea abies and named PARTC were identified in six other conifer species (Picea glauca, Pinus sylvestris, Pinus taeda, Abies sibirica, Taxus baccata and Juniperus communis) covering more than 200 million years of evolution. Here we characterized the structure of this element, assessed its abundance across conifers, studied the modes and timing of its amplification, and evaluated the degree of conservation of its extant copies at nucleotide level over distant species. We demonstrated that the element is ancient, abundant, widespread and its paralogous copies are present in the genera Picea, Pinus and Abies as an LTR-RT family. The amplification leading to the extant copies of PARTC occurred over long evolutionary times spanning 10s of MY and mostly took place after the speciation of the conifers analyzed. The level of conservation of PARTC is striking and may be explained by low substitution rates and limited removal mechanisms for LTR-RTs. These PARTC features and dynamics are representative of a more general scenario for LTR-RTs in gymnosperms quite different from that characterizing the vast majority of LTR-RT elements in angiosperms.

The peculiar landscape of repetitive sequences in the olive (Olea europaea L.) genome

Analyzing genome structure in different species allows to gain an insight into the evolution of plant genome size. Olive (Olea europaea L.) has a medium-sized haploid genome of 1.4 Gb, whose structure is largely uncharacterized, despite the growing importance of this tree as oil crop. Next-generation sequencing technologies and different computational procedures have been used to study the composition of the olive genome and its repetitive fraction. A total of 2.03 and 2.3 genome equivalents of Illumina and 454 reads from genomic DNA, respectively, were assembled following different procedures, which produced more than 200,000 differently redundant contigs, with mean length higher than 1,000 nt. Mapping Illumina reads onto the assembled sequences was used to estimate their redundancy. The genome data set was subdivided into highly and medium redundant and nonredundant contigs. By combining identification and mapping of repeated sequences, it was established that tandem repeats represent a very large portion of the olive genome (∼31% of the whole genome), consisting of six main families of different length, two of which were first discovered in these experiments. The other large redundant class in the olive genome is represented by transposable elements (especially long terminal repeat-retrotransposons). On the whole, the results of our analyses show the peculiar landscape of the olive genome, related to the massive amplification of tandem repeats, more than that reported for any other sequenced plant genome.

Controlled activation of retrotransposition for plant breeding

Plant genomes consist to a large extent of transposable elements (TEs), predominantly retrotransposons. Their accumulation through periodic transposition bursts has shaped the structure and regulatory organization of plant genomes, often contributing to phenotypic traits. Transposon-generated phenotypes selected by humans during plant domestication have been maintained under strict selection during subsequent plant breeding. Our knowledge of the epigenetic, environmental, and developmental regulation of TE activity has advanced considerably in recent years. Here I will consider TEs as an attractive endogenous source of genetic variation that can be liberated in a controlled fashion and thus offer novel phenotypic diversity to be selected for crop improvement.

LTR retrotransposon dynamics in the evolution of the olive (Olea europaea) genome

Improved knowledge of genome composition, especially of its repetitive component, generates important information for both theoretical and applied research. The olive repetitive component is made up of two main classes of sequences: tandem repeats and retrotransposons (REs). In this study, we provide characterization of a sample of 254 unique full-length long terminal repeat (LTR) REs. In the sample, Ty1-Copia elements were more numerous than Ty3-Gypsy elements. Mapping a large set of Illumina whole-genome shotgun reads onto the identified retroelement set revealed that Gypsy elements are more redundant than Copia elements. The insertion time of intact retroelements was estimated based on sister LTR’s divergence. Although some elements inserted relatively recently, the mean insertion age of the isolated retroelements is around 18 million yrs. Gypsy and Copia retroelements showed different waves of transposition, with Gypsy elements especially active between 10 and 25 million yrs ago and nearly inactive in the last 7 million yrs. The occurrence of numerous solo-LTRs related to isolated full-length retroelements was ascertained for two Gypsy elements and one Copia element. Overall, the results reported in this study show that RE activity (both retrotransposition and DNA loss) has impacted the olive genome structure in more ancient times than in other angiosperms.

Rapid diversification of five Oryza AA genomes associated with rice adaptation

Comparative genomic analyses among closely related species can greatly enhance our understanding of plant gene and genome evolution. We report de novo-assembled AA-genome sequences for Oryza nivara, Oryza glaberrima, Oryza barthii, Oryza glumaepatula, and Oryza meridionalis. Our analyses reveal massive levels of genomic structural variation, including segmental duplication and rapid gene family turnover, with particularly high instability in defense-related genes. We show, on a genomic scale, how lineage-specific expansion or contraction of gene families has led to their morphological and reproductive diversification, thus enlightening the evolutionary process of speciation and adaptation. Despite strong purifying selective pressures on most Oryza genes, we documented a large number of positively selected genes, especially those genes involved in flower development, reproduction, and resistance-related processes. These diversifying genes are expected to have played key roles in adaptations to their ecological niches in Asia, South America, Africa and Australia. Extensive variation in noncoding RNA gene numbers, function enrichment, and rates of sequence divergence might also help account for the different genetic adaptations of these rice species. Collectively, these resources provide new opportunities for evolutionary genomics, numerous insights into recent speciation, a valuable database of functional variation for crop improvement, and tools for efficient conservation of wild rice germplasm.

Sequence-Based Analysis of Structural Organization and Composition of the Cultivated Sunflower (Helianthus annuus L.) Genome

Sunflower is an important oilseed crop, as well as a model system for evolutionary studies, but its 3.6 gigabase genome has proven difficult to assemble, in part because of the high repeat content of its genome. Here we report on the sequencing, assembly, and analyses of 96 randomly chosen BACs from sunflower to provide additional information on the repeat content of the sunflower genome, assess how repetitive elements in the sunflower genome are organized relative to genes, and compare the genomic distribution of these repeats to that found in other food crops and model species. We also examine the expression of transposable element-related transcripts in EST databases for sunflower to determine the representation of repeats in the transcriptome and to measure their transcriptional activity. Our data confirm previous reports in suggesting that the sunflower genome is >78% repetitive. Sunflower repeats share very little similarity to other plant repeats such as those of Arabidopsis, rice, maize and wheat; overall 28% of repeats are “novel” to sunflower. The repetitive sequences appear to be randomly distributed within the sequenced BACs. Assuming the 96 BACs are representative of the genome as a whole, then approximately 5.2% of the sunflower genome comprises non TE-related genic sequence, with an average gene density of 18kbp/gene. Expression levels of these transposable elements indicate tissue specificity and differential expression in vegetative and reproductive tissues, suggesting that expressed TEs might contribute to sunflower development. The assembled BACs will also be useful for assessing the quality of several different draft assemblies of the sunflower genome and for annotating the reference sequence.

Could abiotic stress tolerance in wild relatives of rice be used to improve Oryza sativa?

Oryza sativa and Oryza glaberrima have been selected to acquire and partition resources efficiently as part of the process of domestication. However, genetic diversity in cultivated rice is limited compared to wild Oryza species, in spite of 120,000 genotypes being held in gene banks. By contrast, there is untapped diversity in the more than 20 wild species of Oryza, some having been collected from just a few coastal locations (e.g. Oryza schlechteri), while others are widely distributed (e.g. Oryza nivara and Oryza rufipogon). The extent of DNA sequence diversity and phenotypic variation is still being established in wild Oryza, with genetic barriers suggesting a vast range of morphologies and function even within species, such as has been demonstrated for Oryza meridionalis. With increasing climate variability and attempts to make more marginal land arable, abiotic and biotic stresses will be managed over the coming decades by tapping into the genetic diversity of wild relatives of O. sativa. To help create a more targeted approach to sourcing wild rice germplasm for abiotic stress tolerance, we have created a climate distribution map by plotting the natural occurrence of all Oryza species against corresponding temperature and moisture data. We then discuss interspecific variation in phenotype and its significance for rice, followed by a discussion of ways to integrate germplasm from wild relatives into domesticated rice.

The repetitive component of the sunflower genome as shown by different procedures for assembling next generation sequencing reads

Background

Next generation sequencing provides a powerful tool to study genome structure in species whose genomes are far from being completely sequenced. In this work we describe and compare different computational approaches to evaluate the repetitive component of the genome of sunflower, by using medium/low coverage Illumina or 454 libraries.

Results

By varying sequencing technology (Illumina or 454), coverage (0.55 x-1.25 x), assemblers and assembly procedures, six different genomic databases were produced. The annotation of these databases showed that they were composed of different proportions of repetitive DNA families. The final assembly of the sequences belonging to the six databases produced a whole genome set of 283,800 contigs. The redundancy of each contig was estimated by mapping the whole genome set with a large Illumina read set and measuring the number of matched Illumina reads. The repetitive component amounted to 81% of the sunflower genome, that is composed mainly of numerous families of Gypsy and Copia retrotransposons. Also many families of non autonomous retrotransposons and DNA transposons (especially of the Helitron superfamily) were identified.

Conclusions

The results substantially matched those previously obtained by using a Sanger-sequenced shotgun library and a standard 454 whole-genome-shotgun approach, indicating the reliability of the proposed procedures also for other species. The repetitive sequences were collected to produce a database, SUNREP, that will be useful for the annotation of the sunflower genome sequence and for studying the genome evolution in dicotyledons.

The International Oryza Map Alignment Project: development of a genus-wide comparative genomics platform to help solve the 9 billion-people question

The wild relatives of rice contain a virtually untapped reservoir of traits that can be used help drive the 21st century green revolution aimed at solving world food security issues by 2050. To better understand and exploit the 23 species of the Oryza genus the rice research community is developing foundational resources composed of: 1) reference genomes and transcriptomes for all 23 species; 2) advanced mapping populations for functional and breeding studies; and 3) in situ conservation sites for ecological, evolutionary and population genomics. To this end, 16 genome sequencing projects are currently underway, and all completed assemblies have been annotated; and several advanced mapping populations have been developed, and more will be generated, mapped, and phenotyped, to uncover useful alleles. As wild Oryza populations are threatened by human activity and climate change, we also discuss the urgent need for sustainable in situ conservation of the genus.

How important are transposons for plant evolution?

For decades, transposable elements have been known to produce a wide variety of changes in plant gene expression and function. This has led to the idea that transposable element activity has played a key part in adaptive plant evolution. This Review describes the kinds of changes that transposable elements can cause, discusses evidence that those changes have contributed to plant evolution and suggests future strategies for determining the extent to which these changes have in fact contributed to plant adaptation and evolution. Recent advances in genomics and phenomics for a range of plant species, particularly crops, have begun to allow the systematic assessment of these questions.

Comparative sequence analysis of the Ghd7 orthologous regions revealed movement of Ghd7 in the grass genomes

Ghd7 is an important rice gene that has a major effect on several agronomic traits, including yield. To reveal the origin of Ghd7 and sequence evolution of this locus, we performed a comparative sequence analysis of the Ghd7 orthologous regions from ten diploid Oryza species, Brachypodium distachyon, sorghum and maize. Sequence analysis demonstrated high gene collinearity across the genus Oryza and a disruption of collinearity among non-Oryza species. In particular, Ghd7 was not present in orthologous positions except in Oryza species. The Ghd7 regions were found to have low gene densities and high contents of repetitive elements, and that the sizes of orthologous regions varied tremendously. The large transposable element contents resulted in a high frequency of pseudogenization and gene movement events surrounding the Ghd7 loci. Annotation information and cytological experiments have indicated that Ghd7 is a heterochromatic gene. Ghd7 orthologs were identified in B. distachyon, sorghum and maize by phylogenetic analysis; however, the positions of orthologous genes differed dramatically as a consequence of gene movements in grasses. Rather, we identified sequence remnants of gene movement of Ghd7 mediated by illegitimate recombination in the B. distachyon genome.

Functional and structural divergence of an unusual LTR retrotransposon family in plants

Retrotransposons with long terminal repeats (LTRs) more than 3 kb are not frequent in most eukaryotic genomes. Rice LTR retrotransposon, Retrosat2, has LTRs greater than 3.2 kb and two open reading frames (ORF): ORF1 encodes enzymes for retrotransposition whereas no function can be assigned to ORF0 as it is not found in any other organism. A variety of experimental and in silico approaches were used to determine the origin of Retrosat2 and putative function of ORF0. Our data show that not only is Retrosat2 highly abundant in the Oryza genus, it may yet be active in rice. Homologs of Retrosat2 were identified in maize, sorghum, Arabidopsis and other plant genomes suggesting that the Retrosat2 family is of ancient origin. Several putatively cis-acting elements, some multicopy, that regulate retrotransposon replication or responsiveness to environmental factors were found in the LTRs of Retrosat2. Unlike the ORF1, the ORF0 sequences from Retrosat2 and homologs are divergent at the sequence level, 3D-structures and predicted biological functions. In contrast to other retrotransposon families, Retrosat2 and its homologs are dispersed throughout genomes and not concentrated in the specific chromosomal regions, such as centromeres. The genomic distribution of Retrosat2 homologs varies across species which likely reflects the differing evolutionary trajectories of this retrotransposon family across diverse species.

Mutator System Derivatives Isolated from Sugarcane Genome Sequence

Mutator-like transposase is the most represented transposon transcript in the sugarcane transcriptome. Phylogenetic reconstructions derived from sequenced transcripts provided evidence that at least four distinct classes exist (I-IV) and that diversification among these classes occurred early in Angiosperms, prior to the divergence of Monocots/Eudicots. The four previously described classes served as probes to select and further sequence six BAC clones from a genomic library of cultivar R570. A total of 579,352 sugarcane base pairs were produced from these "Mutator system" BAC containing regions for further characterization. The analyzed genomic regions confirmed that the predicted structure and organization of the Mutator system in sugarcane is composed of two true transposon lineages, each containing a specific terminal inverted repeat and two transposase lineages considered to be domesticated. Each Mutator transposase class displayed a particular molecular structure supporting lineage specific evolution. MUSTANG, previously described domesticated genes, are located in syntenic regions across Sacharineae and, as expected for a host functional gene, posses the same gene structure as in other Poaceae. Two sequenced BACs correspond to hom(eo)logous locus with specific retrotransposon insertions that discriminate sugarcane haplotypes. The comparative studies presented, add information to the Mutator systems previously identified in the maize and rice genomes by describing lineage specific molecular structure and genomic distribution pattern in the sugarcane genome. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s12042-012-9104-y) contains supplementary material, which is available to authorized users.

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.

A highly conserved, small LTR retrotransposon that preferentially targets genes in grass genomes

LTR retrotransposons are often the most abundant components of plant genomes and can impact gene and genome evolution. Most reported LTR retrotransposons are large elements (>4 kb) and are most often found in heterochromatic (gene poor) regions. We report the smallest LTR retrotransposon found to date, only 292 bp. The element is found in rice, maize, sorghum and other grass genomes, which indicates that it was present in the ancestor of grass species, at least 50-80 MYA. Estimated insertion times, comparisons between sequenced rice lines, and mRNA data indicate that this element may still be active in some genomes. Unlike other LTR retrotransposons, the small LTR retrotransposons (SMARTs) are distributed throughout the genomes and are often located within or near genes with insertion patterns similar to MITEs (miniature inverted repeat transposable elements). Our data suggests that insertions of SMARTs into or near genes can, in a few instances, alter both gene structures and gene expression. Further evidence for a role in regulating gene expression, SMART-specific small RNAs (sRNAs) were identified that may be involved in gene regulation. Thus, SMARTs may have played an important role in genome evolution and genic innovation and may provide a valuable tool for gene tagging systems in grass.

Isolation and characterization of genome-specific markers in Oryza species with the BB genome

Wild species of rice with many valuable agronomic traits are an important genetic resource for improving cultivated rice by wide hybridization. Genome- or chromosome-specific markers are useful for monitoring genome introgression and for identifying genome components. From 47 random amplified polymorphic DNAs (RAPDs) of nine Oryza species, three bands (Ogla225, Opun225, and Opun246) were found to be genome specific with distinct sizes. Their specificities were further characterized by Southern hybridization, sequence analysis, and fluorescent in situ hybridization (FISH). Ogla225 is specifically amplified from the AA genome but homologous sequences were conserved among Oryza species. Opun225 occurs at a low copy number although is specifically amplified from Oryza punctata. There are estimated 2000-3300 repeats of Opun246 in each haploid genome of Oryza species with the BB or BBCC genome. Clusters of Opun246 repeats were detected at heterochromatic regions on almost all chromosomes of the BB genomes by FISH. Opun246 may be a useful marker for monitoring the introgression of BB genome or for identifying the conserved components of BB genome in genetic resource. The results from this study and our previous study both indicate that numerous unique repeats play role in the differentiation of the BB genome from other Oryza genomes.

Polymorphisms and evolutionary history of retrotransposon insertions in rice promoters

Retrotransposons are ubiquitous in higher plant genomes. The presence or absence of retrotransposons in whole genome and high throughput genomic sequence (HTGS) from cultivated and wild rice was investigated to understand the organization and evolution of retrotransposon insertions in promoter regions. Approximately half of the Oryza sativa subsp. japonica ‘Nipponbare’ promoters with retrotransposons conserved in Oryza sativa subsp. indica ‘93-11’ and four wild rice species showed higher sequence conservation in retrotransposon than nonretrotransposon regions. We further investigated, in detail, the evolutionary dynamics of five retrotransposons in the promoter regions of 95 rice genotypes. Our data suggest that four of five insertions (Rp2–Rp5) occurred in the ancestor of AA genome, while the other insertion (Rp1) predates the ancestral divergence of Oryza officinalis (CC genome). Four retrotransposons (Rp2–Rp5) were present in 52% (Rp2), 29% (Rp3), 53% (Rp4), and 43% (Rp5) of the rice genotypes with AA genome type, and the fifth retrotransposon (Rp1) was present in 95% of the rice genotypes with AA, BBCC, or CC genome types. Furthermore, most of these retrotransposons were found to evolve slower than flanking promoter regions, suggesting a role in promoter function for regulating downstream genes.

Epigenetic contribution to stress adaptation in plants

Plant epigenetics has recently gained unprecedented interest, not only as a subject of basic research but also as a possible new source of beneficial traits for plant breeding. We discuss here mechanisms of epigenetic regulation that should be considered when undertaking the latter. Since these mechanisms are responsible for the formation of heritable epigenetic gene variants (epialleles) and also regulate transposons mobility, both aspects could be exploited to broaden plant phenotypic and genetic variation, which could improve long-term plant adaptation to environmental challenges and, thus, increase productivity.

An siRNA pathway prevents transgenerational retrotransposition in plants subjected to stress

Eukaryotic genomes consist to a significant extent of retrotransposons that are suppressed by host epigenetic mechanisms, preventing their uncontrolled propagation. However, it is not clear how this is achieved. Here we show that in Arabidopsis seedlings subjected to heat stress, a copia-type retrotransposon named ONSEN (Japanese 'hot spring') not only became transcriptionally active but also synthesized extrachromosomal DNA copies. Heat-induced ONSEN accumulation was stimulated in mutants impaired in the biogenesis of small interfering RNAs (siRNAs); however, there was no evidence of transposition occurring in vegetative tissues. After stress, both ONSEN transcripts and extrachromosomal DNA gradually decayed and were no longer detected after 20-30 days. Surprisingly, a high frequency of new ONSEN insertions was observed in the progeny of stressed plants deficient in siRNAs. Insertion patterns revealed that this transgenerational retrotransposition occurred during flower development and before gametogenesis. Therefore in plants with compromised siRNA biogenesis, memory of stress was maintained throughout development, priming ONSEN to transpose during differentiation of generative organs. Retrotransposition was not observed in the progeny of wild-type plants subjected to stress or in non-stressed mutant controls, pointing to a crucial role of the siRNA pathway in restricting retrotransposition triggered by environmental stress. Finally, we found that natural and experimentally induced variants in ONSEN insertions confer heat responsiveness to nearby genes, and therefore mobility bursts may generate novel, stress-responsive regulatory gene networks.

Construction, characterization, and preliminary BAC-end sequence analysis of a bacterial artificial chromosome library of the tea plant (Camellia sinensis)

We describe the construction and characterization of a publicly available BAC library for the tea plant, Camellia sinensis. Using modified methods, the library was constructed with the aim of developing public molecular resources to advance tea plant genomics research. The library consists of a total of 401,280 clones with an average insert size of 135 kb, providing an approximate coverage of 13.5 haploid genome equivalents. No empty vector clones were observed in a random sampling of 576 BAC clones. Further analysis of 182 BAC-end sequences from randomly selected clones revealed a GC content of 40.35% and low chloroplast and mitochondrial contamination. Repetitive sequence analyses indicated that LTR retrotransposons were the most predominant sequence class (86.93%-87.24%), followed by DNA retrotransposons (11.16%-11.69%). Additionally, we found 25 simple sequence repeats (SSRs) that could potentially be used as genetic markers.

Spatio-temporal patterns of genome evolution in allotetraploid species of the genus Oryza

Despite knowledge that polyploidy is widespread and a major evolutionary force in flowering plant diversification, detailed comparative molecular studies on polyploidy have been confined to only a few species and families. The genus Oryza is composed of 23 species that are classified into ten distinct 'genome types' (six diploid and four polyploid), and is emerging as a powerful new model system to study polyploidy. Here we report the identification, sequence and comprehensive comparative annotation of eight homoeologous genomes from a single orthologous region (Adh1-Adh2) from four allopolyploid species representing each of the known Oryza genome types (BC, CD, HJ and KL). Detailed comparative phylogenomic analyses of these regions within and across species and ploidy levels provided several insights into the spatio-temporal dynamics of genome organization and evolution of this region in 'natural' polyploids of Oryza. The major findings of this study are that: (i) homoeologous genomic regions within the same nucleus experience both independent and parallel evolution, (ii) differential lineage-specific selection pressures do not occur between polyploids and their diploid progenitors, (iii) there have been no dramatic structural changes relative to the diploid ancestors, (iv) a variation in the molecular evolutionary rate exists between the two genomes in the BC complex species even though the BC and CD polyploid species appear to have arisen <2 million years ago, and (v) there are no clear distinctions in the patterns of genome evolution in the diploid versus polyploid species.

Orthologous comparisons of the Hd1 region across genera reveal Hd1 gene lability within diploid Oryza species and disruptions to microsynteny in Sorghum

Heading date is one of the most important quantitative traits responsible for the domestication of rice. We compared a 155-kb reference segment of the Oryza sativa ssp. japonica cv. Nipponbare genome surrounding Hd1, a major heading date gene in rice, with orthologous regions from nine diploid Oryza species that diverged over a relatively short time frame (∼16 My) to study sequence evolution around a domestication locus. The orthologous Hd1 region from Sorghum bicolor was included to compare and contrast the evolution in a more distant relative of rice. Consistent with other observations at the adh1/adh2, monoculm1, and sh2/a1 loci in grass species, we found high gene colinearity in the Hd1 region amidst size differences that were lineage specific and long terminal repeat retrotransposon driven. Unexpectedly, the Hd1 gene was deleted in O. glaberrima, whereas the O. rufipogon and O. punctata copies had degenerative mutations, suggesting that other heading date loci might compensate for the loss or nonfunctionality of Hd1 in these species. Compared with the japonica Hd1 region, the orthologous region in sorghum exhibited micro-rearrangements including gene translocations, seven additional genes, and a gene triplication and truncation event predating the divergence from Oryza.

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.

Analysis of collinear regions of Oryza AA and CC genomes

Comparative analyses of genome structure and sequence of closely related species have yielded insights into the evolution and function of plant genomes. A total of 103,844 BAC end sequences delegated approximately 73.8 Mb of O. officinalis that belongs to the CC genome type of the rice genus Oryza were obtained and compared with the genome sequences of rice cultivar, O. sativa ssp. japonica cv. Nipponbare. We found that more than 45% of O. officinalis genome consists of repeat sequences, which is higher than that of Nipponbare cultivar. To further investigate the evolutionary divergence of AA and CC genomes, two BAC-contigs of O. officinalis were compared with the collinear genomic regions of Nipponbare. Of 57 genes predicted in the AA genome orthologous regions, 39 had orthologs in the regions of the CC genome. Alignment of the orthologous regions indicated that the CC genome has undergone expansion in both genic and intergenic regions through primarily retroelement insertion. Particularly, the density of RNA transposable elements was 17.95% and 1.78% in O. officinalis and O. sativa, respectively. This explains why the orthologous region is about 100 kb longer in the CC genome in comparison to the AA genome.