IRAPs in Combination with Highly Informative ISSRs Confer Effective Potentials for Genetic Diversity and Fidelity Assessment in Rhododendron.

The species belonging to the Rhododendron genus are well-known for their colorful corolla. Molecular marker systems have the potential to elucidate genetic diversity as well as to assess genetic fidelity in rhododendrons. In the present study, the reverse transcription domains of long terminal repeat retrotransposons were cloned from rhododendrons and used to develop an inter-retrotransposon amplified polymorphism (IRAP) marker system. Subsequently, 198 polymorphic loci were generated from the IRAP and inter-simple sequence repeat (ISSR) markers, of which 119 were derived from the IRAP markers. It was shown that in rhododendrons, IRAP markers were superior to the ISSRs in some polymorphic parameters, such as the average number of polymorphic loci (14.88 versus 13.17). The combination of the IRAP and ISSR systems was more discriminative for detecting 46 rhododendron accessions than each of the systems on their own. Furthermore, IRAP markers demonstrated more efficiency in genetic fidelity detection of in-vitro-grown R. bailiense Y.P.Ma, C.Q.Zhang and D.F.Chamb, an endangered species recently recorded in Guizhzhou Province, China. The available evidence revealed the distinct properties of IRAP and ISSR markers in the rhododendron-associated applications, and highlighted the availability of highly informative ISSR and IRAP markers in the evaluation of genetic diversity and genetic fidelity of rhododendrons, which may facilitate preservation and genetic breeding of rhododendron plants.

Assessment of genetic diversity among 131 safflower (Carthamus tinctorius L.) accessions using peroxidase gene polymorphism (POGP) markers.

BACKGROUND: Safflower (Carthamus tinctorius L.) is an old oilseed crop with a 1.4 GB genome size and its flowers are used for food coloring, dyes and pharmaceutical industries. It was domesticated from its putative wild ancestor Carthamus palestinus about forty-five hundred years ago in the fertile crescent region.The current study was aimed to determine the genetic diversity, population structure and to check the applicability of iPBS-retrotransposons markers. METHODS AND RESULTS: Eleven POGP primers yielded 70 bands of which 61 were highly polymorphic with 87.14% polymorphism. A great level of genetic variation was examined with higher values of overall gene diversity (0.27), genetic distance (0.53), number of effective alleles (1.46), Shannon's information index (0.41) and polymorphism information contents (0.71). Analysis of molecular variance revealed high genetic variation with 79% within the population. The STRUCTURE, PCoA and Neighbor-joining analysis separated the safflower germplasm into 2 major populations and 1 un-classified population. The accessions which were from Asian countries i.e., China, Afghanistan, Turkey, Iran and Pakistan were genetically similar and clustered together in both populations A and B. The maximum genetic distance was measured 0.88 between Pakistan 26 x Pakistan 24. CONCLUSION: Findings of this research such as maximum diversity indices, higher PIC values showed the effectiveness and utility of POGP markers for the evaluation of genetic relationships among safflower accessions. The results of this study also showed that POGP markers are less effective compared to ISSRs, iPBS-retrotransposons and DArTSeq markers. AMOVA showed high genetic variation (79%) within a population and maximum genetic distance was found between the accessions Pakistan 26- Pakistan 24 and may be suggested as candidate parents for future breeding activities of safflower. The accessions from the fertile crescent region were clustered together and proved the origin of safflower domestication. This study highlights genetic variation among safflower germplasm and could be helpfull for parental selection and planning for future breeding programs.

Applicability of inter-primer binding site iPBS- retrotransposon marker system for the assessment of genetic diversity and population structure of Peruvian rosewood (Aniba rosaeodora Ducke) germplasm.

BACKGROUND: Rosewood (Aniba rosaeodora Ducke), which has a great demand due to its essential oil globally, is an evergreen tree of the Amazon forests. Rosewood natural stands have been depleted through deforestation and the destruction of habitat. Currently, rosewood is included in the ICUN red list of endangered species. METHODS AND RESULTS: The 11 highly polymorphic primers amplified total 305 bands of which 301 (98.69%) were polymorphic. The number of effective alleles (Ne), Shannon's information index (I), overall gene diversity (Ht), gene diversity (h), and polymorphism information content (PIC) were (1.562), (0.505), (0.330), (0.337) and (0.343), respectively. These diversity indices explored high genetic diversity in rosewood germplasm. Among studied germplasm, the Santa Marta population was found most diverse by reflecting higher values of diversity indices while the Zungarococha population was found least diverse. The analysis of molecular variance (AMOVA) revealed that 79% of the genetic variation was within the populations. The STRUCTURE algorithm, unweighted pair group with arithmetic mean (UPGMA), and principal coordinate's analysis (PCoA) separated all germplasms into different population groups according to their geographic locations. Santa Marta population was found more diverse by reflecting higher values of diversity indices. The maximum genetic distance (0.868) was found between the Huajoya-10 and Nanay-3. In this investigation, iPBS- retrotransposon marker system was used to explore the genetic diversity of Peruvian rosewood germplasm. CONCLUSIONS: The results in this study such as higher genetic diversity indices, AMOVA (79%) within population and PIC value (0.343) showed the utility and reproducibility of iPBS-retrotransposons in this species successfully. The STRUCTURE algorithm separated the germplasms into six population groups according to their geographic locations. These results have valuable information for the conservation, management strategies and future breeding activities of rosewood.

Assessment of genetic diversity and variety identification based on developed retrotransposon-based insertion polymorphism (RBIP) markers in sweet potato (Ipomoea batatas (L.) Lam.).

Sweet potato, a dicotyledonous and perennial plant, is the third tuber/root crop species behind potato and cassava in terms of production. Long terminal repeat (LTR) retrotransposons are highly abundant in sweet potato, contributing to genetic diversity. Retrotransposon-based insertion polymorphism (RBIP) is a high-throughput marker system to study the genetic diversity of plant species. To date, there have been no transposon marker-based genetic diversity analyses of sweet potato. Here, we reported a structure-based analysis of the sweet potato genome, a total of 21555 LTR retrotransposons, which belonged to the main LTR-retrotransposon subfamilies Ty3-gypsy and Ty1-copia were identified. After searching and selecting using Hidden Markov Models (HMMs), 1616 LTR retrotransposon sequences containing at least two models were screened. A total of 48 RBIP primers were synthesized based on the high copy numbers of conserved LTR sequences. Fifty-six amplicons with an average polymorphism of 91.07% were generated in 105 sweet potato germplasm resources based on RBIP markers. A Unweighted Pair Group Method with Arithmatic Mean (UPGMA) dendrogram, a model-based genetic structure and principal component analysis divided the sweet potato germplasms into 3 groups containing 8, 53, and 44 germplasms. All the three analyses produced significant groupwise consensus. However, almost all the germplasms contained only one primary locus. The analysis of molecular variance (AMOVA) among the groups indicated higher intergroup genetic variation (53%) than intrapopulation genetic variation. In addition, long-term self-retention may cause some germplasm resources to exhibit variable segregation. These results suggest that these sweet potato germplasms are not well evolutionarily diversified, although geographic speciation could have occurred at a limited level. This study highlights the utility of RBIP markers for determining the intraspecies variability of sweet potato and have the potential to be used as core primer pairs for variety identification, genetic diversity assessment and linkage map construction. The results could provide a good theoretical reference and guidance for germplasm research and breeding.

Retrotransposable Elements: DNA Fingerprinting and the Assessment of Genetic Diversity.

Retrotransposable elements (RTEs) are highly common mobile genetic elements that are composed of several classes and make up the majority of eukaryotic genomes. The "copy-out and paste-in" life cycle of replicative transposition in these dispersive and ubiquitous RTEs leads to new genome insertions without excision of the original element. RTEs are important drivers of species diversity; they exhibit great variety in structure, size, and mechanisms of transposition, making them important putative components in genome evolution. Accordingly, various applications have been developed to explore the polymorphisms in RTE insertion patterns. These applications include conventional or anchored polymerase chain reaction (PCR) and quantitative or digital PCR with primers designed for the 5' or 3' junction. Marker systems exploiting these PCR methods can be easily developed and are inexpensively used in the absence of extensive genome sequence data. The main inter-repeat amplification polymorphism techniques include inter-retrotransposon amplified polymorphism (IRAP), retrotransposon microsatellite amplified polymorphism (REMAP), and Inter-Primer Binding Site (iPBS) for PCR amplification with a single or two primers.

Genomic re-assessment of the transposable element landscape of the potato genome.

KEY MESSAGE: We provide a comprehensive and reliable potato TE landscape, based on a wide variety of identification tools and integrative approaches, producing clear and ready-to-use outputs for the scientific community. Transposable elements (TEs) are DNA sequences with the ability to autoreplicate and move throughout the host genome. TEs are major drivers in stress response and genome evolution. Given their significance, the development of clear and efficient TE annotation pipelines has become essential for many species. The latest de novo TE discovery tools, along with available TEs from Repbase and sRNA-seq data, allowed us to perform a reliable potato TEs detection, classification and annotation through an open-source and freely available pipeline ( https://github.com/DiegoZavallo/TE_Discovery ). Using a variety of tools, approaches and rules, we were able to provide a clearly annotated of characterized TEs landscape. Additionally, we described the distribution of the different types of TEs across the genome, where LTRs and MITEs present a clear clustering pattern in pericentromeric and subtelomeric/telomeric regions respectively. Finally, we analyzed the insertion age and distribution of LTR retrotransposon families which display a distinct pattern between the two major superfamilies. While older Gypsy elements concentrated around heterochromatic regions, younger Copia elements located predominantly on euchromatic regions. Overall, we delivered not only a reliable, ready-to-use potato TE annotation files, but also all the necessary steps to perform de novo detection for other species.

High-throughput retrotransposon-based genetic diversity of maize germplasm assessment and analysis.

Maize is one of the world's most important crops and a model for grass genome research. Long terminal repeat (LTR) retrotransposons comprise most of the maize genome; their ability to produce new copies makes them efficient high-throughput genetic markers. Inter-retrotransposon-amplified polymorphisms (IRAPs) were used to study the genetic diversity of maize germplasm. Five LTR retrotransposons (Huck, Tekay, Opie, Ji, and Grande) were chosen, based on their large number of copies in the maize genome, whereas polymerase chain reaction primers were designed based on consensus LTR sequences. The LTR primers showed high quality and reproducible DNA fingerprints, with a total of 677 bands including 392 polymorphic bands showing 58% polymorphism between maize hybrid lines. These markers were used to identify genetic similarities among all lines of maize. Analysis of genetic similarity was carried out based on polymorphic amplicon profiles and genetic similarity phylogeny analysis. This diversity was expected to display ecogeographical patterns of variation and local adaptation. The clustering method showed that the varieties were grouped into three clusters differing in ecogeographical origin. Each of these clusters comprised divergent hybrids with convergent characters. The clusters reflected the differences among maize hybrids and were in accordance with their pedigree. The IRAP technique is an efficient high-throughput genetic marker-generating method.

On the Search for Retrotransposons: Alternative Protocols to Obtain Sequences to Learn Profile Hidden Markov Models.

Profile hidden Markov models (pHMMs) have been used to search for transposable elements (TEs) in genomes. For the learning of pHMMs aimed to search for TEs of the retrotransposon class, the conventional protocol is to use the whole internal nucleotide portions of these elements as representative sequences. To further explore the potential of pHMMs in such a search, we propose five alternative ways to obtain the sets of representative sequences of TEs other than the conventional protocol. In this study, we are interested in Bel-PAO, Copia, Gypsy, and DIRS superfamilies from the retrotransposon class. We compared the pHMMs of all six protocols. The test results show that, for each TE superfamily, the pHMMs of at least two of the proposed protocols performed better than the conventional one and that the number of correct predictions provided by the latter can be improved by considering together the results of one or more of the alternative protocols.

Establishment of a genome-wide and quantitative protocol for assessment of transcriptional activity at human retrotransposon L1 antisense promoters.

Long interspersed element 1 (L1) retrotransposon sequences are widespread in the human genome, occupying ~500,000 locations. The majority of L1s have lost their retrotransposition capability, although a significant population of human L1s maintains bidirectional transcriptional activity from the internal promoter. While the sense promoter drives transcription of the entire L1 mRNA and leads to L1 retrotransposition, the antisense promoter (ASP) transcribes L1-gene chimeric RNAs that include neighboring exon sequences. Activation mechanisms and functional impacts of L1ASP transcription are thought to vary at every L1ASP location. To explore the locus-specific regulation and function of L1ASP transcription, quantitative methodology is necessary for identifying the genomic positions of highly active L1ASPs on a genome-wide scale. Here, we employed deep-sequencing techniques and built a 3' RACE-based experimental and bioinformatics protocol, named the L1 antisense transcriptome protocol (LATRAP). In LATRAP, the PCR primer and the read mapping scheme were designed to reduce false positives and negatives, which may have been included as hits in previous cloning studies. LATRAP was here applied to the A549 human lung cancer cell line, and 313 L1ASP loci were detected to have transcriptional activity but differed in the number of mapped reads by four orders of magnitude. This indicates that transcriptional activities of the individual L1ASPs can vary greatly and that only a small population of L1ASP loci is active within individual nuclei. LATRAP is the first experimental method for ranking L1ASPs according to their transcriptional activity and will thus open a new avenue to unveiling the locus-specific biology of L1ASPs.

Development and validation of InnoQuant® HY, a system for quantitation and quality assessment of total human and male DNA using high copy targets.

The development and validation of InnoQuant® HY, a real-time PCR system containing four DNA targets-two RE autosomal targets of different sizes, male specific targets, and an internal positive control target-are described herein. The ratio of the two autosomal targets provides a Degradation Index, or a quantitative value of a sample's degradation state. The male specific targets are multi-copy targets located on the Y chromosome, which provides information about a sample's male DNA composition. The experimental results demonstrate InnoQuant HY as a robust qPCR method producing accurate DNA quantitation results even at low dynamic ranges, with reproducibility among population groups. The system is human specific with low level higher primate cross reactivity and is able to consistently and reproducibly detect sub-picogram concentrations of human and human male DNA. The use of high copy number Alu and SVA (>1000 copies per genome) retrotransposable elements as the two autosomal targets significantly enhances both sensitivity and reproducibility of determination of DNA quantitation as well as DNA degradation in forensic samples. The inclusion of a sensitive multi-copy Y-chromosome specific target provides accurate quantitation of DNA from a male in challenging male-female mixtures (i.e. sexual assault samples). Even in the presence of a large excess of DNA from a female, accurate quantitation was achieved with a male to female ratio of 1:128,000. Population database studies reveal an average Short/Y target ratio of the quantification values across all four populations tested was 1.124±0.282, exhibiting the system's reproducibility across multiple populations. The results from InnoQuant HY provide a tool equipping a forensic analyst with crucial data about a sample's DNA quantitation, male:female ratio, degradation state, and the presence or absence of PCR inhibitors. With the information gained from the InnoQuant HY kit, a more streamlined and efficient workflow can be created that minimizes unnecessary sample processing and retesting while maximizing recovery of probative DNA profiles from challenging biological evidence.

Assessment of genetic diversity in Nordic timothy (Phleum pratense L.).

BACKGROUND: Timothy (Phleum pratense L.), a cool-season hexaploid perennial, is the most important forage grass species in Nordic countries. Earlier analyses of genetic diversity in a collection of 96 genebank accessions of timothy with SSR markers demonstrated high levels of diversity but could not resolve population structure. Therefore, we examined a subset of 51 accessions with REMAP markers, which are based on retrotransposons, and compared the diversity results with those obtained with SSR markers. RESULTS: Using four primer combinations, 533 REMAP markers were analyzed, compared with 464 polymorphic alleles in the 13 SSR loci previously. The average marker index, which describes information obtained per experiment (per primer combination or locus) was over six times higher with REMAPs. Most of the variation found was within accessions, with somewhat less, 89 %, for REMAPs, than for SSR, with 93 %. CONCLUSIONS: SSRs revealed differences in the level of diversity slightly better than REMAPs but neither marker type could reveal any clear clustering of accessions based on countries, vegetation zones, or different cultivar types. In our study, reliable evaluation of SSR allele dosages was not possible, so each allele had to be handled as a dominant marker. SSR and REMAP, which report from different mechanisms of generating genetic diversity and from different genomic regions, together indicate a lack of population structure. Taken together, this likely reflects the outcrossing and hexaploid nature of timothy rather than failures of either marker system.

Assessment of genetic and epigenetic changes in virus-free garlic (Allium sativum L.) plants obtained by meristem culture followed by in vitro propagation.

KEY MESSAGE: This is the first report assessing epigenetic variation in garlic. High genetic and epigenetic polymorphism during in vitro culture was detected.Sequencing of MSAP fragments revealed homology with ESTs. Garlic (Allium sativum) is a worldwide crop of economic importance susceptible to viral infections that can cause significant yield losses. Meristem tissue culture is the most employed method to sanitize elite cultivars.Often the virus-free garlic plants obtained are multiplied in vitro (micro propagation). However, it was reported that micro-propagation frequently produces somaclonal variation at the phenotypic level, which is an undesirable trait when breeders are seeking to maintain varietal stability. We employed amplification fragment length polymorphism and methylation sensitive amplified polymorphism (MSAP) methodologies to assess genetic and epigenetic modifications in two culture systems: virus-free plants obtained by meristem culture followed by in vitro multiplication and field culture. Our results suggest that garlic exhibits genetic and epigenetic polymorphism under field growing conditions. However, during in vitro culture system both kinds of polymorphisms intensify indicating that this system induces somaclonal variation. Furthermore, while genetic changes accumulated along the time of in vitro culture, epigenetic polymorphism reached the major variation at 6 months and then stabilize, being demethylation and CG methylation the principal conversions.Cloning and sequencing differentially methylated MSAP fragments allowed us to identify coding and unknown sequences of A. sativum, including sequences belonging to LTR Gypsy retrotransposons. Together, our results highlight that main changes occur in the initial 6 months of micro propagation. For the best of our knowledge, this is the first report on epigenetic assessment in garlic.

Genetic bottlenecks in Turkish okra germplasm and utility of iPBS retrotransposon markers for genetic diversity assessment.

Lack of requisite genetic variation in Turkish okra has necessitated the use of different types of markers for estimating the genetic diversity and identifying the source of variation. Transposable elements, present abundantly in plant genomes, generate genomic diversity through their replication and are thus an excellent source of molecular markers. We hypothesized that inter-primer binding site (iPBS)-retrotransposons could be the source of variation because of their genome plasticity nature. In the present study, genetic diversity of 66 okra landraces was analyzed using iPBS-retrotransposon markers. iPBS-retrotransposons detected 88 bands with 40.2% polymorphism and an average of 6.8 bands per primer. Gene diversity and Shannon's information index ranged from 0.01 to 0.13 and 0.02 to 0.21 for iPBS-retrotransposons and from 0.06 to 0.46 and 0.14 to 0.65 for simple sequence repeat (SSR) markers, respectively. Polymorphism information content value for retrotransposons varied between 0.12 and 0.99, while that for SSR was from 0.52 to 0.81. Neighbor joining analysis based on retrotransposons and SSRs divided all the accessions into four clusters; however, SSR markers were more efficient in clustering the landraces based on their origin. Using the STRUCTURE software for determining population structure, and two populations (at the number of hypothetical subpopulations, K = 2) were identified among the landraces. Low genetic diversity in Turkish okra highlights the need for the introduction of plants from countries with greater genetic diversity for these crops. This study also demonstrates the utility and role of iPBS-retrotransposons, a dominant and ubiquitous part of eukaryotic genomes, for diversity studies in okra.

Terminal-repeat retrotransposons with GAG domain in plant genomes: a new testimony on the complex world of transposable elements.

A novel structure of nonautonomous long terminal repeat (LTR) retrotransposons called terminal repeat with GAG domain (TR-GAG) has been described in plants, both in monocotyledonous, dicotyledonous and basal angiosperm genomes. TR-GAGs are relatively short elements in length (<4 kb) showing the typical features of LTR-retrotransposons. However, they carry only one open reading frame coding for the GAG precursor protein involved for instance in transposition, the assembly, and the packaging of the element into the virus-like particle. GAG precursors show similarities with both Copia and Gypsy GAG proteins, suggesting evolutionary relationships of TR-GAG elements with both families. Despite the lack of the enzymatic machinery required for their mobility, strong evidences suggest that TR-GAGs are still active. TR-GAGs represent ubiquitous nonautonomous structures that could be involved in the molecular diversities of plant genomes.

Prediction of retrotransposons and assessment of genetic variability based on developed retrotransposon-based insertion polymorphism (RBIP) markers in Pyrus L.

Interspecific hybridization has been considered the major mode of evolution in Pyrus (pear), and thus, the genetic relationships within this genus have not been well documented. Retrotransposons are ubiquitous components of plant genomes and 42.4 % of the pear genome was reported to be long terminal repeat (LTR) retrotransposons, implying that retrotransposons might be significant in the evolution of Pyrus. In this study, 1,836 putative full-length LTR retrotransposons were isolated and 196 retrotransposon-based insertion polymorphism (RBIP) primers were developed, of which 24 pairs to the Ppcr1 subfamily of copia retrotransposons were used to analyze genetic diversity among 110 Pyrus accessions from Eurasia. Our results showed that Ppcr1 replicated many times in the development of cultivated Asian pears. The genetic structure analysis and the unweighted pair group method with arithmetic mean (UPGMA) dendrogram indicated that all accessions could be divided into Oriental and Occidental groups. In Oriental pears, wild pea pears clustered separately into independent groups in accordance with their morphological classifications. Cultivars of P. ussuriensis Maxim, P. pyrifolia Nakai, and P. pyrifolia Chinese white pear were mingled together, which inferred that hybridization events occurred during the development of the cultivated Asian pears. In Occidental pears, two clades were obtained in the UPGMA dendrogram in accordance with their geographical distribution; one contained the European species and the other included species from North Africa and West Asia. New findings in this study will be important to further understand the phylogeny of Pyrus and origins of cultivated pears.

Development and validation of InnoQuant™, a sensitive human DNA quantitation and degradation assessment method for forensic samples using high copy number mobile elements Alu and SVA.

There is a constant need in forensic casework laboratories for an improved way to increase the first-pass success rate of forensic samples. The recent advances in mini STR analysis, SNP, and Alu marker systems have now made it possible to analyze highly compromised samples, yet few tools are available that can simultaneously provide an assessment of quantity, inhibition, and degradation in a sample prior to genotyping. Currently there are several different approaches used for fluorescence-based quantification assays which provide a measure of quantity and inhibition. However, a system which can also assess the extent of degradation in a forensic sample will be a useful tool for DNA analysts. Possessing this information prior to genotyping will allow an analyst to more informatively make downstream decisions for the successful typing of a forensic sample without unnecessarily consuming DNA extract. Real-time PCR provides a reliable method for determining the amount and quality of amplifiable DNA in a biological sample. Alu are Short Interspersed Elements (SINE), approximately 300bp insertions which are distributed throughout the human genome in large copy number. The use of an internal primer to amplify a segment of an Alu element allows for human specificity as well as high sensitivity when compared to a single copy target. The advantage of an Alu system is the presence of a large number (>1000) of fixed insertions in every human genome, which minimizes the individual specific variation possible when using a multi-copy target quantification system. This study utilizes two independent retrotransposon genomic targets to obtain quantification of an 80bp "short" DNA fragment and a 207bp "long" DNA fragment in a degraded DNA sample in the multiplex system InnoQuant™. The ratio of the two quantitation values provides a "Degradation Index", or a qualitative measure of a sample's extent of degradation. The Degradation Index was found to be predictive of the observed loss of STR markers and alleles as degradation increases. Use of a synthetic target as an internal positive control (IPC) provides an additional assessment for the presence of PCR inhibitors in the test sample. In conclusion, a DNA based qualitative/quantitative/inhibition assessment system that accurately predicts the status of a biological sample, will be a valuable tool for deciding which DNA test kit to utilize and how much target DNA to use, when processing compromised forensic samples for DNA testing.

Assessment of genetic diversity among Indian potato (Solanum tuberosum L.) collection using microsatellite and retrotransposon based marker systems.

Potato (Solanum tuberosum) is an important non-cereal crop throughout the world and is highly recommended for ensuring global food security. Owing to the complexities in genetics and inheritance pattern of potato, the conventional method of cross breeding for developing improved varieties has been difficult. Identification and tagging of desirable traits with informative molecular markers would aid in the development of improved varieties. Insertional polymorphism of copia-like and gypsy-like long terminal repeat retrotransposons (RTN) were investigated among 47 potato varieties from India using Inter-Retrotransposon Amplified Polymorphism (IRAP) and Retrotransposon Microsatellite Amplified Polymorphism (REMAP) marker techniques and were compared with the DNA profiles obtained with simple sequence repeats (SSRs). The genetic polymorphism, efficiency of polymorphism and effectiveness of marker systems were evaluated to assess the extent of genetic diversity among Indian potato varieties. A total of 139 polymorphic SSR alleles, 270 IRAP and 98 REMAP polymorphic bands, showing polymorphism of 100%, 87.9% and 68.5%, respectively, were used for detailed characterization of the genetic relationships among potato varieties by using cluster analysis and principal coordinate analysis (PCoA). IRAP analysis resulted in the highest number of polymorphic bands with an average of 15 polymorphic bands per assay unit when compared to the other two marker systems. Based on pair-wise comparison, the genetic similarity was calculated using Dice similarity coefficient. The SSRs showed a wide range in genetic similarity values (0.485-0.971) as compared to IRAP (0.69-0.911) and REMAP (0.713-0.947). A Mantel's matrix correspondence test showed a high positive correlation (r=0.6) between IRAP and REMAP, an intermediate value (r=0.58) for IRAP and SSR and the lowest value (r=0.17) for SSR and REMAP. Statistically significant cophenetic correlation coefficient values, of 0.961, 0.941 and 0.905 were observed for REMAP, IRAP and SSR, respectively. The widespread presence and distinct DNA profiles for copia-like and gypsy-like RTNs in the examined genotypes indicate that these elements are active in the genome and may have even contributed to the potato genome organization. Although the three marker systems were capable of distinguishing all the 47 varieties; high reproducibility, low cost and ease of DNA profiling data collection make IRAP and REMAP markers highly efficient whole-genome scanning molecular probes for population genetic studies. Information obtained from the present study regarding the genetic association and distinctiveness provides an useful guide for selection of germplasm for plant breeding and conservation efforts.

Footprint of APOBEC3 on the genome of human retroelements.

Almost half of the human genome is composed of transposable elements. The genomic structures and life cycles of some of these elements suggest they are a result of waves of retroviral infection and transposition over millions of years. The reduction of retrotransposition activity in primates compared to that in nonprimates, such as mice, has been attributed to the positive selection of several antiretroviral factors, such as apolipoprotein B mRNA editing enzymes. Among these, APOBEC3G is known to mutate G to A within the context of GG in the genome of endogenous as well as several exogenous retroelements (the underlining marks the G that is mutated). On the other hand, APOBEC3F and to a lesser extent other APOBEC3 members induce G-to-A changes within the nucleotide GA. It is known that these enzymes can induce deleterious mutations in the genome of retroviral sequences, but the evolution and/or inactivation of retroelements as a result of mutation by these proteins is not clear. Here, we analyze the mutation signatures of these proteins on large populations of long interspersed nuclear element (LINE), short interspersed nuclear element (SINE), and endogenous retrovirus (ERV) families in the human genome to infer possible evolutionary pressure and/or hypermutation events. Sequence context dependency of mutation by APOBEC3 allows investigation of the changes in the genome of retroelements by inspecting the depletion of G and enrichment of A within the APOBEC3 target and product motifs, respectively. Analysis of approximately 22,000 LINE-1 (L1), 24,000 SINE Alu, and 3,000 ERV sequences showed a footprint of GG→AG mutation by APOBEC3G and GA→AA mutation by other members of the APOBEC3 family (e.g., APOBEC3F) on the genome of ERV-K and ERV-1 elements but not on those of ERV-L, LINE, or SINE.

Retrogenes moved out of the z chromosome in the silkworm.

Previous studies on organisms with well-differentiated X and Y chromosomes, such as Drosophila and mammals, consistently detected an excess of genes moving out of the X chromosome and gaining testis-biased expression. Several selective evolutionary mechanisms were shown to be associated with this nonrandom gene traffic, which contributed to the evolution of the X chromosome and autosomes. If selection drives gene traffic, such traffic should also exist in species with Z and W chromosomes, where the females are the heterogametic sex. However, no previous studies on gene traffic in species with female heterogamety have found any nonrandom chromosomal gene movement. Here, we report an excess of retrogenes moving out of the Z chromosome in an organism with the ZW sex determination system, Bombyx mori. In addition, we showed that those "out of Z" retrogenes tended to have ovary-biased expression, which is consistent with the pattern of non-retrogene traffic recently reported in birds and symmetrical to the retrogene movement in mammals and fruit flies out of the X chromosome evolving testis functions. These properties of gene traffic in the ZW system suggest a general role for the heterogamety of sex chromosomes in determining the chromosomal locations and the evolution of sex-biased genes.

MGEScan-non-LTR: computational identification and classification of autonomous non-LTR retrotransposons in eukaryotic genomes.

Computational methods for genome-wide identification of mobile genetic elements (MGEs) have become increasingly necessary for both genome annotation and evolutionary studies. Non-long terminal repeat (non-LTR) retrotransposons are a class of MGEs that have been found in most eukaryotic genomes, sometimes in extremely high numbers. In this article, we present a computational tool, MGEScan-non-LTR, for the identification of non-LTR retrotransposons in genomic sequences, following a computational approach inspired by a generalized hidden Markov model (GHMM). Three different states represent two different protein domains and inter-domain linker regions encoded in the non-LTR retrotransposons, and their scores are evaluated by using profile hidden Markov models (for protein domains) and Gaussian Bayes classifiers (for linker regions), respectively. In order to classify the non-LTR retrotransposons into one of the 12 previously characterized clades using the same model, we defined separate states for different clades. MGEScan-non-LTR was tested on the genome sequences of four eukaryotic organisms, Drosophila melanogaster, Daphnia pulex, Ciona intestinalis and Strongylocentrotus purpuratus. For the D. melanogaster genome, MGEScan-non-LTR found all known 'full-length' elements and simultaneously classified them into the clades CR1, I, Jockey, LOA and R1. Notably, for the D. pulex genome, in which no non-LTR retrotransposon has been annotated, MGEScan-non-LTR found a significantly larger number of elements than did RepeatMasker, using the current version of the RepBase Update library. We also identified novel elements in the other two genomes, which have only been partially studied for non-LTR retrotransposons.