Long Terminal Repeat Retrotransposons Research Articles

Association analysis of repetitive elements and R-loop formation across species

BackgroundAlthough recent studies have revealed the genome-wide distribution of R-loops, our understanding of R-loop formation is still limited. Genomes are known to have a large number of repetitive elements. Emerging evidence suggests that these sequences may play an important regulatory role. However, few studies have investigated the effect of repetitive elements on R-loop formation.ResultsWe found different repetitive elements related to R-loop formation in various species. By controlling length and genomic distributions, we observed that satellite, long interspersed nuclear elements (LINEs), and DNA transposons were each specifically enriched for R-loops in humans, fruit flies, and Arabidopsis thaliana, respectively. R-loops also tended to arise in regions of low-complexity or simple repeats across species. We also found that the repetitive elements associated with R-loop formation differ according to developmental stage. For instance, LINEs and long terminal repeat retrotransposons (LTRs) are more likely to contain R-loops in embryos (fruit fly) and then turn out to be low-complexity and simple repeats in post-developmental S2 cells.ConclusionsOur results indicate that repetitive elements may have species-specific or development-specific regulatory effects on R-loop formation. This work advances our understanding of repetitive elements and R-loop biology.

Saccharomyces cerevisiae RNA lariat debranching enzyme, Dbr1p, is required for completion of reverse transcription by the retrovirus-like element Ty1 and cleaves branched Ty1 RNAs.

RNA debranching enzymes are 2'-5' phosphodiesterases found in all eukaryotes. Their main role is cleavage of intron RNA lariat branch points, promoting RNA turnover via exonucleases. Consistent with this role, cells with reduced RNA debranching enzyme activity accumulate intron RNA lariats. The Saccharomyces cerevisiae RNA debranching enzyme Dbr1p is also a host factor for the yeast long terminal repeat (LTR) retrotransposon Ty1, a model for many aspects of retroviral replication. Fittingly, the human RNA debranching enzyme Dbr1 is a host factor for the human immunodeficiency virus, HIV-1. The yeast and human RNA debranching enzymes act at the reverse transcription stages for Ty1 and HIV-1, respectively. Although efficient production of full-length Ty1 cDNA requires Dbr1p, the findings reported here indicate that production of the earliest distinct cDNA product, minus strand strong stop DNA (-sssDNA), is equivalent in wild type and dbr1∆ mutant cells. Several branched Ty1 RNAs are shown to accumulate in dbr1∆ cells during retrotransposition. These data are consistent with creation of Ty1 RNA branches prior to Ty1 reverse transcription and their removal by Dbr1p to allow efficient extension of early cDNA products. The data support the possibility that RNA branch formation and cleavage play broadly shared, but unknown roles in retroviral and LTR retrotransposon reverse transcription.

The Unique, the Known, and the Unknown of Spumaretrovirus Assembly.

Within the family of Retroviridae, foamy viruses (FVs) are unique and unconventional with respect to many aspects in their molecular biology, including assembly and release of enveloped viral particles. Both components of the minimal assembly and release machinery, Gag and Env, display significant differences in their molecular structures and functions compared to the other retroviruses. This led to the placement of FVs into a separate subfamily, the Spumaretrovirinae. Here, we describe the molecular differences in FV Gag and Env, as well as Pol, which is translated as a separate protein and not in an orthoretroviral manner as a Gag-Pol fusion protein. This feature further complicates FV assembly since a specialized Pol encapsidation strategy via a tripartite Gag-genome–Pol complex is used. We try to relate the different features and specific interaction patterns of the FV Gag, Pol, and Env proteins in order to develop a comprehensive and dynamic picture of particle assembly and release, but also other features that are indirectly affected. Since FVs are at the root of the retrovirus tree, we aim at dissecting the unique/specialized features from those shared among the Spuma- and Orthoretrovirinae. Such analyses may shed light on the evolution and characteristics of virus envelopment since related viruses within the Ortervirales, for instance LTR retrotransposons, are characterized by different levels of envelopment, thus affecting the capacity for intercellular transmission.

Evolution of Tandem Repeats Is Mirroring Post-polyploid Cladogenesis in Heliophila (Brassicaceae).

The unigeneric tribe Heliophileae encompassing more than 100 Heliophila species is morphologically the most diverse Brassicaceae lineage. The tribe is endemic to southern Africa, confined chiefly to the southwestern South Africa, home of two biodiversity hotspots (Cape Floristic Region and Succulent Karoo). The monospecific Chamira (C. circaeoides), the only crucifer species with persistent cotyledons, is traditionally retrieved as the closest relative of Heliophileae. Our transcriptome analysis revealed a whole-genome duplication (WGD) ∼26.15–29.20 million years ago, presumably preceding the Chamira/Heliophila split. The WGD was then followed by genome-wide diploidization, species radiations, and cladogenesis in Heliophila. The expanded phylogeny based on nuclear ribosomal DNA internal transcribed spacer (ITS) uncovered four major infrageneric clades (A–D) in Heliophila and corroborated the sister relationship between Chamira and Heliophila. Herein, we analyzed how the diploidization process impacted the evolution of repetitive sequences through low-coverage whole-genome sequencing of 15 Heliophila species, representing the four clades, and Chamira. Despite the firmly established infrageneric cladogenesis and different ecological life histories (four perennials vs. 11 annual species), repeatome analysis showed overall comparable evolution of genome sizes (288–484 Mb) and repeat content (25.04–38.90%) across Heliophila species and clades. Among Heliophila species, long terminal repeat (LTR) retrotransposons were the predominant components of the analyzed genomes (11.51–22.42%), whereas tandem repeats had lower abundances (1.03–12.10%). In Chamira, the tandem repeat content (17.92%, 16 diverse tandem repeats) equals the abundance of LTR retrotransposons (16.69%). Among the 108 tandem repeats identified in Heliophila, only 16 repeats were found to be shared among two or more species; no tandem repeats were shared by Chamira and Heliophila genomes. Six “relic” tandem repeats were shared between any two different Heliophila clades by a common descent. Four and six clade-specific repeats shared among clade A and C species, respectively, support the monophyly of these two clades. Three repeats shared by all clade A species corroborate the recent diversification of this clade revealed by plastome-based molecular dating. Phylogenetic analysis based on repeat sequence similarities separated the Heliophila species to three clades [A, C, and (B+D)], mirroring the post-polyploid cladogenesis in Heliophila inferred from rDNA ITS and plastome sequences.

Draft genome of the herbaceous bamboo Raddia distichophylla.

Bamboos are important nontimber forest plants widely distributed in the tropical and subtropical regions of Asia, Africa, America, and Pacific islands. They comprise the Bambusoideae in the grass family (Poaceae), including approximately 1700 described species in 127 genera. In spite of the widespread uses of bamboo for food, construction, and bioenergy, the gene repertoire of bamboo still remains largely unexplored. Raddia distichophylla (Schrad. ex Nees) Chase, belonging to the tribe Olyreae (Bambusoideae, Poaceae), is a diploid herbaceous bamboo with only slightly lignified stems. In this study, we report a draft genome assembly of the ∼589 Mb whole-genome sequence of R. distichophylla with a contig N50 length of 86.36 Kb. Repeat sequences account for ∼49.08% of the genome assembly, of which LTR retrotransposons occupy ∼35.99% of the whole genome. A total of 30,763 protein-coding genes were annotated in the R. distichophylla genome with an average transcript size of 2887 bp. Access to this herbaceous bamboo genome sequence will provide novel insights into biochemistry, molecular marker-assisted breeding programs, and germplasm conservation for bamboo species worldwide.

Genome-wide characteristics and evolution analysis of long terminal repeat retrotransposons in Phyllostachys edulis

Objective This study aims to analyze the characteristics of long terminal repeat retrotransposons (LTR-REs) in moso bamboo genome, so as to promote the research on the function of LTR-REs in moso bamboo genome and the genetic diversity of bamboo resources. Method Based on bioinformatics methods, LTR retrotransposons in the second edition of moso bamboo genome were annotated and classified by LTRharvest and RepeatMakser software, and the distribution characteristics, evolution characteristics and insertion time of the obtained LTR retrotransposons were analyzed. Result A total of 1 014 565 LTR retrotransposons and 1 562 families were identified, accounting for 54.97% of moso bamboo genome. Among them, the ratio of solo LTR retrotransposons to intact LTR retrotransposons (S/F) was relatively high (about 1.77∶1.00), indicating that a higher frequency of illegitimate recombination and unbalanced recombination might have occurred in the LTR-REs of moso bamboo genome. LTR retrotransposons were divided into Ty1-copia and Ty3-gypsy superfamilies, and ten lineages included Tork, Reftrofit, Sire, Oryco, Del, Reina, Crm, Tat, Galadriel, and Athila. The preference of Ty1-copia and Ty3-gypsy superfamiles for PBS sites showed an opposite tendency. The longer LTR retrotransposons had longer LTR sequences and more complete structures. The insertion time of LTR retrotransposon in moso bamboo was mainly concentrated in the 0−2.0 Ma region, and it was still in a state of slow growth. Conclusion The high-quality assembly of the second edition of moso bamboo genome can better annotate and analyze the LTR retrotransposons in moso bamboo genome. The LTR harvest method based on structure prediction can more accurately predict the LTR retrotransposons of moso bamboo. The LTR retrotransposons of different lineages have different differentiation and amplification activities during evolution. LTR retrotransposons are generally in a state of continuous amplification, which is one of the reasons for the large genome of moso bamboo. [Ch, 3 fig. 3 tab. 52 ref.]

High-Throughput Profiling of Extrachromosomal Linear DNAs of Long Terminal Repeat Retrotransposons by ALE-seq.

Extrachromosomal linear DNA (eclDNA) is the reverse-transcribed cDNA intermediate derived from long terminal repeat (LTR) transposable elements (TEs) (Cho et al., Nat Plants 5:26-33, 2018). Given that the eclDNAs are the final intermediate of LTR-TE life cycle prior to integration to the host chromosomes, their presence is considered a strong indication of active LTR retrotransposons (Cho et al., Nat Plants 5:26-33, 2018; Lanciano et al., PLoS Genet 13:e1006630, 2017). Here, we describe a method of amplification of LTR extrachromosomal DNA followed by sequencing (ALE-seq) which determines the 5' LTR sequences of eclDNAs. Briefly, ALE-seq consists of two steps of amplification, in vitro transcription of adaptor-ligated eclDNAs and subsequent reverse transcription to cDNAs primed at the conserved primer binding site (PBS) (Cho et al., Nat Plants 5:26-33, 2018). ALE-seq allows the high-throughput identification of novel LTR-TEs which are active in plants that could be potentially useful for crop biotechnology.

Bioinformatics Analysis Guides to LTR Retrotransposon-Derived Extrachromosomal Linear DNAs Identified by ALE-seq.

ALE-seq is a method devised to identify pre-integration intermediates of LTR retrotransposons called extrachromosomal linear DNA, which can be used to predict retrotransposition activity. We describe here a bioinformatic methodology to process reads obtained from the ALE-seq protocol for the effective annotation of novel and active retroelements.

Nanopore RNA Sequencing Revealed Long Non-Coding and LTR Retrotransposon-Related RNAs Expressed at Early Stages of Triticale SEED Development

The intergenic space of plant genomes encodes many functionally important yet unexplored RNAs. The genomic loci encoding these RNAs are often considered “junk”, DNA as they are frequently associated with repeat-rich regions of the genome. The latter makes the annotations of these loci and the assembly of the corresponding transcripts using short RNAseq reads particularly challenging. Here, using long-read Nanopore direct RNA sequencing, we aimed to identify these “junk” RNA molecules, including long non-coding RNAs (lncRNAs) and transposon-derived transcripts expressed during early stages (10 days post anthesis) of seed development of triticale (AABBRR, 2n = 6x = 42), an interspecific hybrid between wheat and rye. Altogether, we found 796 lncRNAs and 20 LTR retrotransposon-related transcripts (RTE-RNAs) expressed at this stage, with most of them being previously unannotated and located in the intergenic as well as intronic regions. Sequence analysis of the lncRNAs provide evidence for the frequent exonization of Class I (retrotransposons) and class II (DNA transposons) transposon sequences and suggest direct influence of “junk” DNA on the structure and origin of lncRNAs. We show that the expression patterns of lncRNAs and RTE-related transcripts have high stage specificity. In turn, almost half of the lncRNAs located in Genomes A and B have the highest expression levels at 10–30 days post anthesis in wheat. Detailed analysis of the protein-coding potential of the RTE-RNAs showed that 75% of them carry open reading frames (ORFs) for a diverse set of GAG proteins, the main component of virus-like particles of LTR retrotransposons. We further experimentally demonstrated that some RTE-RNAs originate from autonomous LTR retrotransposons with ongoing transposition activity during early stages of triticale seed development. Overall, our results provide a framework for further exploration of the newly discovered lncRNAs and RTE-RNAs in functional and genome-wide association studies in triticale and wheat. Our study also demonstrates that Nanopore direct RNA sequencing is an indispensable tool for the elucidation of lncRNA and retrotransposon transcripts.

Genomic and Transcriptomic Survey Provides New Insight into the Organization and Transposition Activity of Highly Expressed LTR Retrotransposons of Sunflower (Helianthus annuus L.).

LTR retrotransposons (RTEs) play a crucial role in plant genome evolution and adaptation. Although RTEs are generally silenced in somatic plant tissues under non-stressed conditions, some expressed RTEs (exRTEs) escape genome defense mechanisms. As our understanding of exRTE organization in plants is rudimentary, we systematically surveyed the genomic and transcriptomic organization and mobilome (transposition) activity of sunflower (Helianthus annuus L.) exRTEs. We identified 44 transcribed RTEs in the sunflower genome and demonstrated their distinct genomic features: more recent insertion time, longer open reading frame (ORF) length, and smaller distance to neighboring genes. We showed that GAG-encoding ORFs are present at significantly higher frequencies in exRTEs, compared with non-expressed RTEs. Most exRTEs exhibit variation in copy number among sunflower cultivars and one exRTE Gagarin produces extrachromosomal circular DNA in seedling, demonstrating recent and ongoing transposition activity. Nanopore direct RNA sequencing of full-length RTE RNA revealed complex patterns of alternative splicing in RTE RNAs, resulting in isoforms that carry ORFs for distinct RTE proteins. Together, our study demonstrates that tens of expressed sunflower RTEs with specific genomic organization shape the hidden layer of the transcriptome, pointing to the evolution of specific strategies that circumvent existing genome defense mechanisms.

Humic acid effects on retrotransposon polymorphisms caused by zinc and iron in the maize (Zea mays L.) genome

Long terminal repeat (LTR) retrotransposons are important pathways to examine various genome reorganizations based on environmental factors. The maize genome is one of the best representative examples. The aim of this study was to examine humic acid (HA) protective effects on genomic template stability (GTS) and LTR retrotransposon polymorphisms in corn seeds subjected to zinc (Zn) and iron (Fe) stress. In this study, maize seedlings were exposed to three doses (20, 40, and 60 mM) of ZnSO4·7H2O and FeSO4 and their combinations at 1500 ppm (0.075 g) HA. Inter-retrotransposon amplified polymorphism and retrotransposon-microsatellite amplified polymorphism techniques were used for genetic analyses. Results indicated that in all doses used, Zn and Fe increased retrotransposon polymorphisms and decreased the percentage of GTS via DNA damage. However, treatment of HA together with Zn and Fe resulted in decreased DNA damage and retrotransposon polymorphisms and increased GTS. These observations suggest that HA can be applied to reduce toxic effects in agricultural areas polluted with heavy metals.

Gene capture by transposable elements leads to epigenetic conflict in maize

Transposable elements (TEs) regularly capture fragments of genes. When the host silences these TEs, siRNAs homologous to the captured regions may also target the genes. This epigenetic crosstalk establishes an intragenomic conflict: silencing the TEs has the cost of silencing the genes. If genes are important, however, natural selection may maintain function by moderating the silencing response, which may also advantage the TEs. In this study, we examined this model by focusing on Helitrons, Pack-MULEs, and Sirevirus LTR retrotransposons in the maize genome. We documented 1263 TEs containing exon fragments from 1629 donor genes. Consistent with epigenetic conflict, donor genes mapped more siRNAs and were more methylated than genes with no evidence of capture. However, these patterns differed between syntelog versus translocated donor genes. Syntelogs appeared to maintain function, as measured by gene expression, consistent with moderation of silencing for functionally important genes. Epigenetic marks did not spread beyond their captured regions and 24nt crosstalk siRNAs were linked with CHH methylation. Translocated genes, in contrast, bore the signature of silencing. They were highly methylated and less expressed, but also overrepresented among donor genes and located away from chromosomal arms, which suggests a link between capture and gene movement. Splitting genes into potential functional categories based on evolutionary constraint supported the synteny-based findings. TE families captured genes in different ways, but the evidence for their advantage was generally less obvious; nevertheless, TEs with captured fragments were older, mapped fewer siRNAs, and were slightly less methylated than TEs without captured fragments. Collectively, our results argue that TE capture triggers an intragenomic conflict that may not affect the function of important genes but may lead to the pseudogenization of less-constrained genes.

Human Endogenous Retrovirus K (HERV-K) can drive gene expression as a promoter in Caenorhabditis elegans

Endogenous retroviruses (ERVs) are retrotransposons present in various metazoan genomes and have been implicated in metazoan evolution as well as in nematodes and humans. The long terminal repeat (LTR) retrotransposons contain several regulatory sequences including promoters and enhancers that regulate endogenous gene expression and thereby control organismal development and response to environmental change. ERVs including the LTR retrotransposons constitute 8% of the human genome and less than 0.6% of the Caenorhabditis elegans (C. elegans) genome, a nematode genetic model system. To investigate the evolutionarily conserved mechanism behind the transcriptional activity of retrotransposons, we generated a transgenic worm model driving green fluorescent protein (GFP) expression using Human endogenous retroviruses (HERV)-K LTR as a promoter. The promoter activity of HERV-K LTR was robust and fluorescence was observed in various tissues throughout the developmental process. Interestingly, persistent GFP expression was specifically detected in the adult vulva muscle. Using deletion constructs, we found that the region from positions 675 to 868 containing the TATA box was necessary for promoter activity driving gene expression in the vulva. Interestingly, we found that the promoter activity of the LTR was dependent on che-1 transcription factor, a sensory neuron driver, and lin-15b, a negative regulator of RNAi and germline gene expression. These results suggest evolutionary conservation of the LTR retrotransposon activity in transcriptional regulation as well as the possibility of che-1 function in non-neuronal tissues.

Genomic analyses reveal evolutionary and geologic context for the plateau fungus Ophiocordyceps sinensis

BackgroundOphiocordyceps sinensis, which is only naturally found in the high-elevation extreme environment of the Tibetan Plateau, has been used in traditional Chinese medicine. Information concerning the evolutionary and geologic context of O. sinensis remains limited, however.MethodsWe constructed the high-quality genome of O. sinensis and provided insight into the evolution and ecology of O. sinensis using comparative genomics.ResultsWe mapped the whole genome of the anamorph/asexual form Hirsutella of O. sinensis using Illumina and PacBio sequencing technologies and obtained a well assembled genome of 119.2 Mbp size. Long-read Single Molecule Real Time (SMRT) sequencing technology generated an assembly with more accurate representation of repeat sequence abundances and placement. Evolutionary analyses indicated that O. sinensis diverged from other fungi 65.9 Mya in the Upper Cretaceous, during the uplift of the Tibetan Plateau. Gene family expansions and contractions in addition to genome inflation via long terminal repeat (LTR) retrotransposon insertions were implicated as an important driver of O. sinensis divergence. The insertion rate of LTR sequences into the O. sinensis genome peaked ~ 30–40 Mya, when the Tibetan Plateau rose rapidly. Gene Ontology (GO) enrichment analysis suggested that O. sinensis contained more genes related to ice binding compared to other closely related fungi, which may aid in their adaptability to the cold Tibetan Plateau. Further, heavy metal resistance genes were in low abundance in the O. sinensis genome, which may help to explain previous observations that O. sinensis tissues contain high levels of heavy metals.ConclusionsOur results reveal the evolutionary, geological, and ecological context for the evolution of the O. sinensis genome and the factors that have contributed to the environmental adaptability of this valuable fungus. These findings suggest that genome inflation via LTR retrotransposon insertions in O. sinensis coincided with the uplift of the Tibetan Plateau. LTRs and the specific genetic mechanisms of O. sinensis contributed to its adaptation to the environment on the plateau.

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

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. 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. 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. 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.

RNA-Seq Transcriptome Analysis Reveals Long Terminal Repeat Retrotransposon Modulation in Human Peripheral Blood Mononuclear Cells after In Vivo Lipopolysaccharide Injection.

Human endogenous retroviruses (HERVs) and mammalian apparent long terminal repeat (LTR) retrotransposons (MaLRs) are retroviral sequences that integrated into germ line cells millions of years ago. Transcripts of these LTR retrotransposons are present in several tissues, and their expression is modulated in pathological conditions, although their function remains often far from being understood. Here, we focused on the HERV/MaLR expression and modulation in a scenario of immune system activation. We used a public data set of human peripheral blood mononuclear cells (PBMCs) RNA-Seq from 15 healthy participants to a clinical trial before and after exposure to lipopolysaccharide (LPS), for which we established an RNA-Seq workflow for the identification of expressed and modulated cellular genes and LTR retrotransposon elements.IMPORTANCE We described the HERV and MaLR transcriptome in PBMCs, finding that about 8.4% of the LTR retrotransposon loci were expressed and identifying the betaretrovirus-like HERVs as those with the highest percentage of expressed loci. We found 4,607 HERV and MaLR loci that were modulated as a result of in vivo stimulation with LPS. The HERV-H group showed the highest number of differentially expressed most intact proviruses. We characterized the HERV and MaLR loci as differentially expressed, checking their genomic context of insertion and observing a general colocalization with genes that are involved and modulated in the immune response, as a consequence of LPS stimulation. The analyses of HERV and MaLR expression and modulation show that these LTR retrotransposons are expressed in PBMCs and regulated in inflammatory settings. The similar regulation of HERVs/MaLRs and genes after LPS stimulation suggests possible interactions of LTR retrotransposons and the immune host response.

The cytosine methylation landscape of spring barley revealed by a new reduced representation bisulfite sequencing pipeline, WellMeth.

Patterns and level of cytosine methylation vary widely among plant species and are associated with genome size as well as the proportion of transposons and other repetitive elements in the genome. We explored epigenetic patterns and diversity in a representative proportion of the spring barley (Hordeum vulgare L.) genome across several commercial and historical cultivars. This study adapted a genotyping-by-sequencing (GBS) approach for the detection of methylated cytosines in genomic DNA. To analyze the data, we developed WellMeth, a complete pipeline for analysis of reduced representation bisulfite sequencing. WellMeth enabled quantification of context-specific DNA methylation at the single-base resolution as well as identification of differentially methylated sites (DMCs) and regions (DMRs). On average, DNA methylation levels were significantly higher than what is commonly observed in many plants species, reaching over 10-fold higher levels than those in Arabidopsis thaliana (L.) Heynh. in the CHH methylation. Preferential methylation was observed within and at the edges of long-terminal repeats (LTR) retrotransposons Gypsy and Copia. From a pairwise comparison of cultivars, numerous DMRs could be identified of which more than 5,000 were conserved within the analyzed set of barley cultivars. The subset of regions overlapping with genes showed enrichment in gene ontology (GO) categories associated with chromatin and cellular structure and organization. A significant correlation between genetic and epigenetic distances suggests that a considerable portion of methylated regions is under strict genetic control in barley. The data presented herein represents the first step in efforts toward a better understanding of genome-level structural and functional aspects of methylation in barley.

TIP_finder: An HPC Software to Detect Transposable Element Insertion Polymorphisms in Large Genomic Datasets.

Transposable elements (TEs) are non-static genomic units capable of moving indistinctly from one chromosomal location to another. Their insertion polymorphisms may cause beneficial mutations, such as the creation of new gene function, or deleterious in eukaryotes, e.g., different types of cancer in humans. A particular type of TE called LTR-retrotransposons comprises almost 8% of the human genome. Among LTR retrotransposons, human endogenous retroviruses (HERVs) bear structural and functional similarities to retroviruses. Several tools allow the detection of transposon insertion polymorphisms (TIPs) but fail to efficiently analyze large genomes or large datasets. Here, we developed a computational tool, named TIP_finder, able to detect mobile element insertions in very large genomes, through high-performance computing (HPC) and parallel programming, using the inference of discordant read pair analysis. TIP_finder inputs are (i) short pair reads such as those obtained by Illumina, (ii) a chromosome-level reference genome sequence, and (iii) a database of consensus TE sequences. The HPC strategy we propose adds scalability and provides a useful tool to analyze huge genomic datasets in a decent running time. TIP_finder accelerates the detection of transposon insertion polymorphisms (TIPs) by up to 55 times in breast cancer datasets and 46 times in cancer-free datasets compared to the fastest available algorithms. TIP_finder applies a validated strategy to find TIPs, accelerates the process through HPC, and addresses the issues of runtime for large-scale analyses in the post-genomic era.

Fundamentally different repetitive element composition of sex chromosomes in Rumex acetosa.

Dioecious species with well-established sex chromosomes are rare in the plant kingdom. Most sex chromosomes increase in size but no comprehensive analysis of the kind of sequences that drive this expansion has been presented. Here we analyse sex chromosome structure in common sorrel (Rumex acetosa), a dioecious plant with XY1Y2 sex determination, and we provide the first chromosome-specific repeatome analysis for a plant species possessing sex chromosomes. We flow-sorted and separately sequenced sex chromosomes and autosomes in R. acetosa using the two-dimensional fluorescence in situ hybridization in suspension (FISHIS) method and Illumina sequencing. We identified and quantified individual repeats using RepeatExplorer, Tandem Repeat Finder and the Tandem Repeats Analysis Program. We employed fluorescence in situ hybridization (FISH) to analyse the chromosomal localization of satellites and transposons. We identified a number of novel satellites, which have, in a fashion similar to previously known satellites, significantly expanded on the Y chromosome but not as much on the X or on autosomes. Additionally, the size increase of Y chromosomes is caused by non-long terminal repeat (LTR) and LTR retrotransposons, while only the latter contribute to the enlargement of the X chromosome. However, the X chromosome is populated by different LTR retrotransposon lineages than those on Y chromosomes. The X and Y chromosomes have significantly diverged in terms of repeat composition. The lack of recombination probably contributed to the expansion of diverse satellites and microsatellites and faster fixation of newly inserted transposable elements (TEs) on the Y chromosomes. In addition, the X and Y chromosomes, despite similar total counts of TEs, differ significantly in the representation of individual TE lineages, which indicates that transposons proliferate preferentially in either the paternal or the maternal lineage.

Improved Brassica oleracea JZS assembly reveals significant changing of LTR-RT dynamics in different morphotypes.

Brassica oleracea is an important vegetable crop that has provided ancestor genomes of the two most important Brassica oil crops, Brassica napus and Brassica carinata. The current B. oleracea reference genome (JZS, also named 02-12) displays problems of large mis-assemblies, low sequence continuity, and low assembly integrity, thus limiting genomic analysis. We reported an updated assembly of the B. oleracea reference genome (JZS v2) obtained through single-molecule sequencing and chromosome conformation capture technologies. We assembled an additional 83.16Mb of genomic sequences, and the updated genome features a contig N50 size of 2.37Mb, representing an ~ 88-fold improvement. We detected a new round of long terminal repeat retrotransposon (LTR-RT) burst in the new assembly. Comparative analysis with the reported genome sequences of two other genomes of B. oleracea (TO1000 and HDEM) identified extensive gene order and gene structural variation. In addition, we found that the genome-specific amplification of Gypsy-like LTR-RTs occurred around 0-1 million years ago (MYA). In particular, the athila, tat, and Del families were extensively amplified in JZS around 0-1 MYA. Moreover, we identified that the syntenic genes were modified due to the insertion of genome-specific LTR-RTs. These results indicated that the genome-specific LTR-RT dynamics were associated with genome diversification in B. oleracea.