Diversity, origin, and distribution of retrotransposons (gypsy and copia) in conifers

The repetitive DNA sequence landscape and DNA methylation in chromosomes of an apomictic tropical forage grass, Cenchrus ciliaris

Cenchrus ciliaris is an apomictic, allotetraploid pasture grass widely distributed in the tropical and subtropical regions of Africa and Asia. In this study, we aimed to investigate the genomic organization and characterize some of the repetitive DNA sequences in this species. Due to the apomictic propagation, various aneuploid genotypes are found, and here, we analyzed a 2n = 4x + 3 = 39 accession. The physical mapping of Ty1-copia and Ty3-gypsy retroelements through fluorescence in situ hybridization with a global assessment of 5-methylcytosine DNA methylation through immunostaining revealed the genome-wide distribution pattern of retroelements and their association with DNA methylation. Approximately one-third of Ty1-copia sites overlapped or spanned centromeric DAPI-positive heterochromatin, while the centromeric regions and arms of some chromosomes were labeled with Ty3-gypsy. Most of the retroelement sites overlapped with 5-methylcytosine signals, except for some Ty3-gypsy on the arms of chromosomes, which did not overlap with anti-5-mC signals. Universal retrotransposon probes did not distinguish genomes of C. ciliaris showing signals in pericentromeric regions of all 39 chromosomes, unlike highly abundant repetitive DNA motifs found in survey genome sequences of C. ciliaris using graph-based clustering. The probes developed from RepeatExplorer clusters gave strong in situ hybridization signals, mostly in pericentromeric regions of about half of the chromosomes, and we suggested that they differentiate the two ancestral genomes in the allotetraploid C. ciliaris, likely having different repeat sequence variants amplified before the genomes came together in the tetraploid.

Genome-Wide Identification and Development of LTR Retrotransposon-Based Molecular Markers for the Melilotus Genus

Melilotus is an important genus of legumes with industrial and medicinal value, partly due to the production of coumarin. To explore the genetic diversity and population structure of Melilotus, 40 accessions were analyzed using long terminal repeat (LTR) retrotransposon-based markers. A total of 585,894,349 bp of LTR retrotransposon sequences, accounting for 55.28% of the Melilotus genome, were identified using bioinformatics tools. A total of 181,040 LTR retrotransposons were identified and classified as Gypsy, Copia, or another type. A total of 350 pairs of primers were designed for assessing polymorphisms in 15 Melilotus albus accessions. Overall, 47 polymorphic primer pairs were screened for their availability and transferability in 18 Melilotus species. All the primer pairs were transferable, and 292 alleles were detected at 47 LTR retrotransposon loci. The average polymorphism information content (PIC) value was 0.66, which indicated that these markers were highly informative. Based on unweighted pair group method with arithmetic mean (UPGMA) dendrogram cluster analysis, the 18 Melilotus species were classified into three clusters. This study provides important data for future breeding programs and for implementing genetic improvements in the Melilotus genus.

Comparative Study of Pine Reference Genomes Reveals Transposable Element Interconnected Gene Networks

Sequencing the giga-genomes of several pine species has enabled comparative genomic analyses of these outcrossing tree species. Previous studies have revealed the wide distribution and extraordinary diversity of transposable elements (TEs) that occupy the large intergenic spaces in conifer genomes. In this study, we analyzed the distribution of TEs in gene regions of the assembled genomes of Pinus taeda and Pinus lambertiana using high-performance computing resources. The quality of draft genomes and the genome annotation have significant consequences for the investigation of TEs and these aspects are discussed. Several TE families frequently inserted into genes or their flanks were identified in both species’ genomes. Potentially important sequence motifs were identified in TEs that could bind additional regulatory factors, promoting gene network formation with faster or enhanced transcription initiation. Node genes that contain many TEs were observed in multiple potential transposable element-associated networks. This study demonstrated the increased accumulation of TEs in the introns of stress-responsive genes of pines and suggests the possibility of rewiring them into responsive networks and sub-networks interconnected with node genes containing multiple TEs. Many such regulatory influences could lead to the adaptive environmental response clines that are characteristic of naturally spread pine populations.

Preferential Subgenome Elimination and Chromosomal Structural Changes Occurring in Newly Formed Tetraploid Wheat-Aegilops ventricosa Amphiploid (AABBDvDvNvNv)

Artificial allopolyploids derived from the genera Triticum and Aegilops have been used as genetic resources for wheat improvement and are a classic example of evolution via allopolyploidization. In this study, we investigated chromosomes and subgenome transmission behavior in the newly formed allopolyploid of wheat group via multicolor Fluorescence in situ hybridization (mc-FISH), using pSc119.2, pTa535, and (GAA)₇ as probe combinations, to enabled us to precisely identify individual chromosomes in 381 S₃ and S₄ generations plants derived from reciprocal crosses between Ae. ventricosa (D^vD^vN^vN^v) and T. turgidum (AABB). A higher rate of aneuploidy, constituting 66.04–86.41% individuals, was observed in these two early generations. Of the four constituent subgenomes, D^v showed the highest frequency of elimination, followed by N^v and B, while A was the most stable. In addition, structural chromosomal changes occurred ubiquitously in the selfed progenies of allopolyploids. Among the constituent subgenomes, B showed the highest number of aberrations. In terms of chromosomal dynamics, there was no significant association between the chromosomal behavior model and the cytoplasm, with the exception of chromosomal loss in the D^v subgenome. The chromosome loss frequency in the D^v subgenome was significantly higher in the T. turgidum × Ae. ventricosa cross than in the Ae. ventricosa × T. turgidum cross. This result indicates that, although the D subgenome showed great instability, allopolyploids containing D subgenome could probably be maintained after a certain hybridization in which the D subgenome donor was used as the maternal parent at its onset stage. Our findings provide valuable information pertaining to the behavior patterns of subgenomes during allopolyploidization. Moreover, the allopolyploids developed here could be used as potential resources for the genetic improvement of wheat.

Functional and morphological evolution in gymnosperms: A portrait of implicated gene families

Gymnosperms diverged from their sister plant clade of flowering plants 300 Mya. Morphological and functional divergence between the two major seed plant clades involved significant changes in their reproductive biology, water-conducting systems, secondary metabolism, stress defense mechanisms, and small RNA-mediated epigenetic silencing. The relatively recent sequencing of several gymnosperm genomes and the development of new genomic resources have enabled whole-genome comparisons within gymnosperms, and between angiosperms and gymnosperms. In this paper, we aim to understand how genes and gene families have contributed to the major functional and morphological differences in gymnosperms, and how this information can be used for applied breeding and biotechnology. In addition, we have analyzed the angiosperm versus gymnosperm evolution of the pleiotropic drug resistance (PDR) gene family with a wide range of functionalities in plants' interaction with their environment including defense mechanisms. Some of the genes reviewed here are newly studied members of gene families that hold potential for biotechnological applications related to commercial and pharmacological value. Some members of conifer gene families can also be exploited for their potential in phytoremediation applications.

ConTEdb: a comprehensive database of transposable elements in conifers

Conifers are the largest and most ubiquitous group of gymnosperms and have significant ecological significance and economic importance. However, the huge and complex genomes have hindered the sequencing and mining of conifer genomes. In this study, we identified 413 423 transposable elements (TEs) from Picea abies, Picea glauca and Pinus taeda using a combination of multiple approaches and classified them into 11 133 families. A comprehensive web-based database, ConTEdb, was constructed and served for researchers. ConTEdb enables users to browse, retrieve and download the TE sequences from the database. Several analysis tools are integrated into ConTEdb to help users mine the TE data easily and effectively. In summary, ConTEdb provides a platform to study TE biology and functional genomics in conifers.

Reverse transcriptase sequences from mulberry LTR retrotransposons: characterization analysis

Copia and Gypsy play important roles in structural, functional and evolutionary dynamics of plant genomes. In this study, a total of 106 and 101, Copia and Gypsy reverse transcriptase (rt) were amplified respectively in the Morus notabilis genome using degenerate primers. All sequences exhibited high levels of heterogeneity, were rich in AT and possessed higher sequence divergence of Copia rt in comparison to Gypsy rt. Two reasons are likely to account for this phenomenon: a) these elements often experience deletions or fragmentation by illegitimate or unequal homologous recombination in the transposition process; b) strong purifying selective pressure drives the evolution of these elements through “selective silencing” with random mutation and eventual deletion from the host genome. Interestingly, mulberry rt clustered with other rt from distantly related taxa according to the phylogenetic analysis. This phenomenon did not result from horizontal transposable element transfer. Results obtained from fluorescence in situ hybridization revealed that most of the hybridization signals were preferentially concentrated in pericentromeric and distal regions of chromosomes, and these elements may play important roles in the regions in which they are found. Results of this study support the continued pursuit of further functional studies of Copia and Gypsy in the mulberry genome.

Obtaining retrotransposon sequences, analysis of their genomic distribution and use of retrotransposon-derived genetic markers in lentil (Lens culinaris Medik.)

Retrotransposons with long terminal repeats (LTR-RTs) are widespread mobile elements in eukaryotic genomes. We obtained a total of 81 partial LTR-RT sequences from lentil corresponding to internal retrotransposon components and LTRs. Sequences were obtained by PCR from genomic DNA. Approximately 37% of the LTR-RT internal sequences presented premature stop codons, pointing out that these elements must be non-autonomous. LTR sequences were obtained using the iPBS technique which amplifies sequences between LTR-RTs. A total of 193 retrotransposon-derived genetic markers, mainly iPBS, were used to obtain a genetic linkage map from 94 F₇ inbred recombinant lines derived from the cross between the cultivar Lupa and the wild ancestor L. culinaris subsp. orientalis. The genetic map included 136 markers located in eight linkage groups. Clusters of tightly linked retrotransposon-derived markers were detected in linkage groups LG1, LG2, and LG6, hence denoting a non-random genomic distribution. Phylogenetic analyses identified the LTR-RT families in which internal and LTR sequences are included. Ty3-gypsy elements were more frequent than Ty1-copia, mainly due to the high Ogre element frequency in lentil, as also occurs in other species of the tribe Vicieae. LTR and internal sequences were used to analyze in silico their distribution among the contigs of the lentil draft genome. Up to 8.8% of the lentil contigs evidenced the presence of at least one LTR-RT similar sequence. A statistical analysis suggested a non-random distribution of these elements within of the lentil genome. In most cases (between 97% and 72%, depending on the LTR-RT type) none of the internal sequences flanked by the LTR sequence pair was detected, suggesting that defective and non-autonomous LTR-RTs are very frequent in lentil. Results support that LTR-RTs are abundant and widespread throughout of the lentil genome and that they are a suitable source of genetic markers useful to carry out further genetic analyses.

Identification and characterization of a LTR retrotransposon from the genome of Cyprinus carpio var. Jian

A Ty3/gypsy-retrotransposon-type transposon was found in the genome of the Jian carp (Cyprinus carpio var. Jian) in a previous study (unpublished), and was designated a JRE retrotransposon (Jian retrotransposon). The full-length JRE retrotransposon is 5126 bp, which includes two long terminal repeats of 470 bp at the 5' end and 453 bp at the 3' end, and two open reading frames between them: 4203 bp encoding the group-specific antigen (GAG) and polyprotein (POL). The pol gene has a typical Ty3/gypsy retrotransposon structure, and the gene order is protease, reverse transcriptase, RNase H, and integrase (PR-RT-RH-IN). A phylogenetic analysis of the pol gene showed that it has similarities of 40.7, 40, and 32.8 %, to retrotransposons of Azumapecten farreri, Mizuhopecten yessoensis, and Xiphophorus maculatus, respectively. Therefore, JRE might belong to the JULE retrotransposon family. The copy number of the JRE transposon in the genome of the Jian carp is 124, determined with real-time quantitative PCR. The mRNA of the JRE retrotransposon is expressed in five Jian carp tissues, the liver, kidney, blood, muscle, and gonad, and slightly higher in the kidney and liver than in the other tissues.

Molecular cytogenetic analysis and genomic organization of major DNA repeats in castor bean (Ricinus communis L.)

This article addresses the bioinformatic, molecular genetic, and cytogenetic study of castor bean (Ricinus communis, 2n = 20), which belongs to the monotypic Ricinus genus within the Euphorbiaceae family. Because castor bean chromosomes are small, karyotypic studies are difficult. However, the use of DNA repeats has yielded new prospects for karyotypic research and genome characterization. In the present study, major DNA repeat sequences were identified, characterized and localized on mitotic metaphase and meiotic pachytene chromosomes. Analyses of the nucleotide composition, curvature models, and FISH localization of the rcsat39 repeat suggest that this repeat plays a key role in building heterochromatic arrays in castor bean. Additionally, the rcsat390 sequences were determined to be chromosome-specific repeats located in the pericentromeric region of mitotic chromosome A (pachytene chromosome 1). The localization of rcsat39, rcsat390, 45S and 5S rDNA genes allowed for the development of cytogenetic landmarks for chromosome identification. General questions linked to heterochromatin formation, DNA repeat distribution, and the evolutionary emergence of the genome are discussed. The article may be of interest to biologists studying small genome organization and short monomer DNA repeats.

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

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

Characterization and chromosomal organization of Ty1-copia retrotransposons in wax gourd

Wax gourd (2n=2x=24) is an important vegetable species in Cucurbitaceae. Because it can be stored for a very long period of time, it plays an important role in ensuring the annual supply and regulating off-season supply of the vegetables. However, the availability of genetic information about wax gourd is limited. This study aimed to identify the useful genetic information for wax gourd. The conserved domains of reverse transcriptase (RT) genes of Ty1-copia retrotransposons were isolated from the genome of wax gourd using degenerate oligonucleotide primers. A total of twenty eight RT sequences were obtained, which showed high heterogeneity with the similarity ranging from 47.5% to 94.3%. Sixteen (57.1%) of them were found to be defective, being disrupted by stop codons and/or frameshift mutations. These 28 sequences were divided into five subfamilies. The comparative phylogenetic analysis with other Cucurbitaceae species from GenBank database showed that most retrotransposons derived from the same genus tended to cluster together, although there were a few exceptions. These results indicate that both vertical transmission and horizontal transmission are the sources of Ty1-copia retrotransposons in wax gourd. Fluorescent in situ hybridization (FISH) with Ty1-copia retrotransposon sequences as probes revealed that this kind of retrotransposons had a dispersed genomic organization, physically distributed among all the chromosomes of wax gourd, with clusters in the heterochromatin regions. This is the first report of Ty1-copia retrotransposons in wax gourd, which would be helpful for our understanding about the organization and evolutions of wax gourd genome and also provide valuable information for our utilization of wax gourd retrotransposons.

Unique features of the loblolly pine (Pinus taeda L.) megagenome revealed through sequence annotation

The largest genus in the conifer family Pinaceae is Pinus, with over 100 species. The size and complexity of their genomes (∼20-40 Gb, 2n = 24) have delayed the arrival of a well-annotated reference sequence. In this study, we present the annotation of the first whole-genome shotgun assembly of loblolly pine (Pinus taeda L.), which comprises 20.1 Gb of sequence. The MAKER-P annotation pipeline combined evidence-based alignments and ab initio predictions to generate 50,172 gene models, of which 15,653 are classified as high confidence. Clustering these gene models with 13 other plant species resulted in 20,646 gene families, of which 1554 are predicted to be unique to conifers. Among the conifer gene families, 159 are composed exclusively of loblolly pine members. The gene models for loblolly pine have the highest median and mean intron lengths of 24 fully sequenced plant genomes. Conifer genomes are full of repetitive DNA, with the most significant contributions from long-terminal-repeat retrotransposons. In depth analysis of the tandem and interspersed repetitive content yielded a combined estimate of 82%.

Isolation and characterization of reverse transcriptase fragments of LTR retrotransposons from the genome of Chenopodium quinoa (Amaranthaceae)

High heterogeneity was observed among conserved domains of reverse transcriptase ( rt ) isolated from quinoa. Only one Ty1- copia rt was highly amplified. Reverse transcriptase sequences were located predominantly in pericentromeric region of quinoa chromosomes. The heterogeneity, genomic abundance, and chromosomal distribution of reverse transcriptase (rt)-coding fragments of Ty1-copia and Ty3-gypsy long terminal repeat retrotransposons were analyzed in the Chenopodium quinoa genome. Conserved domains of the rt gene were amplified and characterized using degenerate oligonucleotide primer pairs. Sequence analyses indicated that half of Ty1-copia rt (51 %) and 39 % of Ty3-gypsy rt fragments contained intact reading frames. High heterogeneity among rt sequences was observed for both Ty1-copia and Ty3-gypsy rt amplicons, with Ty1-copia more heterogeneous than Ty3-gypsy. Most of the isolated rt fragments were present in quinoa genome in low copy numbers, with only one highly amplified Ty1-copia rt sequence family. The gypsy-like RNase H fragments co-amplified with Ty1-copia-degenerate primers were shown to be highly amplified in the quinoa genome indicating either higher abundance of some gypsy families of which rt domains could not be amplified, or independent evolution of this gypsy-region in quinoa. Both Ty1-copia and Ty3-gypsy retrotransposons were preferentially located in pericentromeric heterochromatin of quinoa chromosomes. Phylogenetic analyses of newly amplified rt fragments together with well-characterized retrotransposon families from other organisms allowed identification of major lineages of retroelements in the genome of quinoa and provided preliminary insight into their evolutionary dynamics.

Insights into the loblolly pine genome: characterization of BAC and fosmid sequences

Despite their prevalence and importance, the genome sequences of loblolly pine, Norway spruce, and white spruce, three ecologically and economically important conifer species, are just becoming available to the research community. Following the completion of these large assemblies, annotation efforts will be undertaken to characterize the reference sequences. Accurate annotation of these ancient genomes would be aided by a comprehensive repeat library; however, few studies have generated enough sequence to fully evaluate and catalog their non-genic content. In this paper, two sets of loblolly pine genomic sequence, 103 previously assembled BACs and 90,954 newly sequenced and assembled fosmid scaffolds, were analyzed. Together, this sequence represents 280 Mbp (roughly 1% of the loblolly pine genome) and one of the most comprehensive studies of repetitive elements and genes in a gymnosperm species. A combination of homology and de novo methodologies were applied to identify both conserved and novel repeats. Similarity analysis estimated a repetitive content of 27% that included both full and partial elements. When combined with the de novo investigation, the estimate increased to almost 86%. Over 60% of the repetitive sequence consists of full or partial LTR (long terminal repeat) retrotransposons. Through de novo approaches, 6,270 novel, full-length transposable element families and 9,415 sub-families were identified. Among those 6,270 families, 82% were annotated as single-copy. Several of the novel, high-copy families are described here, with the largest, PtPiedmont, comprising 133 full-length copies. In addition to repeats, analysis of the coding region reported 23 full-length eukaryotic orthologous proteins (KOGS) and another 29 novel or orthologous genes. These discoveries, along with other genomic resources, will be used to annotate conifer genomes and address long-standing questions about gymnosperm evolution.

Karyotype evolution in the Pinaceae: implication with molecular phylogeny

The family Pinaceae is made up mostly of diploid species (2n = 24). Systematization of karyotype analysis was developed to make comparison of intra- and interspecific karyotypes among the Pinaceae more accurate and reliable. Considering all parameters, the genera Pseudotsuga and Pseudolarix have the “most derived” (or advanced) and asymmetric karyotypes in the Pinaceae, followed by Larix, Picea, Abies, and Cedrus. The genus Pinus was the “least derived” (or ancestral) of all the genera of the Pinaceae analyzed. Differences in karyotype formulae and asymmetry indices were found among species within the same genera, suggesting that structural changes may have contributed to the diversification of the genus. This review is a detailed analysis of comparative karyotyping based on similar parameters, including numeric data and cytogenetic information. Telomeric sequence repeats and rDNA distribution in the Pinaceae were surveyed. The role of transposition in rDNA chromosome distribution is analyzed. Cytogenetic implications of hybridization between related species are reported. Likewise, the relationships between molecular phylogenetic and karyotype evolution is discussed in light of several reports. Within many genera, chromosomal organization was conserved despite independent molecular divergence and adaptation through the evolutionary history of the species of the Pinaceae.

Transcriptome profiling and in silico analysis of somatic embryos in Japanese larch (Larix leptolepis)

Japanese larch (Larix leptolepis) is an ecologically and economically important species mainly grown in northeastern China, Japan and Europe. However, erratic flowering and poor germplasm resources caused by high embryo abortion rates have hampered breeding of Larix species. Somatic embryogenesis (SE) is an effective tool for the production of L. leptolepis with desirable characteristics, such as expression of totipotency, preparation of synthetic seeds, and genetic transformation. However, public genomic resources for this species are limited. We sequenced 591,759 raw expressed sequence tags (ESTs) from a 454 sequencing cDNA library of L. leptolepis somatic embryos, resulting in 572,403 high-quality reads. These reads were assembled into 70,927 unique sequences (UniGenes), including 32,321 contigs and 38,606 singletons. After removal of low-quality sequences, 65,115 UniGenes were annotated using the UniProtKB program. Based on their sequence similarity with known proteins, the matched 30,372 sequences from 664 species were estimated to represent approximately 19,000 unique genes. Gene ontology analysis revealed 21,324 UniGenes assigned to 51 categories. By Kyoto Encyclopedia of Genes and Genomes mapping, 25,773 transcripts were associated with 160 biochemical pathways. Further analysis screened four signal transduction pathways represented by 337 enzymes and 17 secondary metabolites. In silico analysis reveals that 207 UniESTs in Larix are homologous to MAPKs genes identified from other model plants, which may be involved in regulating SE development. This study provides an initial insight into the Larix transcriptomes of the pro-embryogenic mass and is a sound basis for future studies.

KEY MESSAGE

We constructed a large, full-length 454 sequencing cDNA library of Larix leptolepis during somatic embryogenesis. More than 590,000 sequences were obtained and a deep-coverage EST database was constructed.

Deciduous trees and the application of universal DNA barcodes: a case study on the circumpolar Fraxinus

The utility of DNA barcoding for identifying representative specimens of the circumpolar tree genus Fraxinus (56 species) was investigated. We examined the genetic variability of several loci suggested in chloroplast DNA barcode protocols such as matK, rpoB, rpoC1 and trnH-psbA in a large worldwide sample of Fraxinus species. The chloroplast intergenic spacer rpl32-trnL was further assessed in search for a potentially variable and useful locus. The results of the study suggest that the proposed cpDNA loci, alone or in combination, cannot fully discriminate among species because of the generally low rates of substitution in the chloroplast genome of Fraxinus. The intergenic spacer trnH-psbA was the best performing locus, but genetic distance-based discrimination was moderately successful and only resulted in the separation of the samples at the subgenus level. Use of the BLAST approach was better than the neighbor-joining tree reconstruction method with pairwise Kimura's two-parameter rates of substitution, but allowed for the correct identification of only less than half of the species sampled. Such rates are substantially lower than the success rate required for a standardised barcoding approach. Consequently, the current cpDNA barcodes are inadequate to fully discriminate Fraxinus species. Given that a low rate of substitution is common among the plastid genomes of trees, the use of the plant cpDNA "universal" barcode may not be suitable for the safe identification of tree species below a generic or sectional level. Supplementary barcoding loci of the nuclear genome and alternative solutions are proposed and discussed.

Slow but not low: genomic comparisons reveal slower evolutionary rate and higher dN/dS in conifers compared to angiosperms

BACKGROUND

Comparative genomics can inform us about the processes of mutation and selection across diverse taxa. Among seed plants, gymnosperms have been lacking in genomic comparisons. Recent EST and full-length cDNA collections for two conifers, Sitka spruce (Picea sitchensis) and loblolly pine (Pinus taeda), together with full genome sequences for two angiosperms, Arabidopsis thaliana and poplar (Populus trichocarpa), offer an opportunity to infer the evolutionary processes underlying thousands of orthologous protein-coding genes in gymnosperms compared with an angiosperm orthologue set.

RESULTS

Based upon pairwise comparisons of 3,723 spruce and pine orthologues, we found an average synonymous genetic distance (dS) of 0.191, and an average dN/dS ratio of 0.314. Using a fossil-established divergence time of 140 million years between spruce and pine, we extrapolated a nucleotide substitution rate of 0.68 × 10(-9) synonymous substitutions per site per year. When compared to angiosperms, this indicates a dramatically slower rate of nucleotide substitution rates in conifers: on average 15-fold. Coincidentally, we found a three-fold higher dN/dS for the spruce-pine lineage compared to the poplar-Arabidopsis lineage. This joint occurrence of a slower evolutionary rate in conifers with higher dN/dS, and possibly positive selection, showcases the uniqueness of conifer genome evolution.

CONCLUSIONS

Our results are in line with documented reduced nucleotide diversity, conservative genome evolution and low rates of diversification in conifers on the one hand and numerous examples of local adaptation in conifers on the other hand. We propose that reduced levels of nucleotide mutation in large and long-lived conifer trees, coupled with large effective population size, were the main factors leading to slow substitution rates but retention of beneficial mutations.

Diversity and evolution of four dispersed repetitive DNA sequences in the genus Secale

We present the first characterization of 360 sequences in six species of the genus Secale of both cultivated and wild accessions. These include four distinct kinds of dispersed repetitive DNA sequences named pSc20H, pSc119.1, pSaO5(411), and pSaD15(940) belonging to the Revolver family. During the evolution of the genus Secale from wild to cultivated accessions, the pSaO5(411)-like sequences became shorter mainly because of the deletion of a trinucleotide tandem repeating unit, the pSc20H-like sequences displayed apparent homogenization in cultivated rye, and the second intron of Revolver became longer. In addition, the pSc20H-, pSc119.1-, and pSaO5(411)-like sequences cloned from wild rye and cultivated rye could be divided into two large clades. No single case of the four kinds of repetitive elements has been inherited by each Secale accession from a lone ancestor. It is reasonable to consider the vertical transmission of the four repetitive elements during the evolution of the genus Secale. The pSc20H- and pSaO5(411)-like sequences showed evolutionary elimination at specific chromosomal locations from wild species to cultivated species. These cases imply that different repetitive DNA sequences have played different roles in the chromosome development and genomic evolution of rye. The present study adds important information to the investigations dealing with characterization of dispersed repetitive elements in wild and cultivated rye.

Isolation and characterization of Ty1-copia group of LTRs in genome of three species of Datura: D. innoxia, D. stramonium and D. metel

Retrotransposons (RT) constitute a major fraction of plant genome. They are implicated in evolution and sequence organization. These elements have been proposed to have major role in evolution and variation in genome size. The sequence information of these RT regions in terms of divergence and conservation could be utilized for determining the interrelationship among various copia retrotransposons within the genome. In order to assess the diversity of Ty1-copia group of retroelements, reverse transcriptase (RT) sequence was amplified from genomes of three medicinally important Datura species: D. innoxia, D. stramonium and D. metel using the primers derived from two conserved domains of RT region. A total of twenty one independent amplicons from RT regions were cloned, sequenced and compared. The intra-family divergence at amino acid level ranged from 4 to 52 %. Though intra-family RT sequences are conserved, no two sequences are identical. Southern blot hybridization suggested that Ty1-copia-like retrotransposons are dispersed throughout the Datura genome. The results indicate a high degree of heterogeneity among the Ty1-copia group of retroelements in Datura species.

Plant centromeric retrotransposons: a structural and cytogenetic perspective

Background

The centromeric and pericentromeric regions of plant chromosomes are colonized by Ty3/gypsy retrotransposons, which, on the basis of their reverse transcriptase sequences, form the chromovirus CRM clade. Despite their potential importance for centromere evolution and function, they have remained poorly characterized. In this work, we aimed to carry out a comprehensive survey of CRM clade elements with an emphasis on their diversity, structure, chromosomal distribution and transcriptional activity.

Results

We have surveyed a set of 190 CRM elements belonging to 81 different retrotransposon families, derived from 33 host species and falling into 12 plant families. The sequences at the C-terminus of their integrases were unexpectedly heterogeneous, despite the understanding that they are responsible for targeting to the centromere. This variation allowed the division of the CRM clade into the three groups A, B and C, and the members of each differed considerably with respect to their chromosomal distribution. The differences in chromosomal distribution coincided with variation in the integrase C-terminus sequences possessing a putative targeting domain (PTD). A majority of the group A elements possess the CR motif and are concentrated in the centromeric region, while members of group C have the type II chromodomain and are dispersed throughout the genome. Although representatives of the group B lack a PTD of any type, they appeared to be localized preferentially in the centromeres of tested species. All tested elements were found to be transcriptionally active.

Conclusions

Comprehensive analysis of the CRM clade elements showed that genuinely centromeric retrotransposons represent only a fraction of the CRM clade (group A). These centromeric retrotransposons represent an active component of centromeres of a wide range of angiosperm species, implying that they play an important role in plant centromere evolution. In addition, their transcriptional activity is consistent with the notion that the transcription of centromeric retrotransposons has a role in normal centromere function.

Comparative analysis of pepper and tomato reveals euchromatin expansion of pepper genome caused by differential accumulation of Ty3/Gypsy-like elements

Background

Among the Solanaceae plants, the pepper genome is three times larger than that of tomato. Although the gene repertoire and gene order of both species are well conserved, the cause of the genome-size difference is not known. To determine the causes for the expansion of pepper euchromatic regions, we compared the pepper genome to that of tomato.

Results

For sequence-level analysis, we generated 35.6 Mb of pepper genomic sequences from euchromatin enriched 1,245 pepper BAC clones. The comparative analysis of orthologous gene-rich regions between both species revealed insertion of transposons exclusively in the pepper sequences, maintaining the gene order and content. The most common type of the transposon found was the LTR retrotransposon. Phylogenetic comparison of the LTR retrotransposons revealed that two groups of Ty3/Gypsy-like elements (Tat and Athila) were overly accumulated in the pepper genome. The FISH analysis of the pepper Tat elements showed a random distribution in heterochromatic and euchromatic regions, whereas the tomato Tat elements showed heterochromatin-preferential accumulation.

Conclusions

Compared to tomato pepper euchromatin doubled its size by differential accumulation of a specific group of Ty3/Gypsy-like elements. Our results could provide an insight on the mechanism of genome evolution in the Solanaceae family.

Diversity and evolution of Ty1-copia retroelements in representative tribes of Bambusoideae subfamily

Ty1-copia retroelements have been found in all major plants and are largely responsible for the huge differences in the genome size. In this study we isolated and sequenced Ty1-copia reverse transcriptase (rt) gene fragments from 44 representative species of bamboo and nine cultivars or forms of Phyllostachys pubescens. Phylogenetic analysis of 72 distinct Ty1-copia rt sequences showed that Ty1-copia retroelements were widespread, diverse and abundant in these species of Bambusoideae subfamily. In addition, a molecular phylogeny of the species of the Bambusoideae subfamily was established by using the internal transcribed spacer sequences of nuclear ribosomal DNA (ITS) sequences. The comparison between ITS- and Ty1-copia rt- based trees is obviously incongruent. The results suggested either the existence of horizontal transfer events between phylogenetically distant species, or an ancestral Ty1-copia retroelement polymorphism followed by different evolution and stochastic losses.

Evolutionary genomics revealed interkingdom distribution of Tcn1-like chromodomain-containing Gypsy LTR retrotransposons among fungi and plants

Background

Chromodomain-containing Gypsy LTR retrotransposons or chromoviruses are widely distributed among eukaryotes and have been found in plants, fungi and vertebrates. The previous comprehensive survey of chromoviruses from mosses (Bryophyta) suggested that genomes of non-seed plants contain the clade which is closely related to the retrotransposons from fungi. The origin, distribution and evolutionary history of this clade remained unclear mainly due to the absence of information concerning the diversity and distribution of LTR retrotransposons in other groups of non-seed plants as well as in fungal genomes.

Results

In present study we preformed in silico analysis of chromodomain-containing LTR retrotransposons in 25 diverse fungi and a number of plant species including spikemoss Selaginella moellendorffii (Lycopodiophyta) coupled with an experimental survey of chromodomain-containing Gypsy LTR retrotransposons from diverse non-seed vascular plants (lycophytes, ferns, and horsetails). Our mining of Gypsy LTR retrotransposons in genomic sequences allowed identification of numerous families which have not been described previously in fungi. Two new well-supported clades, Galahad and Mordred, as well as several other previously unknown lineages of chromodomain-containing Gypsy LTR retrotransposons were described based on the results of PCR-mediated survey of LTR retrotransposon fragments from ferns, horsetails and lycophytes. It appeared that one of the clades, namely Tcn1 clade, was present in basidiomycetes and non-seed plants including mosses (Bryophyta) and lycophytes (genus Selaginella).

Conclusions

The interkingdom distribution is not typical for chromodomain-containing LTR retrotransposons clades which are usually very specific for a particular taxonomic group. Tcn1-like LTR retrotransposons from fungi and non-seed plants demonstrated high similarity to each other which can be explained by strong selective constraints and the 'retained' genes theory or by horizontal transmission.

Expression of a gymnosperm PIN homologous gene correlates with auxin immunolocalization pattern at cotyledon formation and in demarcation of the procambium during Picea abies somatic embryo development and in seedling tissues

In seed plants, the body organization is established during embryogenesis and is uniform across gymnosperms and angiosperms, despite differences during early embryogeny. Evidence from angiosperms implicates the plant hormone auxin and its polar transport, mainly established by the PIN family of auxin efflux transporters, in the patterning of embryos. Here, PaPIN1 from Norway spruce (Picea abies [L.] Karst.), a gene widely expressed in conifer tissues and organs, was characterized and its expression and localization patterns were determined with reverse transcription polymerase chain reaction and in situ hybridization during somatic embryo development and in seedlings. PaPIN1 shares the predicted structure of other PIN proteins, but its central hydrophilic loop is longer than most PINs. In phylogenetic analyses, PaPIN1 clusters with Arabidopsis thaliana (L.) Heynh. PIN3, PIN4 and PIN7, but its expression pattern also suggests similarity to PIN1. The PaPIN1 expression signal was high in the protoderm of pre-cotyledonary embryos, but not if embryos were pre-treated with the auxin transport inhibitor N-1-naphthylphthalamic acid (NPA). This, together with a high auxin immunolocalization signal in this cell layer, suggests a role of PaPIN1 during cotyledon formation. At later stages, high PaPIN1 expression was observed in differentiating procambium, running from the tip of incipient cotyledons down through the embryo axis and to the root apical meristem (RAM), although the mode of RAM specification in conifer embryos differs from that of most angiosperms. Also, the PaPIN1 in situ signal was high in seedling root tips including root cap columella cells. The results thus suggest that PaPIN1 provides an ancient function associated with auxin transport and embryo pattern formation prior to the separation of angiosperms and gymnosperms, in spite of some morphological differences.

HACRE1, a recently inserted copia-like retrotransposon of sunflower (Helianthus annuus L.)

In this paper we report on the isolation and characterization, for the first time, of a complete 6511 bp retrotransposon of sunflower. Considering its protein domain order and sequence similarity to other copia elements of dicotyledons, this retrotransposon was assigned to the copia retrotransposon superfamily and named HACRE1 (Helianthus annuus copia-like retroelement 1). HACRE1 carries 5' and 3' long terminal repeats (LTRs) flanking an internal region of 4661 bp. The LTRs are identical in their sequence except for two deletions of 7 and 5 nucleotides in the 5' LTR. Based on the sequence identity of the LTRs, HACRE1 was estimated to have inserted within the last approximately 84 000 years. The isolated sequence contains a complete open reading frame with only one complete reading frame. The absence of nonsense mutations agrees with the very high sequence identity between LTRs, confirming that HACRE1 insertion is recent. The haploid genome of sunflower (inbred line HCM) contains about 160 copies of HACRE1. This retrotransposon is expressed in leaflets from 7-day-old plantlets under different light conditions, probably in relation to the occurrence of many putative light-related regulatory cis-elements in the LTRs. However, sequenced cDNAs show less variability than HACRE1 genomic sequences, indicating that only a subset of this family is expressed under these conditions.

Possible ancient origin of heterochromatic JNK sequences in chromosomes 2R of Secale vavilovii Grossh

Employing FISH analysis as well as BLAST and CUSTAL W (1.82) programs, we investigated types of DNA nucleotide sequences building an additional heterochromatic band in 2R chromosomes of 3 lines of Secale vavilovii Grossh. The probes used in FISH analysis were designed based on the reverse transcriptase sequence of Ty1-copia and Ty3-gypsy retrotransposons and the 5S rRNA gene sequence. No hybridization signals from the reverse transcriptase probes were observed in the chromosome region where the additional band occurs. On the other hand, signals were observed after hybridization with the 5S rDNA probe, clearly suggesting the presence of that type of sequences in the analyzed heterochromatin band. Using BLAST and CUSTAL W programs, we revealed high similarity of the JNK1 sequence to the 5S rRNA gene from Hordeum chilense (HCH1016, HCH1018, 88%) and to a fragment of the 5S rRNA sequence of H. marinum (HMAR003, 97%). In addition, the same fragment of JNK1 was shown to be very similar to the part of the Angela retrotransposon (92%) as well as to the SNAC 426K20-1 transposon (89%) belonging to CACTA family, both from Triticum monococcum, and to Zingeria biebersteiniana pericentromeric sequences (78%). The similarity of JNK1 to those sequences may be accidental or the JNK1 may represent an ancient mobile genetic element that caught the 5S rRNA sequence. During the evolution those sequences might have been accumulated in the particular region on the 2R chromosome. Our results suggest that the additional heterochromatin band in chromosomes 2R of S. vavilovii is a collection of defective genes and/or mobile genetic elements.

Targeted isolation, sequence assembly and characterization of two white spruce (Picea glauca) BAC clones for terpenoid synthase and cytochrome P450 genes involved in conifer defence reveal insights into a conifer genome

BACKGROUND

Conifers are a large group of gymnosperm trees which are separated from the angiosperms by more than 300 million years of independent evolution. Conifer genomes are extremely large and contain considerable amounts of repetitive DNA. Currently, conifer sequence resources exist predominantly as expressed sequence tags (ESTs) and full-length (FL)cDNAs. There is no genome sequence available for a conifer or any other gymnosperm. Conifer defence-related genes often group into large families with closely related members. The goals of this study are to assess the feasibility of targeted isolation and sequence assembly of conifer BAC clones containing specific genes from two large gene families, and to characterize large segments of genomic DNA sequence for the first time from a conifer.

RESULTS

We used a PCR-based approach to identify BAC clones for two target genes, a terpene synthase (3-carene synthase; 3CAR) and a cytochrome P450 (CYP720B4) from a non-arrayed genomic BAC library of white spruce (Picea glauca). Shotgun genomic fragments isolated from the BAC clones were sequenced to a depth of 15.6- and 16.0-fold coverage, respectively. Assembly and manual curation yielded sequence scaffolds of 172 kbp (3CAR) and 94 kbp (CYP720B4) long. Inspection of the genomic sequences revealed the intron-exon structures, the putative promoter regions and putative cis-regulatory elements of these genes. Sequences related to transposable elements (TEs), high complexity repeats and simple repeats were prevalent and comprised approximately 40% of the sequenced genomic DNA. An in silico simulation of the effect of sequencing depth on the quality of the sequence assembly provides direction for future efforts of conifer genome sequencing.

CONCLUSION

We report the first targeted cloning, sequencing, assembly, and annotation of large segments of genomic DNA from a conifer. We demonstrate that genomic BAC clones for individual members of multi-member gene families can be isolated in a gene-specific fashion. The results of the present work provide important new information about the structure and content of conifer genomic DNA that will guide future efforts to sequence and assemble conifer genomes.

Proliferation of Ty3/gypsy-like retrotransposons in hybrid sunflower taxa inferred from phylogenetic data

Background

Long terminal repeat (LTR) retrotransposons are a class of mobile genetic element capable of autonomous transposition via an RNA intermediate. Their large size and proliferative ability make them important contributors to genome size evolution, especially in plants, where they can reach exceptionally high copy numbers and contribute substantially to variation in genome size even among closely related taxa. Using a phylogenetic approach, we characterize dynamics of proliferation events of Ty3/gypsy-like LTR retrotransposons that led to massive genomic expansion in three Helianthus (sunflower) species of ancient hybrid origin. The three hybrid species are independently derived from the same two parental species, offering a unique opportunity to explore patterns of retrotransposon proliferation in light of reticulate evolutionary events in this species group.

Results

We demonstrate that Ty3/gypsy-like retrotransposons exist as multiple well supported sublineages in both the parental and hybrid derivative species and that the same element sublineage served as the source lineage of proliferation in each hybrid species' genome. This inference is based on patterns of species-specific element numerical abundance within different phylogenetic sublineages as well as through signals of proliferation events present in the distributions of element divergence values. Employing methods to date paralogous sequences within a genome, proliferation events in the hybrid species' genomes are estimated to have occurred approximately 0.5 to 1 million years ago.

Conclusion

Proliferation of the same retrotransposon major sublineage in each hybrid species indicates that similar dynamics of element derepression and amplification likely occurred in each hybrid taxon during their formation. Temporal estimates of these proliferation events suggest an earlier origin for these hybrid species than previously supposed.

Solanum lycopersicum cv. Heinz 1706 chromosome 6: distribution and abundance of genes and retrotransposable elements

We studied the physical and genetic organization of chromosome 6 of tomato (Solanum lycopersicum) cv. Heinz 1706 by combining bacterial artificial chromosome (BAC) sequence analysis, high-information-content fingerprinting, genetic analysis, and BAC-fluorescent in situ hybridization (FISH) mapping data. The chromosome positions of 81 anchored seed and extension BACs corresponded in most cases with the linear marker order on the high-density EXPEN 2000 linkage map. We assembled 25 BAC contigs and eight singleton BACs spanning 2.0 Mb of the short-arm euchromatin, 1.8 Mb of the pericentromeric heterochromatin and 6.9 Mb of the long-arm euchromatin. Sequence data were combined with their corresponding genetic and pachytene chromosome positions into an integrated map that covers approximately a third of the chromosome 6 euchromatin and a small part of the pericentromeric heterochromatin. We then compared physical length (Mb), genetic (cM) and chromosome distances (microm) for determining gap sizes between contigs, revealing relative hot and cold spots of recombination. Through sequence annotation we identified several clusters of functionally related genes and an uneven distribution of both gene and repeat sequences between heterochromatin and euchromatin domains. Although a greater number of the non-transposon genes were located in the euchromatin, the highly repetitive (22.4%) pericentromeric heterochromatin displayed an unexpectedly high gene content of one gene per 36.7 kb. Surprisingly, the short-arm euchromatin was relatively rich in repeats as well, with a repeat content of 13.4%, yet the ratio of Ty3/Gypsy and Ty1/Copia retrotransposable elements across the chromosome clearly distinguished euchromatin (2:3) from heterochromatin (3:2).

Evolution of genome size and complexity in Pinus

BACKGROUND

Genome evolution in the gymnosperm lineage of seed plants has given rise to many of the most complex and largest plant genomes, however the elements involved are poorly understood.

METHODOLOGY/PRINCIPAL FINDINGS

Gymny is a previously undescribed retrotransposon family in Pinus that is related to Athila elements in Arabidopsis. Gymny elements are dispersed throughout the modern Pinus genome and occupy a physical space at least the size of the Arabidopsis thaliana genome. In contrast to previously described retroelements in Pinus, the Gymny family was amplified or introduced after the divergence of pine and spruce (Picea). If retrotransposon expansions are responsible for genome size differences within the Pinaceae, as they are in angiosperms, then they have yet to be identified. In contrast, molecular divergence of Gymny retrotransposons together with other families of retrotransposons can account for the large genome complexity of pines along with protein-coding genic DNA, as revealed by massively parallel DNA sequence analysis of Cot fractionated genomic DNA.

CONCLUSIONS/SIGNIFICANCE

Most of the enormous genome complexity of pines can be explained by divergence of retrotransposons, however the elements responsible for genome size variation are yet to be identified. Genomic resources for Pinus including those reported here should assist in further defining whether and how the roles of retrotransposons differ in the evolution of angiosperm and gymnosperm genomes.

A transcriptionally active copia-like retroelement in Citrus limon

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

Novel clades of chromodomain-containing Gypsy LTR retrotransposons from mosses (Bryophyta)

Retrotransposons are the major component of plant genomes. Chromodomain-containing Gypsy long terminal repeat (LTR) retrotransposons are widely distributed in eukaryotes. Four distinct clades of chromodomain-containing Gypsy retroelements are known from the vascular plants: Reina, CRM, Galadriel and Tekay. At the same time, almost nothing is known about the repertoire of LTR retrotransposons in bryophyte genomes. We have combined a search of chromodomain-containing Gypsy retroelements in Physcomitrella genomic sequences and an experimental investigation of diverse moss species. The computer-based mining of the chromodomain-containing LTR retrotransposons allowed us to describe four different elements from Physcomitrella. Four novel clades were identified that are evolutionarily distinct from the chromodomain-containing Gypsy LTR retrotransposons of other plants.

Envelope-like retrotransposons in the plant kingdom: evidence of their presence in gymnosperms (Pinus pinaster)

Retroviruses differ from retrotransposons due to their infective capacity, which depends critically on the encoded envelope. Some plant retroelements contain domains reminiscent of the env of animal retroviruses but the number of such elements described to date is restricted to angiosperms. We show here the first evidence of the presence of putative env-like gene sequences in a gymnosperm species, Pinus pinaster (maritime pine). Using a degenerate primer approach for conserved domains of RNaseH gene, three clones from putative envelope-like retrotransposons (PpRT2, PpRT3, and PpRT4) were identified. The env-like sequences of P. pinaster clones are predicted to encode proteins with transmembrane domains. These sequences showed identity scores of up to 30% with env-like sequences belonging to different organisms. A phylogenetic analysis based on protein alignment of deduced aminoacid sequences revealed that these clones clustered with env-containing plant retrotransposons, as well as with retrotransposons from invertebrate organisms. The differences found among the sequences of maritime pine clones isolated here suggest the existence of different putative classes of env-like retroelements. The identification for the first time of env-like genes in a gymnosperm species may support the ancestrality of retroviruses among plants shedding light on their role in plant evolution.

A conifer genomics resource of 200,000 spruce (Picea spp.) ESTs and 6,464 high-quality, sequence-finished full-length cDNAs for Sitka spruce (Picea sitchensis)

Background

Members of the pine family (Pinaceae), especially species of spruce (Picea spp.) and pine (Pinus spp.), dominate many of the world's temperate and boreal forests. These conifer forests are of critical importance for global ecosystem stability and biodiversity. They also provide the majority of the world's wood and fiber supply and serve as a renewable resource for other industrial biomaterials. In contrast to angiosperms, functional and comparative genomics research on conifers, or other gymnosperms, is limited by the lack of a relevant reference genome sequence. Sequence-finished full-length (FL)cDNAs and large collections of expressed sequence tags (ESTs) are essential for gene discovery, functional genomics, and for future efforts of conifer genome annotation.

Results

As part of a conifer genomics program to characterize defense against insects and adaptation to local environments, and to discover genes for the production of biomaterials, we developed 20 standard, normalized or full-length enriched cDNA libraries from Sitka spruce (P. sitchensis), white spruce (P. glauca), and interior spruce (P. glauca-engelmannii complex). We sequenced and analyzed 206,875 3'- or 5'-end ESTs from these libraries, and developed a resource of 6,464 high-quality sequence-finished FLcDNAs from Sitka spruce. Clustering and assembly of 147,146 3'-end ESTs resulted in 19,941 contigs and 26,804 singletons, representing 46,745 putative unique transcripts (PUTs). The 6,464 FLcDNAs were all obtained from a single Sitka spruce genotype and represent 5,718 PUTs.

Conclusion

This paper provides detailed annotation and quality assessment of a large EST and FLcDNA resource for spruce. The 6,464 Sitka spruce FLcDNAs represent the third largest sequence-verified FLcDNA resource for any plant species, behind only rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana), and the only substantial FLcDNA resource for a gymnosperm. Our emphasis on capturing FLcDNAs and ESTs from cDNA libraries representing herbivore-, wound- or elicitor-treated induced spruce tissues, along with incorporating normalization to capture rare transcripts, resulted in a rich resource for functional genomics and proteomics studies. Sequence comparisons against five plant genomes and the non-redundant GenBank protein database revealed that a substantial number of spruce transcripts have no obvious similarity to known angiosperm gene sequences. Opportunities for future applications of the sequence and clone resources for comparative and functional genomics are discussed.

Isolation and characterization of genomic retrotransposon sequences from octoploid strawberry (Fragaria x ananassa Duch.)

Strawberry (Fragaria spp.) is a kind of herbaceous perennial plant that propagates vegetatively. The conserved domains of reverse transcriptase (RT) genes of Ty1-copia and Ty3-gypsy groups of LTR retrotransposons were amplified from the cultivated strawberry (Fragaria x ananassa Duch.). Sequence analysis of clones demonstrated that 5 of 19 Ty1-copia group unique sequences and 2 of 10 Ty3-gypsy unique sequences in F. x ananassa genome possessed either stop codon or frameshift. Ty1-copia group sequences are highly heterogeneous (divergence ranged from 1 to 69.8%), but the Ty3-gypsy group sequences are less (divergence ranged from 1 to 10%). Southern dot blot hybridization result suggested that both of the LTR retrotransposons are present in the genome of cultivated strawberry with high copy number (Ty1-copia group 2,875 Ty3-gypsy group 348). RT-PCR amplification from total RNA, which was extracted from leaves of micropropagated strawberry plants, did not yield either of the RT fragments. This is the first report on the presence of RT sequences of Ty1-copia and Ty3-gypsy group retrotransposons in F. x ananassa genome.

Reme1, a Copia retrotransposon in melon, is transcriptionally induced by UV light

For the first time, numerous sequences of Copia and Gypsy retrotransposons from the Cucumis melo genome have been obtained and analyzed. Phylogenetic analyses of sequences of both types of long terminal repeat (LTR) retrotransposons were carried out. The melon genome contains approximately 20,000 Gypsy and 6,800 Copia elements, comprising about 26% of its total size. Starting from a retrotransposon fragment, we have cloned and characterized an entire melon retrotransposon, named Reme1, which is 5,149 bp long. Reme1 belongs to the Superfamily Copia retrotransposons by its protein domain order and sequence similarity to other Copia elements of dicotyledons. The haploid genome of melon (var. "Piel de Sapo") contains about 120 copies of Reme1. Several copies of Reme1 are transcriptionally active, although at low levels, in melon leaves as analyzed by reverse-transcription PCR (RT-PCR) and sequencing. However, the transcript pool is considerably increased when melon leaves are treated with UV light, as has been seen for various retroelements in many organisms. The cDNAs of Reme1 transcripts showed less diversity than do Reme1 genomic sequences, suggesting that a subfamily of these elements is differentially responsive to UV.

Diversity and evolution of Ty1-copia and Ty3-gypsy retroelements in the non-photosynthetic flowering plants Orobanche and Phelipanche (Orobanchaceae)

We present the first study on the diversity and evolution of Ty1-copia and Ty3-gypsy retroelements in a group of non-photosynthetic flowering plants. To this end partial sequences of the reverse transcriptase (rt) gene were obtained from 20 clones for each retroelement type from seven and six accessions of Orobanche and Phelipanche (Orobanchaceae), respectively. Overall sequence similarity is higher in Ty3-gypsy elements than in Ty1-copia elements in agreement with the results from other angiosperm groups. Higher sequence diversity and stronger phylogenetic structure, especially of Ty1-copia sequences, in Orobanche species compared to Phelipanche species support the previously suggested hypothesis (based on karyological and cytological data) that genomes of Orobanche species are more dynamic than those of Phelipanche species. No evidence was found for intraspecific differences of retroelement diversity nor for differences between pest taxa and their putative wild relatives, e.g., O. crenata and O. owerini. The occurrence of a few sequences from Phelipanche species in clades otherwise comprising sequences from Orobanche species might be due to horizontal gene transfer, but the alternative of vertical transmission cannot be rejected unambiguously.

Identification and Characterization of Random Amplified Polymorphic DNA Markers Linked to a Major Gene (Cr2) for Resistance to Cronartium ribicola in Pinus monticola

ABSTRACT DNA markers tightly linked to resistance (R) genes provide a very powerful tool for both marker-assisted selection in plant breeding and positional cloning of R genes. In the present study, a linkage of random amplified polymorphic DNA (RAPD) markers to the single dominant gene (Cr2) for resistance to white pine blister rust fungus (Cronartium ribicola) was investigated in western white pine (Pinus monticola). A mapping population of 128 individual megagametophytes was generated from seeds of a heterozygous resistant tree (Cr2/cr2), and the corresponding seedlings of each megagametophyte were subjected to the test of phenotype segregation by inoculation with C. ribicola. Bulked segregant analysis and haploid segregation analysis identified eight robust RAPD markers linked to Cr2. This constitutes the first Cr2 genetic linkage map spanning 84.7 cM with four markers only 3.2 cM from Cr2. One sequence (U256-1385) of these linked markers was significantly similar to the Ty3/gypsy-like long terminal direct repeats retrotransposons. Another marker, U570-843, had no significant similarity to any entry in either GenBank or the loblolly genomics data bank. As presumed that the average physical distance per centimorgan is about 10 Mb in P. monticola, it is probably unrealistic to use these DNA markers for positional cloning of the Cr2 gene.

Retrotransposon Sequence Variation in Four Asexual Plant Species

Transposable elements (TEs) can be viewed as genetic parasites that persist in populations due to their capacity for increase in copy number and the inefficacy of selection against them. A corollary of this hypothesis is that TEs are more likely to spread within sexual populations and be eliminated or inactivated within asexual populations. While previous work with animals has shown that asexual taxa may contain less TE diversity than sexual taxa, comparable work with plants has been lacking. Here we report the results of a study of Ty1/copia, Ty3/gypsy, and LINE-like retroelement diversity in four asexual plant species. Retroelement-like sequences, with a high degree of conservation both within and between species, were isolated from all four species. The sequences correspond to several previously annotated retroelement subfamilies. They also exhibit a pattern of nucleotide substitution characterized by an excess of synonymous substitutions, suggestive of a history of purifying selection. These findings were compared with retroelement sequence evolution in sexual plant taxa. One likely explanation for the discovery of conserved TE sequences in the genomes of these asexual taxa is simply that asexuality within these taxa evolved relatively recently, such that the loss and breakdown of TEs is not yet detectable through analysis of sequence diversity. This explanation is examined by conducting stochastic simulation of TE evolution and by using published information to infer rough estimates of the ages of asexual taxa.