Genome-wide analysis of long terminal repeat retrotransposons from the cranberry Vaccinium macrocarpon

BACKGROUND: Long terminal repeat (LTR) retrotransposons are widespread in plant genomes and play a large role in the generation of genomic variation. Despite this, their identification and characterization remains challenging, especially for non-model genomes. Hence, LTR retrotransposons remain undercharacterized in Vaccinium genomes, although they may be beneficial for current berry breeding efforts. OBJECTIVE: Exemplarily focusing on the genome of American cranberry (Vaccinium macrocarpon Aiton), we aim to generate an overview of the LTR retrotransposon landscape, highlighting the abundance, transcriptional activity, sequence, and structure of the major retrotransposon lineages. METHODS: Graph-based clustering of whole genome shotgun Illumina reads was performed to identify the most abundant LTR retrotransposons and to reconstruct representative in silico full-length elements. To generate insights into the LTR retrotransposon diversity in V. macrocarpon, we also queried the genome assembly for presence of reverse transcriptases (RTs), the key domain of LTR retrotransposons. Using transcriptomic data, transcriptional activity of retrotransposons corresponding to the consensuses was analyzed. RESULTS: We provide an in-depth characterization of the LTR retrotransposon landscape in the V. macrocarpon genome. Based on 475 RTs harvested from the genome assembly, we detect a high retrotransposon variety, with all major lineages present. To better understand their structural hallmarks, we reconstructed 26 Ty1-copia and 28 Ty3-gypsy in silico consensuses that capture the detected diversity. Accordingly, we frequently identify association with tandemly repeated motifs, extra open reading frames, and specialized, lineage-typical domains. Based on the overall high genomic abundance and transcriptional activity, we suggest that retrotransposons of the Ale and Athila lineages are most promising to monitor retrotransposon-derived polymorphisms across accessions. CONCLUSIONS: We conclude that LTR retrotransposons are major components of the V. macrocarpon genome. The representative consensuses provide an entry point for further Vaccinium genome analyses and may be applied to derive molecular markers for enhancing cranberry selection and breeding.

Taxus yunnanensis genome offers insights into gymnosperm phylogeny and taxol production

Taxol, a natural product derived from Taxus, is one of the most effective natural anticancer drugs and the biosynthetic pathway of Taxol is the basis of heterologous bio-production. Here, we report a high-quality genome assembly and annotation of Taxus yunnanensis based on 10.7 Gb sequences assembled into 12 chromosomes with contig N50 and scaffold N50 of 2.89 Mb and 966.80 Mb, respectively. Phylogenomic analyses show that T. yunnanensis is most closely related to Sequoiadendron giganteum among the sampled taxa, with an estimated divergence time of 133.4−213.0 MYA. As with most gymnosperms, and unlike most angiosperms, there is no evidence of a recent whole-genome duplication in T. yunnanensis. Repetitive sequences, especially long terminal repeat retrotransposons, are prevalent in the T. yunnanensis genome, contributing to its large genome size. We further integrated genomic and transcriptomic data to unveil clusters of genes involved in Taxol synthesis, located on the chromosome 12, while gene families encoding hydroxylase in the Taxol pathway exhibited significant expansion. Our study contributes to the further elucidation of gymnosperm relationships and the Taxol biosynthetic pathway.

Long Terminal Repeat Retrotransposon Afut4 promotes azole resistance of Aspergillus fumigatus by enhancing the expression of sac1 gene

Aspergillus fumigatus causes a series of invasive diseases, including the high-mortality invasive aspergillosis, and has been a serious global health threat because of its increased resistance to the first-line clinical triazoles. We analyzed the whole-genome sequence of 15 A. fumigatus strains from China and found that long terminal repeat retrotransposons (LTR-RTs), including Afut1 , Afut2, Afut3, and Afut4 , are most common and have the largest total nucleotide length among all transposable elements in A. fumigatus . Deleting one of the most enriched Afut4 977-sac1 in azole-resistant strains decreased azole resistance and downregulated its nearby gene, sac1 , but it did not significantly affect the expression of genes of the ergosterol synthesis pathway. We then discovered that 5'LTR of Afut4 977-sac1 had promoter activity and enhanced the adjacent sac1 gene expression. We found that sac1 is important to A. fumigatus , and the upregulated sac1 caused the elevated resistance of A. fumigatus to azoles. Finally, we showed that Afut4 977-sac1 has an evolution pattern similar to that of the whole genome of azole-resistant strains due to azoles; phylogenetic analysis on both the whole genome and Afut4 977-sac1 suggests that the insertion of Afut4 977-sac1 might have preceded the emergence of azole-resistant strains. Taking these data together, we found that LTR-RT Afut4 977-sac1 might be involved in the regulation of azole resistance of A. fumigatus by upregulating its nearby sac1 gene.

Transposons-Based Clonal Diversity in Trematode Involves Parts of CR1 (LINE) in Eu- and Heterochromatin

Trematode parthenitae have long been believed to form clonal populations, but clonal diversity has been discovered in this asexual stage of the lifecycle. Clonal polymorphism in the model species Himasthla elongata has been previously described, but the source of this phenomenon remains unknown. In this work, we traced cercarial clonal diversity using a simplified amplified fragment length polymorphism (SAFLP) method and characterised the nature of fragments in diverse electrophoretic bands. The repetitive elements were identified in both the primary sequence of the H. elongata genome and in the transcriptome data. Long-interspersed nuclear elements (LINEs) and long terminal repeat retrotransposons (LTRs) were found to represent an overwhelming majority of the genome and the transposon transcripts. Most sequenced fragments from SAFLP pattern contained the reverse transcriptase (RT, ORF2) domains of LINEs, and only a few sequences belonged to ORFs of LTRs and ORF1 of LINEs. A fragment corresponding to a CR1-like (LINE) spacer region was discovered and named CR1-renegade (CR1-rng). In addition to RT-containing CR1 transcripts, we found short CR1-rng transcripts in the redia transcriptome and short contigs in the mobilome. Probes against CR1-RT and CR1-rng presented strikingly different pictures in FISH mapping, despite both being fragments of CR1. In silico data and Southern blotting indicated that CR1-rng is not tandemly organised. CR1 involvement in clonal diversity is discussed.

Genome-wide analysis of long terminal repeat retrotransposons from the cranberry Vaccinium macrocarpon

BACKGROUND: Long terminal repeat (LTR) retrotransposons are widespread in plant genomes and play a large role in the generation of genomic variation. Despite this, their identification and characterization remains challenging, especially for non-model genomes. Hence, LTR retrotransposons remain undercharacterized in Vaccinium genomes, although they may be beneficial for current berry breeding efforts. OBJECTIVE: Exemplarily focusing on the genome of American cranberry (Vaccinium macrocarpon Aiton), we aim to generate an overview of the LTR retrotransposon landscape, highlighting the abundance, transcriptional activity, sequence, and structure of the major retrotransposon lineages. METHODS: Graph-based clustering of whole genome shotgun Illumina reads was performed to identify the most abundant LTR retrotransposons and to reconstruct representative in silico full-length elements. To generate insights into the LTR retrotransposon diversity in V. macrocarpon, we also queried the genome assembly for presence of reverse transcriptases (RTs), the key domain of LTR retrotransposons. Using transcriptomic data, transcriptional activity of retrotransposons corresponding to the consensuses was analyzed. RESULTS: We provide an in-depth characterization of the LTR retrotransposon landscape in the V. macrocarpon genome. Based on 475 RTs harvested from the genome assembly, we detect a high retrotransposon variety, with all major lineages present. To better understand their structural hallmarks, we reconstructed 26 Ty1-copia and 28 Ty3-gypsyin silico consensuses that capture the detected diversity. Accordingly, we frequently identify association with tandemly repeated motifs, extra open reading frames, and specialized, lineage-typical domains. Based on the overall high genomic abundance and transcriptional activity, we suggest that retrotransposons of the Ale and Athila lineages are most promising to monitor retrotransposon-derived polymorphisms across accessions. CONCLUSIONS: We conclude that LTR retrotransposons are major components of the V. macrocarpon genome. The representative consensuses provide an entry point for further Vaccinium genome analyses and may be applied to derive molecular markers for enhancing cranberry selection and breeding.

Genome-wide characterization of LTR retrotransposons in the non-model deep-sea annelid Lamellibrachia luymesi

Background Long Terminal Repeat retrotransposons (LTR retrotransposons) are mobile genetic elements composed of a few genes between terminal repeats and, in some cases, can comprise over half of a genome’s content. Available data on LTR retrotransposons have facilitated comparative studies and provided insight on genome evolution. However, data are biased to model systems and marine organisms, including annelids, have been underrepresented in transposable elements studies. Here, we focus on genome of Lamellibrachia luymesi, a vestimentiferan tubeworm from deep-sea hydrocarbon seeps, to gain knowledge of LTR retrotransposons in a deep-sea annelid. Results We characterized LTR retrotransposons present in the genome of L. luymesi using bioinformatic approaches and found that intact LTR retrotransposons makes up about 0.1% of L. luymesi genome. Previous characterization of the genome has shown that this tubeworm hosts several known LTR-retrotransposons. Here we describe and classify LTR retrotransposons in L. luymesi as within the Gypsy, Copia and Bel-pao superfamilies. Although, many elements fell within already recognized families (e.g., Mag, CSRN1), others formed clades distinct from previously recognized families within these superfamilies. However, approximately 19% (41) of recovered elements could not be classified. Gypsy elements were the most abundant while only 2 Copia and 2 Bel-pao elements were present. In addition, analysis of insertion times indicated that several LTR-retrotransposons were recently transposed into the genome of L. luymesi, these elements had identical LTR’s raising possibility of recent or ongoing retrotransposon activity. Conclusions Our analysis contributes to knowledge on diversity of LTR-retrotransposons in marine settings and also serves as an important step to assist our understanding of the potential role of retroelements in marine organisms. We find that many LTR retrotransposons, which have been inserted in the last few million years, are similar to those found in terrestrial model species. However, several new groups of LTR retrotransposons were discovered suggesting that the representation of LTR retrotransposons may be different in marine settings. Further study would improve understanding of the diversity of retrotransposons across animal groups and environments.

LTR-retrotransposon dynamics in common fig (Ficus carica L.) genome

Background Long Terminal Repeat retrotransposons (LTR-REs) are repetitive DNA sequences that constitute a large part of the genome. The improvement of sequencing technologies and sequence assembling strategies has achieved genome sequences with much greater reliability than those of the past, especially in relation to repetitive DNA sequences. Results In this study, we analysed the genome of Ficus carica L., obtained using third generation sequencing technologies and recently released, to characterise the complete complement of full-length LTR-REs to study their dynamics during fig genome evolution. A total of 1867 full-length elements were identified. Those belonging to the Gypsy superfamily were the most abundant; among these, the Chromovirus/Tekay lineage was the most represented. For the Copia superfamily, Ale was the most abundant lineage. Measuring the estimated insertion time of each element showed that, on average, Ivana and Chromovirus/Tekay were the youngest lineages of Copia and Gypsy superfamilies, respectively. Most elements were inactive in transcription, both constitutively and in leaves of plants exposed to an abiotic stress, except for some elements, mostly belonging to the Copia/Ale lineage. A relationship between the inactivity of an element and inactivity of genes lying in close proximity to it was established. Conclusions The data reported in this study provide one of the first sets of information on the genomic dynamics related to LTR-REs in a plant species with highly reliable genome sequence. Fig LTR-REs are highly heterogeneous in abundance and estimated insertion time, and only a few elements are transcriptionally active. In general, the data suggested a direct relationship between estimated insertion time and abundance of an element and an inverse relationship between insertion time (or abundance) and transcription, at least for Copia LTR-REs.

High nucleotide substitution rates associated with retrotransposon proliferation drive dynamic secretome evolution in smut pathogens

Transposable elements (TEs) play a pivotal role in shaping diversity in eukaryotic genomes. The covered smut pathogen on barley, Ustilago hordei, encountered a recent genome expansion. Using long reads, we assembled genomes of 6 U. hordei strains and 3 sister species, to study this genome expansion. We found that larger genome sizes can mainly be attributed to long terminal repeat retrotransposons (LTR-RTs) of the Copia and Gypsy superfamilies. From the studied smuts, LTR-RTs proliferated the most recently and to the furthest extent in the U. hordei genome, in which they make up for 19.5% of the genome. Interestingly, the extent of TE proliferation in different smut species is positively correlated to the mating-type locus size, which is largest in U. hordei with up to ~560 kb. TE transposition within the mating-type loci and their flanking regions are mating-type specific, which is likely due to the very low recombination activity in this region. Furthermore, LTR-RT proliferation was found to be associated with higher nucleotide substitution levels, as genes in genome regions that are rich in dynamic LTR-RTs display higher nucleotide substitution levels. The high nucleotide substitution rate particularly affected the evolution of genes encoding secreted proteins as substitutions more frequently led to amino acid alterations. The mechanism behind this increase in nucleotide substitution rate remains elusive, but seems not to be a consequence of the repeat-induced point mutation (RIP) mechanism, as genes and LTR-RTs did not display typical RIP substitutions.

Genome sequence and genetic diversity analysis of an under-domesticated orphan crop, white fonio (Digitaria exilis)

Background Digitaria exilis, white fonio, is a minor but vital crop of West Africa that is valued for its resilience in hot, dry, and low-fertility environments and for the exceptional quality of its grain for human nutrition. Its success is hindered, however, by a low degree of plant breeding and improvement. Findings We sequenced the fonio genome with long-read SMRT-cell technology, yielding a ∼761 Mb assembly in 3,329 contigs (N50, 1.73 Mb; L50, 126). The assembly approaches a high level of completion, with a BUSCO score of >99%. The fonio genome was found to be a tetraploid, with most of the genome retained as homoeologous duplications that differ overall by ∼4.3%, neglecting indels. The 2 genomes within fonio were found to have begun their independent divergence ∼3.1 million years ago. The repeat content (>49%) is fairly standard for a grass genome of this size, but the ratio of Gypsy to Copia long terminal repeat retrotransposons (∼6.7) was found to be exceptionally high. Several genes related to future improvement of the crop were identified including shattering, plant height, and grain size. Analysis of fonio population genetics, primarily in Mali, indicated that the crop has extensive genetic diversity that is largely partitioned across a north-south gradient coinciding with the Sahel and Sudan grassland domains. Conclusions We provide a high-quality assembly, annotation, and diversity analysis for a vital African crop. The availability of this information should empower future research into further domestication and improvement of fonio.

Evolutionary and genomic comparisons of hybrid uninucleate and nonhybrid Rhizoctonia fungi

The basidiomycetous fungal genus, Rhizoctonia, can cause severe damage to many plants and is composed of multinucleate, binucleate, and uninucleate species differing in pathogenicity. Here we generated chromosome-scale genome assemblies of the three nuclear types of Rhizoctonia isolates. The genomic comparisons revealed that the uninucleate JN strain likely arose by somatic hybridization of two binucleate isolates, and maintained a diploid nucleus. Homeolog gene pairs in the JN genome have experienced both decelerated or accelerated evolution. Homeolog expression dominance occurred between JN subgenomes, in which differentially expressed genes show potentially less evolutionary constraint than the genes without. Analysis of mating-type genes suggested that Rhizoctonia maintains the ancestral tetrapolarity of the Basidiomycota. Long terminal repeat-retrotransposons displayed a reciprocal correlation with the chromosomal GC content in the three chromosome-scale genomes. The more aggressive multinucleate XN strain had more genes encoding enzymes for host cell wall decomposition. These findings demonstrate some evolutionary changes of a recently derived hybrid and in multiple nuclear types of Rhizoctonia.