ITIS, a bioinformatics tool for accurate identification of transposon insertion sites using next-generation sequencing data

Application of next-generation sequencing in the detection of transgenic crop

With the rapid development of transgenic technology and the increasing prevalence of genetically modified (GMO) crops, incidents such as illegal importation, environmental contamination, and safety concerns associated with GMOs have risen significantly in recent years. Consequently, there is a growing demand for more advanced methods of GMO crop detection. Traditional molecular detection techniques, which rely on nucleic acids or proteins, are becoming less effective due to the increasing complexity of GMO crop genomes. In contrast, detection technologies based on second- and third-generation high-throughput sequencing offer promising solutions to these challenges. This review provides a comprehensive overview of the latest advancements in GMO crop detection technologies, categorizing and describing various approaches, and comparing their respective strengths and limitations. The article emphasizes the current state, benefits, challenges, and future prospects of high-throughput sequencing in GMO detection, aiming to guide further research and development in this field.

Homeostasis of flavonoids and triterpenoids most likely modulates starch metabolism for pollen tube penetration in rice

In angiosperms, the timely delivery of sperm cell nuclei by pollen tube (PT) to the ovule is vital for double fertilization. Penetration of PT into maternal stigma tissue is a critical step for sperm cell nuclei delivery, yet little is known about the process. Here, a male-specific and sporophytic mutant xt6, where PTs are able to germinate but unable to penetrate the stigma tissue, is reported in Oryza sativa. Through genetic study, the causative gene was identified as Chalcone synthase (OsCHS1), encoding the first enzyme in flavonoid biosynthesis. Indeed, flavonols were undetected in mutant pollen grains and PTs, indicating that the mutation abolished flavonoid biosynthesis. Nevertheless, the phenotype cannot be rescued by exogenous application of quercetin and kaempferol as reported in maize and petunia, suggesting a different mechanism exists in rice. Further analysis showed that loss of OsCHS1 function disrupted the homeostasis of flavonoid and triterpenoid metabolism and led to the accumulation of triterpenoid, which inhibits significantly α-amylase activity, amyloplast hydrolysis and monosaccharide content in xt6, these ultimately impaired tricarboxylic acid (TCA) cycle, reduced ATP content and lowered the turgor pressure as well. Our findings reveal a new mechanism that OsCHS1 modulates starch hydrolysis and glycometabolism through modulating the metabolic homeostasis of flavonoids and triterpenoids which affects α-amylase activity to maintain PT penetration in rice, which contributes to a better understanding of the function of CHS1 in crop fertility and breeding.

Identification and Genotyping of Transposable Element Insertions From Genome Sequencing Data

Transposable element (TE) mobilization is a significant source of genomic variation and has been associated with various human diseases. The exponential growth of population-scale whole-genome sequencing and rapid innovations in long-read sequencing technologies provide unprecedented opportunities to study TE insertions and their functional impact in human health and disease. Identifying TE insertions, however, is challenging due to the repetitive nature of the TE sequences. Here, we review computational approaches to detecting and genotyping TE insertions using short- and long-read sequencing and discuss the strengths and weaknesses of different approaches. © 2020 Wiley Periodicals LLC.

The F-box protein MIO1/SLB1 regulates organ size and leaf movement in Medicago truncatula

The size of leaf and seed organs, determined by the interplay of cell proliferation and expansion, is closely related to the final yield and quality of forage and crops. Yet the cellular and molecular mechanisms underlying organ size modulation remain poorly understood, especially in legumes. Here, MINI ORGAN1 (MIO1), which encodes an F-box protein SMALL LEAF AND BUSHY1 (SLB1) recently reported to control lateral branching in Medicago truncatula, was identified as a key regulator of organ size. We show that loss-of-function of MIO1/SLB1 severely reduced organ size. Conversely, plants overexpressing MIO1/SLB1 had enlarged organs. Cellular analysis revealed that MIO1/SLB1 controlled organ size mainly by modulating primary cell proliferation during the early stages of leaf development. Biochemical analysis revealed that MIO1/SLB1 could form part of SKP1/Cullin/F-box (SCF) E3 ubiquitin ligase complex, to target BIG SEEDS1 (BS1), a repressor of primary cell division, for degradation. Interestingly, we found that MIO1/SLB1 also played a key role in pulvinus development and leaf movement by modulating cell proliferation of the pulvinus as leaves developed. Our study not only demonstrates a conserved role of MIO1/SLB1 in the control of organ size in legumes, but also sheds light on the novel function of MIO1/SLB1 in leaf movement.

The geometry of the compound leaf plays a significant role in the leaf movement of Medicago truncatula modulated by mtdwarf4a

In most legumes, two typical features found in leaves are diverse compound forms and the pulvinus-driven nyctinastic movement. Many genes have been identified for leaf-shape determination, but the underlying nature of leaf movement as well as its association with the compound form remains largely unknown. Using forward-genetic screening and whole-genome resequencing, we found that two allelic mutants of Medicago truncatula with unclosed leaflets at night were impaired in MtDWARF4A (MtDWF4A), a gene encoding a cytochrome P450 protein orthologous to Arabidopsis DWARF4. The mtdwf4a mutant also had a mild brassinosteroid (BR)-deficient phenotype bearing pulvini without significant deficiency in organ identity. Both mtdwf4a and dwf4 could be fully rescued by MtDWF4A, and mtdwf4a could close their leaflets at night after the application of exogenous 24-epi-BL. Surgical experiments and genetic analysis of double mutants revealed that the failure to exhibit leaf movement in mtdwf4a is a consequence of the physical obstruction of the overlapping leaflet laminae, suggesting a proper geometry of leaflets is important for their movement in M. truncatula. These observations provide a novel insight into the nyctinastic movement of compound leaves, shedding light on the importance of open space for organ movements in plants.

Easymap: A User-Friendly Software Package for Rapid Mapping-by-Sequencing of Point Mutations and Large Insertions

Mapping-by-sequencing strategies combine next-generation sequencing (NGS) with classical linkage analysis, allowing rapid identification of the causal mutations of the phenotypes exhibited by mutants isolated in a genetic screen. Computer programs that analyze NGS data obtained from a mapping population of individuals derived from a mutant of interest to identify a causal mutation are available; however, the installation and usage of such programs requires bioinformatic skills, modifying or combining pieces of existing software, or purchasing licenses. To ease this process, we developed Easymap, an open-source program that simplifies the data analysis workflows from raw NGS reads to candidate mutations. Easymap can perform bulked segregant mapping of point mutations induced by ethyl methanesulfonate (EMS) with DNA-seq or RNA-seq datasets, as well as tagged-sequence mapping for large insertions, such as transposons or T-DNAs. The mapping analyses implemented in Easymap have been validated with experimental and simulated datasets from different plant and animal model species. Easymap was designed to be accessible to all users regardless of their bioinformatics skills by implementing a user-friendly graphical interface, a simple universal installation script, and detailed mapping reports, including informative images and complementary data for assessment of the mapping results. Easymap is available at http://genetics.edu.umh.es/resources/easymap; its Quickstart Installation Guide details the recommended procedure for installation.

Epigenetic conflict on a degenerating Y chromosome increases mutational burden in Drosophila males

Large portions of eukaryotic genomes consist of transposable elements (TEs), and the establishment of transcription-repressing heterochromatin during early development safeguards genome integrity in Drosophila. Repeat-rich Y chromosomes can act as reservoirs for TEs ('toxic' Y effect), and incomplete epigenomic defenses during early development can lead to deleterious TE mobilization. Here, we contrast the dynamics of early TE activation in two Drosophila species with vastly different Y chromosomes of different ages. Zygotic TE expression is elevated in male embryos relative to females in both species, mostly due to expression of Y-linked TEs. Interestingly, male-biased TE expression diminishes across development in D. pseudoobscura, but remains elevated in D. miranda, the species with the younger and larger Y chromosome. The repeat-rich Y of D. miranda still contains many actively transcribed genes, which compromise the formation of silencing heterochromatin. Elevated TE expression results in more de novo insertions of repeats in males compared to females. This lends support to the idea that the 'toxic' Y chromosome can create a mutational burden in males when genome-wide defense mechanisms are compromised, and suggests a previously unappreciated epigenetic conflict on evolving Y chromosomes between transcription of essential genes and silencing of selfish DNA.

A study of male fertility control in Medicago truncatula uncovers an evolutionarily conserved recruitment of two tapetal bHLH subfamilies in plant sexual reproduction

Male sterility is an important tool for plant breeding and hybrid seed production. Male-sterile mutants are largely due to an abnormal development of either the sporophytic or gametophytic anther tissues. Tapetum, a key sporophytic tissue, provides nutrients for pollen development, and its delayed degeneration induces pollen abortion. Numerous bHLH proteins have been documented to participate in the degeneration of the tapetum in angiosperms, but relatively little attention has been given to the evolution of the involved developmental pathways across the phylogeny of land plants. A combination of cellular, molecular, biochemical and evolutionary analyses was used to investigate the male fertility control in Medicago truncatula. We characterized the male-sterile mutant empty anther1 (ean1) and identified EAN1 as a tapetum-specific bHLH transcription factor necessary for tapetum degeneration. Our study uncovered an evolutionarily conserved recruitment of bHLH subfamily II and III(a + c)1 in the regulation of tapetum degeneration. EAN1 belongs to the subfamily II and specifically forms heterodimers with the subfamily III(a + c)1 members, which suggests a heterodimerization mechanism conserved in angiosperms. Our work suggested that the pathway of two tapetal-bHLH subfamilies is conserved in all land plants, and likely was established before the divergence of the spore-producing land plants.

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. TIP_finder version 1.0 is available at https://github.com/simonorozcoarias/TIP_finder.

Insertional mutagenesis of Brachypodium distachyon using the Tnt1 retrotransposable element

Brachypodium distachyon is an annual C3 grass used as a monocot model system in functional genomics research. Insertional mutagenesis is a powerful tool for both forward and reverse genetics studies. In this study, we explored the possibility of using the tobacco retrotransposon Tnt1 to create a transposon-based insertion mutant population in B. distachyon. We developed transgenic B. distachyon plants expressing Tnt1 (R0) and in the subsequent regenerants (R1) we observed that Tnt1 actively transposed during somatic embryogenesis, generating an average of 6.37 insertions per line in a population of 19 independent R1 regenerant plants analyzed. In seed-derived progeny of R1 plants, Tnt1 segregated in a Mendelian ratio of 3:1 and no new Tnt1 transposition was observed. A total of 126 flanking sequence tags (FSTs) were recovered from the analyzed R0 and R1 lines. Analysis of the FSTs showed a uniform pattern of insertion in all the chromosomes (1-5) without any preference for a particular chromosome region. Considering the average length of a gene transcript to be 3.37 kb, we estimated that 29 613 lines are required to achieve a 90% possibility of tagging a given gene in the B. distachyon genome using the Tnt1-based mutagenesis approach. Our results show the possibility of using Tnt1 to achieve near-saturation mutagenesis in B. distachyon, which will aid in functional genomics studies of other C3 grasses.

Cloning and Functional Analysis of Dwarf Gene Mini Plant 1 (MNP1) in Medicago truncatula

Plant height is a vital agronomic trait that greatly determines crop yields because of the close relationship between plant height and lodging resistance. Legumes play a unique role in the worldwide agriculture; however, little attention has been given to the molecular basis of their height. Here, we characterized the first dwarf mutant mini plant 1 (mnp1) of the model legume plant Medicago truncatula. Our study found that both cell length and the cell number of internodes were reduced in a mnp1 mutant. Using the forward genetic screening and subsequent whole-genome resequencing approach, we cloned the MNP1 gene and found that it encodes a putative copalyl diphosphate synthase (CPS) implicated in the first step of gibberellin (GA) biosynthesis. MNP1 was highly homologous to Pisum sativum LS. The subcellular localization showed that MNP1 was located in the chloroplast. Further analysis indicated that GA3 could significantly restore the plant height of mnp1-1, and expression of MNP1 in a cps1 mutant of Arabidopsis partially rescued its mini-plant phenotype, indicating the conservation function of MNP1 in GA biosynthesis. Our results provide valuable information for understanding the genetic regulation of plant height in M. truncatula.

Lateral Leaflet Suppression 1 (LLS1), encoding the MtYUCCA1 protein, regulates lateral leaflet development in Medicago truncatula

In species with compound leaves, the positions of leaflet primordium initiation are associated with local peaks of auxin accumulation. However, the role of auxin during the late developmental stages and outgrowth of compound leaves remains largely unknown. Using genome resequencing approaches, we identified insertion sites at four alleles of the LATERAL LEAFLET SUPPRESSION1 (LLS1) gene, encoding the auxin biosynthetic enzyme YUCCA1 in Medicago truncatula. Linkage analysis and complementation tests showed that the lls1 mutant phenotypes were caused by the Tnt1 insertions that disrupted the LLS1 gene. The transcripts of LLS1 can be detected in primordia at early stages of leaf initiation and later in the basal regions of leaflets, and finally in vein tissues at late leaf developmental stages. Vein numbers and auxin content are reduced in the lls1-1 mutant. Analysis of the lls1 sgl1 and lls1 palm1 double mutants revealed that SGL1 is epistatic to LLS1, and LLS1 works with PALM1 in an independent pathway to regulate the growth of lateral leaflets. Our work demonstrates that the YUCCA1/YUCCA4 subgroup plays very important roles in the outgrowth of lateral leaflets during compound leaf development of M. truncatula, in addition to leaf venation.

ERVcaller: identifying polymorphic endogenous retrovirus and other transposable element insertions using whole-genome sequencing data

MOTIVATION

Approximately 8% of the human genome is derived from endogenous retroviruses (ERVs). In recent years, an increasing number of human diseases have been found to be associated with ERVs. However, it remains challenging to accurately detect the full spectrum of polymorphic (unfixed) ERVs using whole-genome sequencing (WGS) data.

RESULTS

We designed a new tool, ERVcaller, to detect and genotype transposable element (TE) insertions, including ERVs, in the human genome. We evaluated ERVcaller using both simulated and real benchmark WGS datasets. Compared to existing tools, ERVcaller consistently obtained both the highest sensitivity and precision for detecting simulated ERV and other TE insertions derived from real polymorphic TE sequences. For the WGS data from the 1000 Genomes Project, ERVcaller detected the largest number of TE insertions per sample based on consensus TE loci. By analyzing the experimentally verified TE insertions, ERVcaller had 94.0% TE detection sensitivity and 96.6% genotyping accuracy. Polymerase chain reaction and Sanger sequencing in a small sample set verified 86.7% of examined insertion statuses and 100% of examined genotypes. In conclusion, ERVcaller is capable of detecting and genotyping TE insertions using WGS data with both high sensitivity and precision. This tool can be applied broadly to other species.

AVAILABILITY AND IMPLEMENTATION

http://www.uvm.edu/genomics/software/ERVcaller.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Next-generation forward genetic screens: using simulated data to improve the design of mapping-by-sequencing experiments in Arabidopsis

Forward genetic screens have successfully identified many genes and continue to be powerful tools for dissecting biological processes in Arabidopsis and other model species. Next-generation sequencing technologies have revolutionized the time-consuming process of identifying the mutations that cause a phenotype of interest. However, due to the cost of such mapping-by-sequencing experiments, special attention should be paid to experimental design and technical decisions so that the read data allows to map the desired mutation. Here, we simulated different mapping-by-sequencing scenarios. We first evaluated which short-read technology was best suited for analyzing gene-rich genomic regions in Arabidopsis and determined the minimum sequencing depth required to confidently call single nucleotide variants. We also designed ways to discriminate mutagenesis-induced mutations from background Single Nucleotide Polymorphisms in mutants isolated in Arabidopsis non-reference lines. In addition, we simulated bulked segregant mapping populations for identifying point mutations and monitored how the size of the mapping population and the sequencing depth affect mapping precision. Finally, we provide the computational basis of a protocol that we already used to map T-DNA insertions with paired-end Illumina-like reads, using very low sequencing depths and pooling several mutants together; this approach can also be used with single-end reads as well as to map any other insertional mutagen. All these simulations proved useful for designing experiments that allowed us to map several mutations in Arabidopsis.

A benchmark of transposon insertion detection tools using real data

Background

Transposable elements (TEs) are an important source of genomic variability in eukaryotic genomes. Their activity impacts genome architecture and gene expression and can lead to drastic phenotypic changes. Therefore, identifying TE polymorphisms is key to better understand the link between genotype and phenotype. However, most genotype-to-phenotype analyses have concentrated on single nucleotide polymorphisms as they are easier to reliable detect using short-read data. Many bioinformatic tools have been developed to identify transposon insertions from resequencing data using short reads. Nevertheless, the performance of most of these tools has been tested using simulated insertions, which do not accurately reproduce the complexity of natural insertions.

Results

We have overcome this limitation by building a dataset of insertions from the comparison of two high-quality rice genomes, followed by extensive manual curation. This dataset contains validated insertions of two very different types of TEs, LTR-retrotransposons and MITEs. Using this dataset, we have benchmarked the sensitivity and precision of 12 commonly used tools, and our results suggest that in general their sensitivity was previously overestimated when using simulated data. Our results also show that, increasing coverage leads to a better sensitivity but with a cost in precision. Moreover, we found important differences in tool performance, with some tools performing better on a specific type of TEs. We have also used two sets of experimentally validated insertions in Drosophila and humans and show that this trend is maintained in genomes of different size and complexity.

Conclusions

We discuss the possible choice of tools depending on the goals of the study and show that the appropriate combination of tools could be an option for most approaches, increasing the sensitivity while maintaining a good precision.

TDNAscan: A Software to Identify Complete and Truncated T-DNA Insertions

Transfer (T)-DNA insertions in mutants isolated from forward genetic screens are typically identified through thermal asymmetric interlaced polymerase chain reaction (TAIL-PCR), inverse PCR, or plasmid rescue. Despite the popularity and success of these methods, they have limited capabilities, particularly in situations in which the T-DNA is truncated. Here, we present a next generation sequencing (NGS)-based platform to facilitate the identification of complete and truncated T-DNA insertions. Our method enables the detection of the corresponding T-DNA insertion orientation and zygosity as well as insertion annotation. This method, called TDNAscan, was developed to be an open source software. We expect that TDNAscan will be a valuable addition to forward genetics toolkits because it provides a solution to the problem of causal gene identification, particularly genes disrupted by truncated T-DNA insertions. We present a case study in which TDNAscan was used to determine that the recessive Arabidopsis thaliana hypersensitive to latrunculin B (hlb3) mutant isolated in a forward genetic screen of T-DNA mutagenized plants encodes a class II FORMIN.

TranSurVeyor: an improved database-free algorithm for finding non-reference transpositions in high-throughput sequencing data

Transpositions transfer DNA segments between different loci within a genome; in particular, when a transposition is found in a sample but not in a reference genome, it is called a non-reference transposition. They are important structural variations that have clinical impact. Transpositions can be called by analyzing second generation high-throughput sequencing datasets. Current methods follow either a database-based or a database-free approach. Database-based methods require a database of transposable elements. Some of them have good specificity; however this approach cannot detect novel transpositions, and it requires a good database of transposable elements, which is not yet available for many species. Database-free methods perform de novo calling of transpositions, but their accuracy is low. We observe that this is due to the misalignment of the reads; since reads are short and the human genome has many repeats, false alignments create false positive predictions while missing alignments reduce the true positive rate. This paper proposes new techniques to improve database-free non-reference transposition calling: first, we propose a realignment strategy called one-end remapping that corrects the alignments of reads in interspersed repeats; second, we propose a SNV-aware filter that removes some incorrectly aligned reads. By combining these two techniques and other techniques like clustering and positive-to-negative ratio filter, our proposed transposition caller TranSurVeyor shows at least 3.1-fold improvement in terms of F1-score over existing database-free methods. More importantly, even though TranSurVeyor does not use databases of prior information, its performance is at least as good as existing database-based methods such as MELT, Mobster and Retroseq. We also illustrate that TranSurVeyor can discover transpositions that are not known in the current database.

Genome-wide variation in nucleotides and retrotransposons in alpine populations of Arabis alpina (Brassicaceae)

Advances in high-throughput sequencing have promoted the collection of reference genomes and genome-wide diversity. However, the assessment of genomic variation among populations has hitherto mainly been surveyed through single-nucleotide polymorphisms (SNPs) and largely ignored the often major fraction of genomes represented by transposable elements (TEs). Despite accumulating evidence supporting the evolutionary significance of TEs, comprehensive surveys remain scarce. Here, we sequenced the full genomes of 304 individuals of Arabis alpina sampled from four nearby natural populations to genotype SNPs as well as polymorphic long terminal repeat retrotransposons (polymorphic TEs; i.e., presence/absence of TE insertions at specific loci). We identified 291,396 SNPs and 20,548 polymorphic TEs, comparing their contributions to genomic diversity and divergence across populations. Few SNPs were shared among populations and overall showed high population-specific variation, whereas most polymorphic TEs segregated among populations. The genomic context of these two classes of variants further highlighted candidate adaptive loci having a putative impact on functional genes. In particular, 4.96% of the SNPs were identified as nonsynonymous or affecting start/stop codons. In contrast, 43% of the polymorphic TEs were present next to Arabis genes enriched in functional categories related to the regulation of reproduction and responses to biotic as well as abiotic stresses. This unprecedented data set, mapping variation gained from SNPs and complementary polymorphic TEs within and among populations, will serve as a rich resource for addressing microevolutionary processes shaping genome variation.

Genome-wide analysis of flanking sequences reveals that Tnt1 insertion is positively correlated with gene methylation in Medicago truncatula

From a single transgenic line harboring five Tnt1 transposon insertions, we generated a near-saturated insertion population in Medicago truncatula. Using thermal asymmetric interlaced-polymerase chain reaction followed by sequencing, we recovered 388 888 flanking sequence tags (FSTs) from 21 741 insertion lines in this population. FST recovery from 14 Tnt1 lines using the whole-genome sequencing (WGS) and/or Tnt1-capture sequencing approaches suggests an average of 80 insertions per line, which is more than the previous estimation of 25 insertions. Analysis of the distribution pattern and preference of Tnt1 insertions showed that Tnt1 is overall randomly distributed throughout the M. truncatula genome. At the chromosomal level, Tnt1 insertions occurred on both arms of all chromosomes, with insertion frequency negatively correlated with the GC content. Based on 174 546 filtered FSTs that show exact insertion locations in the M. truncatula genome version 4.0 (Mt4.0), 0.44 Tnt1 insertions occurred per kb, and 19 583 genes contained Tnt1 with an average of 3.43 insertions per gene. Pathway and gene ontology analyses revealed that Tnt1-inserted genes are significantly enriched in processes associated with 'stress', 'transport', 'signaling' and 'stimulus response'. Surprisingly, gene groups with higher methylation frequency were more frequently targeted for insertion. Analysis of 19 583 Tnt1-inserted genes revealed that 59% (1265) of 2144 transcription factors, 63% (765) of 1216 receptor kinases and 56% (343) of 616 nucleotide-binding site-leucine-rich repeat genes harbored at least one Tnt1 insertion, compared with the overall 38% of Tnt1-inserted genes out of 50 894 annotated genes in the genome.

On the Population Dynamics of Junk: A Review on the Population Genomics of Transposable Elements

Transposable elements (TEs) play an important role in shaping genomic organization and structure, and may cause dramatic changes in phenotypes. Despite the genetic load they may impose on their host and their importance in microevolutionary processes such as adaptation and speciation, the number of population genetics studies focused on TEs has been rather limited so far compared to single nucleotide polymorphisms (SNPs). Here, we review the current knowledge about the dynamics of transposable elements at recent evolutionary time scales, and discuss the mechanisms that condition their abundance and frequency. We first discuss non-adaptive mechanisms such as purifying selection and the variable rates of transposition and elimination, and then focus on positive and balancing selection, to finally conclude on the potential role of TEs in causing genomic incompatibilities and eventually speciation. We also suggest possible ways to better model TEs dynamics in a population genomics context by incorporating recent advances in TEs into the rich information provided by SNPs about the demography, selection, and intrinsic properties of genomes.

iMGEins: detecting novel mobile genetic elements inserted in individual genomes

Background

Recent advances in sequencing technology have allowed us to investigate personal genomes to find structural variations, which have been studied extensively to identify their association with the physiology of diseases such as cancer. In particular, mobile genetic elements (MGEs) are one of the major constituents of the human genomes, and cause genome instability by insertion, mutation, and rearrangement.

Result

We have developed a new program, iMGEins, to identify such novel MGEs by using sequencing reads of individual genomes, and to explore the breakpoints with the supporting reads and MGEs detected. iMGEins is the first MGE detection program that integrates three algorithmic components: discordant read-pair mapping, split-read mapping, and insertion sequence assembly. Our evaluation results showed its outstanding performance in detecting novel MGEs from simulated genomes, as well as real personal genomes. In detail, the average recall and precision rates of iMGEins are 96.67 and 100%, respectively, which are the highest among the programs compared. In the testing with real human genomes of the NA12878 sample, iMGEins shows the highest accuracy in detecting MGEs within 20 bp proximity of the breakpoints annotated.

Conclusion

In order to study the dynamics of MGEs in individual genomes, iMGEins was developed to accurately detect breakpoints and report inserted MGEs. Compared with other programs, iMGEins has valuable features of identifying novel MGEs and assembling the MGEs inserted.

Electronic supplementary material

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

Benchmarking computational tools for polymorphic transposable element detection

Transposable elements (TEs) are an important source of human genetic variation with demonstrable effects on phenotype. Recently, a number of computational methods for the detection of polymorphic TE (polyTE) insertion sites from next-generation sequence data have been developed. The use of such tools will become increasingly important as the pace of human genome sequencing accelerates. For this report, we performed a comparative benchmarking and validation analysis of polyTE detection tools in an effort to inform their selection and use by the TE research community. We analyzed a core set of seven tools with respect to ease of use and accessibility, polyTE detection performance and runtime parameters. An experimentally validated set of 893 human polyTE insertions was used for this purpose, along with a series of simulated data sets that allowed us to assess the impact of sequence coverage on tool performance. The recently developed tool MELT showed the best overall performance followed by Mobster and then RetroSeq. PolyTE detection tools can best detect Alu insertion events in the human genome with reduced reliability for L1 insertions and substantially lowered performance for SVA insertions. We also show evidence that different polyTE detection tools are complementary with respect to their ability to detect a complete set of insertion events. Accordingly, a combined approach, coupled with manual inspection of individual results, may yield the best overall performance. In addition to the benchmarking results, we also provide notes on tool installation and usage as well as suggestions for future polyTE detection algorithm development.

Transposable element activity, genome regulation and human health

A convergence of novel genome analysis technologies is enabling population genomic studies of human transposable elements (TEs). Population surveys of human genome sequences have uncovered thousands of individual TE insertions that segregate as common genetic variants, i.e. TE polymorphisms. These recent TE insertions provide an important source of naturally occurring human genetic variation. Investigators are beginning to leverage population genomic data sets to execute genome-scale association studies for assessing the phenotypic impact of human TE polymorphisms. For example, the expression quantitative trait loci (eQTL) analytical paradigm has recently been used to uncover hundreds of associations between human TE insertion variants and gene expression levels. These include population-specific gene regulatory effects as well as coordinated changes to gene regulatory networks. In addition, analyses of linkage disequilibrium patterns with previously characterized genome-wide association study (GWAS) trait variants have uncovered TE insertion polymorphisms that are likely causal variants for a variety of common complex diseases. Gene regulatory mechanisms that underlie specific disease phenotypes have been proposed for a number of these trait associated TE polymorphisms. These new population genomic approaches hold great promise for understanding how ongoing TE activity contributes to functionally relevant genetic variation within and between human populations.

Functional Genomics and Seed Development in Medicago truncatula: An Overview

The study of seed development in the model species Medicago truncatula has made a significant contribution to our understanding of this process in crop legumes. Thanks to the availability of comprehensive proteomics and transcriptomics databases, coupled with exhaustive mutant collections, the roles of several regulatory genes in development and maturation are beginning to be deciphered and functionally validated. Advances in next-generation sequencing and the availability of a genomic sequence have made feasible high-density SNP genotyping, allowing the identification of markers tightly linked to traits of agronomic interest. A further major advance is to be expected from the integration of omics resources in functional network construction, which has been used recently to identify "hub" genes central to important traits.

McClintock: An Integrated Pipeline for Detecting Transposable Element Insertions in Whole-Genome Shotgun Sequencing Data

Transposable element (TE) insertions are among the most challenging types of variants to detect in genomic data because of their repetitive nature and complex mechanisms of replication . Nevertheless, the recent availability of large resequencing data sets has spurred the development of many new methods to detect TE insertions in whole-genome shotgun sequences. Here we report an integrated bioinformatics pipeline for the detection of TE insertions in whole-genome shotgun data, called McClintock (https://github.com/bergmanlab/mcclintock), which automatically runs and standardizes output for multiple TE detection methods. We demonstrate the utility of McClintock by evaluating six TE detection methods using simulated and real genome data from the model microbial eukaryote, Saccharomyces cerevisiae We find substantial variation among McClintock component methods in their ability to detect nonreference TEs in the yeast genome, but show that nonreference TEs at nearly all biologically realistic locations can be detected in simulated data by combining multiple methods that use split-read and read-pair evidence. In general, our results reveal that split-read methods detect fewer nonreference TE insertions than read-pair methods, but generally have much higher positional accuracy. Analysis of a large sample of real yeast genomes reveals that most McClintock component methods can recover known aspects of TE biology in yeast such as the transpositional activity status of families, target preferences, and target site duplication structure, albeit with varying levels of accuracy. Our work provides a general framework for integrating and analyzing results from multiple TE detection methods, as well as useful guidance for researchers studying TEs in yeast resequencing data.

RelocaTE2: a high resolution transposable element insertion site mapping tool for population resequencing

Background

Transposable element (TE) polymorphisms are important components of population genetic variation. The functional impacts of TEs in gene regulation and generating genetic diversity have been observed in multiple species, but the frequency and magnitude of TE variation is under appreciated. Inexpensive and deep sequencing technology has made it affordable to apply population genetic methods to whole genomes with methods that identify single nucleotide and insertion/deletion polymorphisms. However, identifying TE polymorphisms, particularly transposition events or non-reference insertion sites can be challenging due to the repetitive nature of these sequences, which hamper both the sensitivity and specificity of analysis tools.

Methods

We have developed the tool RelocaTE2 for identification of TE insertion sites at high sensitivity and specificity. RelocaTE2 searches for known TE sequences in whole genome sequencing reads from second generation sequencing platforms such as Illumina. These sequence reads are used as seeds to pinpoint chromosome locations where TEs have transposed. RelocaTE2 detects target site duplication (TSD) of TE insertions allowing it to report TE polymorphism loci with single base pair precision.

Results and Discussion

The performance of RelocaTE2 is evaluated using both simulated and real sequence data. RelocaTE2 demonstrate high level of sensitivity and specificity, particularly when the sequence coverage is not shallow. In comparison to other tools tested, RelocaTE2 achieves the best balance between sensitivity and specificity. In particular, RelocaTE2 performs best in prediction of TSDs for TE insertions. Even in highly repetitive regions, such as those tested on rice chromosome 4, RelocaTE2 is able to report up to 95% of simulated TE insertions with less than 0.1% false positive rate using 10-fold genome coverage resequencing data. RelocaTE2 provides a robust solution to identify TE insertion sites and can be incorporated into analysis workflows in support of describing the complete genotype from light coverage genome sequencing.

Population and clinical genetics of human transposable elements in the (post) genomic era

Recent technological developments-in genomics, bioinformatics and high-throughput experimental techniques-are providing opportunities to study ongoing human transposable element (TE) activity at an unprecedented level of detail. It is now possible to characterize genome-wide collections of TE insertion sites for multiple human individuals, within and between populations, and for a variety of tissue types. Comparison of TE insertion site profiles between individuals captures the germline activity of TEs and reveals insertion site variants that segregate as polymorphisms among human populations, whereas comparison among tissue types ascertains somatic TE activity that generates cellular heterogeneity. In this review, we provide an overview of these new technologies and explore their implications for population and clinical genetic studies of human TEs. We cover both recent published results on human TE insertion activity as well as the prospects for future TE studies related to human evolution and health.

Profiling Transposable Elements and Their Epigenetic Effects in Non-model Species

Taking transposable elements into consideration in surveys of genetic and epigenetic variation remains challenging in species lacking a high-quality reference genome. Here, molecular techniques reducing genome complexity and specifically targeting restructuring and methylation changes in TE genome fractions are described. In particular, methyl-sensitive transposon display (MSTD) uses isoschizomers and PCR amplifications to assess the methylation environment of TE insertions. MSTD offers reliable insights into genome-wide epigenetic changes associated with TEs, especially when used together with similar techniques tracking random sequences.

Reverse Genetics and High Throughput Sequencing Methodologies for Plant Functional Genomics

In the post-genomic era, increasingly sophisticated genetic tools are being developed with the long-term goal of understanding how the coordinated activity of genes gives rise to a complex organism. With the advent of the next generation sequencing associated with effective computational approaches, wide variety of plant species have been fully sequenced giving a wealth of data sequence information on structure and organization of plant genomes. Since thousands of gene sequences are already known, recently developed functional genomics approaches provide powerful tools to analyze plant gene functions through various gene manipulation technologies. Integration of different omics platforms along with gene annotation and computational analysis may elucidate a complete view in a system biology level. Extensive investigations on reverse genetics methodologies were deployed for assigning biological function to a specific gene or gene product. We provide here an updated overview of these high throughout strategies highlighting recent advances in the knowledge of functional genomics in plants.

Genome-Wide Comparison of Magnaporthe Species Reveals a Host-Specific Pattern of Secretory Proteins and Transposable Elements

Blast disease caused by the Magnaporthe species is a major factor affecting the productivity of rice, wheat and millets. This study was aimed at generating genomic information for rice and non-rice Magnaporthe isolates to understand the extent of genetic variation. We have sequenced the whole genome of the Magnaporthe isolates, infecting rice (leaf and neck), finger millet (leaf and neck), foxtail millet (leaf) and buffel grass (leaf). Rice and finger millet isolates infecting both leaf and neck tissues were sequenced, since the damage and yield loss caused due to neck blast is much higher as compared to leaf blast. The genome-wide comparison was carried out to study the variability in gene content, candidate effectors, repeat element distribution, genes involved in carbohydrate metabolism and SNPs. The analysis of repeat element footprints revealed some genes such as naringenin, 2-oxoglutarate 3-dioxygenase being targeted by Pot2 and Occan, in isolates from different host species. Some repeat insertions were host-specific while other insertions were randomly shared between isolates. The distributions of repeat elements, secretory proteins, CAZymes and SNPs showed significant variation across host-specific lineages of Magnaporthe indicating an independent genome evolution orchestrated by multiple genomic factors.

Exploration of the Drosophila buzzatii transposable element content suggests underestimation of repeats in Drosophila genomes

Background

Many new Drosophila genomes have been sequenced in recent years using new-generation sequencing platforms and assembly methods. Transposable elements (TEs), being repetitive sequences, are often misassembled, especially in the genomes sequenced with short reads. Consequently, the mobile fraction of many of the new genomes has not been analyzed in detail or compared with that of other genomes sequenced with different methods, which could shed light into the understanding of genome and TE evolution. Here we compare the TE content of three genomes: D. buzzatii st-1, j-19, and D. mojavensis.

Results

We have sequenced a new D. buzzatii genome (j-19) that complements the D. buzzatii reference genome (st-1) already published, and compared their TE contents with that of D. mojavensis. We found an underestimation of TE sequences in Drosophila genus NGS-genomes when compared to Sanger-genomes. To be able to compare genomes sequenced with different technologies, we developed a coverage-based method and applied it to the D. buzzatii st-1 and j-19 genome. Between 10.85 and 11.16 % of the D. buzzatii st-1 genome is made up of TEs, between 7 and 7,5 % of D. buzzatii j-19 genome, while TEs represent 15.35 % of the D. mojavensis genome. Helitrons are the most abundant order in the three genomes.

Conclusions

TEs in D. buzzatii are less abundant than in D. mojavensis, as expected according to the genome size and TE content positive correlation. However, TEs alone do not explain the genome size difference. TEs accumulate in the dot chromosomes and proximal regions of D. buzzatii and D. mojavensis chromosomes. We also report a significantly higher TE density in D. buzzatii and D. mojavensis X chromosomes, which is not expected under the current models. Our easy-to-use correction method allowed us to identify recently active families in D. buzzatii st-1 belonging to the LTR-retrotransposon superfamily Gypsy.

Electronic supplementary material

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

Rapid identification of causative insertions underlying Medicago truncatula Tnt1 mutants defective in symbiotic nitrogen fixation from a forward genetic screen by whole genome sequencing

BACKGROUND

In the model legume Medicago truncatula, the near saturation genome-wide Tnt1 insertion mutant population in ecotype R108 is a valuable tool in functional genomics studies. Forward genetic screens have identified many Tnt1 mutants defective in nodule development and symbiotic nitrogen fixation (SNF). However, progress toward identifying the causative mutations of these symbiotic mutants has been slow because of the high copy number of Tnt1 insertions in some mutant plants and inefficient recovery of flanking sequence tags (FSTs) by thermal asymmetric interlaced PCR (TAIL-PCR) and other techniques.

RESULTS

Two Tnt1 symbiotic mutants, NF11217 and NF10547, with defects in nodulation and SNF were isolated during a forward genetic screen. Both TAIL-PCR and whole genome sequencing (WGS) approaches were used in attempts to find the relevant mutant genes in NF11217 and NF10547. Illumina paired-end WGS generated ~16 Gb of sequence data from a 500 bp insert library for each mutant, yielding ~40X genome coverage. Bioinformatics analysis of the sequence data identified 97 and 65 high confidence independent Tnt1 insertion loci in NF11217 and NF10547, respectively. In comparison to TAIL-PCR, WGS recovered more Tnt1 insertions. From the WGS data, we found Tnt1 insertions in the exons of the previously described PHOSPHOLIPASE C (PLC)-like and NODULE INCEPTION (NIN) genes in NF11217 and NF10547 mutants, respectively. Co-segregation analyses confirmed that the symbiotic phenotypes of NF11217 and NF10547 are tightly linked to the Tnt1 insertions in PLC-like and NIN genes, respectively.

CONCLUSIONS

In this work, we demonstrate that WGS is an efficient approach for identification of causative genes underlying SNF defective phenotypes in M. truncatula Tnt1 insertion mutants obtained via forward genetic screens.

Rapid identification of causative insertions underlying Medicago truncatula Tnt1 mutants defective in symbiotic nitrogen fixation from a forward genetic screen by whole genome sequencing

BACKGROUND

In the model legume Medicago truncatula, the near saturation genome-wide Tnt1 insertion mutant population in ecotype R108 is a valuable tool in functional genomics studies. Forward genetic screens have identified many Tnt1 mutants defective in nodule development and symbiotic nitrogen fixation (SNF). However, progress toward identifying the causative mutations of these symbiotic mutants has been slow because of the high copy number of Tnt1 insertions in some mutant plants and inefficient recovery of flanking sequence tags (FSTs) by thermal asymmetric interlaced PCR (TAIL-PCR) and other techniques.

RESULTS

Two Tnt1 symbiotic mutants, NF11217 and NF10547, with defects in nodulation and SNF were isolated during a forward genetic screen. Both TAIL-PCR and whole genome sequencing (WGS) approaches were used in attempts to find the relevant mutant genes in NF11217 and NF10547. Illumina paired-end WGS generated ~16 Gb of sequence data from a 500 bp insert library for each mutant, yielding ~40X genome coverage. Bioinformatics analysis of the sequence data identified 97 and 65 high confidence independent Tnt1 insertion loci in NF11217 and NF10547, respectively. In comparison to TAIL-PCR, WGS recovered more Tnt1 insertions. From the WGS data, we found Tnt1 insertions in the exons of the previously described PHOSPHOLIPASE C (PLC)-like and NODULE INCEPTION (NIN) genes in NF11217 and NF10547 mutants, respectively. Co-segregation analyses confirmed that the symbiotic phenotypes of NF11217 and NF10547 are tightly linked to the Tnt1 insertions in PLC-like and NIN genes, respectively.

CONCLUSIONS

In this work, we demonstrate that WGS is an efficient approach for identification of causative genes underlying SNF defective phenotypes in M. truncatula Tnt1 insertion mutants obtained via forward genetic screens.

Genome ARTIST: a robust, high-accuracy aligner tool for mapping transposon insertions and self-insertions

Background

A critical topic of insertional mutagenesis experiments performed on model organisms is mapping the hits of artificial transposons (ATs) at nucleotide level accuracy. Mapping errors may occur when sequencing artifacts or mutations as single nucleotide polymorphisms (SNPs) and small indels are present very close to the junction between a genomic sequence and a transposon inverted repeat (TIR). Another particular item of insertional mutagenesis is mapping of the transposon self-insertions and, to our best knowledge, there is no publicly available mapping tool designed to analyze such molecular events.

Results

We developed Genome ARTIST, a pairwise gapped aligner tool which works out both issues by means of an original, robust mapping strategy. Genome ARTIST is not designed to use next-generation sequencing (NGS) data but to analyze ATs insertions obtained in small to medium-scale mutagenesis experiments. Genome ARTIST employs a heuristic approach to find DNA sequence similarities and harnesses a multi-step implementation of a Smith-Waterman adapted algorithm to compute the mapping alignments. The experience is enhanced by easily customizable parameters and a user-friendly interface that describes the genomic landscape surrounding the insertion. Genome ARTIST is functional with many genomes of bacteria and eukaryotes available in Ensembl and GenBank repositories. Our tool specifically harnesses the sequence annotation data provided by FlyBase for Drosophila melanogaster (the fruit fly), which enables mapping of insertions relative to various genomic features such as natural transposons. Genome ARTIST was tested against other alignment tools using relevant query sequences derived from the D. melanogaster and Mus musculus (mouse) genomes. Real and simulated query sequences were also comparatively inquired, revealing that Genome ARTIST is a very robust solution for mapping transposon insertions.

Conclusions

Genome ARTIST is a stand-alone user-friendly application, designed for high-accuracy mapping of transposon insertions and self-insertions. The tool is also useful for routine aligning assessments like detection of SNPs or checking the specificity of primers and probes. Genome ARTIST is an open source software and is available for download at www.genomeartist.ro and at GitHub (https://github.com/genomeartist/genomeartist ).

Electronic supplementary material

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

A Snapshot of Functional Genetic Studies in Medicago truncatula

In the current context of food security, increase of plant protein production in a sustainable manner represents one of the major challenges of agronomic research, which could be partially resolved by increased cultivation of legume crops. Medicago truncatula is now a well-established model for legume genomic and genetic studies. With the establishment of genomics tools and mutant populations in M. truncatula, it has become an important resource to answer some of the basic biological questions related to plant development and stress tolerance. This review has an objective to overview a decade of genetic studies in this model plant from generation of mutant populations to nowadays. To date, the three biological fields, which have been extensively studied in M. truncatula, are the symbiotic nitrogen fixation, the seed development, and the abiotic stress tolerance, due to their significant agronomic impacts. In this review, we summarize functional genetic studies related to these three major biological fields. We integrated analyses of a nearly exhaustive list of genes into their biological contexts in order to provide an overview of the forefront research advances in this important legume model plant.

Transposable element detection from whole genome sequence data

The number of software tools available for detecting transposable element insertions from whole genome sequence data has been increasing steadily throughout the last ~5 years. Some of these methods have unique features suiting them for particular use cases, but in general they follow one or more of a common set of approaches. Here, detection and filtering approaches are reviewed in the light of transposable element biology and the current state of whole genome sequencing. We demonstrate that the current state-of-the-art methods still do not produce highly concordant results and provide resources to assist future development in transposable element detection methods.

CONTRAILS: A tool for rapid identification of transgene integration sites in complex, repetitive genomes using low-coverage paired-end sequencing

Transgenic crops have become a staple in modern agriculture, and are typically characterized using a variety of molecular techniques involving proteomics and metabolomics. Characterization of the transgene insertion site is of great interest, as disruptions, deletions, and genomic location can affect product selection and fitness, and identification of these regions and their integrity is required for regulatory agencies. Here, we present CONTRAILS (Characterization of Transgene Insertion Locations with Sequencing), a straightforward, rapid and reproducible method for the identification of transgene insertion sites in highly complex and repetitive genomes using low coverage paired-end Illumina sequencing and traditional PCR. This pipeline requires little to no troubleshooting and is not restricted to any genome type, allowing use for many molecular applications. Using whole genome sequencing of in-house transgenic Glycine max, a legume with a highly repetitive and complex genome, we used CONTRAILS to successfully identify the location of a single T-DNA insertion to single base resolution.

•

We developed a pipeline for transgene identification using paired-end sequencing.

•

This method requires little troubleshooting and is applicable to any genome.

•

Identification of insertion sites is required for deregulation of modified food crops.

•

This assists in identifying potential genomic disruptions in transgenic events.