Differential adaptation of the yeast Candida anglica to fermented food.

A single strain of Candida anglica, isolated from cider, is available in international yeast collections. We present here seven new strains isolated from French PDO cheeses. For one of the cheese strains, we achieved a high-quality genome assembly of 13.7 Mb with eight near-complete telomere-to-telomere chromosomes. The genomes of two additional cheese strains and of the cider strain were also assembled and annotated, resulting in a core genome of 5966 coding sequences. Phylogenetic analysis showed that the seven cheese strains clustered together, away from the cider strain. Mating-type locus analysis revealed the presence of a MATa locus in the cider strain but a MATalpha locus in all cheese strains. The presence of LINE retrotransposons at identical genome position in the cheese strains, and two different karyotypic profiles resulting from chromosomal rearrangements were observed. Together, these findings are consistent with clonal propagation of the cheese strains. Phenotypic trait variations were observed within the cheese population under stress conditions whereas the cider strain was found to have a much greater capacity for growth in all conditions tested.

Conserved and divergent DNA recognition specificities and functions of R2 retrotransposon N-terminal domains.

R2 non-long terminal repeat (non-LTR) retrotransposons are among the most extensively distributed mobile genetic elements in multicellular eukaryotes and show promise for applications in transgene supplementation of the human genome. They insert new gene copies into a conserved site in 28S ribosomal DNA with exquisite specificity. R2 clades are defined by the number of zinc fingers (ZFs) at the N terminus of the retrotransposon-encoded protein, postulated to additively confer DNA site specificity. Here, we illuminate general principles of DNA recognition by R2 N-terminal domains across and between clades, with extensive, specific recognition requiring only one or two compact domains. DNA-binding and protection assays demonstrate broadly shared as well as clade-specific DNA interactions. Gene insertion assays in cells identify the N-terminal domains sufficient for target-site insertion and reveal roles in second-strand cleavage or synthesis for clade-specific ZFs. Our results have implications for understanding evolutionary diversification of non-LTR retrotransposon insertion mechanisms and the design of retrotransposon-based gene therapies.

Structural RNA components supervise the sequential DNA cleavage in R2 retrotransposon.

Retroelements are the widespread jumping elements considered as major drivers for genome evolution, which can also be repurposed as gene-editing tools. Here, we determine the cryo-EM structures of eukaryotic R2 retrotransposon with ribosomal DNA target and regulatory RNAs. Combined with biochemical and sequencing analysis, we reveal two essential DNA regions, Drr and Dcr, required for recognition and cleavage. The association of 3' regulatory RNA with R2 protein accelerates the first-strand cleavage, blocks the second-strand cleavage, and initiates the reverse transcription starting from the 3'-tail. Removing 3' regulatory RNA by reverse transcription allows the association of 5' regulatory RNA and initiates the second-strand cleavage. Taken together, our work explains the DNA recognition and RNA supervised sequential retrotransposition mechanisms by R2 machinery, providing insights into the retrotransposon and application reprogramming.

How to start a LINE: 5' switching rejuvenates LINE retrotransposons in tobacco and related Nicotiana species.

By contrast to their conserved mammalian counterparts, plant long interspersed nuclear elements (LINEs) are highly variable, splitting into many low-copy families. Curiously, LINE families from the retrotransposable element (RTE) clade retain a stronger sequence conservation and hence reach higher copy numbers. The cause of this RTE-typical property is not yet understood, but would help clarify why some transposable elements are removed quickly, whereas others persist in plant genomes. Here, we bring forward a detailed study of RTE LINE structure, diversity and evolution in plants. For this, we argue that the nightshade family is the ideal taxon to follow the evolutionary trajectories of RTE LINEs, given their high abundance, recent activity and partnership to non-autonomous elements. Using bioinformatic, cytogenetic and molecular approaches, we detect 4029 full-length RTE LINEs across the Solanaceae. We finely characterize and manually curate a core group of 458 full-length LINEs in allotetraploid tobacco, show an integration event after polyploidization and trace hybridization by RTE LINE composition of parental genomes. Finally, we reveal the role of the untranslated regions (UTRs) as causes for the unique RTE LINE amplification and evolution pattern in plants. On the one hand, we detected a highly conserved motif at the 3' UTR, suggesting strong selective constraints acting on the RTE terminus. On the other hand, we observed successive rounds of 5' UTR cycling, constantly rejuvenating the promoter sequences. This interplay between exchangeable promoters and conserved LINE bodies and 3' UTR likely allows RTE LINEs to persist and thrive in plant genomes.

Pan-3D genome analysis reveals structural and functional differentiation of soybean genomes.

BACKGROUND: High-order chromatin structure plays important roles in gene regulation. However, the diversity of the three-dimensional (3D) genome across plant accessions are seldom reported. RESULTS: Here, we perform the pan-3D genome analysis using Hi-C sequencing data from 27 soybean accessions and comprehensively investigate the relationships between 3D genomic variations and structural variations (SVs) as well as gene expression. We find that intersection regions between A/B compartments largely contribute to compartment divergence. Topologically associating domain (TAD) boundaries in A compartments exhibit significantly higher density compared to those in B compartments. Pan-3D genome analysis shows that core TAD boundaries have the highest transcription start site (TSS) density and lowest GC content and repeat percentage. Further investigation shows that non-long terminal repeat (non-LTR) retrotransposons play important roles in maintaining TAD boundaries, while Gypsy elements and satellite repeats are associated with private TAD boundaries. Moreover, presence and absence variation (PAV) is found to be the major contributor to 3D genome variations. Nevertheless, approximately 55% of 3D genome variations are not associated with obvious genetic variations, and half of them affect the flanking gene expression. In addition, we find that the 3D genome may also undergo selection during soybean domestication. CONCLUSION: Our study sheds light on the role of 3D genomes in plant genetic diversity and provides a valuable resource for studying gene regulation and genome evolution.

Separable structural requirements for cDNA synthesis, nontemplated extension, and template jumping by a non-LTR retroelement reverse transcriptase.

Broad evolutionary expansion of polymerase families has enabled specialization of their activities for distinct cellular roles. In addition to template-complementary synthesis, many polymerases extend their duplex products by nontemplated nucleotide addition (NTA). This activity is exploited for laboratory strategies of cloning and sequencing nucleic acids and could have important biological function, although the latter has been challenging to test without separation-of-function mutations. Several retroelement and retroviral reverse transcriptases (RTs) support NTA and also template jumping, by which the RT performs continuous complementary DNA (cDNA) synthesis using physically separate templates. Previous studies that aimed to dissect the relationship between NTA and template jumping leave open questions about structural requirements for each activity and their interdependence. Here, we characterize the structural requirements for cDNA synthesis, NTA, template jumping, and the unique terminal transferase activity of Bombyx mori R2 non-long terminal repeat retroelement RT. With sequence alignments and structure modeling to guide mutagenesis, we generated enzyme variants across motifs generally conserved or specific to RT subgroups. Enzyme variants had diverse NTA profiles not correlated with other changes in cDNA synthesis activity or template jumping. Using these enzyme variants and panels of activity assay conditions, we show that template jumping requires NTA. However, template jumping by NTA-deficient enzymes can be rescued using primer duplex with a specific length of 3' overhang. Our findings clarify the relationship between NTA and template jumping as well as additional activities of non-long terminal repeat RTs, with implications for the specialization of RT biological functions and laboratory applications.

The non-LTR retrotransposons of Entamoeba histolytica: genomic organization and biology.

Genome sequence analysis of Entamoeba species revealed various classes of transposable elements. While E. histolytica and E. dispar are rich in non-long terminal repeat (LTR) retrotransposons, E. invadens contains predominantly DNA transposons. Non-LTR retrotransposons of E. histolytica constitute three families of long interspersed nuclear elements (LINEs), and their short, nonautonomous partners, SINEs. They occupy ~ 11% of the genome. The EhLINE1/EhSINE1 family is the most abundant and best studied. EhLINE1 is 4.8 kb, with two ORFs that encode functions needed for retrotransposition. ORF1 codes for the nucleic acid-binding protein, and ORF2 has domains for reverse transcriptase (RT) and endonuclease (EN). Most copies of EhLINEs lack complete ORFs. ORF1p is expressed constitutively, but ORF2p is not detected. Retrotransposition could be demonstrated upon ectopic over expression of ORF2p, showing that retrotransposition machinery is functional. The newly retrotransposed sequences showed a high degree of recombination. In transcriptomic analysis, RNA-Seq reads were mapped to individual EhLINE1 copies. Although full-length copies were transcribed, no full-length 4.8 kb transcripts were seen. Rather, sense transcripts mapped to ORF1, RT and EN domains. Intriguingly, there was strong antisense transcription almost exclusively from the RT domain. These unique features of EhLINE1 could serve to attenuate retrotransposition in E. histolytica.

Population analysis of retrotransposons in giraffe genomes supports RTE decline and widespread LINE1 activity in Giraffidae.

BACKGROUND: The majority of structural variation in genomes is caused by insertions of transposable elements (TEs). In mammalian genomes, the main TE fraction is made up of autonomous and non-autonomous non-LTR retrotransposons commonly known as LINEs and SINEs (Long and Short Interspersed Nuclear Elements). Here we present one of the first population-level analysis of TE insertions in a non-model organism, the giraffe. Giraffes are ruminant artiodactyls, one of the few mammalian groups with genomes that are colonized by putatively active LINEs of two different clades of non-LTR retrotransposons, namely the LINE1 and RTE/BovB LINEs as well as their associated SINEs. We analyzed TE insertions of both types, and their associated SINEs in three giraffe genome assemblies, as well as across a population level sampling of 48 individuals covering all extant giraffe species. RESULTS: The comparative genome screen identified 139,525 recent LINE1 and RTE insertions in the sampled giraffe population. The analysis revealed a drastically reduced RTE activity in giraffes, whereas LINE1 is still actively propagating in the genomes of extant (sub)-species. In concert with the extremely low activity of the giraffe RTE, we also found that RTE-dependent SINEs, namely Bov-tA and Bov-A2, have been virtually immobile in the last 2 million years. Despite the high current activity of the giraffe LINE1, we did not find evidence for the presence of currently active LINE1-dependent SINEs. TE insertion heterozygosity rates differ among the different (sub)-species, likely due to divergent population histories. CONCLUSIONS: The horizontally transferred RTE/BovB and its derived SINEs appear to be close to inactivation and subsequent extinction in the genomes of extant giraffe species. This is the first time that the decline of a TE family has been meticulously analyzed from a population genetics perspective. Our study shows how detailed information about past and present TE activity can be obtained by analyzing large-scale population-level genomic data sets.

Low-bias ncRNA libraries using ordered two-template relay: Serial template jumping by a modified retroelement reverse transcriptase.

Selfish, non-long terminal repeat (non-LTR) retroelements and mobile group II introns encode reverse transcriptases (RTs) that can initiate DNA synthesis without substantial base pairing of primer and template. Biochemical characterization of these enzymes has been limited by recombinant expression challenges, hampering understanding of their properties and the possible exploitation of their properties for research and biotechnology. We investigated the activities of representative RTs using a modified non-LTR RT from Bombyx mori and a group II intron RT from Eubacterium rectale Only the non-LTR RT supported robust and serial template jumping, producing one complementary DNA (cDNA) from several templates each copied end to end. We also discovered an unexpected terminal deoxynucleotidyl transferase activity of the RTs that adds nucleotide(s) of choice to 3' ends of single- and/or double-stranded RNA or DNA. Combining these two types of activity with additional insights about nontemplated nucleotide additions to duplexed cDNA product, we developed a streamlined protocol for fusion of next-generation sequencing adaptors to both cDNA ends in a single RT reaction. When benchmarked using a reference pool of microRNAs (miRNAs), library production by Ordered Two-Template Relay (OTTR) using recombinant non-LTR retroelement RT outperformed all commercially available kits and rivaled the low bias of technically demanding home-brew protocols. We applied OTTR to inventory RNAs purified from extracellular vesicles, identifying miRNAs as well as myriad other noncoding RNAs (ncRNAs) and ncRNA fragments. Our results establish the utility of OTTR for automation-friendly, low-bias, end-to-end RNA sequence inventories of complex ncRNA samples.

TBP and SNAP50 transcription factors bind specifically to the Pr77 promoter sequence from trypanosomatid non-LTR retrotransposons.

BACKGROUND: Trypanosomatid genomes are colonized by active and inactive mobile DNA elements, such as LINE, SINE-like, SIDER and DIRE retrotransposons. These elements all share a 77-nucleotide-long sequence at their 5' ends, known as Pr77, which activates transcription, thereby generating abundant unspliced and translatable transcripts. However, transcription factors that mediates this process have still not been reported. METHODS: TATA-binding protein (TBP) and small nuclear RNA-activating protein 50 kDa (SNAP50) recombinant proteins and specific antibodies raised against them were generated. Protein capture assay, electrophoretic mobility-shift assays (EMSA) and EMSA competition assays carried out using these proteins and nuclear proteins of the parasite together to specific DNA sequences used as probes allowed detecting direct interaction of these transcription factors to Pr77 sequence. RESULTS: This study identified TBP and SNAP50 as part of the DNA-protein complex formed by the Pr77 promoter sequence and nuclear proteins of Trypanosoma cruzi. TBP establishes direct and specific contact with the Pr77 sequence, where the DPE and DPE downstream regions are docking sites with preferential binding. TBP binds cooperatively (Hill coefficient = 1.67) to Pr77 and to both strands of the Pr77 sequence, while the conformation of this highly structured sequence is not involved in TBP binding. Direct binding of SNAP50 to the Pr77 sequence is weak and may be mediated by protein-protein interactions through other trypanosomatid nuclear proteins. CONCLUSIONS: Identification of the transcription factors that mediate Pr77 transcription may help to elucidate how these retrotransposons are mobilized within the trypanosomatid genomes and their roles in gene regulation processes in this human parasite.

TGIRT-seq Protocol for the Comprehensive Profiling of Coding and Non-coding RNA Biotypes in Cellular, Extracellular Vesicle, and Plasma RNAs.

High-throughput RNA sequencing (RNA-seq) has extraordinarily advanced our understanding of gene expression and disease etiology, and is a powerful tool for the identification of biomarkers in a wide range of organisms. However, most RNA-seq methods rely on retroviral reverse transcriptases (RTs), enzymes that have inherently low fidelity and processivity, to convert RNAs into cDNAs for sequencing. Here, we describe an RNA-seq protocol using Thermostable Group II Intron Reverse Transcriptases (TGIRTs), which have high fidelity, processivity, and strand-displacement activity, as well as a proficient template-switching activity that enables efficient and seamless RNA-seq adapter addition. By combining these activities, TGIRT-seq enables the simultaneous profiling of all RNA biotypes from small amounts of starting material, with superior RNA-seq metrics, and unprecedented ability to sequence structured RNAs. The TGIRT-seq protocol for Illumina sequencing consists of three steps: (i) addition of a 3' RNA-seq adapter, coupled to the initiation of cDNA synthesis at the 3' end of a target RNA, via template switching from a synthetic adapter RNA/DNA starter duplex; (ii) addition of a 5' RNA-seq adapter, by using thermostable 5' App DNA/RNA ligase to ligate an adapter oligonucleotide to the 3' end of the completed cDNA; (iii) minimal PCR amplification, to add capture sites and indices for Illumina sequencing. TGIRT-seq for the Illumina sequencing platform has been used for comprehensive profiling of coding and non-coding RNAs in ribodepleted, chemically fragmented cellular RNAs, and for the analysis of intact (non-chemically fragmented) cellular, extracellular vesicle (EV), and plasma RNAs, where it yields continuous full-length end-to-end sequences of structured small non-coding RNAs (sncRNAs), including tRNAs, snoRNAs, snRNAs, pre-miRNAs, and full-length excised linear intron (FLEXI) RNAs. Graphic abstract: Figure 1.Overview of the TGIRT-seq protocol for Illumina sequencing.Major steps are: (1) Template switching from a synthetic R2 RNA/R2R DNA starter duplex with a 1-nt 3' DNA overhang (a mixture of A, C, G, and T residues, denoted N) that base pairs to the 3' nucleotide of a target RNA, and upon initiating reverse transcription by adding dNTPs, seamlessly links an R2R adapter to the 5' end of the resulting cDNA; (2) Ligation of an R1R adapter to the 3' end of the completed cDNA; and (3) Minimal PCR amplification with primers that add Illumina capture sites (P5 and P7) and barcode sequences (indices 5 and 7). The index 7 barcode is required, while the index 5 barcode is optional, to provide unique dual indices (UDIs).

Establishment of Quantitative PCR Assays for Active Long Interspersed Nuclear Element-1 Subfamilies in Mice and Applications to the Analysis of Aging-Associated Retrotransposition.

The retrotransposon long interspersed nuclear element-1 (LINE-1) can autonomously increase its copy number within a host genome through the retrotransposition process. LINE-1 is active in the germline and in neural progenitor cells, and its somatic retrotransposition activity has a broad impact on neural development and susceptibility to neuropsychiatric disorders. The method to quantify the genomic copy number of LINE-1 would be important in unraveling the role of retrotransposition, especially in the brain. However, because of the species-specific evolution of LINE-1 sequences, methods for quantifying the copy number should be independently developed. Here, we developed a quantitative PCR (qPCR) assay to measure the copy number of active LINE-1 subfamilies in mice. Using the assay, we investigated aging-associated alterations of LINE-1 copy number in several brain regions in wild-type mice and Polg+/D257A mice as a model for accelerated aging. We found that aged Polg+/D257A mice showed higher levels of the type GfII LINE-1 in the basal ganglia than the wild-type mice did, highlighting the importance of assays that focus on an individual active LINE-1 subfamily.

Chromosomal Locations of a Non-LTR Retrotransposon, Menolird18, in Cucumis melo and Cucumis sativus, and Its Implication on Genome Evolution of Cucumis Species.

Mobile elements are major regulators of genome evolution through their effects on genome size and chromosome structure in higher organisms. Non-long terminal repeat (non-LTR) retrotransposons, one of the subclasses of transposons, are specifically inserted into repetitive DNA sequences. While studies on the insertion of non-LTR retrotransposons into ribosomal RNA genes and other repetitive DNA sequences have been reported in the animal kingdom, studies in the plant kingdom are limited. Here, using FISH, we confirmed that Menolird18, a member of LINE (long interspersed nuclear element) in non-LTR retrotransposons and found in Cucumis melo, was inserted into ITS and ETS (internal and external transcribed spacers) regions of 18S rDNA in melon and cucumber. Beside the 18S rDNA regions, Menolird18 was also detected in all centromeric regions of melon, while it was located at pericentromeric and sub-telomeric regions in cucumber. The fact that FISH signals of Menolird18 were found in centromeric and rDNA regions of mitotic chromosomes suggests that Menolird18 is a rDNA and centromere-specific non-LTR retrotransposon in melon. Our findings are the first report on a non-LTR retrotransposon that is highly conserved in 2 different plant species, melon and cucumber. The clear distinction of chromosomal localization of Menolird18 in melon and cucumber implies that it might have been involved in the evolutionary processes of the melon (C. melo) and cucumber (C. sativus) genomes.

R2 and Non-Site-Specific R2-Like Retrotransposons of the German Cockroach, Blattella germanica.

The structural and functional organization of the ribosomal RNA gene cluster and the full-length R2 non-LTR retrotransposon (integrated into a specific site of 28S ribosomal RNA genes) of the German cockroach, Blattella germanica, is described. A partial sequence of the R2 retrotransposon of the cockroach Rhyparobia maderae is also analyzed. The analysis of previously published next-generation sequencing data from the B. germanica genome reveals a new type of retrotransposon closely related to R2 retrotransposons but with a random distribution in the genome. Phylogenetic analysis reveals that these newly described retrotransposons form a separate clade. It is shown that proteins corresponding to the open reading frames of newly described retrotransposons exhibit unequal structural domains. Within these retrotransposons, a recombination event is described. New mechanism of transposition activity is discussed. The essential structural features of R2 retrotransposons are conserved in cockroaches and are typical of previously described R2 retrotransposons. However, the investigation of the number and frequency of 5'-truncated R2 retrotransposon insertion variants in eight B. germanica populations suggests recent mobile element activity. It is shown that the pattern of 5'-truncated R2 retrotransposon copies can be an informative molecular genetic marker for revealing genetic distances between insect populations.

The linker region of LINEs modulates DNA cleavage and DNA polymerization.

Long interspersed elements (LINEs) replicate by target primed reverse transcription (TPRT). Insertion involves two half reactions. Each half reaction involves DNA cleavage followed by DNA synthesis. The linker region, located just beyond the reverse transcriptase in the LINE open reading frame, contains a set of predicted helices that may form an α-finger, followed by a gag-like zinc-knuckle. Point mutations of moderately conserved amino-acid residues in the presumptive α-finger severely impair the DNA endonuclease and reverse transcriptase activities of the integration reaction during both half reactions. Mutations in the gag-like zinc-knuckle also impair DNA cleavage and DNA synthesis in some instances. Mutations in core residues that presumably disrupt the protein structure of the presumptive α-finger and the gag-like zinc-knuckle lead to a promiscuous DNA endonuclease and protein-nucleic acid complexes that get stuck in the well during analysis. The linker region appears to function as a protein, DNA, and RNA conformational switching area. The linker is used to properly position nucleic acid substrates into the active sites of the reverse transcriptase and of the DNA endonuclease.

Divergence of 3' ends as a driver of short interspersed nuclear element (SINE) evolution in the Salicaceae.

Short interspersed nuclear elements (SINEs) are small, non-autonomous and heterogeneous retrotransposons that are widespread in plants. To explore the amplification dynamics and evolutionary history of SINE populations in representative deciduous tree species, we analyzed the genomes of the six following Salicaceae species: Populus deltoides, Populus euphratica, Populus tremula, Populus tremuloides, Populus trichocarpa, and Salix purpurea. We identified 11 Salicaceae SINE families (SaliS-I to SaliS-XI), comprising 27 077 full-length copies. Most of these families harbor segmental similarities, providing evidence for SINE emergence by reshuffling or heterodimerization. We observed two SINE groups, differing in phylogenetic distribution pattern, similarity and 3' end structure. These groups probably emerged during the 'salicoid duplication' (~65 million years ago) in the Salix-Populus progenitor and during the separation of the genus Salix (45-65 million years ago), respectively. In contrast to conserved 5' start motifs across species and SINE families, the 3' ends are highly variable in sequence and length. This extraordinary 3'-end variability results from mutations in the poly(A) tail, which were fixed by subsequent amplificational bursts. We show that the dissemination of newly evolved 3' ends is accomplished by a displacement of older motifs, leading to various 3'-end subpopulations within the SaliS families.

A non-LTR retrotransposon activates anthocyanin biosynthesis by regulating a MYB transcription factor in Capsicum annuum.

The flavonoid compound anthocyanin is an important plant metabolite with nutritional and aesthetic value as well as anti-oxidative capacity. MYB transcription factors are key regulators of anthocyanin biosynthesis in plants. In pepper (Capsicum annuum), the CaAn2 gene, encoding an R2R3 MYB transcription factor, regulates anthocyanin biosynthesis. However, no functional study or structural analysis of functional and dysfunctional CaAn2 alleles has been performed. Here, to elucidate the function of CaAn2, we generated transgenic Nicotiana benthamiana and Arabidopsis thaliana plants expressing CaAn2. All of the tissues in these plants were purple. Promoter analysis of CaAn2 in purple C. annuum 'KC00134' plants revealed the insertion of a non-long terminal repeat (LTR) retrotransposon designated Ca-nLTR-A. To determine the promoter activity and functional domain of Ca-nLTR-A, various constructs carrying different domains of Ca-nLTR-A fused with GUS were transformed into N. benthamiana. Promoter analysis showed that the 3' untranslated region (UTR) of the second open reading frame of Ca-nLTR-A is responsible for CaAn2 expression in 'KC00134'. Sequence analysis of Ca-nLTR-A identified transcription factor binding sites known to regulate anthocyanin biosynthesis. This study indicates that insertion of a non-LTR retrotransposon in the promoter may activate expression of CaAn2 by recruiting transcription factors at the 3' UTR and thus provides the first example of exaptation of a non-LTR retrotransposon into a new promoter in plants.

The duck EB66® cell substrate reveals a novel retrotransposon.

During the establishment of the duck EB66® cell line as a new cell substrate for vaccine production in industry, a very low level of reverse transcriptase (RT) activity was detected in the culture supernatant by product-enhanced RT assay but a whole battery of tests failed to evidence infectious particles. Results from extensive biochemical and physical testing demonstrated that RT activity was associated to an intracellular, non-enveloped and dense structure different from an infectious retroviruses. In silico analysis of Anas platyrhynchos genome revealed that the most likely candidates for encoding a ribonucleoprotein (RNP)-associated RT were nine copies of chicken repeat 1 (CR1)-like elements, belonging to the non-long terminal repeat retrotransposons. The presence of the full length Anas platyrhynchos chicken repeat 1-like sequence (APCR1) was confirmed in EB66® cells and the related ribonucleic acid was present in the RT-containing fraction of EB66® cells.

Targeted gene knockin in zebrafish using the 28S rDNA-specific non-LTR-retrotransposon R2Ol.

BACKGROUND: Although most of long interspersed elements (LINEs), one class of non-LTR-retrotransposons, are integrated into the host genome randomely, some elements are retrotransposed into the specific sequences of the genomic regions, such as rRNA gene (rDNA) clusters, telomeric repeats and other repetitive sequenes. Most of the sequence-specific LINEs have been reported mainly among invertebrate species and shown to retrotranspose into the specific sequences in vivo and in vitro systems. Recenlty, 28S rDNA-specific LINE R2 elements are shown to be distributed among widespread vertebrate species, but the sequence-specific retrotransposition of R2 has never been demonstrated in vertebrates. RESULTS: Here we cloned a full length unit of R2 from medaka fish Oryzias latipes, named R2Ol, and engineered it to a targeted gene integration tool in zebrafish. By injecting R2Ol-encoding mRNA into zebrafish embryos, R2Ol retrotransposed precisely into the target site at high efficiency (98%) and was transmitted to the next generation at high frequency (50%). We also generated transgenic zebrafish carrying the enhanced green fluorescent protein (EGFP) reporter gene in 28S rDNA target by the R2Ol retrotransposition system. CONCLUSIONS: Sequence-specific LINE retrotransposes into the precise sequence using target primed reverse transcription (TPRT), possibly providing an alternative and effective targeted gene knockin method in vertebrates.

Impact of non-LTR retrotransposons in the differentiation and evolution of anatomically modern humans.

BACKGROUND: Transposable elements are biologically important components of eukaryote genomes. In particular, non-LTR retrotransposons (N-LTRrs) played a key role in shaping the human genome throughout evolution. In this study, we compared retrotransposon insertions differentially present in the genomes of Anatomically Modern Humans, Neanderthals, Denisovans and Chimpanzees, in order to assess the possible impact of retrotransposition in the differentiation of the human lineage. RESULTS: We first identified species-specific N-LTRrs and established their distribution in present day human populations. These analyses shortlisted a group of N-LTRr insertions that were found exclusively in Anatomically Modern Humans. These insertions are associated with an increase in the number of transcriptional/splicing variants of those genes they inserted in. The analysis of the functionality of genes containing human-specific N-LTRr insertions reflects changes that occurred during human evolution. In particular, the expression of genes containing the most recent N-LTRr insertions is enriched in the brain, especially in undifferentiated neurons, and these genes associate in networks related to neuron maturation and migration. Additionally, we identified candidate N-LTRr insertions that have likely produced new functional variants exclusive to modern humans, whose genomic loci show traces of positive selection. CONCLUSIONS: Our results strongly suggest that N-LTRr impacted our differentiation as a species, most likely inducing an increase in neural complexity, and have been a constant source of genomic variability all throughout the evolution of the human lineage.