Mobster: accurate detection of mobile element insertions in next generation sequencing data

Investigating mobile element variations by statistical genetics

The integration of structural variations (SVs) in statistical genetics provides an opportunity to understand the genetic factors influencing complex human traits and disease. Recent advances in long-read technology and variant calling methods for short reads have improved the accurate discovery and genotyping of SVs, enabling their use in expression quantitative trait loci (eQTL) analysis and genome-wide association studies (GWAS). Mobile elements are DNA sequences that insert themselves into various genome locations. Insertional polymorphisms of mobile elements between humans, called mobile element variations (MEVs), contribute to approximately 25% of human SVs. We recently developed a variant caller that can accurately identify and genotype MEVs from biobank-scale short-read whole-genome sequencing (WGS) datasets and integrate them into statistical genetics. The use of MEVs in eQTL analysis and GWAS has a minimal impact on the discovery of genome loci associated with gene expression and disease; most disease-associated haplotypes can be identified by single nucleotide variations (SNVs). On the other hand, it helps make hypotheses about causal variants or effector variants. Focusing on MEVs, we identified multiple MEVs that contribute to differential gene expression and one of them is a potential cause of skin disease, emphasizing the importance of the integration of MEVs in medical genetics. Here, I will provide an overview of MEVs, MEV calling from WGS, and the integration of MEVs in statistical genetics. Finally, I will discuss the unanswered questions about MEVs, such as rare variants.

Contribution of de novo retroelements to birth defects and childhood cancers

Insertion of active retroelements-L1s, Alus, and SVAs-can disrupt proper genome function and lead to various disorders including cancer. However, the role of de novo retroelements (DNRTs) in birth defects and childhood cancers has not been well characterized due to the lack of adequate data and efficient computational tools. Here, we examine whole-genome sequencing data of 3,244 trios from 12 birth defect and childhood cancer cohorts in the Gabriella Miller Kids First Pediatric Research Program. Using an improved version of our tool xTea (x-Transposable element analyzer) that incorporates a deep-learning module, we identified 162 DNRTs, as well as 2 pseudogene insertions. Several variants are likely to be causal, such as a de novo Alu insertion that led to the ablation of a whole exon in the NF1 gene in a proband with brain tumor. We observe a high de novo SVA insertion burden in both high-intolerance loss-of-function genes and exons as well as more frequent de novo Alu insertions of paternal origin. We also identify potential mosaic DNRTs from embryonic stages. Our study reveals the important roles of DNRTs in causing birth defects and predisposition to childhood cancers.

The role of transposable elements in human evolution and methods for their functional analysis: current status and future perspectives

Transposable elements (TEs) are mobile DNA sequences that can insert themselves into various locations within the genome, causing mutations that may provide advantages or disadvantages to individuals and species. The insertion of TEs can result in genetic variation that may affect a wide range of human traits including genetic disorders. Understanding the role of TEs in human biology is crucial for both evolutionary and medical research. This review discusses the involvement of TEs in human traits and disease susceptibility, as well as methods for functional analysis of TEs.

Mobile element insertions in rare diseases: a comparative benchmark and reanalysis of 60,000 exome samples

Mobile element insertions (MEIs) are a known cause of genetic disease but have been underexplored due to technical limitations of genetic testing methods. Various bioinformatic tools have been developed to identify MEIs in Next Generation Sequencing data. However, most tools have been developed specifically for genome sequencing (GS) data rather than exome sequencing (ES) data, which remains more widely used for routine diagnostic testing. In this study, we benchmarked six MEI detection tools (ERVcaller, MELT, Mobster, SCRAMble, TEMP2 and xTea) on ES data and on GS data from publicly available genomic samples (HG002, NA12878). For all the tools we evaluated sensitivity and precision of different filtering strategies. Results show that there were substantial differences in tool performance between ES and GS data. MELT performed best with ES data and its combination with SCRAMble increased substantially the detection rate of MEIs. By applying both tools to 10,890 ES samples from Solve-RD and 52,624 samples from Radboudumc we were able to diagnose 10 patients who had remained undiagnosed by conventional ES analysis until now. Our study shows that MELT and SCRAMble can be used reliably to identify clinically relevant MEIs in ES data. This may lead to an additional diagnosis for 1 in 3000 to 4000 patients in routine clinical ES.

Towards targeting transposable elements for cancer therapy

Transposable elements (TEs) represent almost half of the human genome. Historically deemed 'junk DNA', recent technological advancements have stimulated a wave of research into the functional impact of TEs on gene-regulatory networks in evolution and development, as well as in diseases including cancer. The genetic and epigenetic evolution of cancer involves the exploitation of TEs, whereby TEs contribute directly to cancer-specific gene activities. This Review provides a perspective on the role of TEs in cancer as being a 'double-edged sword', both promoting cancer evolution and representing a vulnerability that could be exploited in cancer therapy. We discuss how TEs affect transcriptome regulation and other cellular processes in cancer. We highlight the potential of TEs as therapeutic targets for cancer. We also summarize technical hurdles in the characterization of TEs with genomic assays. Last, we outline open questions and exciting future research avenues.

Emerging Opportunities to Study Mobile Element Insertions and Their Source Elements in an Expanding Universe of Sequenced Human Genomes

Three mobile element classes, namely Alu, LINE-1 (L1), and SVA elements, remain actively mobile in human genomes and continue to produce new mobile element insertions (MEIs). Historically, MEIs have been discovered and studied using several methods, including: (1) Southern blots, (2) PCR (including PCR display), and (3) the detection of MEI copies from young subfamilies. We are now entering a new phase of MEI discovery where these methods are being replaced by whole genome sequencing and bioinformatics analysis to discover novel MEIs. We expect that the universe of sequenced human genomes will continue to expand rapidly over the next several years, both with short-read and long-read technologies. These resources will provide unprecedented opportunities to discover MEIs and study their impact on human traits and diseases. They also will allow the MEI community to discover and study the source elements that produce these new MEIs, which will facilitate our ability to study source element regulation in various tissue contexts and disease states. This, in turn, will allow us to better understand MEI mutagenesis in humans and the impact of this mutagenesis on human biology.

A second look at exome sequencing data: detecting mobile elements insertion in a rare disease cohort

About 0.3% of all variants are due to de novo mobile element insertions (MEIs). The massive development of next-generation sequencing has made it possible to identify MEIs on a large scale. We analyzed exome sequencing (ES) data from 3232 individuals (2410 probands) with developmental and/or neurological abnormalities, with MELT, a tool designed to identify MEIs. The results were filtered by frequency, impacted region and gene function. Following phenotype comparison, two candidates were identified in two unrelated probands. The first mobile element (ME) was found in a patient referred for poikilodermia. A homozygous insertion was identified in the FERMT1 gene involved in Kindler syndrome. RNA study confirmed its pathological impact on splicing. The second ME was a de novo Alu insertion in the GRIN2B gene involved in intellectual disability, and detected in a patient with a developmental disorder. The frequency of de novo exonic MEIs in our study is concordant with previous studies on ES data. This project, which aimed to identify pathological MEIs in the coding sequence of genes, confirms that including detection of MEs in the ES pipeline can increase the diagnostic rate. This work provides additional evidence that ES could be used alone as a diagnostic exam.

Detection of trait-associated structural variations using short-read sequencing

Genomic structural variation (SV) affects genetic and phenotypic characteristics in diverse organisms, but the lack of reliable methods to detect SV has hindered genetic analysis. We developed a computational algorithm (MOPline) that includes missing call recovery combined with high-confidence SV call selection and genotyping using short-read whole-genome sequencing (WGS) data. Using 3,672 high-coverage WGS datasets, MOPline stably detected ∼16,000 SVs per individual, which is over ∼1.7-3.3-fold higher than previous large-scale projects while exhibiting a comparable level of statistical quality metrics. We imputed SVs from 181,622 Japanese individuals for 42 diseases and 60 quantitative traits. A genome-wide association study with the imputed SVs revealed 41 top-ranked or nearly top-ranked genome-wide significant SVs, including 8 exonic SVs with 5 novel associations and enriched mobile element insertions. This study demonstrates that short-read WGS data can be used to identify rare and common SVs associated with a variety of traits.

INSurVeyor: improving insertion calling from short read sequencing data

Insertions are one of the major types of structural variations and are defined as the addition of 50 nucleotides or more into a DNA sequence. Several methods exist to detect insertions from next-generation sequencing short read data, but they generally have low sensitivity. Our contribution is two-fold. First, we introduce INSurVeyor, a fast, sensitive and precise method that detects insertions from next-generation sequencing paired-end data. Using publicly available benchmark datasets (both human and non-human), we show that INSurVeyor is not only more sensitive than any individual caller we tested, but also more sensitive than all of them combined. Furthermore, for most types of insertions, INSurVeyor is almost as sensitive as long reads callers. Second, we provide state-of-the-art catalogues of insertions for 1047 Arabidopsis Thaliana genomes from the 1001 Genomes Project and 3202 human genomes from the 1000 Genomes Project, both generated with INSurVeyor. We show that they are more complete and precise than existing resources, and important insertions are missed by existing methods.

WGS Revealed Novel BBS5 Pathogenic Variants, Missed by WES, Causing Ciliary Structure and Function Defects

Bardet-Biedl syndrome (BBS) is an autosomal recessive ciliopathy that affects multiple organs, leading to retinitis pigmentosa, polydactyly, obesity, renal anomalies, cognitive impairment, and hypogonadism. Until now, biallelic pathogenic variants have been identified in at least 24 genes delineating the genetic heterogeneity of BBS. Among those, BBS5 is a minor contributor to the mutation load and is one of the eight subunits forming the BBSome, a protein complex implied in protein trafficking within the cilia. This study reports on a European BBS5 patient with a severe BBS phenotype. Genetic analysis was performed using multiple next-generation sequencing (NGS) tests (targeted exome, TES and whole exome, WES), and biallelic pathogenic variants could only be identified using whole-genome sequencing (WGS), including a previously missed large deletion of the first exons. Despite the absence of family samples, the biallelic status of the variants was confirmed. The BBS5 protein's impact was confirmed on the patient's cells (presence/absence and size of the cilium) and ciliary function (Sonic Hedgehog pathway). This study highlights the importance of WGS and the challenge of reliable structural variant detection in patients' genetic explorations as well as functional tests to assess a variant's pathogenicity.

A Novel Alu Element Insertion in ATM Induces Exon Skipping in Suspected HBOC Patients

The vast majority of patients at risk of hereditary breast and/or ovarian cancer (HBOC) syndrome remain without a molecular diagnosis after routine genetic testing. One type of genomic alteration that is commonly missed by diagnostic pipelines is mobile element insertions (MEIs). Here, we reanalyzed multigene panel data from suspected HBOC patients using the MEI detection tool Mobster. A novel Alu element insertion in ATM intron 54 (ATM:c.8010+30_8010+31insAluYa5) was identified as a potential contributing factor in seven patients. Transcript analysis of patient-derived RNA from three heterozygous carriers revealed exon 54 skipping in 38% of total ATM transcripts. To manifest the direct association between the Alu element insertion and the aberrant splice pattern, HEK293T and MCF7 cells were transfected with wild-type or Alu element-carrying minigene constructs. On average, 77% of plasmid-derived transcripts lacked exon 54 in the presence of the Alu element insertion compared to only 4.7% of transcripts expressed by the wild-type minigene. These results strongly suggest ATM:c.8010+30_8010+31insAluYa5 as the main driver of ATM exon 54 skipping. Since this exon loss is predicted to cause a frameshift and a premature stop codon, mutant transcripts are unlikely to translate into functional proteins. Based on its estimated frequency of up to 0.05% in control populations, we propose to consider ATM:c.8010+30_8010+31insAluYa5 in suspected HBOC patients and to clarify its role in carcinogenesis through future epidemiological and functional analyses. Generally, the implementation of MEI detection tools in diagnostic sequencing pipelines could increase the diagnostic yield, as MEIs are likely underestimated contributors to genetic diseases.

Whole Genome Analysis of Dizygotic Twins With Autism Reveals Prevalent Transposon Insertion Within Neuronal Regulatory Elements: Potential Implications for Disease Etiology and Clinical Assessment

Transposable elements (TEs) have been implicated in autism spectrum disorder (ASD). However, our understanding of their roles is far from complete. Herein, we explored de novo TE insertions (dnTEIs) and de novo variants (DNVs) across the genomes of dizygotic twins with ASD and their parents. The neuronal regulatory elements had a tendency to harbor dnTEIs that were shared between twins, but ASD-risk genes had dnTEIs that were unique to each twin. The dnTEIs were 4.6-fold enriched in enhancers that are active in embryonic stem cell (ESC)-neurons (p < 0.001), but DNVs were 1.5-fold enriched in active enhancers of astrocytes (p = 0.0051). Our findings suggest that dnTEIs and DNVs play a role in ASD etiology by disrupting enhancers of neurons and astrocytes.

Retrotransposon insertion as a novel mutational cause of spinal muscular atrophy

Spinal muscular atrophy (SMA) is an autosomal recessive neuromuscular disorder resulting from biallelic alterations of the SMN1 gene: deletion, gene conversion or, in rare cases, intragenic variants. The disease severity is mainly influenced by the copy number of SMN2, a nearly identical gene, which produces only low amounts of full-length (FL) mRNA. Here we describe the first example of retrotransposon insertion as a pathogenic SMN1 mutational event. The 50-year-old patient is clinically affected by SMA type III with a diagnostic odyssey spanning nearly 30 years. Despite a mild disease course, he carries a single SMN2 copy. Using Exome Sequencing and Sanger sequencing, we characterized a SINE-VNTR-Alu (SVA) type F retrotransposon inserted in SMN1 intron 7. Using RT-PCR and RNASeq experiments on lymphoblastoid cell lines, we documented the dramatic decrease of FL transcript production in the patient compared to subjects with the same SMN1 and SMN2 copy number, thus validating the pathogenicity of this SVA insertion. We described the mutant FL-SMN1-SVA transcript characterized by exon extension and showed that it is subject to degradation by nonsense-mediated mRNA decay. The stability of the SMN-SVA protein may explain the mild course of the disease. This observation exemplifies the role of retrotransposons in human genetic disorders.

Cognitive Impairment in a Complex Family With AAGGG and ACAGG Repeat Expansions in RFC1 Detected by ExpansionHunter Denovo

Background and Objectives

We investigated the genetic basis and brain metabolism and blood flow of a Japanese family with spinocerebellar degeneration (SCD), with multiple affected members for 3 generations.

Methods

After excluding DNA repeat expansion (RE) of common SCD genes by fragment analysis, we performed whole-exome sequencing (WES) and whole-genome sequencing (WGS). Homozygosity mapping was performed using these data. REs were investigated with WGS data using ExpansionHunter Denovo and Expansion Hunter.

Results

WES and WGS were unable to identify likely pathogenic variants, and homozygosity mapping failed to narrow down the locus. However, ExpansionHunter Denovo detected REs in intron 2 of the RFC1 gene and led us to the diagnosis of RFC1-related disorders. Subsequent repeat-primed PCR and Southern blot hybridization analyses revealed that 3 of 6 patients and 1 suspected individual had expansions of AAGGG ((AAGGG)_exp) and (ACAGG)_exp repeats in a compound heterozygous state and 3 had a homozygous (ACAGG)_exp. The patients showed a variety of clinical features, including adult-onset ataxia, sensorimotor neuropathy, head tremor, parkinsonism, dystonia, and cognitive impairment. A comparison of previous reports with those of the family in study suggested that motor neuropathy could be a feature of compound heterozygous patients and biallelic (ACAGG)_exp patients. Cognitive function tests showed cognitive impairment with a predominance of frontal lobe dysfunction. Examination of MRI, SPECT, and ¹⁸F-fluorodeoxyglucose-PET showed clear cortical damage with frontal lobe predominance in 1 case, but no cerebral damage was evident in the other 2 cases.

Discussion

Our report shows the usefulness of WGS and RE detection tools for SCD of unknown cause. The studied family with RFC1-related disorders included patients with (ACAGG)_exp and (AAGGG)_exp in a compound heterozygous state and was characterized by motor neuropathy. Based on the results of cognitive function tests and imaging studies, 1 patient presented with cognitive impairment due to frontal lobe metabolic changes, but there were also patients who presented with cognitive impairment without apparent cerebral metabolic or blood flow, suggesting that other factors are also associated with cognitive impairment.

Methods to Improve Molecular Diagnosis in Genomic Cold Cases in Pediatric Neurology

During the last decade, genetic testing has emerged as an important etiological diagnostic tool for Mendelian diseases, including pediatric neurological conditions. A genetic diagnosis has a considerable impact on disease management and treatment; however, many cases remain undiagnosed after applying standard diagnostic sequencing techniques. This review discusses various methods to improve the molecular diagnostic rates in these genomic cold cases. We discuss extended analysis methods to consider, non-Mendelian inheritance models, mosaicism, dual/multiple diagnoses, periodic re-analysis, artificial intelligence tools, and deep phenotyping, in addition to integrating various omics methods to improve variant prioritization. Last, novel genomic technologies, including long-read sequencing, artificial long-read sequencing, and optical genome mapping are discussed. In conclusion, a more comprehensive molecular analysis and a timely re-analysis of unsolved cases are imperative to improve diagnostic rates. In addition, our current understanding of the human genome is still limited due to restrictions in technologies. Novel technologies are now available that improve upon some of these limitations and can capture all human genomic variation more accurately. Last, we recommend a more routine implementation of high molecular weight DNA extraction methods that is coherent with the ability to use and/or optimally benefit from these novel genomic methods.

The insertion and dysregulation of transposable elements in osteosarcoma and their association with patient event-free survival

The dysregulation of transposable elements (TEs) has been explored in a variety of cancers. However, TE activities in osteosarcoma (OS) have not been extensively studied yet. By integrative analysis of RNA-seq, whole-genome sequencing (WGS), and methylation data, we showed aberrant TE activities associated with dysregulations of TEs in OS tumors. Specifically, expression levels of LINE-1 and Alu of different evolutionary ages, as well as subfamilies of SVA and HERV-K, were significantly up-regulated in OS tumors, accompanied by enhanced DNA repair responses. We verified the characteristics of LINE-1 mediated TE insertions, including target site duplication (TSD) length (centered around 15 bp) and preferential insertions into intergenic and AT-rich regions as well as intronic regions of longer genes. By filtering polymorphic TE insertions reported in 1000 genome project (1KGP), besides 148 tumor-specific somatic TE insertions, we found most OS patient-specific TE insertions (3175 out of 3326) are germline insertions, which are associated with genes involved in neuronal processes or with transcription factors important for cancer development. In addition to 68 TE-affected cancer genes, we found recurrent germline TE insertions in 72 non-cancer genes with high frequencies among patients. We also found that +/− 500 bps flanking regions of transcription start sites (TSS) of LINE-1 (young) and Alu showed lower methylation levels in OS tumor samples than controls. Interestingly, by incorporating patient clinical data and focusing on TE activities in OS tumors, our data analysis suggested that higher TE insertions in OS tumors are associated with a longer event-free survival time.

Targeted massively parallel sequencing of candidate regions on chromosome 22q predisposing to multiple schwannomas: An analysis of 51 individuals in a single-center experience

Constitutional LZTR1 or SMARCB1 pathogenic variants (PVs) have been found in ∼86% of familial and ∼40% of sporadic schwannomatosis cases. Hence, we performed massively parallel sequencing of the entire LZTR1, SMARCB1, and NF2 genomic loci in 35 individuals with schwannomas negative for constitutional first-hit PVs in the LZTR1/SMARCB1/NF2 coding sequences; however, with 22q deletion and/or a different NF2 PV in each tumor, including six cases with only one tumor available. Furthermore, we verified whether any other LZTR1/SMARCB1/NF2 (likely) PVs could be found in 16 cases carrying a SMARCB1 constitutional variant in the 3'-untranslated region (3'-UTR) c.*17C>T, c.*70C>T, or c.*82C>T. As no additional variants were found, functional studies were performed to clarify the effect of these 3'-UTR variants on the transcript. The 3'-UTR variants c.*17C>T and c.*82C>T showed pathogenicity by negatively affecting the SMARCB1 transcript level. Two novel deep intronic SMARCB1 variants, c.500+883T>G and c.500+887G>A, resulting in out-of-frame missplicing of intron 4, were identified in two unrelated individuals. Further resequencing of the entire repeat-masked genomics sequences of chromosome 22q in individuals negative for PVs in the SMARCB1/LZTR1/NF2 coding- and noncoding regions revealed five potential schwannomatosis-predisposing candidate genes, that is, MYO18B, NEFH, SGSM1, SGSM3, and SBF1, pending further verification.

An assessment of bioinformatics tools for the detection of human endogenous retroviral insertions in short-read genome sequencing data

There is a growing interest in the study of human endogenous retroviruses (HERVs) given the substantial body of evidence that implicates them in many human diseases. Although their genomic characterization presents numerous technical challenges, next-generation sequencing (NGS) has shown potential to detect HERV insertions and their polymorphisms in humans. Currently, a number of computational tools to detect them in short-read NGS data exist. In order to design optimal analysis pipelines, an independent evaluation of the available tools is required. We evaluated the performance of a set of such tools using a variety of experimental designs and datasets. These included 50 human short-read whole-genome sequencing samples, matching long and short-read sequencing data, and simulated short-read NGS data. Our results highlight a great performance variability of the tools across the datasets and suggest that different tools might be suitable for different study designs. However, specialized tools designed to detect exclusively human endogenous retroviruses consistently outperformed generalist tools that detect a wider range of transposable elements. We suggest that, if sufficient computing resources are available, using multiple HERV detection tools to obtain a consensus set of insertion loci may be ideal. Furthermore, given that the false positive discovery rate of the tools varied between 8% and 55% across tools and datasets, we recommend the wet lab validation of predicted insertions if DNA samples are available.

Variant calling: Considerations, practices, and developments

The success of many clinical, association, or population genetics studies critically relies on properly performed variant calling step. The variety of modern genomics protocols, techniques, and platforms makes our choices of methods and algorithms difficult and there is no "one size fits all" solution for study design and data analysis. In this review, we discuss considerations that need to be taken into account while designing the study and preparing for the experiments. We outline the variety of variant types that can be detected using sequencing approaches and highlight some specific requirements and basic principles of their detection. Finally, we cover interesting developments that enable variant calling for a broad range of applications in the genomics field. We conclude by discussing technological and algorithmic advances that have the potential to change the ways of calling DNA variants in the nearest future.

Whole-genome analysis reveals the contribution of non-coding de novo transposon insertions to autism spectrum disorder

Background

Retrotransposons have been implicated as causes of Mendelian disease, but their role in autism spectrum disorder (ASD) has not been systematically defined, because they are only called with adequate sensitivity from whole genome sequencing (WGS) data and a large enough cohort for this analysis has only recently become available.

Results

We analyzed WGS data from a cohort of 2288 ASD families from the Simons Simplex Collection by establishing a scalable computational pipeline for retrotransposon insertion detection. We report 86,154 polymorphic retrotransposon insertions—including > 60% not previously reported—and 158 de novo retrotransposition events. The overall burden of de novo events was similar between ASD individuals and unaffected siblings, with 1 de novo insertion per 29, 117, and 206 births for Alu, L1, and SVA respectively, and 1 de novo insertion per 21 births total. However, ASD cases showed more de novo L1 insertions than expected in ASD genes. Additionally, we observed exonic insertions in loss-of-function intolerant genes, including a likely pathogenic exonic insertion in CSDE1, only in ASD individuals.

Conclusions

These findings suggest a modest, but important, impact of intronic and exonic retrotransposon insertions in ASD, show the importance of WGS for their analysis, and highlight the utility of specific bioinformatic tools for high-throughput detection of retrotransposon insertions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-021-00256-w.

Transposable elements that have recently been mobile in the human genome

Background

Transposable elements (TE) comprise nearly half of the human genome and their insertions have profound effects to human genetic diversification and as well as disease. Despite their abovementioned significance, there is no consensus on the TE subfamilies that remain active in the human genome. In this study, we therefore developed a novel statistical test for recently mobile subfamilies (RMSs), based on patterns of overlap with > 100,000 polymorphic indels.

Results

Our analysis produced a catalogue of 20 high-confidence RMSs, which excludes many false positives in public databases. Intriguingly though, it includes HERV-K, an LTR subfamily previously thought to be extinct. The RMS catalogue is strongly enriched for contributions to germline genetic disorders (P = 1.1e-10), and thus constitutes a valuable resource for diagnosing disorders of unknown aetiology using targeted TE-insertion screens. Remarkably, RMSs are also highly enriched for somatic insertions in diverse cancers (P = 2.8e-17), thus indicating strong correlations between germline and somatic TE mobility. Using CRISPR/Cas9 deletion, we show that an RMS-derived polymorphic TE insertion increased the expression of RPL17, a gene associated with lower survival in liver cancer. More broadly, polymorphic TE insertions from RMSs were enriched near genes with allele-specific expression, suggesting widespread effects on gene regulation.

Conclusions

By using a novel statistical test we have defined a catalogue of 20 recently mobile transposable element subfamilies. We illustrate the gene regulatory potential of RMS-derived polymorphic TE insertions, using CRISPR/Cas9 deletion in vitro on a specific candidate, as well as by genome wide analysis of allele-specific expression. Our study presents novel insights into TE mobility and regulatory potential and provides a key resource for human disease genetics and population history studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08085-0.

Optical genome mapping identifies a germline retrotransposon insertion in SMARCB1 in two siblings with atypical teratoid rhabdoid tumors

In a subset of pediatric cancers, a germline cancer predisposition is highly suspected based on clinical and pathological findings, but genetic evidence is lacking, which hampers genetic counseling and predictive testing in the families involved. We describe a family with two siblings born from healthy parents who were both neonatally diagnosed with atypical teratoid rhabdoid tumor (ATRT). This rare and aggressive pediatric tumor is associated with biallelic inactivation of SMARCB1, and in 30% of the cases, a predisposing germline mutation is involved. Whereas the tumors of both siblings showed loss of expression of SMARCB1 and acquired homozygosity of the locus, whole exome and whole genome sequencing failed to identify germline or somatic SMARCB1 pathogenic mutations. We therefore hypothesized that the insertion of a pathogenic repeat‐rich structure might hamper its detection, and we performed optical genome mapping (OGM) as an alternative strategy to identify structural variation in this locus. Using this approach, an insertion of ~2.8 kb within intron 2 of SMARCB1 was detected. Long‐range PCR covering this region remained unsuccessful, but PacBio HiFi genome sequencing identified this insertion to be a SINE‐VNTR‐Alu, subfamily E (SVA‐E) retrotransposon element, which was present in a mosaic state in the mother. This SVA‐E insertion disrupts correct splicing of the gene, resulting in loss of a functional allele. This case demonstrates the power of OGM and long‐read sequencing to identify genomic variations in high‐risk cancer‐predisposing genes that are refractory to detection with standard techniques, thereby completing the clinical and molecular diagnosis of such complex cases and greatly improving counseling and surveillance of the families involved. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

In Silico identification of a common mobile element insertion in exon 4 of RP1

Mobile element insertions (MEIs) typically exceed the read lengths of short-read sequencing technologies and are therefore frequently missed. Recently, a founder Alu insertion in exon 4 of RP1 has been detected in Japanese patients with macular dystrophy by PCR and gel electrophoresis. We aimed to develop a grep search program for the detection of the Alu insertion in exon 4 of RP1 using unprocessed short reads. Among 494 unrelated Korean patients with inherited eye diseases, 273 patients with specific retinal phenotypes who were previously genotyped by targeted panel or whole exome sequencing were selected. Five probands had a single heterozygous truncating RP1 variant, and one of their unaffected parents also carry this variant. To find a hidden genetic variant, whole genome sequencing was performed in two patients, and it revealed AluY c.4052_4053ins328/p.(Tyr1352Alafs*9) insertion in RP1 exon 4. This AluY insertion was additionally identified in other 3 families, which was confirmed by PCR and gel electrophoresis. We developed simplified grep search program to detect this AluY insertion in RP1 exon 4. The simple grep search revealed a median variant allele frequency of 0.282 (interquartile range, 0.232–0.383), with no false-positive results using 120 control samples. The MEI in RP1 exon 4 was a common founder mutation in Korean, occurring in 1.8% of our cohort. The RP1-Alu grep program efficiently detected the AluY insertion, without the preprocessing of raw data or complex installation processes.

Comprehensive identification of transposable element insertions using multiple sequencing technologies

Transposable elements (TEs) help shape the structure and function of the human genome. When inserted into some locations, TEs may disrupt gene regulation and cause diseases. Here, we present xTea (x-Transposable element analyzer), a tool for identifying TE insertions in whole-genome sequencing data. Whereas existing methods are mostly designed for short-read data, xTea can be applied to both short-read and long-read data. Our analysis shows that xTea outperforms other short read-based methods for both germline and somatic TE insertion discovery. With long-read data, we created a catalogue of polymorphic insertions with full assembly and annotation of insertional sequences for various types of retroelements, including pseudogenes and endogenous retroviruses. Notably, we find that individual genomes have an average of nine groups of full-length L1s in centromeres, suggesting that centromeres and other highly repetitive regions such as telomeres are a significant yet unexplored source of active L1s. xTea is available at https://github.com/parklab/xTea .

Cas9 targeted enrichment of mobile elements using nanopore sequencing

Mobile element insertions (MEIs) are repetitive genomic sequences that contribute to genetic variation and can lead to genetic disorders. Targeted and whole-genome approaches using short-read sequencing have been developed to identify reference and non-reference MEIs; however, the read length hampers detection of these elements in complex genomic regions. Here, we pair Cas9-targeted nanopore sequencing with computational methodologies to capture active MEIs in human genomes. We demonstrate parallel enrichment for distinct classes of MEIs, averaging 44% of reads on-targeted signals and exhibiting a 13.4-54x enrichment over whole-genome approaches. We show an individual flow cell can recover most MEIs (97% L1Hs, 93% AluYb, 51% AluYa, 99% SVA_F, and 65% SVA_E). We identify seventeen non-reference MEIs in GM12878 overlooked by modern, long-read analysis pipelines, primarily in repetitive genomic regions. This work introduces the utility of nanopore sequencing for MEI enrichment and lays the foundation for rapid discovery of elusive, repetitive genetic elements.

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.

Whole genome sequencing of 45 Japanese patients with intellectual disability

Intellectual disability (ID) is characterized by significant limitations in both intellectual functioning and adaptive behaviors, originating before the age of 18 years. However, the genetic etiologies of ID are still incompletely elucidated due to the wide range of clinical and genetic heterogeneity. Whole genome sequencing (WGS) has been applied as a single-step clinical diagnostic tool for ID because it detects genetic variations with a wide range of resolution from single nucleotide variants (SNVs) to structural variants (SVs). To explore the causative genes for ID, we employed WGS in 45 patients from 44 unrelated Japanese families and performed a stepwise screening approach focusing on the coding variants in the genes. Here, we report 12 pathogenic and likely pathogenic variants: seven heterozygous variants of ADNP, SATB2, ANKRD11, PTEN, TCF4, SPAST, and KCNA2, three hemizygous variants of SMS, SLC6A8, and IQSEC2, and one homozygous variant in AGTPBP1. Of these, four were considered novel. Furthermore, a novel 76 kb deletion containing exons 1 and 2 in DYRK1A was identified. We confirmed the clinical and genetic heterogeneity and high frequency of de novo causative variants (8/12, 66.7%). This is the first report of WGS analysis in Japanese patients with ID. Our results would provide insight into the correlation between novel variants and expanded phenotypes of the disease.

COG5 variants lead to complex early onset retinal degeneration, upregulation of PERK and DNA damage

Leber congenital amaurosis (LCA), a form of autosomal recessive severe early-onset retinal degeneration, is an important cause of childhood blindness. This may be associated with systemic features or not. Here we identified COG5 compound-heterozygous variants in patients affected with a complex LCA phenotype associated with microcephaly and skeletal dysplasia. COG5 is a component of the COG complex, which facilitates retrograde Golgi trafficking; if disrupted this can result in protein misfolding. To date, variants in COG5 have been associated with a distinct congenital disorder of glycosylation (type IIi) and with a variant of Friedreich’s ataxia. We show that COG5 variants can also result in fragmentation of the Golgi apparatus and upregulation of the UPR modulator, PKR-like endoplasmic reticulum kinase (PERK). In addition, upregulation of PERK induces DNA damage in cultured cells and in murine retina. This study identifies a novel role for COG5 in maintaining ER protein homeostasis and that disruption of that role results in activation of PERK and early-onset retinal degeneration, microcephaly and skeletal dysplasia. These results also highlight the importance of the UPR pathway in early-onset retinal dystrophy and as potential therapeutic targets for patients.

Early-onset breast cancer in a woman with a germline mobile element insertion resulting in BRCA2 disruption: a case report

Mobile element insertions (MEIs) contribute to genomic diversity, but they can be responsible for human disease in some cases. Initial clinical testing (BRCA1, BRCA2 and PALB2) in a 40-year-old female with unilateral breast cancer did not detect any pathogenic variants. Subsequent reanalysis for MEIs detected a novel likely pathogenic insertion of the retrotransposon element (RE) c.7894_7895insSVA in BRCA2. This case highlights the importance of bioinformatic pipeline optimization for the detection of MEIs in genes associated with hereditary cancer, as early detection can significantly impact clinical management.

Phenotype Driven Analysis of Whole Genome Sequencing Identifies Deep Intronic Variants that Cause Retinal Dystrophies by Aberrant Exonization

Purpose

To demonstrate the effectiveness of combining retinal phenotyping and focused variant filtering from genome sequencing (GS) in identifying deep intronic disease causing variants in inherited retinal dystrophies.

Methods

Affected members from three pedigrees with classical enhanced S-cone syndrome (ESCS; Pedigree 1), congenital stationary night blindness (CSNB; Pedigree 2), and achromatopsia (ACHM; Pedigree 3), respectively, underwent detailed ophthalmologic evaluation, optical coherence tomography, and electroretinography. The probands underwent panel-based genetic testing followed by GS analysis. Minigene constructs (NR2E3, GPR179 and CNGB3) and patient-derived cDNA experiments (NR2E3 and GPR179) were performed to assess the functional effect of the deep intronic variants.

Results

The electrophysiological findings confirmed the clinical diagnosis of ESCS, CSNB, and ACHM in the respective pedigrees. Panel-based testing revealed heterozygous pathogenic variants in NR2E3 (NM_014249.3; c.119-2A>C; Pedigree 1) and CNGB3 (NM_019098.4; c.1148delC/p.Thr383Ilefs*13; Pedigree 3). The GS revealed heterozygous deep intronic variants in Pedigrees 1 (NR2E3; c.1100+1124G>A) and 3 (CNGB3; c.852+4751A>T), and a homozygous GPR179 variant in Pedigree 2 (NM_001004334.3; c.903+343G>A). The identified variants segregated with the phenotype in all pedigrees. All deep intronic variants were predicted to generate a splice acceptor gain causing aberrant exonization in NR2E3 [89 base pairs (bp)], GPR179 (197 bp), and CNGB3 (73 bp); splicing defects were validated through patient-derived cDNA experiments and/or minigene constructs and rescued by antisense oligonucleotide treatment.

Conclusions

Deep intronic mutations contribute to missing heritability in retinal dystrophies. Combining results from phenotype-directed gene panel testing, GS, and in silico splice prediction tools can help identify these difficult-to-detect pathogenic deep intronic variants.

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.

rMETL: sensitive mobile element insertion detection with long read realignment

SUMMARY

Mobile element insertion (MEI) is a major category of structure variations (SVs). The rapid development of long read sequencing technologies provides the opportunity to detect MEIs sensitively. However, the signals of MEI implied by noisy long reads are highly complex due to the repetitiveness of mobile elements as well as the high sequencing error rates. Herein, we propose the Realignment-based Mobile Element insertion detection Tool for Long read (rMETL). Benchmarking results of simulated and real datasets demonstrate that rMETL enables to handle the complex signals to discover MEIs sensitively. It is suited to produce high-quality MEI callsets in many genomics studies.

AVAILABILITY AND IMPLEMENTATION

rMETL is available from https://github.com/hitbc/rMETL.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

TypeTE: a tool to genotype mobile element insertions from whole genome resequencing data

Alu retrotransposons account for more than 10% of the human genome, and insertions of these elements create structural variants segregating in human populations. Such polymorphic Alus are powerful markers to understand population structure, and they represent variants that can greatly impact genome function, including gene expression. Accurate genotyping of Alus and other mobile elements has been challenging. Indeed, we found that Alu genotypes previously called for the 1000 Genomes Project are sometimes erroneous, which poses significant problems for phasing these insertions with other variants that comprise the haplotype. To ameliorate this issue, we introduce a new pipeline - TypeTE - which genotypes Alu insertions from whole-genome sequencing data. Starting from a list of polymorphic Alus, TypeTE identifies the hallmarks (poly-A tail and target site duplication) and orientation of Alu insertions using local re-assembly to reconstruct presence and absence alleles. Genotype likelihoods are then computed after re-mapping sequencing reads to the reconstructed alleles. Using a high-quality set of PCR-based genotyping of >200 loci, we show that TypeTE improves genotype accuracy from 83% to 92% in the 1000 Genomes dataset. TypeTE can be readily adapted to other retrotransposon families and brings a valuable toolbox addition for population genomics.

The Medical Genome Reference Bank contains whole genome and phenotype data of 2570 healthy elderly

Population health research is increasingly focused on the genetic determinants of healthy ageing, but there is no public resource of whole genome sequences and phenotype data from healthy elderly individuals. Here we describe the first release of the Medical Genome Reference Bank (MGRB), comprising whole genome sequence and phenotype of 2570 elderly Australians depleted for cancer, cardiovascular disease, and dementia. We analyse the MGRB for single-nucleotide, indel and structural variation in the nuclear and mitochondrial genomes. MGRB individuals have fewer disease-associated common and rare germline variants, relative to both cancer cases and the gnomAD and UK Biobank cohorts, consistent with risk depletion. Age-related somatic changes are correlated with grip strength in men, suggesting blood-derived whole genomes may also provide a biologic measure of age-related functional deterioration. The MGRB provides a broadly applicable reference cohort for clinical genetics and genomic association studies, and for understanding the genetics of healthy ageing.

Healthspan and healthy aging are areas of research with potential socioeconomic impact. Here, the authors present the Medical Genome Reference Bank (MGRB) which consist of over 4,000 individuals aged 70 years and older without a history of the major age-related diseases and report on results from whole-genome sequencing and association analyses.

Mobile element insertion detection in 89,874 clinical exomes

Purpose

Exome sequencing (ES) is increasingly used for the diagnosis of rare genetic disease. However, some pathogenic sequence variants within the exome go undetected due to the technical difficulty of identifying them. Mobile element insertions (MEIs) are a known cause of genetic disease in humans but have been historically difficult to detect via ES and similar targeted sequencing methods.

Methods

We developed and applied a novel MEI detection method prospectively to samples received for clinical ES beginning in November 2017. Positive MEI findings were confirmed by an orthogonal method and reported back to the ordering provider. In this study, we examined 89,874 samples from 38,871 cases.

Results

Diagnostic MEIs were present in 0.03% (95% binomial test confidence interval: 0.02–0.06%) of all cases and account for 0.15% (95% binomial test confidence interval: 0.08–0.25%) of cases with a molecular diagnosis. One diagnostic MEI was a novel founder event. Most patients with pathogenic MEIs had prior genetic testing, three of whom had previous negative DNA sequencing analysis of the diagnostic gene.

Conclusion

MEI detection from ES is a valuable diagnostic tool, reveals molecular findings that may be undetected by other sequencing assays, and increases diagnostic yield by 0.15%.

Transdifferentiation of Mouse Embryonic Fibroblasts into Dopaminergic Neurons Reactivates LINE-1 Repetitive Elements

In mammals, LINE-1 (L1) retrotransposons constitute between 15% and 20% of the genome. Although only a few copies have retained the ability to retrotranspose, evidence in brain and differentiating pluripotent cells indicates that L1 retrotransposition occurs and creates mosaics in normal somatic tissues. The function of de novo insertions remains to be understood. The transdifferentiation of mouse embryonic fibroblasts to dopaminergic neuronal fate provides a suitable model for studying L1 dynamics in a defined genomic and unaltered epigenomic background. We found that L1 elements are specifically re-expressed and mobilized during the initial stages of reprogramming and that their insertions into specific acceptor loci coincides with higher chromatin accessibility and creation of new transcribed units. Those events accompany the maturation of neuronal committed cells. We conclude that L1 retrotransposition is a non-random process correlating with chromatin opening and lncRNA production that accompanies direct somatic cell reprogramming.

•

L1 activation accompanies induced dopaminergic neuron maturation

•

L1 inhibition impairs the transdifferentiation potential of MEFs

•

L1 retrotransposition creates a lineage-specific genetic mosaicism

•

L1 insertions correlates with open chromatin and lncRNA transcription

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.

Comprehensive evaluation of structural variation detection algorithms for whole genome sequencing

BACKGROUND

Structural variations (SVs) or copy number variations (CNVs) greatly impact the functions of the genes encoded in the genome and are responsible for diverse human diseases. Although a number of existing SV detection algorithms can detect many types of SVs using whole genome sequencing (WGS) data, no single algorithm can call every type of SVs with high precision and high recall.

RESULTS

We comprehensively evaluate the performance of 69 existing SV detection algorithms using multiple simulated and real WGS datasets. The results highlight a subset of algorithms that accurately call SVs depending on specific types and size ranges of the SVs and that accurately determine breakpoints, sizes, and genotypes of the SVs. We enumerate potential good algorithms for each SV category, among which GRIDSS, Lumpy, SVseq2, SoftSV, Manta, and Wham are better algorithms in deletion or duplication categories. To improve the accuracy of SV calling, we systematically evaluate the accuracy of overlapping calls between possible combinations of algorithms for every type and size range of SVs. The results demonstrate that both the precision and recall for overlapping calls vary depending on the combinations of specific algorithms rather than the combinations of methods used in the algorithms.

CONCLUSION

These results suggest that careful selection of the algorithms for each type and size range of SVs is required for accurate calling of SVs. The selection of specific pairs of algorithms for overlapping calls promises to effectively improve the SV detection accuracy.

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.

Transposable Elements: Classification, Identification, and Their Use As a Tool For Comparative Genomics

Most genomes are populated by hundreds of thousands of sequences originated from mobile elements. On the one hand, these sequences present a real challenge in the process of genome analysis and annotation. On the other hand, they are very interesting biological subjects involved in many cellular processes. Here we present an overview of transposable elements biodiversity, and we discuss different approaches to transposable elements detection and analyses.

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.

Retrotransposon insertion as a novel mutational event in Bardet-Biedl syndrome

Background

Bardet‐Biedl syndrome (BBS) is an autosomal recessive pleiotropic disorder of the primary cilia that leads to severe visual loss in the teenage years. Approximately 80% of BBS cases are explained by mutations in one of the 21 identified genes. Documented causative mutation types include missense, nonsense, copy number variation (CNV), frameshift deletions or insertions, and splicing variants.

Methods

Whole genome sequencing was performed on a patient affected with BBS for whom no mutations were identified using clinically approved genetic testing of the known genes. Analysis of the WGS was done using internal protocols and publicly available algorithms. The phenotype was defined by retrospective chart review.

Results

We document a female affected with BBS carrying the most common BBS1 mutation (BBS1: Met390Arg) on the maternal allele and an insertion of a ~1.7‐kb retrotransposon in exon 13 on the paternal allele. This retrotransposon insertion was not automatically annotated by the standard variant calling protocols used. This novel variant was identified by visual inspection of the alignment file followed by specific genome analysis with an available algorithm for transposable elements.

Conclusion

This report documents a novel mutation type associated with BBS and highlights the importance of systematically performing transposon detection analysis on WGS data of unsolved cases.

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.

LINE- and Alu-containing genomic instability hotspot at 16q24.1 associated with recurrent and nonrecurrent CNV deletions causative for ACDMPV

Transposable elements modify human genome by inserting into new loci or by mediating homology-, microhomology-, or homeology-driven DNA recombination or repair, resulting in genomic structural variation. Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a rare lethal neonatal developmental lung disorder caused by point mutations or copy-number variant (CNV) deletions of FOXF1 or its distant tissue-specific enhancer. Eighty-five percent of 45 ACDMPV-causative CNV deletions, of which junctions have been sequenced, had at least one of their two breakpoints located in a retrotransposon, with more than half of them being Alu elements. We describe a novel ∼35 kb-large genomic instability hotspot at 16q24.1, involving two evolutionarily young LINE-1 (L1) elements, L1PA2 and L1PA3, flanking AluY, two AluSx, AluSx1, and AluJr elements. The occurrence of L1s at this location coincided with the branching out of the Homo-Pan-Gorilla clade, and was preceded by the insertion of AluSx, AluSx1, and AluJr. Our data show that, in addition to mediating recurrent CNVs, L1 and Alu retrotransposons can predispose the human genome to formation of variably sized CNVs, both of clinical and evolutionary relevance. Nonetheless, epigenetic or other genomic features of this locus might also contribute to its increased instability.

The Role of Transposable Elements in Speciation

Understanding the phenotypic and molecular mechanisms that contribute to genetic diversity between and within species is fundamental in studying the evolution of species. In particular, identifying the interspecific differences that lead to the reduction or even cessation of gene flow between nascent species is one of the main goals of speciation genetic research. Transposable elements (TEs) are DNA sequences with the ability to move within genomes. TEs are ubiquitous throughout eukaryotic genomes and have been shown to alter regulatory networks, gene expression, and to rearrange genomes as a result of their transposition. However, no systematic effort has evaluated the role of TEs in speciation. We compiled the evidence for TEs as potential causes of reproductive isolation across a diversity of taxa. We find that TEs are often associated with hybrid defects that might preclude the fusion between species, but that the involvement of TEs in other barriers to gene flow different from postzygotic isolation is still relatively unknown. Finally, we list a series of guides and research avenues to disentangle the effects of TEs on the origin of new species.

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.

Mapping and phasing of structural variation in patient genomes using nanopore sequencing

Despite improvements in genomics technology, the detection of structural variants (SVs) from short-read sequencing still poses challenges, particularly for complex variation. Here we analyse the genomes of two patients with congenital abnormalities using the MinION nanopore sequencer and a novel computational pipeline-NanoSV. We demonstrate that nanopore long reads are superior to short reads with regard to detection of de novo chromothripsis rearrangements. The long reads also enable efficient phasing of genetic variations, which we leveraged to determine the parental origin of all de novo chromothripsis breakpoints and to resolve the structure of these complex rearrangements. Additionally, genome-wide surveillance of inherited SVs reveals novel variants, missed in short-read data sets, a large proportion of which are retrotransposon insertions. We provide a first exploration of patient genome sequencing with a nanopore sequencer and demonstrate the value of long-read sequencing in mapping and phasing of SVs for both clinical and research applications.

Phylogenetic Conflict in Bears Identified by Automated Discovery of Transposable Element Insertions in Low-Coverage Genomes

Phylogenetic reconstruction from transposable elements (TEs) offers an additional perspective to study evolutionary processes. However, detecting phylogenetically informative TE insertions requires tedious experimental work, limiting the power of phylogenetic inference. Here, we analyzed the genomes of seven bear species using high-throughput sequencing data to detect thousands of TE insertions. The newly developed pipeline for TE detection called TeddyPi (TE detection and discovery for Phylogenetic Inference) identified 150,513 high-quality TE insertions in the genomes of ursine and tremarctine bears. By integrating different TE insertion callers and using a stringent filtering approach, the TeddyPi pipeline produced highly reliable TE insertion calls, which were confirmed by extensive in vitro validation experiments. Analysis of single nucleotide substitutions in the flanking regions of the TEs shows that these substitutions correlate with the phylogenetic signal from the TE insertions. Our phylogenomic analyses show that TEs are a major driver of genomic variation in bears and enabled phylogenetic reconstruction of a well-resolved species tree, despite strong signals for incomplete lineage sorting and introgression. The analyses show that the Asiatic black, sun, and sloth bear form a monophyletic clade, in which phylogenetic incongruence originates from incomplete lineage sorting. TeddyPi is open source and can be adapted to various TE and structural variation callers. The pipeline makes it possible to confidently extract thousands of TE insertions even from low-coverage genomes (∼10×) of nonmodel organisms. This opens new possibilities for biologists to study phylogenies and evolutionary processes as well as rates and patterns of (retro-)transposition and structural variation.