Rice transposable elements: a survey of 73,000 sequence-tagged-connectors

New insights into Oryza genome evolution: high gene colinearity and differential retrotransposon amplification

An approximately 247-kb genomic region from FF genome of wild rice Oryza brachyantha, possessing the smallest Oryza genome, was compared to the orthologous approximately 450-kb region from AA genome, O. sativa L. ssp. japonica. 37 of 38 genes in the orthologous regions are shared between japonica and O. brachyantha. Analyses of nucleotide substitution in coding regions suggest the two genomes diverged approximately 10 million years ago. Comparisons of transposable elements (TEs) reveal that the density of DNA TEs in O. brachyantha is comparable to O. sativa; however, the density of RNA TEs is dramatically lower. The genomic fraction of RNA TEs in japonica is two times greater than in O. brachyantha. Differences, particularly in RNA TEs, in this region and in BAC end sequences from five wild and two cultivated Oryza species explain major genome size differences between sativa and brachyantha. Gene expression analyses of three ObDREB1 genes in the sequenced region indicate orthologous genes retain similar expression patterns following cold stress. Our results demonstrate that size and number of RNA TEs play a major role in genomic differentiation and evolution in Oryza. Additionally, distantly related O. brachyantha shares colinearity with O. sativa, offering opportunities to use comparative genomics to explore the genetic diversity of wild species to improve cultivated rice.

Construction and characterization of two BAC libraries from Brachypodium distachyon, a new model for grass genomics

Brachypodium is well suited as a model system for temperate grasses because of its compact genome and a range of biological features. In an effort to develop resources for genome research in this emerging model species, we constructed 2 bacterial artificial chromosome (BAC) libraries from an inbred diploid Brachypodium distachyon line, Bd21, using restriction enzymes HindIII and BamHI. A total of 73,728 clones (36,864 per BAC library) were picked and arrayed in 192,384-well plates. The average insert size for the BamHI and HindIII libraries is estimated to be 100 and 105 kb, respectively, and inserts of chloroplast origin account for 4.4% and 2.4%, respectively. The libraries individually represent 9.4- and 9.9-fold haploid genome equivalents with combined 19.3-fold genome coverage, based on a genome size of 355 Mb reported for the diploid Brachypodium, implying a 99.99% probability that any given specific sequence will be present in each library. Hybridization of the libraries with 8 starch biosynthesis genes was used to empirically evaluate this theoretical genome coverage; the frequency at which these genes were present in the library clones gave an estimated coverage of 11.6- and 19.6-fold genome equivalents. To obtain a first view of the sequence composition of the Brachypodium genome, 2185 BAC end sequences (BES) representing 1.3 Mb of random genomic sequence were compared with the NCBI GenBank database and the GIRI repeat database. Using a cutoff expectation value of E<10-10, only 3.3% of the BESs showed similarity to repetitive sequences in the existing database, whereas 40.0% had matches to the sequences in the EST database, suggesting that a considerable portion of the Brachypodium genome is likely transcribed. When the BESs were compared with individual EST databases, more matches hit wheat than maize, although their EST collections are of a similar size, further supporting the close relationship between Brachypodium and the Triticeae. Moreover, 122 BESs have significant matches to wheat ESTs mapped to individual chromosome bin positions. These BACs represent colinear regions containing the mapped wheat ESTs and would be useful in identifying additional markers for specific wheat chromosome regions.

A family of human microRNA genes from miniature inverted-repeat transposable elements

While hundreds of novel microRNA (miRNA) genes have been discovered in the last few years alone, the origin and evolution of these non-coding regulatory sequences remain largely obscure. In this report, we demonstrate that members of a recently discovered family of human miRNA genes, hsa-mir-548, are derived from Made1 transposable elements. Made1 elements are short miniature inverted-repeat transposable elements (MITEs), which consist of two 37 base pair (bp) terminal inverted repeats that flank 6 bp of internal sequence. Thus, Made1 elements are nearly perfect palindromes, and when expressed as RNA they form highly stable hairpin loops. Apparently, these Made1-related structures are recognized by the RNA interference enzymatic machinery and processed to form 22 bp mature miRNA sequences. Consistent with their origin from MITEs, hsa-mir-548 genes are primate-specific and have many potential paralogs in the human genome. There are more than 3,500 putative hsa-mir-548 target genes; analysis of their expression profiles and functional affinities suggests cancer-related regulatory roles for hsa-mir-548. Taken together, the characteristics of Made1 elements, and MITEs in general, point to a specific mechanism for the generation of numerous small regulatory RNAs and target sites throughout the genome. The evolutionary lineage-specific nature of MITEs could also provide for the generation of novel regulatory phenotypes related to species diversification. Finally, we propose that MITEs may represent an evolutionary link between siRNAs and miRNAs.

Mobilization and evolutionary history of miniature inverted-repeat transposable elements (MITEs) in Beta vulgaris L

We have identified three families of miniature inverted-repeat transposable elements (VulMITEs) in the genome of sugar beet (Beta vulgaris L.), evidently derived from a member of the Vulmar family of mariner transposons. While VulMITEs I are typical stowaway-like MITEs, VulMITEs II and VulMITEs III are rearranged stowaway elements of increased size. The integration of divergent moderately and highly repetitive sequences into VulMITEs II and, in particular in VulMITEs III, respectively, shows that amplification of repetitive DNA by MITEs contribute to the increase of genome size with possible implications for plant genome evolution. Fluorescent in-situ hybridization (FISH), for the first time visualizing stowaway MITE distribution on plant chromosomes, revealed a dispersed localization of VulMITEs along all B. vulgaris chromosomes. Analysis of the flanking sequences identified a dispersed repeat as target site for the integration of the stowaway element VulMITE I. Recent transposition of VulMITE I, which most likely occurred during the domestication of cultivated beets, was concluded from insertional polymorphisms between different B. vulgaris cultivars and species.

Transposition of a fungal miniature inverted-repeat transposable element through the action of a Tc1-like transposase

The mimp1 element previously identified in the ascomycete fungus Fusarium oxysporum has hallmarks of miniature inverted-repeat transposable elements (MITEs): short size, terminal inverted repeats (TIRs), structural homogeneity, and a stable secondary structure. Since mimp1 has no coding capacity, its mobilization requires a transposase-encoding element. On the basis of the similarity of TIRs and target-site preference with the autonomous Tc1-like element impala, together with a correlated distribution of both elements among the Fusarium genus, we investigated the ability of mimp1 to jump upon expression of the impala transposase provided in trans. Under these conditions, we present evidence that mimp1 transposes by a cut-and-paste mechanism into TA dinucleotides, which are duplicated upon insertion. Our results also show that mimp1 reinserts very frequently in genic regions for at least one-third of the cases. We also show that the mimp1/impala double-component system is fully functional in the heterologous species F. graminearum, allowing the development of a highly efficient tool for gene tagging in filamentous fungi.

Improved detection of episomal Banana streak viruses by multiplex immunocapture PCR

Banana streak viruses (BSV) are currently the main viral constraint to Musa germplasm movement, genetic improvement and mass propagation. Therefore, it is necessary to develop and implement BSV detection strategies that are both reliable and sensitive, such as PCR-based techniques. Unfortunately, BSV endogenous pararetrovirus sequences (BSV EPRVs) are present in the genome of Musa balbisiana. They interfere with PCR-based detection of episomal BSV in infected banana and plantain, such as immunocapture PCR. Therefore, a multiplex, immunocapture PCR (M-IC-PCR) was developed for the detection of BSV. Musa sequence tagged microsatellite site (STMS) primers were selected and used in combination with BSV species-specific primers in order to monitor possible contamination by Musa genomic DNA, using multiplex PCR. Furthermore, immunocapture conditions were optimized in order to prevent Musa DNA from interfering with episomal BSV DNA during the PCR step. This improved detection method successfully allowed the accurate, specific and sensitive detection of episomal DNA only from distinct BSV species. Its implementation should benefit PCR-based detection of viruses for which homologous sequences are present in the genome of their hosts, including transgenic plants expressing viral sequences.

Analysis of papaya BAC end sequences reveals first insights into the organization of a fruit tree genome

Papaya (Carica papaya L.) is a major tree fruit crop of tropical and subtropical regions with an estimated genome size of 372 Mbp. We present the analysis of 4.7% of the papaya genome based on BAC end sequences (BESs) representing 17 million high-quality bases. Microsatellites discovered in 5,452 BESs and flanking primer sequences are available to papaya breeding programs at http://www.genomics.hawaii.edu/papaya/BES . Sixteen percent of BESs contain plant repeat elements, the vast majority (83.3%) of which are class I retrotransposons. Several novel papaya-specific repeats were identified. Approximately 19.1% of the BESs have homology to Arabidopsis cDNA. Increasing numbers of completely sequenced plant genomes and BES projects enable novel approaches to comparative plant genomics. Paired BESs of Carica, Arabidopsis, Populus, Brassica and Lycopersicon were mapped onto the completed genomes of Arabidopsis and Populus. In general the level of microsynteny was highest between closely related organisms. However, papaya revealed a higher degree of apparent synteny with the more distantly related poplar than with the more closely related Arabidopsis. This, as well as significant colinearity observed between peach and poplar genome sequences, support recent observations of frequent genome rearrangements in the Arabidopsis lineage and suggest that the poplar genome sequence may be more useful for elucidating the papaya and other rosid genomes. These insights will play a critical role in selecting species and sequencing strategies that will optimally represent crop genomes in sequence databases.

Ac-like transposons in populations of wild diploid Triticeae species: comparative analysis of chromosomal distribution

Data are presented on the intra- and interspecific differences/similarities in chromosomal patterns of Ac-like elements (hAT family) in ecologically contrasted populations of three Triticeae species - Aegilops speltoides, Triticum urartu, and Hordeum spontaneum. Application of original computer software made it possible to precisely map transposon clusters and to link them to known chromosomal markers (rDNA sites, centromeres, and heterochromatin regions). From our data we can specify the most visible features of Ac-like elements chromosomal distribution: preferential concentration in chromosomal proximal regions; high percentage of clusters on the border between euchromatin and heterochromatin; complementary chromosomal arrangement towards En/Spm transposons (CACTA); population-specific insertions into centromeres; more differences in total cluster numbers between populations of self-pollinated species than between populations of cross-pollinated species. The application of statistical simulation (Resampling) method to analysis of data indicates that ecology may play a certain role in dynamics of Ac-like elements. Comparison of real Ayala distances, as well as real chromosomal distribution of Ac-like elements in populations of two species with different mating systems with the same but randomly simulated parameters, revealed that non-random population structure in the Mediterranean floral zone suffers and becomes chaotic in the Irano-Turanian zone.

Horizontal transfer of a plant transposon

The majority of well-documented cases of horizontal transfer between higher eukaryotes involve the movement of transposable elements between animals. Surprisingly, although plant genomes often contain vast numbers of these mobile genetic elements, no evidence of horizontal transfer of a nuclear-encoded transposon between plant species has been detected to date. The most mutagenic known plant transposable element system is the Mutator system in maize. Mu-like elements (MULEs) are widespread among plants, and previous analysis has suggested that the distribution of various subgroups of MULEs is patchy, consistent with horizontal transfer. We have sequenced portions of MULE transposons from a number of species of the genus Setaria and compared them to each other and to publicly available databases. A subset of these elements is remarkably similar to a small family of MULEs in rice. A comparison of noncoding and synonymous sequences revealed that the observed similarity is not due to selection at the amino acid level. Given the amount of time separating Setaria and rice, the degree of similarity between these elements excludes the possibility of simple vertical transmission of this class of MULEs. This is the first well-documented example of horizontal transfer of any nuclear-encoded genes between higher plants.

Sequencing and analysis of MULE transposons and their surrounding genomic regions from closely related grass species and rice provides evidence of horizontal transfer in plants.

Cotton genome mapping with new microsatellites from Acala ‘Maxxa’ BAC-ends

Fine mapping and positional cloning will eventually improve with the anchoring of additional markers derived from genomic clones such as BACs. From 2,603 new BAC-end genomic sequences from Gossypium hirsutum Acala 'Maxxa', 1,316 PCR primer pairs (designated as MUSB) were designed to flank microsatellite or simple sequence repeat motif sequences. Most (1164 or 88%) MUSB primer pairs successfully amplified DNA from three species of cotton with an average of three amplicons per marker and 365 markers (21%) were polymorphic between G. hirsutum and G. barbadense. An interspecific RIL population developed from the above two entries was used to map 433 marker loci and 46 linkage groups with a genetic distance of 2,126.3 cM covering approximately 45% of the cotton genome and an average distance between two loci of 4.9 cM. Based on genome-specific chromosomes identified in G. hirsutum tetraploid (A and D), 56.9% of the coverage was located on the A subgenome while 39.7% was assigned to the D subgenome in the genetic map, suggesting that the A subgenome may be more polymorphic and recombinationally active than originally thought. The linkage groups were assigned to 23 of the 26 chromosomes. This is the first genetic map in which the linkage groups A01 and A02/D03 have been assigned to specific chromosomes. In addition the MUSB-derived markers from BAC-end sequences markers allows fine genetic and QTL mapping of important traits and for the first time provides reconciliation of the genetic and physical maps. Limited QTL analyses suggested that loci on chromosomes 2, 3, 12, 15 and 18 may affect variation in fiber quality traits. The original BAC clones containing the newly mapped MUSB that tag the QTLs provide critical DNA regions for the discovery of gene sequences involved in biological processes such as fiber development and pest resistance in cotton.

Decoding the rice genome

Rice cultivation is one of the most important agricultural activities on earth, with nearly 90% of it being produced in Asia. It belongs to the family of crops that includes wheat, maize and barley, and it supplies more than 50% of calories consumed by the world population. Its immense economic value and a relatively small genome size makes it a focal point for scientific investigations, so much so that four whole genome sequence drafts with varying qualities have been generated by both public and privately funded ventures. The availability of a complete and high-quality map-based sequence has provided the opportunity to study genome organization and evolution. Most importantly, the order and identity of 37,544 genes of rice have been unraveled. The sequence provides the required ingredients for functional genomics and molecular breeding programs aimed at unraveling intricate cellular processes and improving rice productivity.

The Oryza bacterial artificial chromosome library resource: construction and analysis of 12 deep-coverage large-insert BAC libraries that represent the 10 genome types of the genus Oryza

Rice (Oryza sativa L.) is the most important food crop in the world and a model system for plant biology. With the completion of a finished genome sequence we must now functionally characterize the rice genome by a variety of methods, including comparative genomic analysis between cereal species and within the genus Oryza. Oryza contains two cultivated and 22 wild species that represent 10 distinct genome types. The wild species contain an essentially untapped reservoir of agriculturally important genes that must be harnessed if we are to maintain a safe and secure food supply for the 21st century. As a first step to functionally characterize the rice genome from a comparative standpoint, we report the construction and analysis of a comprehensive set of 12 BAC libraries that represent the 10 genome types of Oryza. To estimate the number of clones required to generate 10 genome equivalent BAC libraries we determined the genome sizes of nine of the 12 species using flow cytometry. Each library represents a minimum of 10 genome equivalents, has an average insert size range between 123 and 161 kb, an average organellar content of 0.4%-4.1% and nonrecombinant content between 0% and 5%. Genome coverage was estimated mathematically and empirically by hybridization and extensive contig and BAC end sequence analysis. A preliminary analysis of BAC end sequences of clones from these libraries indicated that LTR retrotransposons are the predominant class of repeat elements in Oryza and a roughly linear relationship of these elements with genome size was observed.

Random sheared fosmid library as a new genomic tool to accelerate complete finishing of rice (Oryza sativa spp. Nipponbare) genome sequence: sequencing of gap-specific fosmid clones uncovers new euchromatic portions of the genome

The International Rice Genome Sequencing Project has recently announced the high-quality finished sequence that covers nearly 95% of the japonica rice genome representing 370 Mbp. Nevertheless, the current physical map of japonica rice contains 62 physical gaps corresponding to approximately 5% of the genome, that have not been identified/represented in the comprehensive array of publicly available BAC, PAC and other genomic library resources. Without finishing these gaps, it is impossible to identify the complete complement of genes encoded by rice genome and will also leave us ignorant of some 5% of the genome and its unknown functions. In this article, we report the construction and characterization of a tenfold redundant, 40 kbp insert fosmid library generated by random mechanical shearing. We demonstrated its utility in refining the physical map of rice by identifying and in silico mapping 22 gap-specific fosmid clones with particular emphasis on chromosomes 1, 2, 6, 7, 8, 9 and 10. Further sequencing of 12 of the gap-specific fosmid clones uncovered unique rice genome sequence that was not previously reported in the finished IRGSP sequence and emphasizes the need to complete finishing of the rice genome.

Characterisation of pararetrovirus-like sequences in the genome of potato (Solanum tuberosum)

Three families of pararetrovirus-like sequences were isolated from the genome of potato using PCR of a characteristic fragment extending from the end of the transactivator domain. The potato pararetrovirus-like sequences are abundant in the nuclear genome of potato as demonstrated by their hybridisation to high-molecular weight DNA in Southern transfers and by fluorescence in situ hybridisation. Sequencing of cloned PCR products demonstrated that the potato pararetrovirus-like sequences were similar to other pararetroviral sequences and also to some expressed sequences from tobacco and tomato, notably from callus and Agrobacterium-infected tissue. It is possible that the potato pararetroviral sequences defend against viral genes via silencing mechanisms, although, as in Petunia or banana, their transcription and recombination may cause infection under stress conditions.

Updating of transposable element annotations from large wheat genomic sequences reveals diverse activities and gene associations

Triticeae species (including wheat, barley and rye) have huge and complex genomes due to polyploidization and a high content of transposable elements (TEs). TEs are known to play a major role in the structure and evolutionary dynamics of Triticeae genomes. During the last 5 years, substantial stretches of contiguous genomic sequence from various species of Triticeae have been generated, making it necessary to update and standardize TE annotations and nomenclature. In this study we propose standard procedures for these tasks, based on structure, nucleic acid and protein sequence homologies. We report statistical analyses of TE composition and distribution in large blocks of genomic sequences from wheat and barley. Altogether, 3.8 Mb of wheat sequence available in the databases was analyzed or re-analyzed, and compared with 1.3 Mb of re-annotated genomic sequences from barley. The wheat sequences were relatively gene-rich (one gene per 23.9 kb), although wheat gene-derived sequences represented only 7.8% (159 elements) of the total, while the remainder mainly comprised coding sequences found in TEs (54.7%, 751 elements). Class I elements [mainly long terminal repeat (LTR) retrotransposons] accounted for the major proportion of TEs, in terms of sequence length as well as element number (83.6% and 498, respectively). In addition, we show that the gene-rich sequences of wheat genome A seem to have a higher TE content than those of genomes B and D, or of barley gene-rich sequences. Moreover, among the various TE groups, MITEs were most often associated with genes: 43.1% of MITEs fell into this category. Finally, the TRIM and copia elements were shown to be the most active TEs in the wheat genome. The implications of these results for the evolution of diploid and polyploid wheat species are discussed.

Characterization of Stowaway MITEs in pea (Pisum sativum L.) and identification of their potential master elements

We have investigated miniature inverted-repeat transposable elements (MITEs) of the Stowaway family and corresponding Mariner-like master elements that could potentially facilitate their mobilization in the genome of the garden pea (Pisum sativum L.). The population of pea Stowaway MITEs consists of 103-104 copies dispersed in the genome. Judging from a sequence analysis of 17 isolated Stowaway elements and their flanking genomic regions, the elements are relatively uniform in size and sequence and occur in the vicinity of genes as well as within repetitive sequences. Insertional polymorphism of several elements was detected among various Pisum accessions, suggesting they were still transpositionally active during diversification of these taxa. The identification of several Mariner-like elements (MLEs) harboring intact open reading frames, capable of encoding a transposase, further supports a recent mobilization of the Stowaway elements. Using transposase-coding sequences as a hybridization probe, we estimated that there are about 50 MLE sequences in the pea genome. Among the 5 elements sequenced, 3 distinct subfamilies showing mutual similarities within their transposase-coding regions, but otherwise diverged in sequence, were distinguished and designated as Psmar-1 to Psmar-3. The terminal inverted repeats (TIRs) of these MLE subfamilies differed in their homology to the TIRs of Stowaway MITEs. The homlogy ranged from 9 bp in Psmar-3 to 30 bp in Psmar-1, which corresponds to the complete Stowaway TIR sequence. Based on this feature, the Psmar-1 elements are believed to be the most likely candidates for the master elements of the Stowaway MITEs in pea.

Pack-MULEs: theft on a massive scale

It has been known for some time that plant transposons can capture and mobilize cellular genes. Recent work by Jiang and coworkers((1)) has revealed that this process has happened on a massive scale. They found that portions of more than 1000 genes in rice have been captured and mobilized by members of the MULE family of transposons. In rice, and perhaps other plants as well, it appears that thousands of genes and portions of genes have been duplicated, transposed and rearranged. These results have fascinating implications for our understanding of the mode and tempo of gene evolution in plants.

The map-based sequence of the rice genome

Rice, one of the world's most important food plants, has important syntenic relationships with the other cereal species and is a model plant for the grasses. Here we present a map-based, finished quality sequence that covers 95% of the 389 Mb genome, including virtually all of the euchromatin and two complete centromeres. A total of 37,544 non-transposable-element-related protein-coding genes were identified, of which 71% had a putative homologue in Arabidopsis. In a reciprocal analysis, 90% of the Arabidopsis proteins had a putative homologue in the predicted rice proteome. Twenty-nine per cent of the 37,544 predicted genes appear in clustered gene families. The number and classes of transposable elements found in the rice genome are consistent with the expansion of syntenic regions in the maize and sorghum genomes. We find evidence for widespread and recurrent gene transfer from the organelles to the nuclear chromosomes. The map-based sequence has proven useful for the identification of genes underlying agronomic traits. The additional single-nucleotide polymorphisms and simple sequence repeats identified in our study should accelerate improvements in rice production.

The oryza map alignment project: the golden path to unlocking the genetic potential of wild rice species

The wild species of the genus Oryza offer enormous potential to make a significant impact on agricultural productivity of the cultivated rice species Oryza sativa and Oryza glaberrima. To unlock the genetic potential of wild rice we have initiated a project entitled the 'Oryza Map Alignment Project' (OMAP) with the ultimate goal of constructing and aligning BAC/STC based physical maps of 11 wild and one cultivated rice species to the International Rice Genome Sequencing Project's finished reference genome--O. sativa ssp. japonica c. v. Nipponbare. The 11 wild rice species comprise nine different genome types and include six diploid genomes (AA, BB, CC, EE, FF and GG) and four tetrapliod genomes (BBCC, CCDD, HHKK and HHJJ) with broad geographical distribution and ecological adaptation. In this paper we describe our strategy to construct robust physical maps of all 12 rice species with an emphasis on the AA diploid O. nivara--thought to be the progenitor of modern cultivated rice.

Abundance, distribution and dynamics of retrotransposable elements and transposons: similarities and differences

Retrotransposable elements and transposons are generally both found in most eukaryotes. These two classes of elements are usually distinguished on the basis of their differing mechanisms of transposition. However, their respective frequencies, their intragenomic dynamics and distributions, and the frequencies of their horizontal transfer from one species to another can also differ. The main objective of this review is to compare these two types of elements from a new perspective, using data provided by genome sequencing projects and relating this to the theoretical and observed dynamics. It is shown that the traditional division into two classes, based on the transposition mechanisms, becomes less obvious when other factors are taken into consideration. A great diversity in distribution and dynamics within each class is observed. In contrast, the impact on and the interactions with the genome can show striking similarities between families of the two classes.

A tiling microarray expression analysis of rice chromosome 4 suggests a chromosome-level regulation of transcription

The complete genome sequence of cultivated rice (Oryza sativa) provides an unprecedented opportunity to understand the biology of this model cereal. An essential and necessary step in this effort is the determination of the coding information and expression patterns of each sequenced chromosome. Here, we report an analysis of the transcriptional activity of rice chromosome 4 using a tiling path microarray based on PCR-generated genomic DNA fragments. Six representative rice organ types were examined using this microarray to catalog the transcribed regions of rice chromosome 4 and to reveal organ- and developmental stage-specific transcription patterns. This analysis provided expression support for 82% of the gene models in the chromosome. Transcriptional activities in 1643 nonannotated regions were also detected. Comparison with cytologically defined chromatin features indicated that in juvenile-stage rice the euchromatic region is more actively transcribed than is the transposon-rich heterochromatic portion of the chromosome. Interestingly, increased transcription of transposon-related gene models in certain heterochromatic regions was observed in mature-stage rice organs and in suspension-cultured cells. These results suggest a close correlation between transcriptional activity and chromosome organization and the developmental regulation of transcription activity at the chromosome level.

Identification of an active transposon in intact rice plants

A transposable element that is active in intact plants has been identified in rice (Oryza sativa L.). The 607-bp element itself, termed nonautonomous DNA-based active rice transposon (nDart), has no coding capacity. It was found inserted in the gene encoding Mg-protoporphyrin IX methyltransferase in a chlorophyll-deficient albino mutant isolated from backcross progeny derived from a cross between wild-type japonica varieties. The nDart has 19-bp terminal inverted repeats (TIRs) and, when mobilized, generates an 8-bp target-site duplication (TSD). At least 13 nDart elements were identified in the genome sequence of the japonica cultivar Nipponbare. Database searches identified larger elements, termed DNA-based active rice transposon (Dart) that contained one ORF for a protein that contains a region with high similarity to the hAT dimerization motif. Dart shares several features with nDart, including identical TIRs, similar subterminal sequences and the generation of an 8-bp TSD. These shared features indicate that the nonautonomous element nDart is an internal deletion derivative of the autonomous element Dart. We conclude that these active transposon systems belong to the hAT superfamily of class II transposons. Because the transposons are active in intact rice plants, they should be useful tools for tagging genes in studies of functional genomics.

Rim 2/Hipa CACTA transposon display: a new genetic marker technique in Oryza species

Background

Transposons constitute the major fractions of repetitive sequences in eukaryotes, and have been crucial in the shaping of current genomes. Transposons are generally divided into two classes according to the mechanism underlying their transposition: RNA intermediate class 1 and DNA intermediate class 2. CACTA is a class 2 transposon superfamily, which is found exclusively in plants. As some transposons, including the CACTA superfamily, are highly abundant in plant species, and their nucleotide sequences are highly conserved within a family, they can be utilized as genetic markers, using a slightly modified version of the conventional AFLP protocol. Rim2 /Hipa is a CACTA transposon family having 16 bp consensus TIR sequences to be present in high copy numbers in rice genome. This research was carried out in order to develop a Rim2/Hipa CACTA-AFLP or Rim2/Hipa CACTA-TD (transposon display, hereafter Rim2/Hipa-TD) protocol for the study of genetic markers in map construction and the study of genetic diversity in rice.

Results

Rim2/Hipa-TD generated ample polymorphic profiles among the different rice accessions, and the amplification profiles were highly reproducible between different thermocyclers and Taq polymerases. These amplification profiles allowed for clear distinction between two different ecotypes, Japonica and Indica, of Oryza sativa. In the analysis of RIL populations, the Rim2/Hipa-TD markers were found to be segregated largely in a dominant manner, although in a few cases, non-parental bands were observed in the segregating populations. Upon linkage analysis, the Rim2/Hipa-TD markers were found to be distributed in the regions proximal to the centromeres of the chromosomes. The distribution of the Rim2/Hipa CACTA elements was surveyed in 15 different Oryza species via Rim2/Hipa-TD. While Rim2/Hipa-TD yielded ample amplification profiles between 100 to 700 bp in the AA diploid Oryza species, other species having BB, CC, EE, BBCC and CCDD, profiles demonstrated that most of the amplified fragments were larger than 400 bp, and that our methods were insufficient to clearly distinguish between these fragments. However, the overall amplification profiles between species in the Oryza genus were fully distinct. Phenetic relationships among the AA diploid Oryza species, as evidenced by the Rim2/Hipa-TD markers, were matched with their geographical distributions.

Conclusion

The abundance of the Rim2/Hipa TIR sequences is very informative since the Rim2/Hipa-TD produced high polymorphic profiles with ample reproducibility within a species as well as between species in the Oryza genus. Therefore, Rim2/Hipa-TD markers can be useful in the development of high-density of genetic map around the centromeric regions. Rim2/Hipa-TD may also prove useful in evaluations of genetic variation and species relationships in the Oryza species.

Plant MITEs: useful tools for plant genetics and genomics

MITEs (Miniature inverted-repeat transposable elements) are reminiscence of non-autonomous DNA (class II) elements, which are distinguished from other transposable elements by their small size, short terminal inverted repeats (TIRs), high copy numbers, genic preference, and DNA sequence identity among family members. Although MITEs were first discovered in plants and still actively reshaping genomes, they have been isolated from a wide range of eukaryotic organisms. MITEs can be divided into Tourist-like, Stowaway-like, and pogo-like groups, according to similarities of their TIRs and TSDs (target site duplications). In despite of several models to explain the origin and amplification of MITEs, their mechanisms of transposition and accumulation in eukaryotic genomes remain poorly understood owing to insufficient experimental data. The unique properties of MITEs have been exploited as useful genetic tools for plant genome analysis. Utilization of MITEs as effective and informative genomic markers and potential application of MITEs in plants systematic, phylogenetic, and genetic studies are discussed.

Mutator-like element in the yeast Yarrowia lipolytica displays multiple alternative splicings

A new type of DNA transposon, Mutyl, has been identified in the sequenced genome of the yeast Yarrowia lipolytica. This transposon is 7,413 bp long and carries two open reading frames (ORFs) which potentially encode proteins of 459 and 1,178 amino acids, respectively. Whereas the first ORF shows no significant homology to previously described proteins, the second ORF shows sequence similarities with various Mutator-like element (MULE)-encoded transposases, including the bacterial transposase signature sequence. Other MULE features shared by Mutyl include a zinc finger motif in the putative transposase, a 22-bp-long imperfect inverted repeat at each end, and a 9- to 10-bp duplication of its target site in the chromosome. Of the five copies of Mutyl present in the genome, one has a deletion of the first 8 bases, and the others are full length with a single base change in one element. The first potential gene of Mutyl, mutB, was shown to be expressed in exponentially growing cells. Its sequence contains a predicted intron with two 5' splice sites, a single branch point, and two 3' splice sites. Its mRNA is alternatively spliced, as judged by reverse transcription-PCR, and generates four mRNAs corresponding to protein-coding sequences of 128, 156, 161, and 190 amino acids. Of the three distinct lineages characterized in Y. lipolytica, strains from the German lineage and the French lineage do not carry Mutyl. A study of the distribution of Mutyl in strains of the French lineage evidenced a recent transposition event. Taken together, these results indicate that Mutyl is still active.

The Genomes of Oryza sativa: a history of duplications

We report improved whole-genome shotgun sequences for the genomes of indica and japonica rice, both with multimegabase contiguity, or almost 1,000-fold improvement over the drafts of 2002. Tested against a nonredundant collection of 19,079 full-length cDNAs, 97.7% of the genes are aligned, without fragmentation, to the mapped super-scaffolds of one or the other genome. We introduce a gene identification procedure for plants that does not rely on similarity to known genes to remove erroneous predictions resulting from transposable elements. Using the available EST data to adjust for residual errors in the predictions, the estimated gene count is at least 38,000–40,000. Only 2%–3% of the genes are unique to any one subspecies, comparable to the amount of sequence that might still be missing. Despite this lack of variation in gene content, there is enormous variation in the intergenic regions. At least a quarter of the two sequences could not be aligned, and where they could be aligned, single nucleotide polymorphism (SNP) rates varied from as little as 3.0 SNP/kb in the coding regions to 27.6 SNP/kb in the transposable elements. A more inclusive new approach for analyzing duplication history is introduced here. It reveals an ancient whole-genome duplication, a recent segmental duplication on Chromosomes 11 and 12, and massive ongoing individual gene duplications. We find 18 distinct pairs of duplicated segments that cover 65.7% of the genome; 17 of these pairs date back to a common time before the divergence of the grasses. More important, ongoing individual gene duplications provide a never-ending source of raw material for gene genesis and are major contributors to the differences between members of the grass family.

Comparative genome sequencing of indica and japonica rice reveals that duplication of genes and genomic regions has played a major part in the evolution of grass genomes

From mapping to sequencing, post-sequencing and beyond

The Rice Genome Research Program (RGP) in Japan has been collaborating with the international community in elucidating a complete high-quality sequence of the rice genome. As the pioneer in large-scale analysis of the rice genome, the RGP has successfully established the fundamental tools for genome research such as a genetic map, a yeast artificial chromosome (YAC)-based physical map, a transcript map and a phage P1 artificial chromosome (PAC)/bacterial artificial chromosome (BAC) sequence-ready physical map, which serve as common resources for genome sequencing. Among the 12 rice chromosomes, the RGP is in charge of sequencing six chromosomes covering 52% of the 390 Mb total length of the genome. The contribution of the RGP to the realization of decoding the rice genome sequence with high accuracy and deciphering the genetic information in the genome will have a great impact in understanding the biology of the rice plant that provides a major food source for almost half of the world's population. A high-quality draft sequence (phase 2) was completed in December 2002. Since then, much of the finished quality sequence (phase 3) has become available in public databases. With the completion of sequencing in December 2004, it is expected that the genome sequence would facilitate innovative research in functional and applied genomics. A map-based genome sequence is indispensable for further improvement of current rice varieties and for development of novel varieties carrying agronomically important traits such as high yield potential and tolerance to both biotic and abiotic stresses. In addition to genome sequencing, various related projects have been initiated to generate valuable resources, which could serve as indispensable tools in clarifying the structure and function of the rice genome. These resources have been made available to the scientific community through the Rice Genome Resource Center (RGRC) of the National Institute of Agrobiological Sciences (NIAS) to enable rapid progress in research that will lead to thorough understanding of the rice plant. As the next trend in rice genome research will focus on determining the function of about 40,000-50,000 genes predicted in the genome as well as applying various genomics tools in rice breeding, an unlimited access to rice DNA and seed stocks will provide a broad community of scientists with the necessary materials for formulating new concepts, developing innovative research and making new scientific discoveries in rice genomics.

Expression of a xyloglucan endotransglucosylase/hydrolase gene, Mt-XTH1, from Medicago truncatula is induced systemically in mycorrhizal roots

Xyloglucan endotransglucosylase/hydrolases (XTH) are enzymes that catalyze the hydrolysis and transglycosylation of xyloglucan polymers in plant cell walls. Previously, we isolated a cDNA from mycorrhizal roots of Medicago truncatula that is predicted to encode an XTH [van Buuren, M.L., Maldonado-Mendoza, I.E., Trieu, A.T., Blaylock, L.A., Harrison, M.J., 1999. Novel genes induced during an arbuscular mycorrhizal (AM) symbiosis between M. truncatula and G. versiforme. Mol. Plant-Microb. Interact. 12, 171-181.]. Here, we identified the corresponding XTH gene, designated Mt-XTH1. The Mt-XTH1 gene contains four exons separated by three introns and resides on a 15-kb Xba1 fragment adjacent to a second XTH gene designated Mt-XTH2. Mt-XTH2 shares the same exon-intron structure as Mt-XTH1. Exons 2, 3 and 4 and introns 1 and 2 are identical to Mt-XTH1, while exon 1 and intron 3 are divergent, both in sequence and in length. Mt-XTH1 is induced following colonization of the roots by AM fungi but does not respond to changes in phosphate status. Analysis of transgenic roots expressing an Mt-XTH1 promoterColon, two colonsuidA fusion revealed that the Mt-XTH1 promoter directs expression in cells throughout the root system with significantly higher levels of activity in mycorrhizal roots. Mt-XTH1 expression is elevated not only in the regions of the roots colonized by the fungus, but also at sites distal to the infected regions. These expression patterns are consistent with activation in response to a systemic signal.

Sequence composition, organization, and evolution of the core Triticeae genome

We investigated the composition and the basis of genome expansion in the core Triticeae genome using Aegilops tauschii, the D-genome donor of bread wheat. We sequenced an unfiltered genomic shotgun (trs) and a methylation-filtration (tmf) library of A. tauschii, and analyzed wheat expressed sequence tags (ESTs) to estimate the expression of genes and transposable elements (TEs). The sampled D-genome sequences consisted of 91.6% repetitive elements, 2.5% known genes, and 5.9% low-copy sequences of unknown function. TEs constituted 68.2% of the D-genome compared with 50% in maize and 14% in rice. The DNA transposons constituted 13% of the D-genome compared with 2% in maize. TEs were methylated unevenly within and among elements and families, and most were transcribed which contributed to genome expansion in the core Triticeae genome. The copy number of a majority of repeat families increased gradually following polyploidization. Certain TE families occupied discrete chromosome territories. Nested insertions and illegitimate recombination occurred extensively between the TE families, and a majority of the TEs contained internal deletions. The GC content varied significantly among the three sequence sets examined ranging from 42% in tmf to 46% in trs and 52% in the EST. Based on enrichment of genic sequences, methylation-filtration offers one option, although not as efficient as in maize, for isolating gene-rich regions from the large genome of wheat.

Reconstruction of putative DNA virus from endogenous rice tungro bacilliform virus-like sequences in the rice genome: implications for integration and evolution

Background

Plant genomes contain various kinds of repetitive sequences such as transposable elements, microsatellites, tandem repeats and virus-like sequences. Most of them, with the exception of virus-like sequences, do not allow us to trace their origins nor to follow the process of their integration into the host genome. Recent discoveries of virus-like sequences in plant genomes led us to set the objective of elucidating the origin of the repetitive sequences. Endogenous rice tungro bacilliform virus (RTBV)-like sequences (ERTBVs) have been found throughout the rice genome. Here, we reconstructed putative virus structures from RTBV-like sequences in the rice genome and characterized to understand evolutionary implication, integration manner and involvements of endogenous virus segments in the corresponding disease response.

Results

We have collected ERTBVs from the rice genomes. They contain rearranged structures and no intact ORFs. The identified ERTBV segments were shown to be phylogenetically divided into three clusters. For each phylogenetic cluster, we were able to make a consensus alignment for a circular virus-like structure carrying two complete ORFs. Comparisons of DNA and amino acid sequences suggested the closely relationship between ERTBV and RTBV. The Oryza AA-genome species vary in the ERTBV copy number. The species carrying low-copy-number of ERTBV segments have been reported to be extremely susceptible to RTBV. The DNA methylation state of the ERTBV sequences was correlated with their copy number in the genome.

Conclusions

These ERTBV segments are unlikely to have functional potential as a virus. However, these sequences facilitate to establish putative virus that provided information underlying virus integration and evolutionary relationship with existing virus. Comparison of ERTBV among the Oryza AA-genome species allowed us to speculate a possible role of endogenous virus segments against its related disease.

Genomic organization of the complex alpha-gliadin gene loci in wheat

To better understand the molecular evolution of the large alpha-gliadin gene family, a half-million bacterial artificial chromosome (BAC) library clones from tetraploid durum wheat, Triticum turgidum ssp. durum (2n = 4x = 28, genome AB), were screened for large genomic segments carrying the alpha-gliadin genes of the Gli-2 loci on the group 6 homoeologous chromosomes. The resulting 220 positive BAC clones--each containing between one and four copies of alpha-gliadin sequences--were fingerprinted for contig assembly to produce contiguous chromosomal regions covering the Gli-2 loci. While contigs consisting of as many as 21 BAC clones and containing up to 17 alpha-gliadin genes were formed, many BAC clones remained as singletons. The accuracy of the order of BAC clones in the contigs was verified by Southern hybridization analysis of the BAC fingerprints using an alpha-gliadin probe. These results indicate that alpha-gliadin genes are not evenly dispersed in the Gli-2 locus regions. Hybridization of these BACs with probes for long terminal repeat retrotransposons was used to determine the abundance and distribution of repetitive DNA in this region. Sequencing of BAC ends indicated that 70% of the sequences were significantly similar to different classes of retrotransposons, suggesting that these elements are abundant in this region. Several mechanisms underlying the dynamic evolution of the Gli-2 loci are discussed.

BAC end sequences and a physical map reveal transposable element content and clustering patterns in the genome of Magnaporthe grisea

Transposable elements (TEs) are viewed as major contributors to the evolution of fungal genomes. Genomic resources such as BAC libraries are an underutilized resource for studying genome-wide TE distribution. Using the BAC end sequences and physical map that are available for the rice blast fungus, Magnaporthe grisea, we describe a likelihood ratio test designed to identify clustering of TEs in the genome. A significant variation in the distribution of three TEs, MAGGY, MGL, and Pot2 was observed among the fingerprint contigs of the physical map. We utilized a draft sequence of M. grisea chromosome 7 to validate our results and found a similar pattern of clustering. By examining individual BAC end sequences, we found evidence for 11 unique integrations of MAGGY or MGL into Pot2 but no evidence for the reciprocal integration of Pot2 into another TE. This suggests that: (a) the presence of Pot2 in the genome predates that of the other TEs, (b) Pot2 was less transpositionally active than other TEs, or (c) that MAGGY and MGL have integration site preference for Pot2. High transition/transversion mutation ratios as well as bias in transition site context was observed in MAGGY and MGL elements, but not in Pot2 elements. These features are consistent with the effects of a Repeat-Induced Point (RIP) mutation-like process occurring in MAGGY and MGL elements. This study illustrates the general utility of a physical map and BAC end sequences for the study of genome-wide repetitive DNA content and organization.

Construction of a BAC library of Korean ginseng and initial analysis of BAC-end sequences

We estimated the genome size of Korean ginseng (Panax ginseng C.A. Meyer), a medicinal herb, constructed a HindIII BAC library, and analyzed BAC-end sequences to provide an initial characterization of the library. The 1C nuclear DNA content of Korean ginseng was estimated to be 3.33 pg (3.12 x 10(3) Mb). The BAC library consists of 106,368 clones with an average size of 98.61 kb, amounting to 3.34 genome equivalents. Sequencing of 2167 BAC clones generated 2492 BAC-end sequences with an average length of 400 bp. Analysis using BLAST and motif searches revealed that 10.2%, 20.9% and 3.8% of the BAC-end sequences contained protein-coding regions, transposable elements and microsatellites, respectively. A comparison of the functional categories represented by the protein-coding regions found in BAC-end sequences with those of Arabidopsis revealed that proteins pertaining to energy metabolism, subcellular localization, cofactor requirement and transport facilitation were more highly represented in the P. ginseng sample. In addition, a sequence encoding a glucosyltransferase-like protein implicated in the ginsenoside biosynthesis pathway was also found. The majority of the transposable element sequences found belonged to the gypsy type (67.6%), followed by copia (11.7%) and LINE (8.0%) retrotransposons, whereas DNA transposons accounted for only 2.1% of the total in our sequence sample. Higher levels of transposable elements than protein-coding regions suggest that mobile elements have played an important role in the evolution of the genome of Korean ginseng, and contributed significantly to its complexity. We also identified 103 microsatellites with 3-38 repeats in their motifs. The BAC library and BAC-end sequences will serve as a useful resource for physical mapping, positional cloning and genome sequencing of P. ginseng.

Jittery, a Mutator distant relative with a paradoxical mobile behavior: excision without reinsertion

The unstable mutation bz-m039 arose in a maize (Zea mays) stock that originated from a plant infected with barley stripe mosaic virus. The instability of the mutation is caused by a 3.9-kb mobile element that has been named Jittery (Jit). Jit has terminal inverted repeats (TIRs) of 181 bp, causes a 9-bp direct duplication of the target site, and appears to excise autonomously. It is predicted to encode a single 709-amino acid protein, JITA, which is distantly related to the MURA transposase protein of the Mutator system but is more closely related to the MURA protein of Mutator-like elements (MULEs) from Arabidopsis thaliana and rice (Oryza sativa). Like MULEs, Jit resembles Mutator in the length of the element's TIRs, the size of the target site duplication, and in the makeup of its transposase but differs from the autonomous element Mutator-Don Robertson in that it encodes a single protein. Jit also differs from Mutator elements in the high frequency with which it excises to produce germinal revertants and in its copy number in the maize genome: Jit-like TIRs are present at low copy number in all maize lines and teosinte accessions examined, and JITA sequences occur in only a few maize inbreds. However, Jit cannot be considered a bona fide transposon in its present host line because it does not leave footprints upon excision and does not reinsert in the genome. These unusual mobile element properties are discussed in light of the structure and gene organization of Jit and related elements.

Genome-wide comparative analysis of the transposable elements in the related species Arabidopsis thaliana and Brassica oleracea

Transposable elements (TEs) are the major component of plant genomes where they contribute significantly to the >1,000-fold genome size variation. To understand the dynamics of TE-mediated genome expansion, we have undertaken a comparative analysis of the TEs in two related organisms: the weed Arabidopsis thaliana (125 megabases) and Brassica oleracea ( approximately 600 megabases), a species with many crop plants. Comparison of the whole genome sequence of A. thaliana with a partial draft of B. oleracea has permitted an estimation of the patterns of TE amplification, diversification, and loss that has occurred in related species since their divergence from a common ancestor. Although we find that nearly all TE lineages are shared, the number of elements in each lineage is almost always greater in B. oleracea. Class 1 (retro) elements are the most abundant TE class in both species with LTR and non-LTR elements comprising the largest fraction of each genome. However, several families of class 2 (DNA) elements have amplified to very high copy number in B. oleracea where they have contributed significantly to genome expansion. Taken together, the results of this analysis indicate that amplification of both class 1 and class 2 TEs is responsible, in part, for B. oleracea genome expansion since divergence from a common ancestor with A. thaliana. In addition, the observation that B. oleracea and A. thaliana share virtually all TE lineages makes it unlikely that wholesale removal of TEs is responsible for the compact genome of A. thaliana.

Sequencing the maize genome

Sequencing of complex genomes can be accomplished by enriching shotgun libraries for genes. In maize, gene-enrichment by copy-number normalization (high C(0)t) and methylation filtration (MF) have been used to generate up to two-fold coverage of the gene-space with less than 1 million sequencing reads. Simulations using sequenced bacterial artificial chromosome (BAC) clones predict that 5x coverage of gene-rich regions, accompanied by less than 1x coverage of subclones from BAC contigs, will generate high-quality mapped sequence that meets the needs of geneticists while accommodating unusually high levels of structural polymorphism. By sequencing several inbred strains, we propose a strategy for capturing this polymorphism to investigate hybrid vigor or heterosis.

Structural and functional analysis of rice genome

Rice is an excellent system for plant genomics as it represents a modest size genome of 430 Mb. It feeds more than half the population of the world. Draft sequences of the rice genome, derived by whole-genome shotgun approach at relatively low coverage (4-6 X), were published and the International Rice Genome Sequencing Project (IRGSP) declared high quality (>10 X), genetically anchored, phase 2 level sequence in 2002. In addition, phase 3 level finished sequence of chromosomes 1, 4 and 10 (out of 12 chromosomes of rice) has already been reported by scientists from IRGSP consortium. Various estimates of genes in rice place the number at >50,000. Already, over 28,000 full-length cDNAs have been sequenced, most of which map to genetically anchored genome sequence. Such information is very useful in revealing novel features of macro- and micro-level synteny of rice genome with other cereals. Microarray analysis is unraveling the identity of rice genes expressing in temporal and spatial manner and should help target candidate genes useful for improving traits of agronomic importance. Simultaneously, functional analysis of rice genome has been initiated by marker-based characterization of useful genes and employing functional knock-outs created by mutation or gene tagging. Integration of this enormous information is expected to catalyze tremendous activity on basic and applied aspects of rice genomics.

Using rice to understand the origin and amplification of miniature inverted repeat transposable elements (MITEs)

Recent studies of rice miniature inverted repeat transposable elements (MITEs), largely fueled by the availability of genomic sequence, have provided answers to many of the outstanding questions regarding the existence of active MITEs, their source of transposases (TPases) and their chromosomal distribution. Although many questions remain about MITE origins and mode of amplification, data accumulated over the past two years have led to the formulation of testable models.

Evolutionary history of Oryza sativa LTR retrotransposons: a preliminary survey of the rice genome sequences

Background

LTR Retrotransposons transpose through reverse transcription of an RNA intermediate and are ubiquitous components of all eukaryotic genomes thus far examined. Plant genomes, in particular, have been found to be comprised of a remarkably high number of LTR retrotransposons. There is a significant body of direct and indirect evidence that LTR retrotransposons have contributed to gene and genome evolution in plants.

Results

To explore the evolutionary history of long terminal repeat (LTR) retrotransposons and their impact on the genome of Oryza sativa, we have extended an earlier computer-based survey to include all identifiable full-length, fragmented and solo LTR elements in the rice genome database as of April 2002. A total of 1,219 retroelement sequences were identified, including 217 full-length elements, 822 fragmented elements, and 180 solo LTRs. In order to gain insight into the chromosomal distribution of LTR-retrotransposons in the rice genome, a detailed examination of LTR-retrotransposon sequences on Chromosome 10 was carried out. An average of 22.3 LTR-retrotransposons per Mb were detected in Chromosome 10.

Conclusions

Gypsy-like elements were found to be >4 × more abundant than copia-like elements. Eleven of the thirty-eight investigated LTR-retrotransposon families displayed significant subfamily structure. We estimate that at least 46.5% of LTR-retrotransposons in the rice genome are older than the age of the species (< 680,000 years). LTR-retrotransposons present in the rice genome range in age from those just recently inserted up to nearly 10 million years old. Approximately 20% of LTR retrotransposon sequences lie within putative genes. The distribution of elements across chromosome 10 is non-random with the highest density (48 elements per Mb) being present in the pericentric region.

Coevolution of viruses with hosts and vectors and possible paleontology

The coevolution of viruses with their hosts and vectors depends on the evolution of the hosts and vectors coupled with factors involved in virus evolution. The long-term perspective involves the origin of life forms, the evolution of host and vector (especially arthropods) kingdoms and families, and changes in biological diversity induced mainly by the last five great extinctions. In the medium term, the diversification of hosts and vectors is important, and in the short term, recent events, especially humans, have had a great impact on virus coevolution. As there are few, if any, examples of conventional fossils of viruses, evidence for their evolution related to host and vector evolution is being found from other sources, especially virus-induced cellular structures and recent developments in molecular biology. Recognizing these other sources is becoming important for paleontologists gaining an understanding of the influence that viruses have had on the development of higher organisms.

Sources and predictors of resolvable indel polymorphism assessed using rice as a model

The principal sources of genetic variation that can be assayed with restriction enzymes are base substitutions and insertions/deletions (indels). The likelihood of detecting indels as restriction fragment length polymorphisms (RFLPs) is determined by the size and frequency of the indels, and the ability to resolve small indels as RFLPs is limited by the distribution of restriction fragment sizes. In this study, we use aligned sequences from the indica and japonica subspecies of rice ( Oryza sativaL.) to quantify and compare the ability of restriction enzymes to detect indels. We look specifically at two abundant transposable element-derived indel sources: miniature inverted repeat transposable elements (MITEs) and long terminal repeat (LTR) retroelements. From this analysis we conclude that indels rather than base substitutions are the prevailing source of the polymorphism detected in rice. We show that, although MITE derived indels are more abundant than LTR-retroelement derived indels, LTR-retroelements have a greater capacity to generate visible restriction fragment length polymorphism because of their larger size. We find that the variation in the detectability of indels among restriction enzymes can be explained by differences in the frequency and dispersion of their restriction sites in the genome. The parameters that describe the fragment size distributions obtained with the restriction enzymes are highly correlated across the sequenced genomes of rice, Arabidopsis and human, with the exception of some extreme deviations in frequency for particular recognition sequences corresponding to variations in the levels and modes of DNA methylation in the three disparate organisms. Thus, we can predict the relative ability of a restriction enzyme to detect indels derived from a specific source based on the distribution of restriction fragment sizes, even when this is estimated for a distantly related genome.

Maize Genome Sequencing by Methylation Filtration

Gene enrichment strategies offer an alternative to sequencing large and repetitive genomes such as that of maize. We report the generation and analysis of nearly 100,000 undermethylated (or methylation filtration) maize sequences. Comparison with the rice genome reveals that methylation filtration results in a more comprehensive representation of maize genes than those that result from expressed sequence tags or transposon insertion sites sequences. About 7% of the repetitive DNA is unmethylated and thus selected in our libraries, but potentially active transposons and unmethylated organelle genomes can be identified. Reverse transcription polymerase chain reaction can be used to finish the maize transcriptome.

Physical maps and recombination frequency of six rice chromosomes

We constructed physical maps of rice chromosomes 1, 2, and 6-9 with P1-derived artificial chromosome (PAC) and bacterial artificial chromosome (BAC) clones. These maps, with only 20 gaps, cover more than 97% of the predicted length of the six chromosomes. We submitted a total of 193 Mbp of non-overlapping sequences to public databases. We analyzed the DNA sequences of 1316 genetic markers and six centromere-specific repeats to facilitate characterization of chromosomal recombination frequency and of the genomic composition and structure of the centromeric regions. We found marked changes in the relative recombination rate along the length of each chromosome. Chromosomal recombination at the centromere core and surrounding regions on the six chromosomes was completely suppressed. These regions have a total physical length of about 23 Mbp, corresponding to 11.4% of the entire size of the six chromosomes. Chromosome 6 has the longest quiescent region, with about 5.6 Mbp, followed by chromosome 8, with quiescent region about half this size. Repetitive sequences accounted for at least 40% of the total genomic sequence on the partly sequenced centromeric region of chromosome 1. Rice CentO satellite DNA is arrayed in clusters and is closely associated with the presence of Centromeric Retrotransposon of Rice (CRR)- and RIce RetroElement 7 (RIRE7)-like retroelement sequences. We also detected relatively small coldspot regions outside the centromeric region; their repetitive content and gene density were similar to those of regions with normal recombination rates. Sequence analysis of these regions suggests that either the amount or the organization patterns of repetitive sequences may play a role in the inactivation of recombination.

Genomic structure of the NtPDR1 gene, harboring the two miniature inverted-repeat transposable elements, NtToya1 and NtStowaway101

Here we report the genomic structure including the promoter sequence and coding region of NtPDR1 (Nicotiana tabacum Pleiotropic Drug Resistance 1), which is an elicitor-responsive gene encoding an ATP binding cassette (ABC) transporter that might be involved in the defense response in tobacco, as we reported recently. The NtPDR1 gene consists of 20 exons and 19 introns. Among the introns, the first and fifth are much larger than the others and harbor typical miniature inverted-repeat transposable elements (MITEs). One of the MITE elements in the first intron, termed NtToya1, belongs to the Toya family that was recently described in rice, while the other element in the fifth intron, termed NtStowaway101, shows high homology with the Stowaway elements of the IS630-Tc1-mariner family. Many of the genes we found to harbor Toya and Stowaway elements in Nicotiana species by BLAST search are also involved in stress responses or plant-pathogen interactions. The existence of putative cis-elements (a GCC box, three W boxes, and several JA-responsive elements) in the promoter region supports our previous finding that this gene is strongly inducible by elicitation and methyljasmonate, and that this ABC transporter might be essential for plant defense responses. Furthermore, Southern blot analysis and PCR amplification of the introns harboring the MITE-like elements from genomic DNA of three Nicotiana species suggests that NtPDR1 originated from N. sylvestris.

Genomic characterization of Rim2/Hipa elements reveals a CACTA-like transposon superfamily with unique features in the rice genome

The availability of huge amounts of rice genome sequence now permits large-scale analysis of the structure and molecular characteristics of the previously identified transposase-encoding Rim2 (also called Hipa) element, which is transcriptionally activated by infection with the fungal pathogen Magnaporthe grisea and by treatment with the corresponding fungal elicitor. Based on genomic cloning and data mining from 230 Mb of rice genome sequence, 347 Rim2 elements, with an average size of 5.8 kb, were identified. This indicates that an estimated total of 600-700 Rim2 elements are present in the whole genome. Rim2 insertions occur non-randomly on the chromosomes, as visualized by fluorescence in situ hybridization. The elements harbor 16-bp terminal inverted repeats with the core sequence CACTG, 16-bp sub-terminal repeats, internal variable regions, 3-bp target sequence duplications in the flanking regions, and genes coding for Rim2 proteins (the putative transposase) and hydroxyproline-rich glycoproteins. High levels of insertion into genic regions are observed for members of this family, and the transposition history of the family can be deduced from the high level of shared sequences and analysis of repeat target sites of the elements. Phylogenetic analysis indicates that the putative RIM2 proteins fall into a subgroup distinct from the TNP2-like subgroup of transposases. Southern hybridization with genomic DNA from monocotyledonous and dicotyledonous plants demonstrates that the RIM2-coding sequence is unique to the Oryza genome. Our results demonstrate that the Rim2 elements from rice belong to a distinct superfamily of CACTA-like elements with evolutionary diversity.

Locating sequence on FPC maps and selecting a minimal tiling path

This study discusses three software tools, the first two aid in integrating sequence with an FPC physical map and the third automatically selects a minimal tiling path given genomic draft sequence and BAC end sequences. The first tool, FSD (FPC Simulated Digest), takes a sequenced clone and adds it back to the map based on a fingerprint generated by an in silico digest of the clone. This allows verification of sequenced clone positions and the integration of sequenced clones that were not originally part of the FPC map. The second tool, BSS (Blast Some Sequence), takes a query sequence and positions it on the map based on sequence associated with the clones in the map. BSS has multiple uses as follows: (1) When the query is a file of marker sequences, they can be added as electronic markers. (2) When the query is draft sequence, the results of BSS can be used to close gaps in a sequenced clone or the physical map. (3) When the query is a sequenced clone and the target is BAC end sequences, one may select the next clone for sequencing using both sequence comparison results and map location. (4) When the query is whole-genome draft sequence and the target is BAC end sequences, the results can be used to select many clones for a minimal tiling path at once. The third tool, pickMTP, automates the majority of this last usage of BSS. Results are presented using the rice FPC map, BAC end sequences, and whole-genome shotgun from Syngenta.

Zaba: a novel miniature transposable element present in genomes of legume plants

A novel family of miniature transposable elements, named Zaba, was identified in pea (Pisum sativum) and subsequently also in other legume species using computer analysis of their DNA sequences. Zaba elements are 141-190 bp long, generate 10-bp target site duplications, and their terminal inverted repeats make up most of the sequence. Zaba elements thus resemble class 3 foldback transposons. The elements are only moderately repetitive in pea (tens to hundreds copies per haploid genome), but they are present in up to thousands of copies in the genomes of several Medicago and Vicia species. More detailed analysis of the elements from pea, including isolation of new sequences from a genomic library, revealed that a fraction of these elements are truncated, and that their last transposition probably did not occur recently. A search for Zaba sequences in EST databases showed that at least some elements are transcribed, most probably due to their association with genic regions.

Plant LTR-retrotransposons and MITEs: control of transposition and impact on the evolution of plant genes and genomes

Transposons are genetic elements that can move, and sometimes spread, within genomes, and that constitute an important fraction of eukaryote genomes. Two types of transposons, long terminal repeat (LTR)-retrotransposons and miniature inverted-repeat transposable elements (MITEs), are highly represented in plant genomes, and can account for as much as 50-80% of the total DNA content. In the last few years it has been shown that, in spite of their mutagenic capacity, both LTR-retrotransposons and MITEs can be found associated to genes, suggesting that their activity has influenced the evolution of plant genes. In this review we will summarise recent data on the control of the activity and the impact of both LTR-retrotransposons and MITEs on the evolution of plant genes and genomes.

Internet Contig Explorer (iCE)--a tool for visualizing clone fingerprint maps

Fingerprinted clone physical maps have proven useful in various applications, supporting both whole-genome and region-specific DNA sequencing as well as gene cloning studies. Fingerprint maps have been generated for several genomes, including those of human, mouse, rat, the nematodes Caenorhabditis elegans and Caenorhabditis briggsae, Arabidopsis thaliana and rice. Fingerprint maps of other genomes, including those of fungi, bacteria, poplar, and the cow, are being generated. The increasing use of fingerprint maps in genomic research has spawned a need in the research community for intuitive computer tools that facilitate viewing of the maps and the underlying fingerprint data. In this report we describe a new Java-based application called iCE (Internet Contig Explorer) that has been designed to provide views of fingerprint maps and associated data. Users can search for and display individual clones, contigs, clone fingerprints, clone insert sizes and markers. Users can also load into the software lists of particular clones of interest and view their fingerprints. iCE is being used at our Genome Centre to offer up to the research community views of the mouse, rat, bovine, C. briggsae, and several fungal genome bacterial artificial chromosome (BAC) fingerprint maps we have either completed or are currently constructing. We are also using iCE as part of the Rat Genome Sequencing Project to manage our provision of rat BAC clones for sequencing at the Human Genome Sequencing Center at the Baylor College of Medicine.