Analysis of plant LTR-retrotransposons at the fine-scale family level reveals individual molecular patterns

Frequency of cytosine methylation in the adjacent regions of soybean retrotransposon SORE-1 depends on chromosomal location

Mobilization of transposable elements (TEs) is suppressed by epigenetic mechanisms involving cytosine methylation. However, few studies have focused on clarifying relationships between epigenetic influences of TEs on the adjacent DNA regions and time after insertion of TEs into the genome and/or their chromosomal location. Here we addressed these issues using soybean retrotransposon SORE-1. We analyzed SORE-1, inserted in exon 1 of the GmphyA2 gene, one of the newest insertions in this family so far identified. Cytosine methylation was detected in this element but was barely present in the adjacent regions. These results were correlated, respectively, with the presence and absence of the production of short interfering RNAs. Cytosine methylation profiles of 74 SORE-1 elements in the Williams 82 reference genome indicated that methylation frequency in the adjacent regions of SORE-1 was profoundly higher in pericentromeric regions than in euchromatic chromosome arms and was only weakly correlated with the length of time after insertion into the genome. Notably, the higher level of methylation in the 5' adjacent regions of SORE-1 coincided with the presence of repetitive elements in pericentromeric regions. Together, these results suggest that epigenetic influence of SORE-1 on the adjacent regions is influenced by its location on the chromosome.

A target cultivar-specific identification system based on the chromatographic printed array strip method for eight prominent Japanese citrus cultivars

Citrus is a major cultivated crop in Japan, and new cultivars are of great interest in the Japanese and global market. Recently, the infringement of breeders' rights to citrus cultivars bred in Japan has become a problem related to the agricultural product export strategy promoted by the Japanese government. Cultivar identification systems using DNA markers are an effective tool for protecting breeders' rights. Here, a novel target cultivar-specific identification system using the chromatographic printed array strip method was developed for eight prominent Japanese citrus cultivars. A polymorphic InDel fragment specific to each cultivar was explored through the screening of published citrus InDel markers and next-generation sequencing of retrotransposon libraries. The cultivar-specific DNA marker set for each cultivar comprised 1-3 polymorphic InDel fragments in combination with a PCR-positive DNA marker for the ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit gene. The DNA markers were detected within 3 hours from DNA extraction to the detection by the C-PAS4 membrane stick following multiplex PCR. The developed system is superior as a convenient, rapid, and cost-effective DNA diagnostic method during inspection. The proposed target cultivar-specific identification system is expected to serve as an efficient tool for the injunction of suspicious registered cultivars, contributing to the protection of breeders' rights.

Isolation and sequencing of a single copy of an introgressed chromosome from a complex genome for gene and SNP identification

This manuscript describes the identification, isolation and sequencing of a single chromosome containing high value resistance genes from a complex polyploid where sequencing the whole genome is too costly. The large complex genomes of many crops constrain the use of new technologies for genome-assisted selection and genetic improvement. One method to simplify a genome is to break it into individual chromosomes by flow cytometry; however, in many crop species most chromosomes cannot be isolated individually. Flow sorting of a single copy of a chromosome has been developed in wheat, and here we demonstrate its use to identify markers of interest in an Erianthus/Sacchurum hybrid. Erianthus/Saccharum hybrids are of interest because Erianthus is known to be highly resistant to soil borne diseases which cause extensive sugarcane yield losses in Australia. Sugarcane (Saccharum) cultivars are autopolyploids with a highly complex genome and over 100 chromosomes. Flow cytometry for sugarcane, as in most crops, does not resolve individual chromosomes to a karyotype peak for sorting. To isolate a single chromosome, we used genomic in situ hybridization (GISH) to identify the flow karyotype region containing the Erianthus chromosomes, flow sorted single chromosomes from this region, PCR screened for the Erianthus chromosomes and sequenced them. One Erianthus chromosome amplified and sequenced well, and from this data we could identify 57 resistant type genes and SNPs in nearly half of these genes. We developed KASP SNP assays and demonstrated that the identified SNP markers segregated as expected in a small introgression population. The pipeline we developed here to flow sort and sequence single chromosomes could be used in any crop with a large complex genome to rapidly discover and develop markers to important loci.

Primer Binding Site (PBS) Profiling of Genetic Diversity of Natural Populations of Endemic Species Allium ledebourianum Schult

Endemic species are especially vulnerable to biodiversity loss caused by isolation or habitat specificity, small population size, and anthropogenic factors. Endemic species biodiversity analysis has a critically important global value for the development of conservation strategies. The rare onion Allium ledebourianum is a narrow-lined endemic species, with natural populations located in the extreme climatic conditions of the Kazakh Altai. A. ledebourianum populations are decreasing everywhere due to anthropogenic impact, and therefore, this species requires preservation and protection. Conservation of this rare species is associated with monitoring studies to investigate the genetic diversity of natural populations. Fundamental components of eukaryote genome include multiple classes of interspersed repeats. Various PCR-based DNA fingerprinting methods are used to detect chromosomal changes related to recombination processes of these interspersed elements. These methods are based on interspersed repeat sequences and are an effective approach for assessing the biological diversity of plants and their variability. We applied DNA profiling approaches based on conservative sequences of interspersed repeats to assess the genetic diversity of natural A. ledebourianum populations located in the territory of Kazakhstan Altai. The analysis of natural A. ledebourianum populations, carried out using the DNA profiling approach, allowed the effective differentiation of the populations and assessment of their genetic diversity. We used conservative sequences of tRNA primer binding sites (PBS) of the long-terminal repeat (LTR) retrotransposons as PCR primers. Amplification using the three most effective PBS primers generated 628 PCR amplicons, with an average of 209 amplicons. The average polymorphism level varied from 34% to 40% for all studied samples. Resolution analysis of the PBS primers showed all of them to have high or medium polymorphism levels, which varied from 0.763 to 0.965. Results of the molecular analysis of variance showed that the general biodiversity of A. ledebourianum populations is due to interpopulation (67%) and intrapopulation (33%) differences. The revealed genetic diversity was higher in the most distant population of A. ledebourianum LD64, located on the Sarymsakty ridge of Southern Altai. This is the first genetic diversity study of the endemic species A. ledebourianum using DNA profiling approaches. This work allowed us to collect new genetic data on the structure of A. ledebourianum populations in the Altai for subsequent development of preservation strategies to enhance the reproduction of this relict species. The results will be useful for the conservation and exploitation of this species, serving as the basis for further studies of its evolution and ecology.

Full-length LTR retroelements in Capsicum annuum revealed a few species-specific family bursts with insertional preferences

Capsicum annuum is a species that has undergone an expansion of the size of its genome caused mainly by the amplification of repetitive DNA sequences, including mobile genetic elements. Based on information obtained from sequencing the genome of pepper, the estimated fraction of retroelements is approximately 81%, and previous results revealed an important contribution of lineages derived from Gypsy superfamily. However, the dynamics of the retroelements in the C. annuum genome is poorly understood. In this way, the present work seeks to investigate the phylogenetic diversity and genomic abundance of the families of autonomous (complete and intact) LTR retroelements from C. annuum and inspect their distribution along its chromosomes. In total, we identified 1151 structurally full-length retroelements (340 Copia; 811 Gypsy) grouped in 124 phylogenetic families in the base of their retrotranscriptase. All the evolutive lineages of LTR retroelements identified in plants were present in pepper; however, three of them comprise 83% of the entire LTR retroelements population, the lineages Athila, Del/Tekay, and Ale/Retrofit. From them, only three families represent 70.8% of the total number of the identified retroelements. A massive family-specific wave of amplification of two of them occurred in the last 0.5 Mya (GypsyCa_16; CopiaCa_01), whereas the third is more ancient and occurred 3.0 Mya (GypsyCa_13). Fluorescent in situ hybridization performed with family and lineage-specific probes revealed contrasting patterns of chromosomal affinity. Our results provide a database of the populations LTR retroelements specific to C. annuum genome. The most abundant families were analyzed according to chromosome insertional preferences, suppling useful tools to the design of retroelement-based markers specific to the species.

InpactorDB: A Classified Lineage-Level Plant LTR Retrotransposon Reference Library for Free-Alignment Methods Based on Machine Learning

Long terminal repeat (LTR) retrotransposons are mobile elements that constitute the major fraction of most plant genomes. The identification and annotation of these elements via bioinformatics approaches represent a major challenge in the era of massive plant genome sequencing. In addition to their involvement in genome size variation, LTR retrotransposons are also associated with the function and structure of different chromosomal regions and can alter the function of coding regions, among others. Several sequence databases of plant LTR retrotransposons are available for public access, such as PGSB and RepetDB, or restricted access such as Repbase. Although these databases are useful to identify LTR-RTs in new genomes by similarity, the elements of these databases are not fully classified to the lineage (also called family) level. Here, we present InpactorDB, a semi-curated dataset composed of 130,439 elements from 195 plant genomes (belonging to 108 plant species) classified to the lineage level. This dataset has been used to train two deep neural networks (i.e., one fully connected and one convolutional) for the rapid classification of these elements. In lineage-level classification approaches, we obtain up to 98% performance, indicated by the F1-score, precision and recall scores.

The Cassandra retrotransposon landscape in sugar beet (Beta vulgaris) and related Amaranthaceae: recombination and re-shuffling lead to a high structural variability

BACKGROUND AND AIMS

Plant genomes contain many retrotransposons and their derivatives, which are subject to rapid sequence turnover. As non-autonomous retrotransposons do not encode any proteins, they experience reduced selective constraints leading to their diversification into multiple families, usually limited to a few closely related species. In contrast, the non-coding Cassandra terminal repeat retrotransposons in miniature (TRIMs) are widespread in many plants. Their hallmark is a conserved 5S rDNA-derived promoter in their long terminal repeats (LTRs). As sugar beet (Beta vulgaris) has a well-described LTR retrotransposon landscape, we aim to characterize TRIMs in beet and related genomes.

METHODS

We identified Cassandra retrotransposons in the sugar beet reference genome and characterized their structural relationships. Genomic organization, chromosomal localization, and distribution of Cassandra-TRIMs across the Amaranthaceae were verified by Southern and fluorescent in situ hybridization.

KEY RESULTS

All 638 Cassandra sequences in the sugar beet genome contain conserved LTRs and thus constitute a single family. Nevertheless, variable internal regions required a subdivision into two Cassandra subfamilies within B. vulgaris. The related Chenopodium quinoa harbours a third subfamily. These subfamilies vary in their distribution within Amaranthaceae genomes, their insertion times and the degree of silencing by small RNAs. Cassandra retrotransposons gave rise to many structural variants, such as solo LTRs or tandemly arranged Cassandra retrotransposons. These Cassandra derivatives point to an interplay of template switch and recombination processes - mechanisms that likely caused Cassandra's subfamily formation and diversification.

CONCLUSIONS

We traced the evolution of Cassandra in the Amaranthaceae and detected a considerable variability within the short internal regions, whereas the LTRs are strongly conserved in sequence and length. Presumably these hallmarks make Cassandra a prime target for unequal recombination, resulting in the observed structural diversity, an example of the impact of LTR-mediated evolutionary mechanisms on the host genome.

Structural and Functional Annotation of Transposable Elements Revealed a Potential Regulation of Genes Involved in Rubber Biosynthesis by TE-Derived siRNA Interference in Hevea brasiliensis

The natural rubber biosynthetic pathway is well described in Hevea, although the final stages of rubber elongation are still poorly understood. Small Rubber Particle Proteins and Rubber Elongation Factors (SRPPs and REFs) are proteins with major function in rubber particle formation and stabilization. Their corresponding genes are clustered on a scaffold1222 of the reference genomic sequence of the Hevea brasiliensis genome. Apart from gene expression by transcriptomic analyses, to date, no deep analyses have been carried out for the genomic environment of SRPPs and REFs loci. By integrative analyses on transposable element annotation, small RNAs production and gene expression, we analysed their role in the control of the transcription of rubber biosynthetic genes. The first in-depth annotation of TEs (Transposable Elements) and their capacity to produce TE-derived siRNAs (small interfering RNAs) is presented, only possible in the Hevea brasiliensis clone PB 260 for which all data are available. We observed that 11% of genes are located near TEs and their presence may interfere in their transcription at both genetic and epigenetic level. We hypothesized that the genomic environment of rubber biosynthesis genes has been shaped by TE and TE-derived siRNAs with possible transcriptional interference on their gene expression. We discussed possible functionalization of TEs as enhancers and as donors of alternative transcription start sites in promoter sequences, possibly through the modelling of genetic and epigenetic landscapes.

Inter-retrotransposon amplified polymorphism markers revealed long terminal repeat retrotransposon insertion polymorphism in flax cultivated on the experimental fields around Chernobyl

Ionizing radiation in environment comes from various natural and anthropogenic sources. The effect of radioactivity released after the CNPP (Chernobyl Nuclear Power Plant) on plant systems remains of great interest. Even now, more than three decades after the nuclear accident, the long-lived radionuclides represent a strong stress factor. Herein, the emphasis has been placed on analysis of genetic variability represented by activation of LTR (Long Terminal Repeat)-retrotransposons. Polymorphism in LTR-retrotransposon insertions has been investigated throughout the genome of two flax varieties, Kyivskyi and Bethune. For this purpose, two retrotransposon-based marker techniques, IRAP (Inter-Retrotransposon Amplified Polymorphism) and iPBS (inter-Primer Binding Site), have been employed. The hypothesis that chronic radioactive stress may induce mechanism of retransposition has been supported by the activation of FL9, FL11 and FL12 LTR-retrotransposons in flax seeds harvested from radioactive environment. Out of two retrotransposon-based approaches, IRAP appears to be more suitable for identification of LTR-retrotransposon polymorphism. Even though the LTR-retrotransposon polymorphism was identified, the results suggest the high level of plant adaptation in the radioactive Chernobyl area. However, it is not really surprising that plants developed an effective strategy to survive in radio-contaminated environment over the past 30 years.

Assembly of the 373k gene space of the polyploid sugarcane genome reveals reservoirs of functional diversity in the world's leading biomass crop

BACKGROUND

Sugarcane cultivars are polyploid interspecific hybrids of giant genomes, typically with 10-13 sets of chromosomes from 2 Saccharum species. The ploidy, hybridity, and size of the genome, estimated to have >10 Gb, pose a challenge for sequencing.

RESULTS

Here we present a gene space assembly of SP80-3280, including 373,869 putative genes and their potential regulatory regions. The alignment of single-copy genes in diploid grasses to the putative genes indicates that we could resolve 2-6 (up to 15) putative homo(eo)logs that are 99.1% identical within their coding sequences. Dissimilarities increase in their regulatory regions, and gene promoter analysis shows differences in regulatory elements within gene families that are expressed in a species-specific manner. We exemplify these differences for sucrose synthase (SuSy) and phenylalanine ammonia-lyase (PAL), 2 gene families central to carbon partitioning. SP80-3280 has particular regulatory elements involved in sucrose synthesis not found in the ancestor Saccharum spontaneum. PAL regulatory elements are found in co-expressed genes related to fiber synthesis within gene networks defined during plant growth and maturation. Comparison with sorghum reveals predominantly bi-allelic variations in sugarcane, consistent with the formation of 2 "subgenomes" after their divergence ∼3.8-4.6 million years ago and reveals single-nucleotide variants that may underlie their differences.

CONCLUSIONS

This assembly represents a large step towards a whole-genome assembly of a commercial sugarcane cultivar. It includes a rich diversity of genes and homo(eo)logous resolution for a representative fraction of the gene space, relevant to improve biomass and food production.

The population genetic structure approach adds new insights into the evolution of plant LTR retrotransposon lineages

Long terminal repeat retrotransposons (LTR-RTs) in plant genomes differ in abundance, structure and genomic distribution, reflecting the large number of evolutionary lineages. Elements within lineages can be considered populations, in which each element is an individual in its genomic environment. In this way, it would be reasonable to apply microevolutionary analyses to understand transposable element (TE) evolution, such as those used to study the genetic structure of natural populations. Here, we applied a Bayesian method to infer genetic structure of populations together with classical phylogenetic and dating tools to analyze LTR-RT evolution using the monocot Setaria italica as a model species. In contrast to a phylogeny, the Bayesian clusterization method identifies populations by assigning individuals to one or more clusters according to the most probabilistic scenario of admixture, based on genetic diversity patterns. In this work, each LTR-RT insertion was considered to be one individual and each LTR-RT lineage was considered to be a single species. Nine evolutionary lineages of LTR-RTs were identified in the S. italica genome that had different genetic structures with variable numbers of clusters and levels of admixture. Comprehensive analysis of the phylogenetic, clusterization and time of insertion data allowed us to hypothesize that admixed elements represent sequences that harbor ancestral polymorphic sequence signatures. In conclusion, application of microevolutionary concepts in genome evolution studies is suitable as a complementary approach to phylogenetic analyses to address the evolutionary history and functional features of TEs.

Comparative analysis of miniature inverted-repeat transposable elements (MITEs) and long terminal repeat (LTR) retrotransposons in six Citrus species

BACKGROUND

Miniature inverted-repeat transposable elements (MITEs) and long terminal repeat (LTR) retrotransposons are ubiquitous in plants genomes, and highly important in their evolution and diversity. However, their mechanisms of insertion/amplification and roles in Citrus genome's evolution/diversity are still poorly understood.

RESULTS

To address this knowledge gap, we developed different computational pipelines to analyze, annotate and classify MITEs and LTR retrotransposons in six different sequenced Citrus species. We identified 62,010 full-length MITEs from 110 distinguished families. We observed MITEs tend to insert in gene related regions and enriched in promoters. We found that DTM63 is possibly an active Mutator-like MITE family in the traceable past and may still be active in Citrus. The insertion of MITEs resulted in massive polymorphisms and played an important role in Citrus genome diversity and gene structure variations. In addition, 6630 complete LTR retrotransposons and 13,371 solo-LTRs were identified. Among them, 12 LTR lineages separated before the differentiation of mono- and dicotyledonous plants. We observed insertion and deletion of LTR retrotransposons was accomplished with a dynamic balance, and their half-life in Citrus was ~ 1.8 million years.

CONCLUSIONS

These findings provide insights into MITEs and LTR retrotransposons and their roles in genome diversity in different Citrus genomes.

The control of endopolygalacturonase expression by the sugarcane RAV transcription factor during aerenchyma formation

The development of lysigenous aerenchyma starts with cell expansion and degradation of pectin from the middle lamella, leading to cell wall modification, and culminating with cell separation. Here we report that nutritional starvation of sugarcane induced gene expression along sections of the first 5 cm of the root and between treatments. We selected two candidate genes: a RAV transcription factor, from the ethylene response factors superfamily, and an endopolygalacturonase (EPG), a glycosyl hydrolase related to homogalacturonan hydrolysis from the middle lamella. epg1 and rav1 transcriptional patterns suggest they are essential genes at the initial steps of pectin degradation during aerenchyma development in sugarcane. Due to the high complexity of the sugarcane genome, rav1 and epg1 were sequenced from 17 bacterial artificial chromosome clones containing hom(e)ologous genomic regions, and the sequences were compared with those of Sorghum bicolor. We used one hom(e)olog sequence from each gene for transactivation assays in tobacco. rav1 was shown to bind to the epg1 promoter, repressing β-glucuronidase activity. RAV repression upon epg1 transcription is the first reported link between ethylene regulation and pectin hydrolysis during aerenchyma formation. Our findings may help to elucidate cell wall degradation in sugarcane and therefore contribute to second-generation bioethanol production.

Comparative Analysis of Homologous Sequences of Saccharum officinarum and Saccharum spontaneum Reveals Independent Polyploidization Events

Sugarcane (Saccharum spp. hybrids) is an economically important crop widely grown in tropical and subtropical regions for sugar and ethanol production. However, the large genome size, high ploidy level, interspecific hybridization and aneuploidy make sugarcane one of the most complex genomes and have long hampered genome research in sugarcane. Modern sugarcane cultivars are derived from interspecific hybridization between S. officinarum and S. spontaneum with 80-90% of the genome from S. officinarum and 10-20% of the genome from S. spontaneum. We constructed bacterial artificial chromosome (BAC) libraries of S. officinarum variety LA Purple (2n = 8x = 80) and S. spontaneum haploid clone AP85-441 (2n = 4x = 32), and selected and sequenced 97 BAC clones from the two Saccharum BAC libraries. A total of 5,847,280 bp sequence from S. officinarum and 5,011,570 bp from S. spontaneum were assembled and 749 gene models were annotated in these BACs. A relatively higher gene density and lower repeat content were observed in S. spontaneum BACs than in S. officinarum BACs. Comparative analysis of syntenic regions revealed a high degree of collinearity in genic regions between Saccharum and Sorghum bicolor and between S. officinarum and S. spontaneum. In the syntenic regions, S. spontaneum showed expansion relative to S. officinarum, and both S. officinarum and S. spontaneum showed expansion relative to sorghum. Among the 75 full-length LTR retrotransposons identified in the Saccharum BACs, none of them are older than 2.6 mys and no full-length LTR elements are shared between S. officinarum and S. spontaneum. In addition, divergence time estimated using a LTR junction marker and a syntenic gene shared by 3 S. officinarum and 1 S. spontaneum BACs revealed that the S. spontaneum intergenic region was distant to those from the 3 homologous regions in S. officinarum. Our results suggested that S. officinarum and S. spontaneum experienced at least two rounds of independent polyploidization in each lineage after their divergence from a common ancestor.

Next generation sequencing reveals packaging of host RNAs by brome mosaic virus

Although RNA viruses evolved the mechanisms of specific encapsidation, miss-packaging of cellular RNAs has been reported in such RNA virus systems as flock house virus or cucumber necrosis virus. To find out if brome mosaic virus (BMV), a tripartite RNA virus, can package cellular RNAs, BMV was propagated in barley and in Nicotiana benthamiana hosts, purified by cesium chloride (CsCl) gradient ultracentrifugation followed by nuclease treatment to remove any contaminating cellular (host) RNAs. The extracted virion RNA was then sequenced by using next-generation sequencing (NGS RNA-Seq) with the Illumina protocol. Bioinformatic analysis revealed the content of host RNAs ranging from 0.07% for BMV extracted from barley to 0.10% for the virus extracted from N. benthamiana. The viruses from two sources appeared to co-encapsidate different patterns of host-RNAs, including ribosomal RNAs (rRNAs), messenger RNAs (mRNAs) but also mitochondrial and plastid RNAs and, interestingly, transposable elements, both transposons and retrotransposons. Our data reveal that BMV virions can carry host RNAs, having a potential to mediate horizontal gene transfer (HGT) in plants.

Genomic abundance and transcriptional activity of diverse gypsy and copia long terminal repeat retrotransposons in three wild sunflower species

BACKGROUND

Long terminal repeat (LTR) retrotransposons are highly abundant in plant genomes and require transcriptional activity for their proliferative mode of replication. These sequences exist in plant genomes as diverse sublineages within the main element superfamilies (i.e., gypsy and copia). While transcriptional activity of these elements is increasingly recognized as a regular attribute of plant transcriptomes, it is currently unknown the extent to which different sublineages of these elements are transcriptionally active both within and across species. In the current report, we utilize next generation sequencing methods to examine genomic copy number abundance of diverse LTR retrotransposon sublineages and their corresponding levels of transcriptional activity in three diploid wild sunflower species, Helianthus agrestis, H. carnosus and H. porteri.

RESULTS

The diploid sunflower species under investigation differ in genome size 2.75-fold, with 2C values of 22.93 for H. agrestis, 12.31 for H. carnosus and 8.33 for H. porteri. The same diverse gypsy and copia sublineages of LTR retrotransposons were identified across species, but with gypsy sequences consistently more abundant than copia and with global gypsy sequence abundance positively correlated with nuclear genome size. Transcriptional activity was detected for multiple copia and gypsy sequences, with significantly higher activity levels detected for copia versus gypsy. Interestingly, of 11 elements identified as transcriptionally active, 5 exhibited detectable expression in all three species and 3 exhibited detectable expression in two species.

CONCLUSIONS

Combined analyses of LTR retrotransposon genomic abundance and transcriptional activity across three sunflower species provides novel insights into genome size evolution and transposable element dynamics in this group. Despite considerable variation in nuclear genome size among species, relatively conserved patterns of LTR retrotransposon transcriptional activity were observed, with a highly overlapping set of copia and gypsy sequences observed to be transcriptionally active across species. A higher proportion of copia versus gypsy elements were found to be transcriptionally active and these sequences also were expressed at higher levels.

Diversity, distribution and dynamics of full-length Copia and Gypsy LTR retroelements in Solanum lycopersicum

Transposable elements are the most abundant components of plant genomes and can dramatically induce genetic changes and impact genome evolution. In the recently sequenced genome of tomato (Solanum lycopersicum), the estimated fraction of elements corresponding to retrotransposons is nearly 62%. Given that tomato is one of the most important vegetable crop cultivated and consumed worldwide, understanding retrotransposon dynamics can provide insight into its evolution and domestication processes. In this study, we performed a genome-wide in silico search of full-length LTR retroelements in the tomato nuclear genome and annotated 736 full-length Gypsy and Copia retroelements. The dispersion level across the 12 chromosomes, the diversity and tissue-specific expression of those elements were estimated. Phylogenetic analysis based on the retrotranscriptase region revealed the presence of 12 major lineages of LTR retroelements in the tomato genome. We identified 97 families, of which 77 and 20 belong to the superfamilies Copia and Gypsy, respectively. Each retroelement family was characterized according to their element size, relative frequencies and insertion time. These analyses represent a valuable resource for comparative genomics within the Solanaceae, transposon-tagging and for the design of cultivar-specific molecular markers in tomato.

Characterization, Genomic Organization, Abundance, and Chromosomal Distribution of Ty1-copia Retrotransposons in Erianthus arundinaceus

Erianthus arundinaceus is an important wild species of the genus Saccharum with many valuable traits. However, the composition and structure of its genome are largely unknown, which have hindered its utilization in sugarcane breeding and evolutionary research. Retrotransposons constitute an appreciable fraction of plant genomes and may have played a significant role in the evolution and sequence organization of genomes. In the current study, we investigate the phylogenetic diversity and genomic abundance of Ty1-copia retrotransposons for the first time and inspect their chromosomal distribution patterns in E. arundinaceus. In total, 70 Ty1-copia reverse transcriptase (RT) sequences with significant levels of heterogeneity were obtained. The phylogenetic analysis revealed these Ty1-copia retrotransposons were classified into four distinct evolutionary lineages (Tork/TAR, Tork/Angela, Retrofit/Ale, and Sire/Maximus). Dot-blot analysis showed estimated the total copy number of Ty1-copia retrotransposons to be about 4.5 × 10³ in the E. arundinaceus genome, indicating they were a significant component. Fluorescence in situ hybridization revealed that Ty1-copia retrotransposons from the four lineages had strikingly similar patterns of chromosomal enrichment, being exclusively enriched in the subterminal heterochromatic regions of most E. arundinaceus chromosomes. This is the first clear evidence of the presence of Ty1-copia retrotransposons in the subterminal heterochromatin of E. arundinaceus. Altogether, these results promote the understanding of the diversification of Ty1-copia retrotransposons and shed light on their chromosomal distribution patterns in E. arundinaceus.

sRNAs as possible regulators of retrotransposon activity in Cryptococcus gattii VGII

Background

The absence of Argonaute genes in the fungal pathogen Cryptococcus gattii R265 and other VGII strains indicates that yeasts of this genotype cannot have a functional RNAi pathway, an evolutionarily conserved gene silencing mechanism performed by small RNAs. The success of the R265 strain as a pathogen that caused the Pacific Northwest and Vancouver Island outbreaks may imply that RNAi machinery loss could be beneficial under certain circumstances during evolution. As a result, a hypermutant phenotype would be created with high rates of genome retrotransposition, for instance. This study therefore aimed to evaluate in silicio the effect of retrotransposons and their control mechanisms by small RNAs on genomic stability and synteny loss of C. gattii R265 through retrotransposons sequence comparison and orthology analysis with other 16 C. gattii genomic sequences available.

Results

Retrotransposon mining identified a higher sequence count to VGI genotype compared to VGII, VGIII, and VGIV. However, despite the lower retrotransposon number, VGII exhibited increased synteny loss and genome rearrangement events. RNA-Seq analysis indicated highly expressed retrotransposons as well as sRNA production.

Conclusions

Genome rearrangement and synteny loss may suggest a greater retrotransposon mobilization caused by RNAi pathway absence, but the effective presence of sRNAs that matches retrotransposon sequences means that an alternative retrotransposon silencing mechanism could be active in genomic integrity maintenance of C. gattii VGII strains.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3688-4) contains supplementary material, which is available to authorized users.

Genomic skimming for identification of medium/highly abundant transposable elements in Arundo donax and Arundo plinii

Transposable elements (TEs) are the most abundant genetic material for almost all eukaryotic genomes. Their effects on the host genomes range from an extensive size variation to the regulation of gene expression, altering gene function and creating new genes. Because of TEs pivotal contribute to the host genome structure and regulation, their identification and characterization provide a wealth of useful data for gaining an in-depth understanding of host genome functioning. The giant reed (Arundo donax) is a perennial rhizomatous C₃ grass, octadecaploid, with an estimated nuclear genome size of 2744 Mbp. It is a promising feedstock for second-generation biofuels and biomethane production. To identify and characterize the most repetitive TEs in the genomes of A. donax and its ancestral A. plinii species, we carried out low-coverage whole genome shotgun sequencing for both species. Using a de novo repeat identification approach, 33,041 and 28,237 non-redundant repetitive sequences were identified and characterized in A. donax and A. plinii genomes, representing 37.55 and 31.68% of each genome, respectively. Comparative phylogenetic analyses, including the major TE classes identified in A. donax and A. plinii, together with rice and maize TE paralogs, were carried out to understand the evolutionary relationship of the most abundant TE classes. Highly conserved copies of RIRE1-like Ty1-Copia elements were discovered in two Arundo spp. in which they represented nearly 3% of each genomic sequence. We identified and characterized the medium/highly repetitive TEs in two unexplored polyploid genomes, thus generating useful information for the study of the genomic structure, composition, and functioning of these two non-model species. We provided a valuable resource that could be exploited in any effort aimed at sequencing and assembling these two genomes.

Virus-like attachment sites as structural landmarks of plants retrotransposons

BACKGROUND

The genomic data available nowadays has enabled the study of repetitive sequences and their relationship to viruses. Among them, long terminal repeat retrotransposons (LTR-RTs) are the largest component of most plant genomes, the Gypsy and Copia superfamilies being the most common. Recently it has been found that Del lineage, an LTR-RT of Gypsy superfamily, has putative virus-like attachment (vl-att) sites. This signature, originally described for retroviruses, is recognized by retroviral integrase conferring specificity to the integration process.

RESULTS

Here we retrieved 26,092 putative complete LTR-RTs from 10 lineages found in 10 fully sequenced angiosperm genomes and found putative vl-att sites that are a conserved structural landmark across these genomes. Furthermore, we reveal that each plant genome has a distinguishable LTR-RT lineage amplification pattern that could be related to the vl-att sites diversity. We used these patterns to generate a specific quick-response (QR) code for each genome that could be used as a barcode of identification of plants in the future.

CONCLUSIONS

The universal distribution of vl-att sites represents a new structural feature common to plant LTR-RTs and retroviruses. This is an important finding that expands the information about the structural similarity between LTR-RT and retroviruses. We speculate that the sequence diversity of vl-att sites could be important for the life cycle of retrotransposons, as it was shown for retroviruses. All the structural vl-att site signatures are strong candidates for further functional studies. Moreover, this is the first identification of specific LTR-RT content and their amplification patterns in a large dataset of LTR-RT lineages and angiosperm genomes. These distribution patterns could be used in the future with biotechnological identification purposes.

The impact of the Tekay chromoviral elements on genome organisation and evolution of Anemone s.l. (Ranunculaceae)

We studied the highly abundant chromoviral Tekay clade in species from three sister genera - Anemone, Pulsatilla and Hepatica (Ranunculaceae). With this clade, we performed a concomitant survey of its phylogenetic diversity, chromosomal organisation and transcriptional activity in Anemone s.l. in order to investigate dynamics of the Tekay elements at a finer scale than previously achieved in this or any other flowering clade. The phylogenetic tree built from Tekay sequences conformed to expected evolutionary relationships of the species; exceptions being A. nemorosa and A. sylvestris, which appeared more closely related that expected, and we invoke hybridisation events to explain the observed topology. The separation of elements into six clusters could be explained by episodic bursts of activity since divergence from a common ancestor at different points in their respective evolutionary histories. In Anemone s.l. the Tekay elements do not have a preferential position on chromosomes, i.e. they can have a: (i) centromeric/pericentromeric position; (ii) interstitial position in DAPI-positive AT-rich heterochromatic regions; can be (iii) dispersed throughout chromosomes; or even (iv) be absent from large heterochromatic blocks. Widespread transcriptional activity of the Tekay elements in Anemone s.l. taxa indicate that some copies of Tekay elements could still be active in this plant group, contributing to genome evolution and speciation within Anemone s.l. Identification of Tekay elements in Anemone s.l. provides valuable information for understanding how different localisation patterns might help to facilitate plant genome organisation in a structural and functional manner.

A First Glimpse of Wild Lupin Karyotype Variation As Revealed by Comparative Cytogenetic Mapping

Insight into plant genomes at the cytomolecular level provides useful information about their karyotype structure, enabling inferences about taxonomic relationships and evolutionary origins. The Old World lupins (OWL) demonstrate a high level of genomic diversification involving variation in chromosome numbers (2n = 32-52), basic chromosome numbers (x = 5-7, 9, 13) and in nuclear genome size (2C DNA = 0.97-2.68 pg). Lupins comprise both crop and wild species and provide an intriguing system to study karyotype evolution. In order to investigate lupin chromosome structure, heterologous FISH was used. Sixteen BACs that had been generated as chromosome markers for the reference species, Lupinus angustifolius, were used to identify chromosomes in the wild species and explore karyotype variation. While all "single-locus" in L. angustifolius, in the wild lupins these clones proved to be "single-locus," "single-locus" with additional signals, "repetitive" or had no detectable BAC-FISH signal. The diverse distribution of the clones in the targeted genomes suggests a complex evolution history, which possibly involved multiple chromosomal changes such as fusions/fissions and repetitive sequence amplification. Twelve BACs were sequenced and we found numerous transposable elements including DNA transposons as well as LTR and non-LTR retrotransposons with varying quantity and composition among the different lupin species. However, at this preliminary stage, no correlation was observed between the pattern of BAC-FISH signals and the repeat content in particular BACs. Here, we describe the first BAC-based chromosome-specific markers for the wild species: L. cosentinii, L. cryptanthus, L. pilosus, L. micranthus and one New World lupin, L. multiflorus. These BACs could constitute the basis for an assignment of the chromosomal and genetic maps of other lupins, e.g., L. albus and L. luteus. Moreover, we identified karyotype variation that helps illustrate the relationships between the lupins and the extensive cytological diversity within this group. In this study we premise that lupin genomes underwent at least two rounds of fusion and fission events resulting in the reduction in chromosome number from 2n = 52 through 2n = 40 to 2n = 32, followed by chromosome number increment to 2n = 42.

Genome-wide Annotation and Comparative Analysis of Long Terminal Repeat Retrotransposons between Pear Species of P. bretschneideri and P. Communis

Recent sequencing of the Oriental pear (P. bretschneideri Rehd.) genome and the availability of the draft genome sequence of Occidental pear (P. communis L.), has provided a good opportunity to characterize the abundance, distribution, timing, and evolution of long terminal repeat retrotransposons (LTR-RTs) in these two important fruit plants. Here, a total of 7247 LTR-RTs, which can be classified into 148 families, have been identified in the assembled Oriental pear genome. Unlike in other plant genomes, approximately 90% of these elements were found to be randomly distributed along the pear chromosomes. Further analysis revealed that the amplification timeframe of elements varies dramatically in different families, super-families and lineages, and the Copia-like elements have highest activity in the recent 0.5 million years (Mys). The data also showed that two genomes evolved with similar evolutionary rates after their split from the common ancestor ~0.77–1.66 million years ago (Mya). Overall, the data provided here will be a valuable resource for further investigating the impact of transposable elements on gene structure, expression, and epigenetic modification in the pear genomes.

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.

Encapsidation of Host RNAs by Cucumber Necrosis Virus Coat Protein during both Agroinfiltration and Infection

Next-generation sequence analysis of virus-like particles (VLPs) produced during agroinfiltration of cucumber necrosis virus (CNV) coat protein (CP) and of authentic CNV virions was conducted to assess if host RNAs can be encapsidated by CNV CP. VLPs containing host RNAs were found to be produced during agroinfiltration, accumulating to approximately 1/60 the level that CNV virions accumulated during infection. VLPs contained a variety of host RNA species, including the major rRNAs as well as cytoplasmic, chloroplast, and mitochondrial mRNAs. The most predominant host RNA species encapsidated in VLPs were chloroplast encoded, consistent with the efficient targeting of CNV CP to chloroplasts during agroinfiltration. Interestingly, droplet digital PCR analysis showed that the CNV CP mRNA expressed during agroinfiltration was the most efficiently encapsidated mRNA, suggesting that the CNV CP open reading frame may contain a high-affinity site or sites for CP binding and thus contribute to the specificity of CNV RNA encapsidation. Approximately 0.09% to 0.7% of the RNA derived from authentic CNV virions contained host RNA, with chloroplast RNA again being the most prominent species. This is consistent with our previous finding that a small proportion of CNV CP enters chloroplasts during the infection process and highlights the possibility that chloroplast targeting is a significant aspect of CNV infection. Remarkably, 6 to 8 of the top 10 most efficiently encapsidated nucleus-encoded RNAs in CNV virions correspond to retrotransposon or retrotransposon-like RNA sequences. Thus, CNV could potentially serve as a vehicle for horizontal transmission of retrotransposons to new hosts and thereby significantly influence genome evolution.

IMPORTANCE

Viruses predominantly encapsidate their own virus-related RNA species due to the possession of specific sequences and/or structures on viral RNA which serve as high-affinity binding sites for the coat protein. In this study, we show, using next-generation sequence analysis, that CNV also encapsidates host RNA species, which account for ∼0.1% of the RNA packaged in CNV particles. The encapsidated host RNAs predominantly include chloroplast RNAs, reinforcing previous observations that CNV CP enters chloroplasts during infection. Remarkably, the most abundantly encapsidated cytoplasmic mRNAs consisted of retrotransposon-like RNA sequences, similar to findings recently reported for flock house virus (A. Routh, T. Domitrovic, and J. E. Johnson, Proc Natl Acad Sci U S A 109:1907-1912, 2012). Encapsidation of retrotransposon sequences may contribute to their horizontal transmission should CNV virions carrying retrotransposons infect a new host. Such an event could lead to large-scale genomic changes in a naive plant host, thus facilitating host evolutionary novelty.

Using quantitative PCR with retrotransposon-based insertion polymorphisms as markers in sugarcane

Sugarcane is the main source of the world's sugar and is becoming increasingly important as a source of biofuel. The highly polyploid and heterozygous nature of the sugarcane genome has meant that characterization of the genome has lagged behind that of other important crops. Here we developed a method using a combination of quantitative PCR with a transposable marker system to score the relative number of alleles with a transposable element (TE) present at a particular locus. We screened two genera closely related to Saccharum (Miscanthus and Erianthus), wild Saccharum, traditional cultivars, and 127 modern cultivars from Brazilian and Australian breeding programmes. We showed how this method could be used in various ways. First, we showed that the method could be extended to be used as part of a genotyping system. Secondly, the history of insertion and timing of the three TEs examined supports our current understanding of the evolution of the Saccharum complex. Thirdly, all three TEs were found in only one of the two main lineages leading to the modern sugarcane cultivars and are therefore the first TEs identified that could potentially be used as markers for Saccharum spontaneum.

Building the sugarcane genome for biotechnology and identifying evolutionary trends

Background

Sugarcane is the source of sugar in all tropical and subtropical countries and is becoming increasingly important for bio-based fuels. However, its large (10 Gb), polyploid, complex genome has hindered genome based breeding efforts. Here we release the largest and most diverse set of sugarcane genome sequences to date, as part of an on-going initiative to provide a sugarcane genomic information resource, with the ultimate goal of producing a gold standard genome.

Results

Three hundred and seventeen chiefly euchromatic BACs were sequenced. A reference set of one thousand four hundred manually-annotated protein-coding genes was generated. A small RNA collection and a RNA-seq library were used to explore expression patterns and the sRNA landscape. In the sucrose and starch metabolism pathway, 16 non-redundant enzyme-encoding genes were identified. One of the sucrose pathway genes, sucrose-6-phosphate phosphohydrolase, is duplicated in sugarcane and sorghum, but not in rice and maize. A diversity analysis of the s6pp duplication region revealed haplotype-structured sequence composition. Examination of hom(e)ologous loci indicate both sequence structural and sRNA landscape variation. A synteny analysis shows that the sugarcane genome has expanded relative to the sorghum genome, largely due to the presence of transposable elements and uncharacterized intergenic and intronic sequences.

Conclusion

This release of sugarcane genomic sequences will advance our understanding of sugarcane genetics and contribute to the development of molecular tools for breeding purposes and gene discovery.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-540) contains supplementary material, which is available to authorized users.

Virus-like attachment sites and plastic CpG islands:landmarks of diversity in plant Del retrotransposons

Full-length Del elements from ten angiosperm genomes, 5 monocot and 5 dicot, were retrieved and putative attachment (att) sites were identified. In the 2432 Del elements, two types of U5 att sites and a single conserved type of U3 att site were identified. Retroviral att sites confer specificity to the integration process, different att sites types therefore implies lineage specificity. While some features are common to all Del elements, CpG island patterns within the LTRs were particular to lineage specific clusters. All eudicot copies grouped into one single clade while the monocots harbour a more diverse collection of elements. Furthermore, full-length Del elements and truncated copies were unevenly distributed amongst chromosomes. Elements of Del lineage are organized in plants into three clusters and each cluster is composed of elements with distinct LTR features. Our results suggest that the Del lineage efficiently amplified in the monocots and that one branch is probably a newly emerging sub-lineage. Finally, sequences in all groups are under purifying selection. These results show the LTR region is dynamic and important in the evolution of LTR-retrotransposons, we speculate that it is a trigger for retrotransposon diversification.

De novo assembly and transcriptome analysis of contrasting sugarcane varieties

Sugarcane is an important crop and a major source of sugar and alcohol. In this study, we performed de novo assembly and transcriptome annotation for six sugarcane genotypes involved in bi-parental crosses. The de novo assembly of the sugarcane transcriptome was performed using short reads generated using the Illumina RNA-Seq platform. We produced more than 400 million reads, which were assembled into 72,269 unigenes. Based on a similarity search, the unigenes showed significant similarity to more than 28,788 sorghum proteins, including a set of 5,272 unigenes that are not present in the public sugarcane EST databases; many of these unigenes are likely putative undescribed sugarcane genes. From this collection of unigenes, a large number of molecular markers were identified, including 5,106 simple sequence repeats (SSRs) and 708,125 single-nucleotide polymorphisms (SNPs). This new dataset will be a useful resource for future genetic and genomic studies in this species.

Differential detection of transposable elements between Saccharum species

Cultivars of sugarcane (Saccharum) are hybrids between species S. officinarum (x = 10, 2n = 8x = 80) and S. spontaneum (x = 8, 2n = 5 – 16x = 40 – 128). These accessions have 100 to 130 chromosomes, 80–85% of which are derived from S. officinarum, 10–15% from S. spontaneum, and 5–10% are possible recombinants between the two genomes. The aim of this study was to analyze the repetition of DNA sequences in S. officinarum and S. spontaneum. For this purpose, genomic DNA from S. officinarum was digested with restriction enzymes and the fragments cloned. Sixty-eight fragments, approximately 500 bp, were cloned, sequenced and had their identity analyzed in NCBI, and in the rice, maize, and sorghum genome databases using BLAST. Twelve clones containing partial transposable elements, one single-copy control, one DNA repetitive clone control and two genome controls were analyzed by DNA hybridization on membrane, using genomic probes from S. officinarum and S. spontaneum. The hybridization experiment revealed that six TEs had a similar repetitive DNA pattern in the genomes of S. officinarum and S. spontaneum, while six TEs were more abundant in the genome of S. officinarum. We concluded that the species S. officinarum and S. spontaneum have differential accumulation LTR retrotransposon families, suggesting distinct insertion or modification patterns.

Isolation and characterization of novel Ty1-copia-like retrotransposons from lily

Species of the genus Lilium are well known for their large genomes. Although expansion of noncoding repeated DNA is believed to account for this genome size, retroelement del Ty3-gypsy is the only one described so far in the genus Lilium. We isolated Ty1-copia elements from Lilium longiflorum and named them LIREs (lily retrotransposons). The long terminal repeats, primer binding site, and polypurine tract sequences are highly similar among the LIRE elements, indicating that they are in the same lineage. Although the protein-coding regions were highly decayed, the sequence motifs of the integrase, reverse transcriptase, and RNase H domains were identifiable as belonging to the order of Ty1-copia elements. Phylogenetic analysis and primer binding site sequences revealed that these elements belonged to the Ale lineage among the six lineages of plant Ty1-copia elements. Base substitutions in the long terminal repeats estimated that the integration times of the LIRE Ty1-copia elements were between 0.7 and 5.5 mya. In situ hybridization showed that the LIRE elements were present in all the chromosomes of L. longiflorum and L. lancifolium, but absent in centromeres, telomeres, and 45S rRNA sites in both species. The LIRE elements were present very abundantly in species of the genus Lilium, but absent in other genera of the family Liliaceae, implying that the LIRE elements might have contributed to the expansion of the genome in the genus Lilium.

Characterization of the LTR retrotransposon repertoire of a plant clade of six diploid and one tetraploid species

Comparisons of closely related species are needed to understand the fine-scale dynamics of retrotransposon evolution in flowering plants. Towards this goal, we classified the long terminal repeat (LTR) retrotransposons from six diploid and one tetraploid species of Orobanchaceae. The study species are the autotrophic, non-parasitic Lindenbergia philippensis (as an out-group) and six closely related holoparasitic species of Orobanche [O. crenata, O. cumana, O. gracilis (tetraploid) and O. pancicii] and Phelipanche (P. lavandulacea and P. ramosa). All major plant LTR retrotransposon clades could be identified, and appear to be inherited from a common ancestor. Species of Orobanche, but not Phelipanche, are enriched in Ty3/Gypsy retrotransposons due to a diversification of elements, especially chromoviruses. This is particularly striking in O. gracilis, where tetraploidization seems to have contributed to the Ty3/Gypsy enrichment and led to the emergence of seven large species-specific families of chromoviruses. The preferential insertion of chromoviruses in heterochromatin via their chromodomains might have favored their diversification and enrichment. Our phylogenetic analyses of LTR retrotransposons from Orobanchaceae also revealed that the Bianca clade of Ty1/Copia and the SMART-related elements are much more widely distributed among angiosperms than previously known.

Biofuel and energy crops: high-yield Saccharinae take center stage in the post-genomics era

The Saccharinae, especially sugarcane, Miscanthus and sorghum, present remarkable characteristics for bioenergy production. Biotechnology of these plants will be important for a sustainable feedstock supply. Herein, we review knowledge useful for their improvement and synergies gained by their parallel study.

Global analysis of the sugarcane microtranscriptome reveals a unique composition of small RNAs associated with axillary bud outgrowth

Axillary bud outgrowth determines shoot architecture and is under the control of endogenous hormones and a fine-tuned gene-expression network, which probably includes small RNAs (sRNAs). Although it is well known that sRNAs act broadly in plant development, our understanding about their roles in vegetative bud outgrowth remains limited. Moreover, the expression profiles of microRNAs (miRNAs) and their targets within axillary buds are largely unknown. Here, we employed sRNA next-generation sequencing as well as computational and gene-expression analysis to identify and quantify sRNAs and their targets in vegetative axillary buds of the biofuel crop sugarcane (Saccharum spp.). Computational analysis allowed the identification of 26 conserved miRNA families and two putative novel miRNAs, as well as a number of trans-acting small interfering RNAs. sRNAs associated with transposable elements and protein-encoding genes were similarly represented in both inactive and developing bud libraries. Conversely, sequencing and quantitative reverse transcription-PCR results revealed that specific miRNAs were differentially expressed in developing buds, and some correlated negatively with the expression of their targets at specific stages of axillary bud development. For instance, the expression patterns of miR159 and its target GAMYB suggested that they may play roles in regulating abscisic acid-signalling pathways during sugarcane bud outgrowth. Our work reveals, for the first time, differences in the composition and expression profiles of diverse sRNAs and targets between inactive and developing vegetative buds that, together with the endogenous balance of specific hormones, may be important in regulating axillary bud outgrowth.

[Identification and analysis methods of plant LTR retrotransposon sequences]

LTR retrotransposons are an important class of eukaryotic transposable elements, which are ubiquitous and highly heterogeneous in plant and play a major role in genome evolution of eukaryote. They are now extensively employed in gene function and genetic diversity analyses. Identification of LTR retrotransposons is the precondition for its application. Therefore, it has important theoretical significance and practical application value in studying identification and analysis methods LTR retrotransposon sequences. Bioinformatic software of the sequence analysis, according to the work principle, can be classified roughly into two types: sequence alignment and sequence identification of conserved domains. Alignment software, such as BLAST and DNAstar, produce the corresponding sequence information through comparison of sequence similarity; however, this kind of software cannot be applied for full length sequences. According to the principle, LTR retro-transposon sequence identification software can be roughly sorted into four types: de novo repeat discovery method, com-parative genomic method, homology-based method, and structure-based method. For example, LTR_Finder based on de novo repeat discovery method can accurately predict and annotate LTR retrotransposons for full length sequences; Repeat-Masker, which is based on homology-based method, can discover LTR retrotransposons by comparing the similarity with known sequences in the database. In this article, different methods of identification and analysis of retrotransposon se-quences were compared and analyzed, and a set of flow of LTR retrotransposons sequence analysis was summarized in order to provide the reference for LTR retrotransposons sequence analysis.

A novel linkage map of sugarcane with evidence for clustering of retrotransposon-based markers

BACKGROUND

The development of sugarcane as a sustainable crop has unlimited applications. The crop is one of the most economically viable for renewable energy production, and CO2 balance. Linkage maps are valuable tools for understanding genetic and genomic organization, particularly in sugarcane due to its complex polyploid genome of multispecific origins. The overall objective of our study was to construct a novel sugarcane linkage map, compiling AFLP and EST-SSR markers, and to generate data on the distribution of markers anchored to sequences of scIvana_1, a complete sugarcane transposable element, and member of the Copia superfamily.

RESULTS

The mapping population parents ('IAC66-6' and 'TUC71-7') contributed equally to polymorphisms, independent of marker type, and generated markers that were distributed into nearly the same number of co-segregation groups (or CGs). Bi-parentally inherited alleles provided the integration of 19 CGs. The marker number per CG ranged from two to 39. The total map length was 4,843.19 cM, with a marker density of 8.87 cM. Markers were assembled into 92 CGs that ranged in length from 1.14 to 404.72 cM, with an estimated average length of 52.64 cM. The greatest distance between two adjacent markers was 48.25 cM. The scIvana_1-based markers (56) were positioned on 21 CGs, but were not regularly distributed. Interestingly, the distance between adjacent scIvana_1-based markers was less than 5 cM, and was observed on five CGs, suggesting a clustered organization.

CONCLUSIONS

Results indicated the use of a NBS-profiling technique was efficient to develop retrotransposon-based markers in sugarcane. The simultaneous maximum-likelihood estimates of linkage and linkage phase based strategies confirmed the suitability of its approach to estimate linkage, and construct the linkage map. Interestingly, using our genetic data it was possible to calculate the number of retrotransposon scIvana_1 (~60) copies in the sugarcane genome, confirming previously reported molecular results. In addition, this research possibly will have indirect implications in crop economics e.g., productivity enhancement via QTL studies, as the mapping population parents differ in response to an important fungal disease.