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

Retrotransposons: How the continuous evolutionary front shapes plant genomes for response to heat stress

Long terminal repeat retrotransposons (LTR retrotransposons) are the most abundant group of mobile genetic elements in eukaryotic genomes and are essential in organizing genomic architecture and phenotypic variations. The diverse families of retrotransposons are related to retroviruses. As retrotransposable elements are dispersed and ubiquitous, their "copy-out and paste-in" life cycle of replicative transposition leads to new genome insertions without the excision of the original element. The overall structure of retrotransposons and the domains responsible for the various phases of their replication is highly conserved in all eukaryotes. The two major superfamilies of LTR retrotransposons, Ty1/Copia and Ty3/Gypsy, are distinguished and dispersed across the chromosomes of higher plants. Members of these superfamilies can increase in copy number and are often activated by various biotic and abiotic stresses due to retrotransposition bursts. LTR retrotransposons are important drivers of species diversity and exhibit great variety in structure, size, and mechanisms of transposition, making them important putative actors in genome evolution. Additionally, LTR retrotransposons influence the gene expression patterns of adjacent genes by modulating potential small interfering RNA (siRNA) and RNA-directed DNA methylation (RdDM) pathways. Furthermore, comparative and evolutionary analysis of the most important crop genome sequences and advanced technologies have elucidated the epigenetics and structural and functional modifications driven by LTR retrotransposon during speciation. However, mechanistic insights into LTR retrotransposons remain obscure in plant development due to a lack of advancement in high throughput technologies. In this review, we focus on the key role of LTR retrotransposons response in plants during heat stress, the role of centromeric LTR retrotransposons, and the role of LTR retrotransposon markers in genome expression and evolution.

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.

Genome-Wide Identification and Analysis of High-Copy-Number LTR Retrotransposons in Asian Pears

A large proportion of the genome of 'Suli' pear (Pyrus pyrifolia) contains long terminal repeat retrotransposons (LTR-RTs), which suggests that LTR-RTs have played important roles in the evolution of Pyrus. Further analysis of retrotransposons, particularly of high-copy-number LTR-RTs in different species, will provide new insights into the evolutionary history of Pyrus. A total of 4912 putative LTR-RTs classified into 198 subfamilies were identified in the 'Suli' pear genome. Six Asian pear accessions, including cultivars and wild species, were resequenced. The comparison of copy number for each LTR-RT subfamily was evaluated in Pyrus accessions, and data showed up to four-fold differences for some subfamilies. This contrast suggests different fates for retrotransposon families in the evolution of Pyrus. Fourteen high-copy-number subfamilies were identified in Asian pears, and more than 50% of the LTR-RTs in the genomes of all Pyrus accessions were from these 14 identified LTR-RT subfamilies. Their average insertion time was 3.42 million years ago, which suggests that these subfamilies were recently inserted into the genome. Many homologous and specific retrotransposon insertion sites were identified in oriental and occidental pears, suggesting that the duplication of retrotransposons has occurred throughout almost the entire origin and evolution of Pyrus species. The LTR-RTs show high heterogeneity, and their copy numbers vary in different Pyrus species. Thus, our findings suggest that LTR-RTs are an important source of genetic variation among Pyrus species.

Chromosome level high-density integrated genetic maps improve the Pyrus bretschneideri 'DangshanSuli' v1.0 genome

Background

Chromosomal level reference genomes provide a crucial foundation for genomics research such as genome-wide association studies (GWAS) and whole genome selection. The chromosomal-level sequences of both the European (Pyrus communis) and Chinese (P. bretschneideri) pear genomes have not been published in public databases so far.

Results

To anchor the scaffolds of P. bretschneideri ‘DangshanSuli’ (DS) v1.0 genome into pseudo-chromosomes, two genetic maps (MH and YM maps) were constructed using half sibling populations of Chinese pear crosses, ‘Mantianhong’ (MTH) × ‘Hongxiangsu’ (HXS) and ‘Yuluxiang’ (YLX) × MTH, from 345 and 162 seedlings, respectively, which were prepared for SNP discovery using genotyping-by-sequencing (GBS) technology. The MH and YM maps, each with 17 linkage groups (LGs), were constructed from 2606 and 2489 SNP markers and spanned 1847 and 1668 cM, respectively, with average marker intervals of 0.7. The two maps were further merged with a previously published genetic map (BD) based on the cross ‘Bayuehong’ (BYH) × ‘Dangshansuli’ (DS) to build a new integrated MH-YM-BD map. By using 7757 markers located on the integrated MH-YM-BD map, 898 scaffolds (400.57 Mb) of the DS v1.0 assembly were successfully anchored into 17 pseudo-chromosomes, accounting for 78.8% of the assembled genome size. About 88.31% of them (793 scaffolds) were directionally anchored with two or more markers on the pseudo-chromosomes. Furthermore, the errors in each pseudo-chromosome (especially 1, 5, 7 and 11) were manually corrected and pseudo-chromosomes 1, 5 and 7 were extended by adding 19, 12 and 14 scaffolds respectively in the newly constructed DS v1.1 genome. Synteny analyses revealed that the DS v1.1 genome had high collinearity with the apple genome, and the homologous fragments between pseudo-chromosomes were similar to those found in previous studies. Moreover, the red-skin trait of Asian pear was mapped to an identical locus as identified previously.

Conclusions

The accuracy of DS v1.1 genome was improved by using larger mapping populations and merged genetic map. With more than 400 MB anchored to 17 pseudo-chromosomes, the new DS v1.1 genome provides a critical tool that is essential for studies of pear genetics, genomics and molecular breeding.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5224-6) contains supplementary material, which is available to authorized users.

Use of retrotransposon-derived genetic markers to analyse genomic variability in plants

Transposable elements (TEs) are common mobile genetic elements comprising several classes and making up the majority of eukaryotic genomes. The movement and accumulation of TEs has been a major force shaping the genes and genomes of most organisms. Most eukaryotic genomes are dominated by retrotransposons and minimal DNA transposon accumulation. The 'copy and paste' lifecycle of replicative transposition produces new genome insertions without excising the original element. Horizontal TE transfer among lineages is rare. TEs represent a reservoir of potential genomic instability and RNA-level toxicity. Many TEs appear static and nonfunctional, but some are capable of replicating and mobilising to new positions, and somatic transposition events have been observed. The overall structure of retrotransposons and the domains responsible for the phases of their replication are highly conserved in all eukaryotes. TEs are important drivers of species diversity and exhibit great variety in their structure, size and transposition mechanisms, making them important putative actors in evolution. Because TEs are abundant in plant genomes, various applications have been developed to exploit polymorphisms in TE insertion patterns, including conventional or anchored PCR, and quantitative or digital PCR with primers for the 5' or 3' junction. Alternatively, the retrotransposon junction can be mapped using high-throughput next-generation sequencing and bioinformatics. With these applications, TE insertions can be rapidly, easily and accurately identified, or new TE insertions can be found. This review provides an overview of the TE-based applications developed for plant species and assesses the contributions of TEs to the analysis of plants' genetic diversity.

Combined Analyses of Chloroplast DNA Haplotypes and Microsatellite Markers Reveal New Insights Into the Origin and Dissemination Route of Cultivated Pears Native to East Asia

Asian pear plays an important role in the world pear industry, accounting for over 70% of world total production volume. Commercial Asian pear production relies on four major pear cultivar groups, Japanese pear (JP), Chinese white pear (CWP), Chinese sand pear (CSP), and Ussurian pear (UP), but their origins remain controversial. We estimated the genetic diversity levels and structures in a large sample of existing local cultivars to investigate the origins of Asian pears using twenty-five genome-covering nuclear microsatellite (simple sequence repeats, nSSR) markers and two non-coding chloroplast DNA (cpDNA) regions (trnL-trnF and accD-psaI). High levels of genetic diversity were detected for both nSSRs (H_E = 0.744) and cpDNAs (Hd = 0.792). The major variation was found within geographic populations of cultivated pear groups, demonstrating a close relationship among cultivar groups. CSPs showed a greater genetic diversity than CWPs and JPs, and lowest levels of genetic differentiation were detected among them. Phylogeographical analyses indicated that the CSP, CWP, and JP were derived from the same progenitor of Pyrus pyrifolia in China. A dissemination route of cultivated P. pyrifolia estimated by approximate Bayesian computation suggested that cultivated P. pyrifolia from the Middle Yangtze River Valley area contributed the major genetic resources to the cultivars, excluding those of southwestern China. Three major genetic groups of cultivated Pyrus pyrifolia were revealed using nSSRs and a Bayesian statistical inference: (a) JPs; (b) cultivars from South-Central China northward to northeastern China, covering the main pear production area in China; (c) cultivars from southwestern China to southeastern China, including Yunnan, Guizhou, Guangdong, Guangxi, and Fujian Provinces. This reflected the synergistic effects of ecogeographical factors and human selection during cultivar spread and improvement. The analyses indicated that UP cultivars might be originated from the interspecific hybridization of wild Pyrus ussuriensis with cultivated Pyrus pyrifolia. The combination of uniparental DNA sequences and nuclear markers give us a better understanding of origins and genetic relationships for Asian pear groups and will be beneficial for the future improvement of Asian pear cultivars.

Primitive Genepools of Asian Pears and Their Complex Hybrid Origins Inferred from Fluorescent Sequence-Specific Amplification Polymorphism (SSAP) Markers Based on LTR Retrotransposons

Recent evidence indicated that interspecific hybridization was the major mode of evolution in Pyrus. The genetic relationships and origins of the Asian pear are still unclear because of frequent hybrid events, fast radial evolution, and lack of informative data. Here, we developed fluorescent sequence-specific amplification polymorphism (SSAP) markers with lots of informative sites and high polymorphism to analyze the population structure among 93 pear accessions, including nearly all species native to Asia. Results of a population structure analysis indicated that nearly all Asian pear species experienced hybridization, and originated from five primitive genepools. Four genepools corresponded to four primary Asian species: P. betulaefolia, P. pashia, P. pyrifolia, and P. ussuriensis. However, cultivars of P. ussuriensis were not monophyletic and introgression occurred from P. pyrifolia. The specific genepool detected in putative hybrids between occidental and oriental pears might be from occidental pears. The remaining species, including P. calleryana, P. xerophila, P. sinkiangensis, P. phaeocarpa, P. hondoensis, and P. hopeiensis in Asia, were inferred to be of hybrid origins and their possible genepools were identified. This study will be of great help for understanding the origin and evolution of Asian pears.

Isolation and characterization of putative functional long terminal repeat retrotransposons in the Pyrus genome

Background

Long terminal repeat (LTR)-retrotransposons constitute 42.4 % of the genome of the ‘Suli’ pear (Pyrus pyrifolia white pear group), implying that retrotransposons have played important roles in Pyrus evolution. Therefore, further analysis of retrotransposons will enhance our understanding of the evolutionary history of Pyrus.

Results

We identified 1836 LTR-retrotransposons in the ‘Suli’ pear genome, of which 440 LTR-retrotransposons were predicted to contain at least two of three gene models (gag, integrase and reverse transcriptase). Because these were most likely to be functional transposons, we focused our analyses on this set of 440. Most of the LTR-retrotransposons were estimated to have inserted into the genome less than 2.5 million years ago. Sequence analysis showed that the reverse transcriptase component of the identified LTR-retrotransposons was highly heterogeneous. Analyses of transcripts assembled from RNA-Seq databases of two cultivars of Pyrus species showed that LTR-retrotransposons were expressed in the buds and fruit of Pyrus. A total of 734 coding sequences in the ‘Suli’ genome were disrupted by the identified LTR-retrotransposons. Five high-copy-number LTR-retrotransposon families were identified in Pyrus. These families were rarely found in the genomes of Malus and Prunus, but were distributed extensively in Pyrus and abundance varied between species.

Conclusions

We identified potentially functional, full-length LTR-retrotransposons with three gene models in the ‘Suli’ genome. The analysis of RNA-seq data demonstrated that these retrotransposons are expressed in the organs of pears. The differential copy number of LTR-retrotransposon families between Pyrus species suggests that the transposition of retrotransposons is an important evolutionary force driving the genetic divergence of species within the genus.

Electronic supplementary material

The online version of this article (doi:10.1186/s13100-016-0058-8) contains supplementary material, which is available to authorized users.