Non-random distribution of transposable elements in the nuclear genome of plants

Impact of transposable elements on the evolution of complex living systems and their epigenetic control

Transposable elements (TEs) contribute to genomic innovations, as well as genome instability, across a wide variety of species. Popular designations such as 'selfish DNA' and 'junk DNA,' common in the 1980s, may be either inaccurate or misleading, while a more enlightened view of the TE-host relationship covers a range from parasitism to mutualism. Both plant and animal hosts have evolved epigenetic mechanisms to reduce the impact of TEs, both by directly silencing them and by reducing their ability to transpose in the genome. However, TEs have also been co-opted by both plant and animal genomes to perform a variety of physiological functions, ranging from TE-derived proteins acting directly in normal biological functions to innovations in transcription factor activity and also influencing gene expression. Their presence, in fact, can affect a range of features at genome, phenotype, and population levels. The impact TEs have had on evolution is multifaceted, and many aspects still remain unexplored. In this review, the epigenetic control of TEs is contextualized according to the evolution of complex living systems.

Integrated cytogenetics and genomics analysis of transposable elements in the Nile tilapia, Oreochromis niloticus

Integration of cytogenetics and genomics has become essential to a better view of architecture and function of genomes. Although the advances on genomic sequencing have contributed to study genes and genomes, the repetitive DNA fraction of the genome is still enigmatic and poorly understood. Among repeated DNAs, transposable elements (TEs) are major components of eukaryotic chromatin and their investigation has been hindered even after the availability of whole sequenced genomes. The cytogenetic mapping of TEs in chromosomes has proved to be of high value to integrate information from the micro level of nucleotide sequence to a cytological view of chromosomes. Different TEs have been cytogenetically mapped in cichlids; however, neither details about their genomic arrangement nor appropriated copy number are well defined by these approaches. The current study integrates TEs distribution in Nile tilapia Oreochromis niloticus genome based on cytogenetic and genomics/bioinformatics approach. The results showed that some elements are not randomly distributed and that some are genomic dependent on each other. Moreover, we found extensive overlap between genomics and cytogenetics data and that tandem duplication may be the major mechanism responsible for the genomic dynamics of TEs here analyzed. This paper provides insights in the genomic organization of TEs under an integrated view based on cytogenetics and genomics.

Patterns of nucleotide asymmetries in plant and animal genomes

Symmetry in biology provides many intriguing puzzles to the scientist's mind. Chargaff's second parity rule states a symmetric distribution of oligonucleotides within a single strand of double-stranded DNA. While this rule has been verified in a wide range of microbial genomes, it still awaits explanation. In our study, we inquired into patterns of mono- and trinucleotide intra-strand parity in complex plant genomic sequences that became available during the last few years, and compared these to equally complex animal genomes. The degree and patterns of deviation from Chargaff's second rule were different between plant and animal species. We observed a universal inter-chromosomal homogeneity of mononucleotide skews in coding sequences of plant chromosomes, while the base composition of animal coding sequences differed between chromosomes even within a single species. We also found differences in the base composition of dicot introns in comparison to those of monocots. These genome-wide patterns were limited to genic regions and were not encountered in inter-genic sequences. We discuss the implications of our findings in relation to hypotheses about functional correlations of intra-strand parity which have hitherto been put forward. Furthermore, we propose more recent polyploidization and subsequent homogenization of homoeologues as a possible reason for more homogeneous skew patterns in plants.

Adventures in the enormous: a 1.8 million clone BAC library for the 21.7 Gb genome of loblolly pine

Loblolly pine (LP; Pinus taeda L.) is the most economically important tree in the U.S. and a cornerstone species in southeastern forests. However, genomics research on LP and other conifers has lagged behind studies on flowering plants due, in part, to the large size of conifer genomes. As a means to accelerate conifer genome research, we constructed a BAC library for the LP genotype 7-56. The LP BAC library consists of 1,824,768 individually-archived clones making it the largest single BAC library constructed to date, has a mean insert size of 96 kb, and affords 7.6X coverage of the 21.7 Gb LP genome. To demonstrate the efficacy of the library in gene isolation, we screened macroarrays with overgos designed from a pine EST anchored on LP chromosome 10. A positive BAC was sequenced and found to contain the expected full-length target gene, several gene-like regions, and both known and novel repeats. Macroarray analysis using the retrotransposon IFG-7 (the most abundant repeat in the sequenced BAC) as a probe indicates that IFG-7 is found in roughly 210,557 copies and constitutes about 5.8% or 1.26 Gb of LP nuclear DNA; this DNA quantity is eight times the Arabidopsis genome. In addition to its use in genome characterization and gene isolation as demonstrated herein, the BAC library should hasten whole genome sequencing of LP via next-generation sequencing strategies/technologies and facilitate improvement of trees through molecular breeding and genetic engineering. The library and associated products are distributed by the Clemson University Genomics Institute (www.genome.clemson.edu).

Toward the analysis of the petunia MADS box gene family by reverse and forward transposon insertion mutagenesis approaches: B, C, and D floral organ identity functions require SEPALLATA-like MADS box genes in petunia

We have initiated a systematic functional analysis of the MADS box, intervening region, K domain, C domain-type MADS box gene family in petunia. The starting point for this has been a reverse-genetics approach, aiming to select for transposon insertions into any MADS box gene. We have developed and applied a family signature insertion screening protocol that is highly suited for this purpose, resulting in the isolation of 32 insertion mutants in 20 different MADS box genes. In addition, we identified three more MADS box gene insertion mutants using a candidate-gene approach. The defined insertion lines provide a sound foundation for a systematic functional analysis of the MADS box gene family in petunia. Here, we focus on the analysis of Floral Binding Protein2 (FBP2) and FBP5 genes that encode the E-function, which in Arabidopsis has been shown to be required for B and C floral organ identity functions. fbp2 mutants display sepaloid petals and ectopic inflorescences originating from the third floral whorl, whereas fbp5 mutants appear as wild type. In fbp2 fbp5 double mutants, reversion of floral organs to leaf-like organs is increased further. Strikingly, ovules are replaced by leaf-like structures in the carpel, indicating that in addition to the B- and C-functions, the D-function, which specifies ovule development, requires E-function activity. Finally, we compare our data with results obtained using cosuppression approaches and conclude that the latter might be less suited for assigning functions to individual members of the MADS box gene family.

The maize LAG1-O mutant suggests that reproductive cell lineages show unique gene expression patterns early in vegetative development

Patterns of transposable element activity often provide useful information about how and when organisms regulate gene expression. The maize lowered Ac/Ds germinal reversion 1 (LAG1)-O mutation causes unusually low rates of germinal reversion by Ac/Ds-induced alleles even though these same alleles revert frequently and early in somatic development. LAG1-O suppresses Ds transposition at multiple, unlinked loci, and does not affect Spm elements, indicating that the mutation acts in trans and may be specific to Ac/Ds elements. Our data suggest that LAG1-O suppression gradually reduces Ac/Ds activity in the meristem and newly formed leaves until, by the floral transition, transposition is undetectable even with PCR-based assays. This suppression persists during tassel development and does not appear to be released until some point after meiosis. Competitive RT-PCR results show no difference in Ac transposase mRNA levels between LAG1-O and lag1(+) tassels, suggesting that suppression is post-transcriptional. The pattern of LAG1-O expression is consistent with a model in which at least some gene expression specific to those meristem cells that will ultimately give rise to floral tissue and therefore gametes begins very early in plant development, and then persists throughout development.

The distribution of T-DNA in the genomes of transgenic Arabidopsis and rice

Almost all the nuclear genes of four Gramineae (maize, wheat, barley, rice) and pea are located in DNA fractions covering only a 1-2% GC range and representing between 10 and 25% of the different genomes. These DNA fractions comprise large gene-rich regions (collectively called the 'gene space') separated by vast gene-empty, repeated sequences. In contrast, in Arabidopsis thaliana, genes are distributed in DNA fractions covering an 8% GC range and representing 85% of the genome. Here, we investigated the integration of a transferred DNA (T-DNA) in the genomes of Arabidopsis and rice and found different patterns of integration, which are correlated with the different gene distributions. While T-DNA integrates essentially everywhere in the Arabidopsis genome, integration was detected only in the gene space, namely in the gene-rich, transcriptionally active, regions of the rice genome. The implications of these results for the integration of foreign DNA are discussed.

Recent, extensive, and preferential insertion of members of the miniature inverted-repeat transposable element family Heartbreaker into genic regions of maize

A 314-bp DNA element called Heartbreaker-hm1 (Hbr-hm1) was previously identified in the 3' untranslated region of a mutant allele of the maize disease resistance gene HM1. This element has structural features of miniature inverted-repeat transposable elements (MITEs) and is a member of a large family of approximately 4,000 copies in the maize genome. Unlike previously described MITEs, most members of the Hbr family display over 90% sequence identity. This, coupled with the insertion of an Hbr element into an allele of the HM1 gene, suggested that this family might have spread recently throughout the genome. Consistent with this view is the finding that Hbr insertion sites are remarkably polymorphic. Ten of ten loci containing Hbr elements were found to be polymorphic for the presence or absence of Hbr among a collection of maize inbred lines and teosinte strains. Despite the fact that over 80% of the maize genome contain moderate to highly repetitive DNA, we find that randomly chosen Hbr elements are predominantly in single or low copy regions. Furthermore, when used to query both the public and private databases of plant genes, over 50% of the sequences flanking these Hbr elements resulted in significant "hits." Taken together, these data indicate that the presence or absence of Hbr elements is a significant contributory factor to the high level of polymorphism associated with maize genic regions.

Distribution of genes in the genome of Arabidopsis thaliana and its implications for the genome organization of plants

Previous work has shown that, in the large genomes of three Gramineae [rice, maize, and barley: 415, 2,500, and 5,300 megabases (Mb), respectively] most genes are clustered in long DNA segments (collectively called the "gene space") that represent a small fraction (12-24%) of nuclear DNA, cover a very narrow (0.8-1.6%) GC range, and are separated by vast expanses of gene-empty sequences. In the present work, we have analyzed the small (ca. 120 Mb) nuclear genome of Arabidopsis thaliana and shown that its organization is drastically different from that of the genomes of Gramineae. Indeed, (i) genes are distributed over about 85% of the main band of DNA in CsCl and cover an 8% GC range; (ii) ORFs are fairly evenly distributed in long (>50 kb) sequences from GenBank that amount to about 10 Mb; and (iii) the GC levels of protein-coding sequences (and of their third codon positions) are correlated with the GC levels of their flanking sequences. The different pattern of gene distribution of Arabidopsis compared with Gramineae appears to be because the genomes of the latter comprise (i) many large gene-empty regions separating gene clusters and (ii) abundant transposons in the intergenic sequences of gene clusters. Both sequences are absent or very scarce in the Arabidopsis genome. These observations provide a comparative view of angiosperm genome organization.

Position effects and epigenetic silencing of plant transgenes

Nuclear processes that silence plant transgenes are being revealed by analyses of natural triggers of epigenetic modifications, particularly cytosine methylation, and by comparisons of the genomic environments of differentially expressed transgene loci. It is increasingly apparent that plant genomes can sense and respond to the presence of foreign DNA in certain sequence contexts and at multiple dispersed sites. Determining the basis of this sensitivity and how nuclear defense systems are activated poses major challenges for the future.

Transposon Display identifies individual transposable elements in high copy number lines

The dTph1 transposable element family of Petunia hybrida line W138 consists of between 100 and 200 members. A strategy that allows simultaneous detection of individual elements is described. Sequences flanking dTph1 elements are amplified by means of a ligation-mediated PCR. The resulting fragments are locus-specific and can be analysed by polyacrylamide gel electrophoresis. One of the applications of Transposon Display is the isolation of dTph1-tagged genes. Fragments that co-segregate with a mutant phenotype can be extracted from the gel and reamplified, providing access to tagged genes, as demonstrated in a reconstruction experiment. Data on the molecular identification of a phenotypic mutant, isolated in a random tagging experiment is also presented. Upon sequencing, the obtained candidate fragment was found to be identical to part of the previously identified Fbp1 gene.

The distribution of genes in the genomes of Gramineae

Recent investigations showed that most maize genes are present in compositional fractions of nuclear DNA that cover only a 1-2% GC (molar fraction of guanosine plus cytosine in DNA) range and represent only 10-20% of the genome. These fractions, which correspond to compositional genome compartments that are distributed on all chromosomes, were collectively called the "gene space." Outside the gene space, the maize genome appears to contain no genes, except for some zein genes and for ribosomal genes. Here, we investigated the distribution of genes in the genomes of two other Gramineae, rice and barley, and used a new set of probes to study further the gene distribution of maize. We found that the distribution of genes in these three genomes is basically similar in that all genes, except for ribosomal genes and some storage protein genes, were located in gene spaces that (i) cover GC ranges of 0.8%, 1.0%, and 1.6% and represent 12%, 17%, and 24% of the genomes of barley, maize, and rice, respectively; (ii) are due to a remarkably uniform base composition in the sequences surrounding the genes, which are now known to consist mainly of transposons; (iii) have sizes approximately proportional to genome sizes, suggesting that expansion-contraction phenomena proceed in parallel in the gene space and in the gene-empty regions of the genome; and (iv) only hybridize on the gene spaces (and not on the other DNA fractions) of other Gramineae.

Comparative analysis of the chromosomal and genomic organization of Ty1-copia-like retrotransposons in pteridophytes, gymnosperms and angiosperms

We have investigated the physical distribution of the reverse transcriptase genes of Ty1-copia-like retrotransposable elements from 12 plant species belonging to different subdivisions by hybridization in situ on chromosome preparations. Ty1-copia-like elements showed different and non-random hybridization patterns. A dispersed distribution throughout most of the chromosomes with reduced hybridization at some regions or with some weak clustering at other regions was found in Allium cepa, Beta vulgaris, Brassica campestris, Brassica oleracea, Pennisetum glaucum, Pinus elliottii, Selaginella apoda, Vicia faba and Vicia narbonensis. Reduced hybridization occurred mainly at centromeric regions, nucleolus-organizing regions and regions known to be mainly composed of tandemly repeated sequences. In the fern Pteris cretica the retroelements showed a dispersed genomic organization with clustering at some chromosomal regions and whole chromosomes showing little signal. In Arabidopsis thaliana and Cicer arietinum Ty1-copia-like elements were found in clusters at the paracentromeric heterochromatin, a novel organization for a repetitive element in A. thaliana. New retroelement families were isolated from A. thaliana and from Beta vulgaris. Alignment of the deduced peptide sequences with Ty1-copia-like elements from other plants showed considerable divergence which was used to calculate their relationships, indicating the value of reverse transcriptase gene analysis in phylogenetic and biodiversity studies.

The gene distribution of the maize genome

Previous investigations from our laboratory showed that the genomes of plants, like those of vertebrates, are mosaics of isochores, i.e., of very long DNA segments that are compositionally homogeneous and that can be subdivided into a small number of families characterized by different GC levels (GC is the mole fraction of guanine+cytosine). Compositional DNA fractions corresponding to different isochore families were used to investigate, by hybridization with appropriate probes, the gene distribution in vertebrate genomes. Here we report such a study on the genome of a plant, maize. The gene distribution that we found is most striking, in that almost all genes are present in isochores covering an extremely narrow (1-2%) GC range and only representing 10-20% of the genome. This gene distribution, which seems to characterize other Gramineae as well, is remarkably different from the gene distribution previously found in vertebrate genomes.

Interplasmid transposition of Drosophila hobo elements in non-drosophilid insects

A modified hobo element from Drosophila melanogaster was introduced into embryos of the housefly, Musca domestica (family Muscidae) and the Queensland fruitfly, Bactrocera tryoni (family Tephritidae) to assess its ability to transpose. Hobo was capable of transposition in these species and transposition products had all of the hallmarks of hobo transposition products recovered from D. melanogaster, including the movement only of sequences precisely delimited by the inverted terminal repeats of hobo, the creation of an 8 bp duplication of the insertion site and an absolute requirement for hobo-encoded transposase. Transposition of hobo into the target gene resulted in a non-random distribution of insertion sites, with 10 of 38 independent insertions into the same nucleotide position. The results indicate that hobo can transpose in heterologous species, further demonstrating the similarity of hobo to Ac (Activator) of Zea mays and Tam3 of Antirrhinum majus. Hobo has excellent potential to act as a gene vector or gene tagging agent in nondrosophilid insects.

Invasion of the hobo transposable element studied by in situ hybridization on polytene chromosomes of Drosophila melanogaster

The invasion kinetics of hobo transposable element in the Drosophila melanogaster genome was studied by in situ hybridization on the polytene chromosomes. Six independent lines of Drosophila melanogaster flies that had been previously transformed by microinjection of the pHFL1 plasmid containing a complete hobo element were followed over 50 generations. We observed that hobo elements were scattered on each of the chromosome arms, with more insertion sites on the 3R arm. The total number of insertion sites remains quite small, between four and six, at generation 52. On the 2R arm, a short inversion appeared once at generation 52. Most of the integration sites reported here were already described for several transposons but some of them appear to be hotspots for hobo elements.