LIONS: analysis suite for detecting and quantifying transposable element initiated transcription from RNA-seq.

SUMMARY: Transposable elements (TEs) influence the evolution of novel transcriptional networks yet the specific and meaningful interpretation of how TE-derived transcriptional initiation contributes to the transcriptome has been marred by computational and methodological deficiencies. We developed LIONS for the analysis of RNA-seq data to specifically detect and quantify TE-initiated transcripts. AVAILABILITY AND IMPLEMENTATION: Source code, container, test data and instruction manual are freely available at www.github.com/ababaian/LIONS. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Distinct isoform of FABP7 revealed by screening for retroelement-activated genes in diffuse large B-cell lymphoma.

Remnants of ancient transposable elements (TEs) are abundant in mammalian genomes. These sequences harbor multiple regulatory motifs and hence are capable of influencing expression of host genes. In response to environmental changes, TEs are known to be released from epigenetic repression and to become transcriptionally active. Such activation could also lead to lineage-inappropriate activation of oncogenes, as one study described in Hodgkin lymphoma. However, little further evidence for this mechanism in other cancers has been reported. Here, we reanalyzed whole transcriptome data from a large cohort of patients with diffuse large B-cell lymphoma (DLBCL) compared with normal B-cell centroblasts to detect genes ectopically expressed through activation of TE promoters. We have identified 98 such TE-gene chimeric transcripts that were exclusively expressed in primary DLBCL cases and confirmed several in DLBCL-derived cell lines. We further characterized a TE-gene chimeric transcript involving a fatty acid-binding protein gene (LTR2-FABP7), normally expressed in brain, that was ectopically expressed in a subset of DLBCL patients through the use of an endogenous retroviral LTR promoter of the LTR2 family. The LTR2-FABP7 chimeric transcript encodes a novel chimeric isoform of the protein with characteristics distinct from native FABP7. In vitro studies reveal a dependency for DLBCL cell line proliferation and growth on LTR2-FABP7 chimeric protein expression. Taken together, these data demonstrate the significance of TEs as regulators of aberrant gene expression in cancer and suggest that LTR2-FABP7 may contribute to the pathogenesis of DLBCL in a subgroup of patients.

Epigenetic mechanism causes Wnt9b deficiency and nonsyndromic cleft lip and palate in the A/WySn mouse strain.

BACKGROUND: The heritable multifactorial etiology of human nonsyndromic cleft lip with or without cleft palate (CL ± P) is not understood. CL ± P occurs in 15% of neonates in the homozygous A/WySn mouse strain, with a multifactorial genetic etiology, the clf1 and clf2 variant genes. Clf1 acts as a mutant allele of Wnt9b but its coding sequence is normal. An IAP (intracisternal A particle) retrotransposon inserted near the Wnt9b gene is associated with clf1. METHODS: Transcription of noncoding sequence between the IAP and the Wnt9b gene was examined in A/WySn embryos. The levels of Wnt9b transcript and of an "IAP antisense" transcript initiated in the IAP and extending into the noncoding interval were assayed in A/WySn and C57BL/6J whole embryos or heads across embryonic days 8 to 12. Methylation of the 5' LTR of the IAP was examined in E12 A/WySn embryo heads. RESULTS: Mean Wnt9b transcript levels were lower in A/WySn than in C57BL/6J at all ages examined and lower in CL ± P embryos than in their normal littermates. The "IAP antisense" transcript was found in all A/WySn embryos and was highest in CL ± P embryos. The IAP at Wnt9b was generally unmethylated in CL ± P embryos and approximately 50% methylated in normal littermates. CONCLUSION: The clf1 mutation in A/WySn is a "metastable epiallele", in which stochastic deficiency in some individuals of DNA methylation of a retrotransposon uniquely inserted near the Wnt9b gene allows transcriptional activity of the retrotransposon and interference with transcription from Wnt9b. Methylation of metastable epialleles should be investigated in human nonsyndromic CL ± P.

Retrotransposon-induced heterochromatin spreading in the mouse revealed by insertional polymorphisms.

The "arms race" relationship between transposable elements (TEs) and their host has promoted a series of epigenetic silencing mechanisms directed against TEs. Retrotransposons, a class of TEs, are often located in repressed regions and are thought to induce heterochromatin formation and spreading. However, direct evidence for TE-induced local heterochromatin in mammals is surprisingly scarce. To examine this phenomenon, we chose two mouse embryonic stem (ES) cell lines that possess insertionally polymorphic retrotransposons (IAP, ETn/MusD, and LINE elements) at specific loci in one cell line but not the other. Employing ChIP-seq data for these cell lines, we show that IAP elements robustly induce H3K9me3 and H4K20me3 marks in flanking genomic DNA. In contrast, such heterochromatin is not induced by LINE copies and only by a minority of polymorphic ETn/MusD copies. DNA methylation is independent of the presence of IAP copies, since it is present in flanking regions of both full and empty sites. Finally, such spreading into genes appears to be rare, since the transcriptional start sites of very few genes are less than one Kb from an IAP. However, the B3galtl gene is subject to transcriptional silencing via IAP-induced heterochromatin. Hence, although rare, IAP-induced local heterochromatin spreading into nearby genes may influence expression and, in turn, host fitness.

The clf2 gene has an epigenetic role in the multifactorial etiology of cleft lip and palate in the A/WySn mouse strain.

BACKGROUND: The A/WySn mouse strain with 15 to 20% penetrance of cleft lip and palate (CLP) is an animal model for human multifactorial CLP. The CLP is due to two unlinked genes that interact epistatically, Wnt9b(clf1) and clf2, plus a maternal effect. The Wnt9b(clf1) mutation is an IAP transposon insertion. The clf2 gene, with unknown function, was located in a 13.6 Mb region of chromosome 13 containing 145 genes. METHODS: To reduce the clf2 candidate region, 1146 mice segregating for A/WySn and C57BL/6J alleles at clf2 were screened for recombinants by simple sequence-length polymorphism haplotypes; recombinants' testcross progeny were typed for CLP and simple-sequence length polymorphisms. To identify the function of clf2, the effect of clf2 genotype on risk of CLP was tested in Wnt9b(null/null) knockouts and in compound mutants (Wnt9b(clf1/null) ), and the methylation of the IAP at Wnt9b was assayed in the Wnt9b(clf1/null) mutants by combined bisulfite restriction analysis. RESULTS: The location of clf2 was redefined to 3.0 Mb between Cntnap3 and AK029746 containing 48 genes, of which 30 are Zfp genes. The clf2 genotype had no detectable effect on Wnt9b(null/null) embryos, but strongly affected risk of CLP and methylation of the IAP in Wnt9b(clf1/null) embryos. CLP was associated with low levels of methylation of the IAP. CONCLUSIONS: The clf2 gene is the first identified polymorphism that affects the epigenetic methylation and silencing of IAP retrotransposons. This CLP model raises the question of whether parallel epigenetic factors are involved in risk and environmental sensitivity of human CLP.

Genome-wide assessments reveal extremely high levels of polymorphism of two active families of mouse endogenous retroviral elements.

Endogenous retroviral elements (ERVs) in mice are significant genomic mutagens, causing approximately 10% of all reported spontaneous germ line mutations in laboratory strains. The majority of these mutations are due to insertions of two high copy ERV families, the IAP and ETn/MusD elements. This significant level of ongoing retrotranspositional activity suggests that inbred mice are highly variable in content of these two ERV groups. However, no comprehensive genome-wide studies have been performed to assess their level of polymorphism. Here we compared three test strains, for which sufficient genomic sequence is available, to each other and to the reference C57BL/6J genome and detected very high levels of insertional polymorphism for both ERV families, with an estimated false discovery rate of only 0.4%. Specifically, we found that at least 60% of IAP and 25% of ETn/MusD elements detected in any strain are absent in one or more of the other three strains. The polymorphic nature of a set of 40 ETn/MusD elements found within gene introns was confirmed using genomic PCR on DNA from a panel of mouse strains. For some cases, we detected gene-splicing abnormalities involving the ERV and obtained additional evidence for decreased gene expression in strains carrying the insertion. In total, we identified nearly 700 polymorphic IAP or ETn/MusD ERVs or solitary LTRs that reside in gene introns, providing potential candidates that may contribute to gene expression differences among strains. These extreme levels of polymorphism suggest that ERV insertions play a significant role in genetic drift of mouse lines.

Inter-Strain Epigenomic Profiling Reveals a Candidate IAP Master Copy in C3H Mice.

Insertions of endogenous retroviruses cause a significant fraction of mutations in inbred mice but not all strains are equally susceptible. Notably, most new Intracisternal A particle (IAP) ERV mutagenic insertions have occurred in C3H mice. We show here that strain-specific insertional polymorphic IAPs accumulate faster in C3H/HeJ mice, relative to other sequenced strains, and that IAP transcript levels are higher in C3H/HeJ embryonic stem (ES) cells compared to other ES cells. To investigate the mechanism for high IAP activity in C3H mice, we identified 61 IAP copies in C3H/HeJ ES cells enriched with H3K4me3 (a mark of active promoters) and, among those tested, all are unmethylated in C3H/HeJ ES cells. Notably, 13 of the 61 are specific to C3H/HeJ and are members of the non-autonomous 1Δ1 IAP subfamily that is responsible for nearly all new insertions in C3H. One copy is full length with intact open reading frames and hence potentially capable of providing proteins in trans to other 1Δ1 elements. This potential "master copy" is present in other strains, including 129, but its 5' long terminal repeat (LTR) is methylated in 129 ES cells. Thus, the unusual IAP activity in C3H may be due to reduced epigenetic repression coupled with the presence of a master copy.

Retroviral elements and their hosts: insertional mutagenesis in the mouse germ line.

The inbred mouse is an invaluable model for human biology and disease. Nevertheless, when considering genetic mechanisms of variation and disease, it is important to appreciate the significant differences in the spectra of spontaneous mutations that distinguish these species. While insertions of transposable elements are responsible for only approximately 0.1% of de novo mutations in humans, the figure is 100-fold higher in the laboratory mouse. This striking difference is largely due to the ongoing activity of mouse endogenous retroviral elements. Here we briefly review mouse endogenous retroviruses (ERVs) and their influence on gene expression, analyze mechanisms of interaction between ERVs and the host cell, and summarize the variety of mutations caused by ERV insertions. The prevalence of mouse ERV activity indicates that the genome of the laboratory mouse is presently behind in the "arms race" against invasion.

Mouse germ line mutations due to retrotransposon insertions.

Transposable element (TE) insertions are responsible for a significant fraction of spontaneous germ line mutations reported in inbred mouse strains. This major contribution of TEs to the mutational landscape in mouse contrasts with the situation in human, where their relative contribution as germ line insertional mutagens is much lower. In this focussed review, we provide comprehensive lists of TE-induced mouse mutations, discuss the different TE types involved in these insertional mutations and elaborate on particularly interesting cases. We also discuss differences and similarities between the mutational role of TEs in mice and humans.

A novel isoform of IL-33 revealed by screening for transposable element promoted genes in human colorectal cancer.

Remnants of ancient transposable elements (TEs) are abundant in mammalian genomes. These sequences contain multiple regulatory motifs and hence are capable of influencing expression of host genes. TEs are known to be released from epigenetic repression and can become transcriptionally active in cancer. Such activation could also lead to lineage-inappropriate activation of oncogenes, as previously described in lymphomas. However, there are few reports of this mechanism occurring in non-blood cancers. Here, we re-analyzed whole transcriptome data from a large cohort of patients with colon cancer, compared to matched normal colon control samples, to detect genes or transcripts ectopically expressed through activation of TE promoters. Among many such transcripts, we identified six where the affected gene has described role in cancer and where the TE-driven gene mRNA is expressed in primary colon cancer, but not normal matched tissue, and confirmed expression in colon cancer-derived cell lines. We further characterized a TE-gene chimeric transcript involving the Interleukin 33 (IL-33) gene (termed LTR-IL-33), that is ectopically expressed in a subset of colon cancer samples through the use of an endogenous retroviral long terminal repeat (LTR) promoter of the MSTD family. The LTR-IL-33 chimeric transcript encodes a novel shorter isoform of the protein, which is missing the initial N-terminus (including many conserved residues) of Native IL-33. In vitro studies showed that LTR-IL-33 expression is required for optimal CRC cell line growth as 3D colonospheres. Taken together, these data demonstrate the significance of TEs as regulators of aberrant gene expression in colon cancer.

Isolated peripheral neuropathy in atypical metachromatic leukodystrophy: a recurrent mutation.

BACKGROUND: Metachromatic leukodystrophy (MLD) is a genetic neurodegenerative disorder resulting from a deficiency of arylsulfatase A. Late onset forms are relatively rare. Central nervous system (CNS) involvement is characteristic at all ages. METHODS: A patient in her late 40s with peripheral neuropathy was assessed by EEG, evoked potentials, CT and nerve conduction studies. Nerve and muscle biopsy samples were investigated by electron microscopy. Arylsulfatase A activity in leukocytes and excreted cerebroside sulfate were determined. The arylsulfatase A gene was investigated for mutations using polymerase chain reaction (PCR) and DNA sequencing. The identified mutation was expressed transiently in African green monkey kidney (COS) cells to determine the effect of the mutation on arylsulfatase A activity. RESULTS: Central nervous system functions were normal. Nerve conduction velocities were decreased. Sural nerve biopsy showed inclusions typical of MLD. Arylsulfatase A was less than 5% of normal. A homozygous mutation thr286pro was identified in the arylsulfatase A gene and demonstrated to be deleterious through transient expression studies. CONCLUSIONS: Our patient has a progressive peripheral neuropathy but has apparently intact CNS function at her present age of 57 years. Biochemical, physiological and pathological findings are consistent with a diagnosis of MLD. A homozygous mutation, thr286pro, found in her arylsulfatase A gene, decreased enzyme activity to a level consistent with a late onset form of MLD.

Visualized computational predictions of transcriptional effects by intronic endogenous retroviruses.

When endogenous retroviruses (ERVs) or other transposable elements (TEs) insert into an intron, the consequence on gene transcription can range from negligible to a complete ablation of normal transcripts. With the advance of sequencing technology, more and more insertionally polymorphic or private TE insertions are being identified in humans and mice, of which some could have a significant impact on host gene expression. Nevertheless, an efficient and low cost approach to prioritize their potential effect on gene transcription has been lacking. By building a computational model based on artificial neural networks (ANN), we demonstrate the feasibility of using machine-learning approaches to predict the likelihood that intronic ERV insertions will have major effects on gene transcription, focusing on the two ERV families, namely Intracisternal A-type Particle (IAP) and Early Transposon (ETn)/MusD elements, which are responsible for the majority of ERV-induced mutations in mice. We trained the ANN model using properties associated with these ERVs known to cause germ-line mutations (positive cases) and properties associated with likely neutral ERVs of the same families (negative cases), and derived a set of prediction plots that can visualize the likelihood of affecting gene transcription by ERV insertions. Our results show a highly reliable prediction power of our model, and offer a potential approach to computationally screen for other types of TE insertions that may affect gene transcription or even cause disease.

Epigenetic interplay between mouse endogenous retroviruses and host genes.

BACKGROUND: Transposable elements are often the targets of repressive epigenetic modifications such as DNA methylation that, in theory, have the potential to spread toward nearby genes and induce epigenetic silencing. To better understand the role of DNA methylation in the relationship between transposable elements and genes, we assessed the methylation state of mouse endogenous retroviruses (ERVs) located near genes. RESULTS: We found that ERVs of the ETn/MusD family show decreased DNA methylation when near transcription start sites in tissues where the nearby gene is expressed. ERVs belonging to the IAP family, however, are generally heavily methylated, regardless of the genomic environment and the tissue studied. Furthermore, we found full-length ETn and IAP copies that display differential DNA methylation between their two long terminal repeats (LTRs), suggesting that the environment surrounding gene promoters can prevent methylation of the nearby LTR. Spreading from methylated ERV copies to nearby genes was rarely observed, with the regions between the ERVs and genes apparently acting as a boundary, enriched in H3K4me3 and CTCF, which possibly protects the unmethylated gene promoter. Furthermore, the flanking regions of unmethylated ERV copies harbor H3K4me3, consistent with spreading of euchromatin from the host gene toward ERV insertions. CONCLUSIONS: We have shown that spreading of DNA methylation from ERV copies toward active gene promoters is rare. We provide evidence that genes can be protected from ERV-induced heterochromatin spreading by either blocking the invasion of repressive marks or by spreading euchromatin toward the ERV copy.

Evidence for epigenetic maintenance of Ly49a monoallelic gene expression.

Although structurally unrelated, the human killer cell Ig-like (KIR) genes and the rodent lectin-like Ly49 genes serve similar functional roles in NK cells. Moreover, both gene families display variegated, monoallelic expression patterns established at the transcriptional level. DNA methylation has been shown to play an important role in maintenance of expression patterns of KIR genes, which have CpG island promoters. The potential role of DNA methylation in expression of Ly49 genes, which have CpG-poor promoters, is unknown. In this study, we show that hypomethylation of the region encompassing the Pro-2 promoter of Ly49a and Ly49c in primary C57BL/6 NK cells correlates with expression of the gene. Using C57BL/6 x BALB/c F1 hybrid mice, we demonstrate that the expressed allele of Ly49a is hypomethylated while the nonexpressed allele is heavily methylated, indicating a role for epigenetics in maintaining monoallelic Ly49 gene expression. Furthermore, the Ly49a Pro-2 region is heavily methylated in fetal NK cells but variably methylated in nonlymphoid tissues. Finally, in apparent contrast to the KIR genes, we show that DNA methylation and the histone acetylation state of the Pro-2 region are strictly linked with Ly49a expression status.

Genomic deletions and precise removal of transposable elements mediated by short identical DNA segments in primates.

Insertion of transposable elements is a major cause of genomic expansion in eukaryotes. Less is understood, however, about mechanisms underlying contraction of genomes. In this study, we show that retroelements can, in rare cases, be precisely deleted from primate genomes, most likely via recombination between 10- to 20-bp target site duplications (TSDs) flanking the retroelement. The deleted loci are indistinguishable from pre-integration sites, effectively reversing the insertion. Through human-chimpanzee-Rhesus monkey genomic comparisons, we estimate that 0.5%-1% of apparent retroelement "insertions" distinguishing humans and chimpanzees actually represent deletions. Furthermore, we demonstrate that 19% of genomic deletions of 200-500 bp that have occurred since the human-chimpanzee divergence are associated with flanking identical repeats of at least 10 bp. A large number of deletions internal to Alu elements were also found flanked by homologies. These results suggest that illegitimate recombination between short direct repeats has played a significant role in human genome evolution. Moreover, this study lends perspective to the view that insertions of retroelements represent unidirectional genetic events.

Investigations of the genomic region that contains the clf1 mutation, a causal gene in multifactorial cleft lip and palate in mice.

BACKGROUND: Human nonsyndromic cleft lip and palate, CL(P), is genetically complex, with one contributing gene on chromosome 17q. A potentially homologous gene, clf1 on distal chromosome 11, is part of the digenic cause of the 10-30% CL(P) in the A/WySn mouse strain. Here we report our progress toward identifying the clf1 mutation. METHODS: Transcription from all of the known and predicted genes in the 1.5-Mb candidate region was examined in A/WySn and control (AXB-4/Pgn) ED10-11 embryo heads. The marker haplotype for 28 inbred strains across the clf1 region was obtained. The entire transcripts of Wnt9b and Wnt3 in A/WySn were sequenced. Using long PCR, the genomic region from Wnt3 throughWnt9b was screened in A/WySn for an inserted retrotransposon. RESULTS: Gosr2, Wnt9b, Wnt3, Nsf, Arf2, Crhr1, Mapt, Cdc27, Myl4, Itgb3, chr11_20.152, chr11_20.154, chr11_20.155, and chr11_20.156 are expressed in ED10-11 heads. None is absent or detectably reduced in A/WySn. The ancestral pre-clf1 mutation haplotype was found in CBA/J mice. By a test-cross, CBA/J was confirmed to lack the clf1 mutation. Three single-nucleotide variants in A/WySn (vs. C57BL/6J) were found in each of the 3' untranslated regions (3'UTRs) of Wnt3 and of Wnt9b, respectively; their presence in CBA/J shows that none are the clf1 mutation. An inserted intracisternal A particle (IAP) retrotransposon located 6.6 kb from the 3' end of Wnt9b was found in A/WySn and in all clf1 strains tested. This IAP is absent in C57BL/6J and CBA/J. CONCLUSIONS: The clf1 mutation is a genomic alteration present in A/WySn and absent in the ancestral chromosomal segment in CBA/J. The IAP retrotransposon insertion near Wnt9b in A/WySn fits this criterion; we predict that interference with Wnt9b function by this IAP is the clf1 mutation.

Spectrum of mutations in the arylsulfatase A gene in a Canadian DNA collection including two novel frameshift mutations, a new missense mutation (C488R) and an MLD mutation (R84Q) in cis with a pseudodeficiency allele.

We describe three novel mutations in the human arylsulfatase A gene in three patients with MLD, an autosomal recessive lysosomal storage disorder. An insertion, 2590_2591insCCCC in exon 8 and a deletion, 752_758delGCCGGCC, in exon 3 will both result in frameshifts. A mutation in exon 8, 2566T-->C, results in a missense mutation C488R, disrupting an unusual cysteine-knot at the C-terminal end of the protein. All three mutations are heterozygous with previously documented mutations. A previously reported mutation, R84Q was identified on a pseudodeficiency allele. These mutations are part of a heterogeneous spectrum of mutations found in a collection of DNA samples from MLD patients from across Canada and the USA.

Ly49 genes in non-rodent mammals.

The Ly49 family of natural killer (NK) receptors is encoded by a highly polymorphic multigene family in the mouse and is also present in multiple copies in the rat. However, this gene exists as a single copy in primates and is mutated to non-function in humans. We recently showed that the cow also likely has only one Ly49 gene, but it is unclear what the Ly49 gene content is for other mammals. We have now isolated Ly49 cDNAs from the domestic cat, dog and pig and show that the corresponding gene appears to be single copy in these three species. The open reading frame is intact in all the genes and the putative proteins contain an immune tyrosine-based inhibition motif (ITIM), suggesting a role as an inhibitory receptor. In contrast to the other mammals, several Ly49-like genes appear to exist in the horse, indicating that amplification of this locus has occurred in a non-rodent lineage. Finally, phylogenetic analysis suggests that the rodent Ly49 genes have evolved more rapidly than their counterparts in mammals where the gene has remained as a single copy.

Sequence analysis of the ly49 cluster in C57BL/6 mice: a rapidly evolving multigene family in the immune system.

The cytotoxic activity of murine natural killer cells is controlled in part through the action of genes belonging to the Ly49 family. Members of this multigene family are found in a region on mouse chromosome 6 termed the natural killer gene complex. Using data available through public databases, we performed sequence analysis of a 620-kb region in C57Bl/6 (B6) mice that contains the Ly49 genes. The contiguous genomic sequence has allowed us to describe the complete B6 Ly49 gene repertoire, which includes two recently described genes as well as three partial genes. We have shown that the genes in the cluster have evolved through a series of large duplication events involving units of one or more genes and we have attempted to characterize the nature of the duplication end points. Finally, we have used information regarding gene sequence relationships and insertion of repetitive elements to construct a model for the evolution of the gene cluster. Our study illustrates that the Ly49 cluster represents an example of a rapidly evolving gene family, and continued analysis of this region in other strains will undoubtedly provide further insight into mechanisms for generating genomic diversity.

Structure and expression of mobile ETnII retroelements and their coding-competent MusD relatives in the mouse.

ETnII elements are mobile members of the repetitive early transposon family of mouse long terminal repeat (LTR) retroelements and have caused a number of mutations by inserting into genes. ETnII sequences lack retroviral genes, but the recent discovery of related MusD retroviral elements with regions similar to gag, pro, and pol suggests that MusD provides the proteins necessary for ETnII transposition in trans. For this study, we analyzed all ETnII elements in the draft sequence of the C57BL/6J genome and classified them into three subtypes (alpha, beta, and gamma) based on structural differences. We then used database searches and quantitative real-time PCR to determine the copy number and expression of ETnII and MusD elements in various mouse strains. In 7.5-day-old embryos of a mouse strain in which two mutations due to ETnII-beta insertions have been identified (SELH/Bc), we detected a three- to sixfold higher level of ETnII-beta and MusD transcripts than in control strains (C57BL/6J and LM/Bc). The increased ETnII transcription level can in part be attributed to a higher number of ETnII-beta elements, but 70% of the MusD transcripts appear to have been derived from one or a few MusD elements that are not detectable in C57BL/6J mice. This element belongs to a young MusD subgroup with intact open reading frames and identical LTRs, suggesting that the overexpressed element(s) in SELH/Bc mice might provide the proteins for the retrotransposition of ETnII and MusD elements. We also show that ETnII is expressed up to 30-fold more than MusD, which could explain why only ETnII, but not MusD, elements have been positively identified as new insertions.