Full-length Element Research Articles

HiTE: a fast and accurate dynamic boundary adjustment approach for full-length transposable element detection and annotation

Recent advancements in genome assembly have greatly improved the prospects for comprehensive annotation of Transposable Elements (TEs). However, existing methods for TE annotation using genome assemblies suffer from limited accuracy and robustness, requiring extensive manual editing. In addition, the currently available gold-standard TE databases are not comprehensive, even for extensively studied species, highlighting the critical need for an automated TE detection method to supplement existing repositories. In this study, we introduce HiTE, a fast and accurate dynamic boundary adjustment approach designed to detect full-length TEs. The experimental results demonstrate that HiTE outperforms RepeatModeler2, the state-of-the-art tool, across various species. Furthermore, HiTE has identified numerous novel transposons with well-defined structures containing protein-coding domains, some of which are directly inserted within crucial genes, leading to direct alterations in gene expression. A Nextflow version of HiTE is also available, with enhanced parallelism, reproducibility, and portability.

The hAT family hopper transposon exists as highly similar yet discontinuous elements in the Bactrocera tephritid fly genus.

The hAT family transposable element, hopper, was originally discovered as a defective 3120-bp full-length element in a wild-type strain of the oriental fruit fly, Bactrocera dorsalis (Hendel) (Diptera: Tephritidae), and subsequently a functional 3131-bp element, hopperBdwe, was isolated from a white eye mutant strain. The latter study showed that closely related elements exist in melonfly, Zeugodacus cucurbitae (Coquillett) (Diptera: Tephritidae), a closely related subgenus, suggesting that hopper could have a widespread presence in the Bactrocera genus. To further understand the distribution of hopper within and beyond the B. dorsalis species complex, primer pairs from hopperBdwe and its adjacent genomic insertion site were used to survey the presence and relatedness of hopper in five species within the complex and four species beyond the complex. Based on sequence identity of a 1.94 kb internal nucleotide sequence, the closest relationships were with mutated elements from B. dorsalis s.s. and species synonymized with B. dorsalis including B. papayae, B. philippinensis and B. invadens, ranging in identity between 88.4% and 99.5%. Notably, Bactrocera carambolae (Drew & Hancock) (Diptera: Tephritidae), which is most closely related to B. dorsalis beyond the synonymized species, shared hopper identities of 97.3%-99.5%. Beyond the B. dorsalis complex, Z. cucurbitae,Bactrocera tryoni (Froggatt) (Diptera: Tephritidae) and Bactrocera zonata (Saunders) (Diptera: Tephritidae) shared identities of 83.1%-97.1%, while hopper was absent from the Bactrocera oleae (Gmelin) (Diptera: Tephritidae) strain tested. While the functional autonomous hopperBdwe element was not detected in these species, another closely related hopper element isolated from a B. dorsalis genetic sexing strain has an uninterrupted transposase open reading frame. The discontinuous presence of hopper in the Bactrocera genus has implications for its use for genomic manipulation and understanding the phylogenetic relationship of these species.

The discovery of multiple active mys-related LTR-retroelements within the Neotominae subfamily of cricetid rodents.

Retrotransposon families in the rodent family Cricetidae have been understudied in contrast to Muridae, both taxa classified within the superfamily Muroidea. Therefore, we carried out a study to advance our knowledge of the unique mys LTR-retroelement identified in Peromyscus leucopus, by incorporating intra-ORF PCR, quantitative dot blots, DNA and protein library screens, the generation of molecular phylogenies, and analyses of orthologous LTR-retroelement loci. These analyses led to the discovery of three additional related families of LTR-retroelements, which include a 2900bp full-length element of mys-related sequences (mysRS), an 8000bp element containing the mys ORF1 sequence (mORF1) with ERV-related sequences downstream in the reverse orientation, as well as an 1800bp element primarily consisting of mys ORF2 (mORF2) related sequences flanked by LTRs. Our data revealed only a few full-length mys elements among genera of the Neotominae subfamily of cricetid rodents, most existing as partial copies. The mysRS and mORF1 elements are also limited to the genomes of the Neotominae subfamily, whereas mORF2 appears to be restricted to the Peromyscus genus. Molecular phylogenies demonstrating concerted evolution along with an assessment of orthologous loci in Peromyscus for the presence or absence of elements are consistent with activity of these novel LTR-retroelement families within this genus. Together with known activity of various families of non-LTR retroelements in Peromyscus species, we propose that retrotransposons have been continually contributing to the dynamics of Peromyscus genomes promoting genomic diversity and may be correlated with the evolution of more than 50 identified Peromyscus species.

Abstract 1462: Oncogenic activation of FOXR2 driven by somatic acquisition of a LINE-1 promoter in pediatric high-grade glioma

Abstract Long Interspersed Element-1 (LINE-1 or L1) is the only autonomously active mobile genetic element in the human genome. It has been shown to mobilize (i.e., retrotranspose) in cancers, causing de novo somatic L1 insertions. Somatic L1 retrotransposition can cause aberrant splicing of tumor suppressor genes by utilizing cryptic splice donor (SD) sites such as the nucleotide 97 SD in the L1 5’UTR or disrupting repressive regulatory regions of proto-oncogenes, inducing aberrant activation, leading to cancer. In contrast to these reported mechanisms, here we report the first evidence of a driver oncogene activation due to somatic acquisition of an L1 promoter in a pediatric high-grade glioma (HGG). While investigating activation of FOXR2 (forkhead box R2), a known oncogene in CNS embryonal tumors, in a recurrent pediatric diffuse HGG that showed a methylation profile of CNS FOXR2-activated tumors, we identified overexpression of a non-canonical FOXR2 isoform. Transcription of FOXR2 initiated from ~11kb upstream of exon 2 at a site harboring a somatic structural variation (SV) reminiscent of an L1 insertion. Targeted PCR amplification and subsequent PacBio sequencing revealed a 5’ inverted/deleted L1 insertion (~3kb sequence) containing an intact L1 5’UTR, presence of a 3’ transduction sequence, two poly-A tails, and target site duplications. The 3’ transduction sequence matches the sequence downstream an intact polymorphic full-length L1 element at 6p24.1, the likely source L1 element that initiated the FOXR2 somatic insertion. Strikingly, tumor RNA-seq data reveals FOXR2 transcription initiated from the intact sense pol-II promoter of the inserted L1 5’UTR, generating a fusion transcript joining the L1 5’UTR 97 SD to a canonical splice acceptor within exon 2 of FOXR2. Targeted bisulfite sequencing of the FOXR2 L1 5’UTR revealed a hypomethylated state, further supporting its role in FOXR2 transcription initiation. FOXR2 activation is likely a tumor initiating event as the following was observed in the primary tumor that occurred two years prior: a matching recurrent tumor methylation profile, the presence of L1/FOXR2 fusion transcripts, and overexpression of FOXR2. Additional driver mutations, clonal TP53 R175H mutation and PDGFRA amplification, identified in the recurrent tumor but absent in the primary tumor, support the somatic L1 insertion as a tumor initiating event occurring prior to later acquired alterations. Our data represent the first example of oncogenic activation by L1 “promoter donation” via somatic retrotransposition as a novel tumor initiating mechanism. Future large-scale research and development of new computational methods for detecting rearranged L1 insertions are required to determine the prevalence of L1 “promoter donation” in cancer or other diseases. Citation Format: Xiaolong Chen, Diane A. Flasch, Bensheng Ju, Heather L. Mulder, John Easton, Lu Wang, Suzanne J. Baker, Jason Chiang, Jinghui Zhang. Oncogenic activation of FOXR2 driven by somatic acquisition of a LINE-1 promoter in pediatric high-grade glioma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 1462.

Identification of LINE retrotransposons and long non-coding RNAs expressed in the octopus brain

BackgroundTransposable elements (TEs) widely contribute to the evolution of genomes allowing genomic innovations, generating germinal and somatic heterogeneity, and giving birth to long non-coding RNAs (lncRNAs). These features have been associated to the evolution, functioning, and complexity of the nervous system at such a level that somatic retrotransposition of long interspersed element (LINE) L1 has been proposed to be associated to human cognition. Among invertebrates, octopuses are fascinating animals whose nervous system reaches a high level of complexity achieving sophisticated cognitive abilities. The sequencing of the genome of the Octopus bimaculoides revealed a striking expansion of TEs which were proposed to have contributed to the evolution of its complex nervous system. We recently found a similar expansion also in the genome of Octopus vulgaris. However, a specific search for the existence and the transcription of full-length transpositionally competent TEs has not been performed in this genus.ResultsHere, we report the identification of LINE elements competent for retrotransposition in Octopus vulgaris and Octopus bimaculoides and show evidence suggesting that they might be transcribed and determine germline and somatic polymorphisms especially in the brain. Transcription and translation measured for one of these elements resulted in specific signals in neurons belonging to areas associated with behavioral plasticity. We also report the transcription of thousands of lncRNAs and the pervasive inclusion of TE fragments in the transcriptomes of both Octopus species, further testifying the crucial activity of TEs in the evolution of the octopus genomes.ConclusionsThe neural transcriptome of the octopus shows the transcription of thousands of putative lncRNAs and of a full-length LINE element belonging to the RTE class. We speculate that a convergent evolutionary process involving retrotransposons activity in the brain has been important for the evolution of sophisticated cognitive abilities in this genus.

Excision and reinsertion of Ac macrotransposons in maize.

Eukaryotic Macrotransposons (MTns) can be formed by 2 nearby elements flanking a segment of host DNA. The maize Ac transposon can form Ac::MTns, but little is known about Ac::MTn transposition activities. Here, we studied 3 Ac::MTns at the maize p1 locus, each of which is composed of a segment of maize p1 genomic DNA (up to 15 kb) bounded by a fractured Ac element (fAc, 2039 bp), and a full-length Ac element in direct orientation. The resulting Ac::MTns are of 16, 16.5, and 22 kb total length. From these 3 Ac::MTns, we identified 10 independent cases of macrotransposition, and observed similar features of transposition between Ac::MTn and standard Ac/Ds, including characteristic excision footprints and insertion target site duplications. Nine out of the 10 Ac::MTn reinsertion targets were genetically linked to the donor sites, another similarity with Ac/Ds standard transposition. We also identified a MTn-like structure in the maize B73 reference genome and 5 NAM founder lines. The MTn in diverse lines is flanked by target site duplications, confirming the historic occurrence of MTn transposition during genome evolution. Our results show that Ac::MTns are capable of mobilizing segments of DNA long enough to include a typical full-length plant gene and in theory could erode gene colinearity in syntenic regions during plant genome evolution.

Reactivation of transposable elements following hybridization in fission yeast.

Hybridization is thought to reactivate transposable elements (TEs) that were efficiently suppressed in the genomes of the parental hosts. Here, we provide evidence for this “genomic shock hypothesis” in the fission yeast Schizosaccharomyces pombe. In this species, two divergent lineages (Sp and Sk) have experienced recent, likely human-induced, hybridization. We used long-read sequencing data to assemble genomes of 37 samples derived from 31 S. pombe strains spanning a wide range of ancestral admixture proportions. A comprehensive TE inventory revealed exclusive presence of long terminal repeat (LTR) retrotransposons. Sequence analysis of active full-length elements, as well as solo LTRs, revealed a complex history of homologous recombination. Population genetic analyses of syntenic sequences placed insertion of many solo LTRs before the split of the Sp and Sk lineages. Most full-length elements were inserted more recently, after hybridization. With the exception of a single full-length element with signs of positive selection, both solo LTRs and, in particular, full-length elements carry signatures of purifying selection indicating effective removal by the host. Consistent with reactivation upon hybridization, the number of full-length LTR retrotransposons, varying extensively from zero to 87 among strains, significantly increases with the degree of genomic admixture. This study gives a detailed account of global TE diversity in S. pombe, documents complex recombination histories within TE elements, and provides evidence for the “genomic shock hypothesis.”

DIRS retrotransposons amplify via linear, single-stranded cDNA intermediates.

The Dictyostelium Intermediate Repeat Sequence 1 (DIRS-1) is the name-giving member of the DIRS order of tyrosine recombinase retrotransposons. In Dictyostelium discoideum, DIRS-1 is highly amplified and enriched in heterochromatic centromers of the D. discoideum genome. We show here that DIRS-1 it tightly controlled by the D. discoideum RNA interference machinery and is only mobilized in mutants lacking either the RNA dependent RNA polymerase RrpC or the Argonaute protein AgnA. DIRS retrotransposons contain an internal complementary region (ICR) that is thought to be required to reconstitute a full-length element from incomplete RNA transcripts. Using different versions of D. discoideum DIRS-1 equipped with retrotransposition marker genes, we show experimentally that the ICR is in fact essential to complete retrotransposition. We further show that DIRS-1 produces a mixture of single-stranded, mostly linear extrachromosomal cDNA intermediates. If this cDNA is isolated and transformed into D. discoideum cells, it can be used by DIRS-1 proteins to complete productive retrotransposition. This work provides the first experimental evidence to propose a general retrotransposition mechanism of the class of DIRS like tyrosine recombinase retrotransposons.

Transcription of soybean retrotransposon SORE-1 is temporally upregulated in developing ovules.

Transcription of soybean retrotransposon SORE-1 was temporally upregulated during ovule development. This transcriptional pattern was under intrinsic control conferred by the long terminal repeat of SORE-1. Transcriptionally active retrotransposons are capable of inducing random disruption of genes, providing a powerful tool for mutagenesis. Activation of retrotransposons in reproductive cells, in particular, can lead to heritable changes. Here, we examined developmental control of transcription of soybean retrotransposon SORE-1. Transgenic Arabidopsis plants that contain β-glucuronidase (GUS) reporter gene fused with the SORE-1 long terminal repeat (LTR) had GUS staining in the ovule. Quantitative analysis of transcripts in plants with this DNA construct and those with the full-length SORE-1 element indicated a temporal upregulation of SORE-1 transcription during ovule development. A comparable phenomenon was also observed in soybean plants that had a recent insertion of this element in the GmphyA2 gene. These results provide evidence that the temporal upregulation of SORE-1 in the reproductive organ is sufficiently controlled by its LTR and indicate that the intrinsic expression pattern of SORE-1 is consistent with its mutagenic property.

Diverse AR-V7 cistromes in castration-resistant prostate cancer are governed by HoxB13

The constitutively active androgen receptor (AR) splice variant 7 (AR-V7) plays an important role in the progression of castration-resistant prostate cancer (CRPC). Although biomarker studies established the role of AR-V7 in resistance to AR-targeting therapies, how AR-V7 mediates genomic functions in CRPC remains largely unknown. Using a ChIP-exo approach, we show AR-V7 binds to distinct genomic regions and recognizes a full-length androgen-responsive element in CRPC cells and patient tissues. Remarkably, we find dramatic differences in AR-V7 cistromes across diverse CRPC cells and patient tissues, regulating different target gene sets involved in CRPC progression. Surprisingly, we discover that HoxB13 is universally required for and colocalizes with AR-V7 binding to open chromatin across CRPC genomes. HoxB13 pioneers AR-V7 binding through direct physical interaction, and collaborates with AR-V7 to up-regulate target oncogenes. Transcriptional coregulation by HoxB13 and AR-V7 was further supported by their coexpression in tumors and circulating tumor cells from CRPC patients. Importantly, HoxB13 silencing significantly decreases CRPC growth through inhibition of AR-V7 oncogenic function. These results identify HoxB13 as a pivotal upstream regulator of AR-V7-driven transcriptomes that are often cell context-dependent in CRPC, suggesting that HoxB13 may serve as a therapeutic target for AR-V7-driven prostate tumors.

L1 retrotransposition is a common feature of mammalian hepatocarcinogenesis.

The retrotransposon Long Interspersed Element 1 (LINE-1 or L1) is a continuing source of germline and somatic mutagenesis in mammals. Deregulated L1 activity is a hallmark of cancer, and L1 mutagenesis has been described in numerous human malignancies. We previously employed retrotransposon capture sequencing (RC-seq) to analyze hepatocellular carcinoma (HCC) samples from patients infected with hepatitis B or hepatitis C virus and identified L1 variants responsible for activating oncogenic pathways. Here, we have applied RC-seq and whole-genome sequencing (WGS) to an Abcb4 (Mdr2)−/− mouse model of hepatic carcinogenesis and demonstrated for the first time that L1 mobilization occurs in murine tumors. In 12 HCC nodules obtained from 10 animals, we validated four somatic L1 insertions by PCR and capillary sequencing, including TF subfamily elements, and one GF subfamily example. One of the TF insertions carried a 3′ transduction, allowing us to identify its donor L1 and to demonstrate that this full-length TF element retained retrotransposition capacity in cultured cancer cells. Using RC-seq, we also identified eight tumor-specific L1 insertions from 25 HCC patients with a history of alcohol abuse. Finally, we used RC-seq and WGS to identify three tumor-specific L1 insertions among 10 intra-hepatic cholangiocarcinoma (ICC) patients, including one insertion traced to a donor L1 on Chromosome 22 known to be highly active in other cancers. This study reveals L1 mobilization as a common feature of hepatocarcinogenesis in mammals, demonstrating that the phenomenon is not restricted to human viral HCC etiologies and is encountered in murine liver tumors.

Transposable elements in the Anopheles funestus transcriptome.

Transposable elements (TEs) are present in most of the eukaryotic genomes and their impact on genome evolution is increasingly recognized. Although there is extensive information on the TEs present in several eukaryotic genomes, less is known about the expression of these elements at the transcriptome level. Here we present a detailed analysis regarding the expression of TEs in Anopheles funestus, the second most important vector of human malaria in Africa. Several transcriptionally active TE families belonging both to Class I and II were identified and characterized. Interestingly, we have identified a full-length putative active element (including the presence of full length TIRs in the genomic sequence) belonging to the hAT superfamily, which presents active members in other insect genomes. This work contributes to a comprehensive understanding of the landscape of transposable elements in A. funestus transcriptome. Our results reveal that TEs are abundant and diverse in the mosquito and that most of the TE families found in the genome are represented in the mosquito transcriptome, a fact that could indicate activity of these elements.The vast diversity of TEs expressed in A. funestus suggests that there is ongoing amplification of several families in this organism.

Lentiviral Vectors Incorporating Ubiquitous Chromatin Opening Element Driving Canine CD18 Expression

Leukocyte adhesion deficiency type 1 (LAD1) in humans is caused due to mutations in the ITGB2 gene encoding the leukocyte CD18 subunit (b2 integrin). This results in defective leukocyte adhesion and migration leading to recurrent episodes of life-threatening bacterial infection. Canine leukocyte adhesion deficiency (CLAD) represents a disease-specific large animal model of LAD1 in which new therapeutic approaches could be tested. Our previous studies have demonstrated variable efficiency of CD18 expression under the control of several promoters. These include cellular promoters such as those of human elongation factor 1a (hEF1a): long (1169bp) and short (248bp) fragments, human phosphoglycerate kinase (hPGK), human CD11b and human CD18 genes. In addition, murine stem cell virus (MSCV) promoter has also been demonstrated to lead to very high levels of CD18 expression in CLAD CD34+ cells thereby reversing the CLAD phenotype in dogs previously treated with both foamy and lentiviral vectors. But, due to potential genotoxicity associated with the use of viral promoters, we continued our efforts in search of novel cellular promoters. One such promoter is the ubiquitous chromatin opening element (UCOE) from the human heterogeneous ribonucleoprotein A2/B1 and chromobox homolog 3 (HNRPA2B1-CBX3) loci. UCOE has been previously shown to display reproducible and stable transgene expression within the context of a self-inactivating (SIN) lentiviral vector in the absence of classical enhancer activity (Zhang et al., Blood 2007).It has also been shown to confer resistance to DNA methylation-mediated transgene silencing even upon integration into the heterochromatin regions of the host chromosome (Zhang et al., Mol Ther. 2010). Since the full-length element is about 2.6 kb, we cloned and tested different fragment lengths of the UCOE promoter in a SIN lentiviral vector (pCL20) in CLAD CD34+ cells in vitro. Efficiency of expression of CD18 obtained with the six promoter fragments of UCOE (in bp), namely U3'631, U3'1262, U3'652, U5'1357, U5'723 and U5'655 were compared to those obtained with an MSCV promoter. Functional viral titers were first determined using a human LAD EBV-transformed B-cell line that lacks endogenous human CD18. When comparable titers of each vector were used in an overnight transduction of CLAD CD34+ cells after a 24h cytokine prestimulation in vitro, the percentage of CD18+ cells 5 days after transduction were as follows: U3'631 - 8.49%, U3'1262 - 15.9%, U3'652 - 21.3% (tested at MOI 100), U5'1357 - 2.05% (tested at MOI 30), U5'723 - 2.44% (tested at MOI 20), U5'655 - 3.01% (tested at MOI 50) and MSCV - 35.3% (tested at MOI 100). The CD18 expression levels driven by some of these promoter fragments were comparable to those driven by cellular promoters mentioned previously. The UCOE is promising in that it could overcome possible gene silencing effects when used in vivo, unlike promoters such as EF1a and PGK which were largely subjected to post-transcriptional gene silencing with sub-therapeutic levels of CD18 as previously tested in the dog model. Hence, functional correction of the CD18 defect could be achieved with candidate UCOE-incorporating SIN lentiviral vector(s) when used in the treatment of CLAD dogs. DisclosuresNo relevant conflicts of interest to declare.

The Flow of the Gibbon LAVA Element Is Facilitated by the LINE-1 Retrotransposition Machinery.

LINE-Alu-VNTR-Alu-like (LAVA) elements comprise a family of non-autonomous, composite, non-LTR retrotransposons specific to gibbons and may have played a role in the evolution of this lineage. A full-length LAVA element consists of portions of repeats found in most primate genomes: CT-rich, Alu-like, and VNTR regions from the SVA retrotransposon, and portions of the AluSz and L1ME5 elements. To evaluate whether the gibbon genome currently harbors functional LAVA elements capable of mobilization by the endogenous LINE-1 (L1) protein machinery and which LAVA components are important for retrotransposition, we established a trans-mobilization assay in HeLa cells. Specifically, we tested if a full-length member of the older LAVA subfamily C that was isolated from the gibbon genome and named LAVAC, or its components, can be mobilized in the presence of the human L1 protein machinery. We show that L1 proteins mobilize the LAVAC element at frequencies exceeding processed pseudogene formation and human SVAE retrotransposition by > 100-fold and ≥3-fold, respectively. We find that only the SVA-derived portions confer activity, and truncation of the 3′ L1ME5 portion increases retrotransposition rates by at least 100%. Tagged de novo insertions integrated into intronic regions in cell culture, recapitulating findings in the gibbon genome. Finally, we present alternative models for the rise of the LAVA retrotransposon in the gibbon lineage.

Discovery of a new repeat family in the Callithrix jacchus genome.

We identified a novel repeat family, termed Platy-1, in the Callithrix jacchus (common marmoset) genome that arose around the time of the divergence of platyrrhines and catarrhines and established itself as a repeat family in New World monkeys (NWMs). A full-length Platy-1 element is ∼100 bp in length, making it the shortest known short interspersed element (SINE) in primates, and harbors features characteristic of non-LTR retrotransposons. We identified 2268 full-length Platy-1 elements across 62 subfamilies in the common marmoset genome. Our subfamily reconstruction and phylogenetic analyses support Platy-1 propagation throughout the evolution of NWMs in the lineage leading to C. jacchus. Platy-1 appears to have reached its amplification peak in the common ancestor of current day marmosets and has since moderately declined. However, identification of more than 200 Platy-1 elements identical to their respective consensus sequence, and the presence of polymorphic elements within common marmoset populations, suggests ongoing retrotransposition activity. Platy-1, a SINE, appears to have originated from an Alu element, and hence is likely derived from 7SL RNA. Our analyses illustrate the birth of a new repeat family and its propagation dynamics in the lineage leading to the common marmoset over the last 40 million years.

Dead Element Replicating: Degenerate R2 Element Replication and rDNA Genomic Turnover in the Bacillus rossius Stick Insect (Insecta: Phasmida)

R2 is an extensively investigated non-LTR retrotransposon that specifically inserts into the 28S rRNA gene sequences of a wide range of metazoans, disrupting its functionality. During R2 integration, first strand synthesis can be incomplete so that 5’ end deleted copies are occasionally inserted. While active R2 copies repopulate the locus by retrotransposing, the non-functional truncated elements should frequently be eliminated by molecular drive processes leading to the concerted evolution of the rDNA array(s). Although, multiple R2 lineages have been discovered in the genome of many animals, the rDNA of the stick insect Bacillus rossius exhibits a peculiar situation: it harbors both a canonical, functional R2 element (R2Brfun) as well as a full-length but degenerate element (R2Brdeg). An intensive sequencing survey in the present study reveals that all truncated variants in stick insects are present in multiple copies suggesting they were duplicated by unequal recombination. Sequencing results also demonstrate that all R2Brdeg copies are full-length, i. e. they have no associated 5' end deletions, and functional assays indicate they have lost the active ribozyme necessary for R2 RNA maturation. Although it cannot be completely ruled out, it seems unlikely that the degenerate elements replicate via reverse transcription, exploiting the R2Brfun element enzymatic machinery, but rather via genomic amplification of inserted 28S by unequal recombination. That inactive copies (both R2Brdeg or 5'-truncated elements) are not eliminated in a short term in stick insects contrasts with findings for the Drosophila R2, suggesting a widely different management of rDNA loci and a lower efficiency of the molecular drive while achieving the concerted evolution.

Endogenous retrovirus insertion in the KIT oncogene determines white and white spotting in domestic cats.

The Dominant White locus (W) in the domestic cat demonstrates pleiotropic effects exhibiting complete penetrance for absence of coat pigmentation and incomplete penetrance for deafness and iris hypopigmentation. We performed linkage analysis using a pedigree segregating White to identify KIT (Chr. B1) as the feline W locus. Segregation and sequence analysis of the KIT gene in two pedigrees (P1 and P2) revealed the remarkable retrotransposition and evolution of a feline endogenous retrovirus (FERV1) as responsible for two distinct phenotypes of the W locus, Dominant White, and white spotting. A full-length (7125 bp) FERV1 element is associated with white spotting, whereas a FERV1 long terminal repeat (LTR) is associated with all Dominant White individuals. For purposes of statistical analysis, the alternatives of wild-type sequence, FERV1 element, and LTR-only define a triallelic marker. Taking into account pedigree relationships, deafness is genetically linked and associated with this marker; estimated P values for association are in the range of 0.007 to 0.10. The retrotransposition interrupts a DNAase I hypersensitive site in KIT intron 1 that is highly conserved across mammals and was previously demonstrated to regulate temporal and tissue-specific expression of KIT in murine hematopoietic and melanocytic cells. A large-population genetic survey of cats (n = 270), representing 30 cat breeds, supports our findings and demonstrates statistical significance of the FERV1 LTR and full-length element with Dominant White/blue iris (P < 0.0001) and white spotting (P < 0.0001), respectively.

Discovery and analysis of an active long terminal repeat-retrotransposable element in Aspergillus oryzae

Wild-type Aspergillus oryzae RIB40 contains two copies of the AO090005001597 gene. We previously constructed A. oryzae RIB40 strain, RKuAF8B, with multiple chromosomal deletions, in which the AO090005001597 copy number was found to be increased significantly. Sequence analysis indicated that AO090005001597 is part of a putative 6,000-bp retrotransposable element, flanked by two long terminal repeats (LTRs) of 669 bp, with characteristics of retroviruses and retrotransposons, and thus designated AoLTR (A. oryzae LTR-retrotransposable element). AoLTR comprised putative reverse transcriptase, RNase H, and integrase domains. The deduced amino acid sequence alignment of AoLTR showed 94% overall identity with AFLAV, an A. flavus Tf1/sushi retrotransposon. Quantitative real-time RT-PCR showed that AoLTR gene expression was significantly increased in the RKuAF8B, in accordance with the increased copy number. Inverse PCR indicated that the full-length retrotransposable element was randomly integrated into multiple genomic locations. However, no obvious phenotypic changes were associated with the increased AoLTR gene copy number.

Insights into the Loblolly Pine Genome: Characterization of BAC and Fosmid Sequences

Despite their prevalence and importance, the genome sequences of loblolly pine, Norway spruce, and white spruce, three ecologically and economically important conifer species, are just becoming available to the research community. Following the completion of these large assemblies, annotation efforts will be undertaken to characterize the reference sequences. Accurate annotation of these ancient genomes would be aided by a comprehensive repeat library; however, few studies have generated enough sequence to fully evaluate and catalog their non-genic content. In this paper, two sets of loblolly pine genomic sequence, 103 previously assembled BACs and 90,954 newly sequenced and assembled fosmid scaffolds, were analyzed. Together, this sequence represents 280 Mbp (roughly 1% of the loblolly pine genome) and one of the most comprehensive studies of repetitive elements and genes in a gymnosperm species. A combination of homology and de novo methodologies were applied to identify both conserved and novel repeats. Similarity analysis estimated a repetitive content of 27% that included both full and partial elements. When combined with the de novo investigation, the estimate increased to almost 86%. Over 60% of the repetitive sequence consists of full or partial LTR (long terminal repeat) retrotransposons. Through de novo approaches, 6,270 novel, full-length transposable element families and 9,415 sub-families were identified. Among those 6,270 families, 82% were annotated as single-copy. Several of the novel, high-copy families are described here, with the largest, PtPiedmont, comprising 133 full-length copies. In addition to repeats, analysis of the coding region reported 23 full-length eukaryotic orthologous proteins (KOGS) and another 29 novel or orthologous genes. These discoveries, along with other genomic resources, will be used to annotate conifer genomes and address long-standing questions about gymnosperm evolution.

β-Thalassemia Due to Intronic LINE-1 Insertion in theβ-GlobinGene (HBB): Molecular Mechanisms Underlying Reduced Transcript Levels of theβ-GlobinL1Allele

We describe the molecular etiology of β(+)-thalassemia that is caused by the insertion of the full-length transposable element LINE-1 (L1) into the intron-2 of the β-globin gene (HBB). The transcript level of the affected β-globin gene was severely reduced. The remaining transcripts consisted of full-length, correctly processed β-globin mRNA and a minute amount of three aberrantly spliced transcripts with a decreased half-life due to activation of the nonsense-mediated decay pathway. The lower steady-state amount of mRNA produced by the β-globin(L1) allele also resulted from a reduced rate of transcription and decreased production of full-length β-globin primary transcripts. The promoter and enhancer sequences of the β-globin(L1) allele were hypermethylated; however, treatment with a demethylating agent did not restore the impaired transcription. A histone deacetylase inhibitor partially reactivated the β-globin(L1) transcription despite permanent β-globin(L1) promoter CpG methylation. This result indicates that the decreased rate of transcription from the β-globin(L1) allele is associated with an altered chromatin structure. Therefore, the molecular defect caused by intronic L1 insertion in the β-globin gene represents a novel etiology of β-thalassemia.