Development of an efficient retrotransposon-based fingerprinting method for rapid pea variety identification

Analysis of genetic diversity in pigeon pea germplasm using retrotransposon-based molecular markers

Pigeon pea (Cajanus cajan), an important legume crop is predominantly cultivated in tropical and subtropical regions of Asia and Africa. It is normally considered to have a low degree of genetic diversity, an impediment in undertaking crop improvement programmes.We have analysed genetic polymorphism of domesticated pigeon pea germplasm (47 accessions) across the world using earlier characterized panzee retrotransposon-based molecularmarkers. Itwas conjectured that since retrotransposons are interspersed throughout the genome, retroelements-based markers would be able to uncover polymorphism possibly inherent in the diversity of retroelement sequences. Two PCR-based techniques, sequence-specific amplified polymorphism (SSAP) and retrotransposon microsatellite amplified polymorphism (REMAP) were utilized for the analyses.We show that a considerable degree of polymorphism could be detected using these techniques. Three primer combinations in SSAP generated 297 amplified products across 47 accessions with an average of 99 amplicons per assay. Degree of polymorphism varied from 84-95%. In the REMAP assays, the number of amplicons was much less but up to 73% polymorphism could be detected. On the basis of similarity coefficients, dendrograms were constructed. The results demonstrate that the retrotransposon-based markers could serve as a better alternative for the assessment of genetic diversity in crops with apparent low genetic base.

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.

Useful parasites: the evolutionary biology and biotechnology applications of transposable elements

Transposable elements usually comprise the most abundant nongenic fraction of eukaryotic genomes. Because of their capacity to selfreplicate and to induce a wide range of mutations, transposable elements have long been considered as 'parasitic' or 'selfish'. Today, we recognize that the findings about genomic changes affected by transposable elements have considerably altered our view of the ways in which genomes evolve and work. Numerous studies have provided evidences that mobile elements have the potential to act as agents of evolution by increasing, rearranging and diversifying the genetic repertoire of their hosts. With large-scale sequencing becoming increasingly available, more and more scientists come across transposable element sequences in their data. I will provide examples that transposable elements, although having signatures of 'selfish' DNA, play a significant biological role in the maintainance of genome integrity and providing novel regulatoty networks. These features, along with the transpositional and mutagenic capacity to produce a raw genetic diversity, make the genome mobile fraction, a key player in species adaptation and microevolution. The last but not least, transposable elements stand as informative DNA markers that may complement other conventional DNA markers. Altogether, transposable elements represent a promising, but still largely unexplored research niche and deserve to be included into the agenda of molecular ecologists, evolutionary geneticists, conservation biologists and plant breeders.

Ty1-copia elements reveal diverse insertion sites linked to polymorphisms among flax (Linum usitatissimum L.) accessions

BACKGROUND

Initial characterization of the flax genome showed that Ty1-copia retrotransposons are abundant, with several members being recently inserted, and in close association with genes. Recent insertions indicate a potential for ongoing transpositional activity that can create genomic diversity among accessions, cultivars or varieties. The polymorphisms generated constitute a good source of molecular markers that may be associated with phenotype if the insertions alter gene activity. Flax, where accessions are bred mainly for seed nutritional properties or for fibers, constitutes a good model for studying the relationship of transpositional activity with diversification and breeding. In this study, we estimated copy number and used a type of transposon display known as Sequence-Specific Amplification Polymorphisms (SSAPs), to characterize six families of Ty1-copia elements across 14 flax accessions. Polymorphic insertion sites were sequenced to find insertions that could potentially alter gene expression, and a preliminary test was performed with selected genes bearing transposable element (TE) insertions.

RESULTS

Quantification of six families of Ty1-copia elements indicated different abundances among TE families and between flax accessions, which suggested diverse transpositional histories. SSAPs showed a high level of polymorphism in most of the evaluated retrotransposon families, with a trend towards higher levels of polymorphism in low-copy number families. Ty1-copia insertion polymorphisms among cultivars allowed a general distinction between oil and fiber types, and between spring and winter types, demonstrating their utility in diversity studies. Characterization of polymorphic insertions revealed an overwhelming association with genes, with insertions disrupting exons, introns or within 1 kb of coding regions. A preliminary test on the potential transcriptional disruption by TEs of four selected genes evaluated in three different tissues, showed one case of significant impact of the insertion on gene expression.

CONCLUSIONS

We demonstrated that specific Ty1-copia families have been active since breeding commenced in flax. The retrotransposon-derived polymorphism can be used to separate flax types, and the close association of many insertions with genes defines a good source of potential mutations that could be associated with phenotypic changes, resulting in diversification processes.

Transcriptionally active LTR retrotransposons in Eucalyptus genus are differentially expressed and insertionally polymorphic

BACKGROUND

In Eucalyptus genus, studies on genome composition and transposable elements (TEs) are particularly scarce. Nearly half of the recently released Eucalyptus grandis genome is composed by retrotransposons and this data provides an important opportunity to understand TE dynamics in Eucalyptus genome and transcriptome.

RESULTS

We characterized nine families of transcriptionally active LTR retrotransposons from Copia and Gypsy superfamilies in Eucalyptus grandis genome and we depicted genomic distribution and copy number in two Eucalyptus species. We also evaluated genomic polymorphism and transcriptional profile in three organs of five Eucalyptus species. We observed contrasting genomic and transcriptional behavior in the same family among different species. RLC_egMax_1 was the most prevalent family and RLC_egAngela_1 was the family with the lowest copy number. Most families of both superfamilies have their insertions occurring <3 million years, except one Copia family, RLC_egBianca_1. Protein theoretical models suggest different properties between Copia and Gypsy domains. IRAP and REMAP markers suggested genomic polymorphisms among Eucalyptus species. Using EST analysis and qRT-PCRs, we observed transcriptional activity in several tissues and in all evaluated species. In some families, osmotic stress increases transcript values.

CONCLUSION

Our strategy was successful in isolating transcriptionally active retrotransposons in Eucalyptus, and each family has a particular genomic and transcriptional pattern. Overall, our results show that retrotransposon activity have differentially affected genome and transcriptome among Eucalyptus species.

Genetic diversity of maize germplasm assessed by retrotransposon-based markers

Maize is one of the most important crops and also a model for grass genome research. Transposable elements comprise over 78% of the maize genome and their ability to generate new copies makes them good potential markers. Interretrotransposon-amplified polymorphism (IRAP) and retrotransposon microsatellite amplified polymorphism (REMAP) protocols were used for the first time in maize to study the genetic variability between maize cultivars. Ten PCR primers were selected based on a systematic analysis of the sequence conservation in the extremities of different high copy number transposable elements, whereas one primer was chosen based on a microsatellite sequence. Of the 16 primer combinations tested, 14 produced polymorphic bands. These markers were used to identify genetic similarity among 20 maize cultivars selected by their different kernel oil content. Genetic similarity analysis was performed based on the polymorphic band profiles and dendrograms were developed by the unweighted pair-group method with arithmetic averages. Clustering technique revealed that samples were grouped into three clusters that differed in their kernel oil content and size, and in their relative embryo size. In the current investigation, there is evidence that IRAP/REMAP may be useful as markers in maize.

Retrotransposon insertional polymorphism in Iranian bread wheat cultivars and breeding lines revealed by IRAP and REMAP markers

Inter-retrotransposon amplified polymorphisms (IRAPs) and retrotransposon-microsatellite amplified polymorphisms (REMAPs) were used to detect retrotransposon integration events and genetic diversity in 101 Iranian bread wheat (Triticum aestivum L.) cultivars and breeding lines. The 9 IRAP primers amplified 128 loci, and 20 REMAP primers amplified 263 loci. Percentage of polymorphic loci, average expected heterozygosity, number of effective alleles, and Shannon's information index for the REMAP markers were slightly higher than those for the IRAP markers. The same estimated parameters calculated for native and nonnative retrotransposons were not considerably different. A Mantel test between IRAP and REMAP cophenetic matrices evidenced no significant correlation. Cluster analysis based on the Dice similarity coefficient and complete linkage algorithm using IRAP+REMAP loci identified five groups among the genotypes studied that could be applied as crossing parents in T. aestivum breeding programs.

The application of LTR retrotransposons as molecular markers in plants

Retrotransposons are a major agent of genome evolution. Various molecular marker systems have been developed that exploit the ubiquitous nature of these genetic elements and their property of stable integration into dispersed chromosomal loci that are polymorphic within species. The key methods, SSAP, IRAP, REMAP, RBIP, and ISBP, all detect the sites at which the retrotransposon DNA, which is conserved between families of elements, is integrated into the genome. Marker systems exploiting these methods can be easily developed and inexpensively deployed in the absence of extensive genome sequence data. They offer access to the dynamic and polymorphic, nongenic portion of the genome and thereby complement methods, such as gene-derived SNPs, that target primarily the genic fraction.

Genetic diversity of cultivated flax (Linum usitatissimum L.) germplasm assessed by retrotransposon-based markers

Retrotransposon segments were characterized and inter-retrotransposon amplified polymorphism (IRAP) markers developed for cultivated flax (Linum usitatissimum L.) and the Linum genus. Over 75 distinct long terminal repeat retrotransposon segments were cloned, the first set for Linum, and specific primers designed for them. IRAP was then used to evaluate genetic diversity among 708 accessions of cultivated flax comprising 143 landraces, 387 varieties, and 178 breeding lines. These included both traditional and modern, oil (86), fiber (351), and combined-use (271) accessions, originating from 36 countries, and 10 wild Linum species. The set of 10 most polymorphic primers yielded 141 reproducible informative data points per accession, with 52% polymorphism and a 0.34 Shannon diversity index. The maximal genetic diversity was detected among wild Linum species (100% IRAP polymorphism and 0.57 Jaccard similarity), while diversity within cultivated germplasm decreased from landraces (58%, 0.63) to breeding lines (48%, 0.85) and cultivars (50%, 0.81). Application of Bayesian methods for clustering resulted in the robust identification of 20 clusters of accessions, which were unstratified according to origin or user type. This indicates an overlap in genetic diversity despite disruptive selection for fiber versus oil types. Nevertheless, eight clusters contained high proportions (70-100%) of commercial cultivars, whereas two clusters were rich (60%) in landraces. These findings provide a basis for better flax germplasm management, core collection establishment, and exploration of diversity in breeding, as well as for exploration of the role of retrotransposons in flax genome dynamics.

Analysis of plant diversity with retrotransposon-based molecular markers

Retrotransposons are both major generators of genetic diversity and tools for detecting the genomic changes associated with their activity because they create large and stable insertions in the genome. After the demonstration that retrotransposons are ubiquitous, active and abundant in plant genomes, various marker systems were developed to exploit polymorphisms in retrotransposon insertion patterns. These have found applications ranging from the mapping of genes responsible for particular traits and the management of backcrossing programs to analysis of population structure and diversity of wild species. This review provides an insight into the spectrum of retrotransposon-based marker systems developed for plant species and evaluates the contributions of retrotransposon markers to the analysis of population diversity in plants.

Genome merger: from sequence rearrangements in triticale to their elimination in wheat-rye addition lines

Genetic and epigenetic modifications resulting from different genomes adjusting to a common nuclear environment have been observed in polyploids. Sequence restructuring within genomes involving retrotransposon/microsatellite-rich regions has been reported in triticale. The present study uses inter-retrotransposon amplified polymorphisms (IRAP) and retrotransposon microsatellite amplified polymorphisms (REMAP) to assess genome rearrangements in wheat-rye addition lines obtained by the controlled backcrossing of octoploid triticale to hexaploid wheat followed by self-fertilization. The comparative analysis of IRAP and REMAP banding profiles, involving a complete set of wheat-rye addition lines, and their parental species revealed in those lines the presence of wheat-origin bands absent in triticale, and the absence of rye-origin and triticale-specific bands. The presence in triticale x wheat backcrosses (BC) of rye-origin bands that were absent in the addition lines demonstrated that genomic rearrangement events were not a direct consequence of backcrossing, but resulted from further genome structural rearrangements in the BC plant progeny. PCR experiments using primers designed from different rye-origin sequences showed that the absence of a rye-origin band in wheat-rye addition lines results from sequence elimination rather than restrict changes on primer annealing sites, as noted in triticale. The level of genome restructuring events evaluated in all seven wheat-rye addition lines, compared to triticale, indicated that the unbalanced genome merger situation observed in the addition lines induced a new round of genome rearrangement, suggesting that the lesser the amount of rye chromatin introgressed into wheat the larger the outcome of genome reshuffling.

Evolutionary conserved lineage of Angela-family retrotransposons as a genome-wide microsatellite repeat dispersal agent

A detailed examination of 45 pea (Pisum sativum L.) simple sequence repeat (SSR) loci revealed that 21 of them included homologous sequences corresponding to the long terminal repeat (LTR) of a novel retrotransposon. Further investigation, including full-length sequencing, led to its classification as an RLC-Angela-family-FJ434420 element. The LTR contained a variable region ranging from a simple TC repeat (TC)(11) to more complex repeats of TC/CA, (TC)(12-30), (CA)(18-22) and was up to 146 bp in length. These elements are the most abundant Ty1/copia retrotransposons identified in the pea genome and also occur in other legume species. It is interesting that analysis of 63 LTR-derived sequences originating from 30 legume species showed high phylogenetic conservation in their sequence, including the position of the variable SSR region. This extraordinary conservancy led us to the proposition of a new lineage, named MARTIANS, within the Angela family. Similar LTR structures and partial sequence similarities were detected in more distant members of this Angela family, the barley BARE-1 and rice RIRE-1 elements. Comparison of the LTR sequences from pea and Medicago truncatula elements indicated that microsatellites arise through the expansion of a pre-existing repeat motif. Thus, the presence of an SSR region within the LTR seems to be a typical feature of this MARTIANS lineage, and the evidence gathered from a wide range of species suggests that these elements may facilitate amplification and genome-wide dispersal of associated SSR sequences. The implications of this finding regarding the evolution of SSRs within the genome, as well as their utilization as molecular markers, are discussed.

Epigenetic changes detected in micropropagated hop plants

Micropropagation is a widely used technique in hops (Humulus lupulus L.). However, to the best of our knowledge, the genetic and epigenetic stability of the microplants has never been tested before. In the present study, two hop accessions were established in vitro and micropropagated for 2 years. The genetic and epigenetic stability of the in vitro plants was analyzed with several molecular techniques: random amplified DNA polymorphism (RAPD), retrotransposon microsatellite amplified polymorphism (REMAP), and methylation-sensitive amplification polymorphism (MSAP). No genetic variation among control and treated plants was found, even after 12 cycles of micropropagation. Epigenetic variation was detected, first, when field and in vitro samples were compared. Nearly a 30% of the detected fragments presented the same pattern of alterations in all the vitroplants. Second, lower levels of epigenetic variation were detected among plants from the different subcultures. Part of this detected variation seemed to be accumulated along the 12 sequential subcultures tested.

Genetic diversity and population structure of pea (Pisum sativum L.) varieties derived from combined retrotransposon, microsatellite and morphological marker analysis

One hundred and sixty-four accessions representing Czech and Slovak pea (Pisum sativum L.) varieties bred over the last 50 years were evaluated for genetic diversity using morphological, simple sequence repeat (SSR) and retrotransposon-based insertion polymorphism (RBIP) markers. Polymorphic information content (PIC) values of 10 SSR loci and 31 RBIP markers were on average high at 0.89 and 0.73, respectively. The silhouette method after the Ward clustering produced the most probable cluster estimate, identifying nine clusters from molecular data and five to seven clusters from morphological characters. Principal component analysis of nine qualitative and eight quantitative morphological parameters explain over 90 and 93% of total variability, respectively, in the first three axes. Multidimensional scaling of molecular data revealed a continuous structure for the set. To enable integration and evaluation of all data types, a Bayesian method for clustering was applied. Three clusters identified using morphology data, with clear separation of fodder, dry seed and afila types, were resolved by DNA data into 17, 12 and five sub-clusters, respectively. A core collection of 34 samples was derived from the complete collection by BAPS Bayesian analysis. Values for average gene diversity and allelic richness for molecular marker loci and diversity indexes of phenotypic data were found to be similar between the two collections, showing that this is a useful approach for representative core selection.

Polyploidization as a retraction force in plant genome evolution: sequence rearrangements in triticale

BACKGROUND

Polyploidization is a major evolutionary process in plants where hybridization and chromosome doubling induce enormous genomic stress and can generate genetic and epigenetic modifications. However, proper evaluation of DNA sequence restructuring events and the precise characterization of sequences involved are still sparse.

METHODOLOGY/PRINCIPAL FINDINGS

Inter Retrotransposons Amplified Polymorphism (IRAP), Retrotransposons Microsatellite Amplified Polymorphism (REMAP) and Inter Simple Sequence Repeat (ISSR) largely confirmed the absence of any intraspecific variation in wheat, rye and triticale. The comparative analysis of banding profiles between wheat and rye inbred lines revealed 34% of monomorphic (common to both parental species) bands for the ten different primer combinations used. The analysis of triticale plants uncovered nearly 51% of rearranged bands in the polyploid, being the majority of these modifications, due to the loss of rye bands (83%). Sequence analysis of rye fragments absent in triticale revealed for instance homology with hydroxyproline-rich glycoproteins (HRGP), a protein that belongs to a major family of inducible defence response proteins. Conversely, a wheat-specific band absent in triticale comprises a nested structure of copia-like retrotransposons elements, namely Claudia and Barbara. Sequencing of a polyploid-specific band (absent in both parents) revealed a microsatellite related sequence. Cytological studies using Fluorescent In Situ Hybridization (FISH) with REMAP products revealed a widespread distribution of retrotransposon and/or microsatellite flanking sequences on rye chromosomes, with a preferential accumulation in heterochromatic sub-telomeric domains.

CONCLUSIONS/SIGNIFICANCE

Here, we used PCR-based molecular marker techniques involving retrotransposons and microsatellites to uncover polyploidization induced genetic restructuring in triticale. Sequence analysis of rearranged genomic fragments either from rye or wheat origin showed these to be retrotransposon-related as well as coding sequences. Further FISH analysis revealed possible chromosome hotspots for sequence rearrangements. The role of chromatin condensation on the origin of genomic rearrangements mediated by polyploidization in triticale is also discussed.

Variety discrimination in pea (Pisum sativum L.) by molecular, biochemical and morphological markers

The distinctness, uniformity and stability (DUS) requirements involve expensive, space- and time-consuming measurements of morphological traits. Moreover, for a majority of traits, interactions between genotype and environment complicate the evaluation. Molecular markers have a potential to facilitate this procedure, increase the reliability of decisions, and substantially save the time and space needed for experiments. We chose 25 varieties of pea (Pisum sativum L.) from the list of recommended varieties for cultivation in the Czech Republic, and made both a standard classification by 12 morphological descriptors and a classification by biochemical-molecular markers. Two isozyme systems, 10 microsatellite loci, 2 retrotransposons for multilocus inter-retrotransposon amplified polymorphism (IRAP), and 12 retrotransposon-based insertion polymorphism (RBIP) DNA markers were analysed. The main objective of the study was to examine the potential of each method for discrimination between pea varieties. The results demonstrate a high potential and resolving power of DNA-based methods. Superior in terms of high information content and discrimination power were SSR markers, owing to high allelic variation, which was the only biochemical-molecular method allowing clear identification of all varieties. Retrotransposon markers in RBIP format proved to be the most robust and easy to score method, while multilocus IRAP produced informative fingerprint already in a single analysis. Isozyme analysis offered a fast and less expensive alternative. The results showed that molecular identification could be used to assess distinctness and complement morphological assessment, especially in cases where the time frame plays an important role. Currently developed pea marker systems might serve also for germplasm management and genetic diversity studies.

Assessment of genetic and epigenetic stability in long-term in vitro shoot culture of pea (Pisum sativum L.)

In vitro clonal propagation of plants should generate identical copies of the selected genotype. However, associated stress might result in a breakdown of control mechanisms and consequent instability of the genome. We have used several molecular methods to assess the genetic stability of long-term propagated (24 years) multiple shoot in vitro culture of pea (Pisum sativum L.). We focused on assessing the stability of repetitive sequences, such as simple sequence repeats (SSR) and retrotransposons, both comprising a large part of genome. No differences were found when seedlings (Co-2004) or original seed (Co-1982) controls and long-term or newly established in vitro (one subculture cycle) samples were investigated by the SSR, inter-repeats (ISSR) or inter-retrotransposon amplified polymorphism (IRAP) method. However, the more global amplified fragment length polymorphism (AFLP) and particularly the methylation sensitive MSAP methods detected 11 and 18% polymorphism among samples, respectively. Interestingly, investigation of the global cytosine methylation status by HPCE measurement revealed no statistically significant differences. Some evidence of retrotransposon re-arrangement was observed by sequence-specific amplification polymorphism. This occurred mostly in the abundant Ty3-gypsy type Cyclop element and to a smaller extent in the Ogre element. Alternatively, no polymorphism was detected among the PDR-1 element of the Ty1-copia type retrotransposon. Based on these results, multiple shoot culture of pea maintained over a long period may be considered as a true to type multiplication method of the original genotype.

IRAP and REMAP for retrotransposon-based genotyping and fingerprinting

Retrotransposons can be used as markers because their integration creates new joints between genomic DNA and their conserved ends. To detect polymorphisms for retrotransposon insertion, marker systems generally rely on PCR amplification between these ends and some component of flanking genomic DNA. We have developed two methods, retrotransposon-microsatellite amplified polymorphism (REMAP) analysis and inter-retrotransposon amplified polymorphism (IRAP) analysis, that require neither restriction enzyme digestion nor ligation to generate the marker bands. The IRAP products are generated from two nearby retrotransposons using outward-facing primers. In REMAP, amplification between retrotransposons proximal to simple sequence repeats (microsatellites) produces the marker bands. Here, we describe protocols for the IRAP and REMAP techniques, including methods for PCR amplification with a single primer or with two primers and for agarose gel electrophoresis of the product using optimal electrophoresis buffers and conditions. This protocol can be completed in 1-2 d.