MIRs are classic, tRNA-derived SINEs that amplified before the mammalian radiation

TAD evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find 14% of all human TAD boundaries to be shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons compared to species-specific boundaries. CRISPR-Cas9 knockouts of an ultraconserved boundary in a mouse model lead to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in the upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations and showcases the functional importance of TAD evolution.

TAD Evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species, and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find that only 14% of all human TAD boundaries are shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons, compared to species-specific boundaries. CRISPR-Cas9 knockouts of two ultraconserved boundaries in mouse models leads to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations, and showcase the functional importance of TAD evolution.

Utility of Circulating Cell-Free DNA in Assessing Microsatellite Instability and Loss of Heterozygosity in Breast Cancer Using Human Identification Approach

The diagnostic and prognostic utility of circulating cell-free DNA (cfDNA) in breast cancer (BC) patients was recently reported. Here, we investigated the use of cfDNA to examine microsatellite instability (MSI) and loss of heterozygosity (LOH) for early BC diagnosis. cfDNA and genomic DNA from 41 female BC patients and 40 healthy controls were quantified using NanoDrop spectrophotometry and real-time PCR. The stability of genomic and cfDNA was assessed using a high-resolution AmpFlSTR MiniFiler human identification kit. Significant increases in cfDNA plasma concentrations were observed in BC patients compared to controls. The genotype distribution of the eight autosomal short tandem repeat (STR) loci D7S820, D13S317, D21S11, D2S1338, D18S51, D16S539, FGA, and CSF1PO were in Hardy-Weinberg equilibrium. Significant differences in the allele frequencies of D7S820 allele-8, D21S11 allele-29, allele-30.2, allele-32.2, and CSF1PO allele-11 were seen between BC patients and controls. LOH and MSI were detected in 36.6% of the cfDNA of patients compared to genomic DNA. This study highlights the utility of plasma-derived cfDNA for earlier, less invasive, and cost-effective cancer diagnosis and molecular stratification. It also highlights the potential value of cfDNA in molecular profiling and biomarkers discovery in precision and forensic medicine.

Suggested Absence of Horizontal Transfer of Retrotransposons between Humans and Domestic Mammal Species

Transposable element sequences are usually vertically inherited but have also spread across taxa via horizontal transfer. Previous investigations of ancient horizontal transfer of transposons have compared consensus sequences, but this method resists detection of recent single or low copy number transfer events. The relationship between humans and domesticated animals represents an opportunity for potential horizontal transfer due to the consistent shared proximity and exposure to parasitic insects, which have been identified as plausible transfer vectors. The relatively short period of extended human-animal contact (tens of thousands of years or less) makes horizontal transfer of transposons between them unlikely. However, the availability of high-quality reference genomes allows individual element comparisons to detect low copy number events. Using pairwise all-versus-all megablast searches of the complete suite of retrotransposons of thirteen domestic animals against human, we searched a total of 27,949,823 individual TEs. Based on manual comparisons of stringently filtered BLAST search results for evidence of vertical inheritance, no plausible instances of HTT were identified. These results indicate that significant recent HTT between humans and domesticated animals has not occurred despite the close proximity, either due to the short timescale, inhospitable recipient genomes, a failure of vector activity, or other factors.

MIR-NATs repress MAPT translation and aid proteostasis in neurodegeneration

The human genome expresses thousands of natural antisense transcripts (NAT) that can regulate epigenetic state, transcription, RNA stability or translation of their overlapping genes1,2. Here we describe MAPT-AS1, a brain-enriched NAT that is conserved in primates and contains an embedded mammalian-wide interspersed repeat (MIR), which represses tau translation by competing for ribosomal RNA pairing with the MAPT mRNA internal ribosome entry site3. MAPT encodes tau, a neuronal intrinsically disordered protein (IDP) that stabilizes axonal microtubules. Hyperphosphorylated, aggregation-prone tau forms the hallmark inclusions of tauopathies4. Mutations in MAPT cause familial frontotemporal dementia, and common variations forming the MAPT H1 haplotype are a significant risk factor in many tauopathies5 and Parkinson's disease. Notably, expression of MAPT-AS1 or minimal essential sequences from MAPT-AS1 (including MIR) reduces-whereas silencing MAPT-AS1 expression increases-neuronal tau levels, and correlate with tau pathology in human brain. Moreover, we identified many additional NATs with embedded MIRs (MIR-NATs), which are overrepresented at coding genes linked to neurodegeneration and/or encoding IDPs, and confirmed MIR-NAT-mediated translational control of one such gene, PLCG1. These results demonstrate a key role for MAPT-AS1 in tauopathies and reveal a potentially broad contribution of MIR-NATs to the tightly controlled translation of IDPs6, with particular relevance for proteostasis in neurodegeneration.

Single-cell transcriptomic analyses of dairy cattle ruminal epithelial cells during weaning

Using the 10× Genomics Chromium Controller, we obtained scRNA-seq data of 5064 and 1372 individual cells from two Holstein calf ruminal epithelial tissues before and after weaning, respectively. We detected six distinct cell clusters, designated their cell types, and reported their marker genes. We then examined these clusters' underlining cell types and relationships by performing cell cycle, pseudotime trajectory, regulatory network, weighted gene co-expression network and gene ontology analyses. By integrating these cell marker genes with Holstein GWAS signals, we found they were enriched for animal production and body conformation traits. Finally, we confirmed their cell identities by comparing them with human and mouse stomach epithelial cells. This study presents an initial effort to implement single-cell transcriptomic analysis in cattle, and demonstrates ruminal tissue epithelial cell types and their developments during weaning, opening the door for new discoveries about tissue/cell type roles in complex traits at single-cell resolution.

Retrotransposons as Drivers of Mammalian Brain Evolution

Retrotransposons, a large and diverse class of transposable elements that are still active in humans, represent a remarkable force of genomic innovation underlying mammalian evolution. Among the features distinguishing mammals from all other vertebrates, the presence of a neocortex with a peculiar neuronal organization, composition and connectivity is perhaps the one that, by affecting the cognitive abilities of mammals, contributed mostly to their evolutionary success. Among mammals, hominids and especially humans display an extraordinarily expanded cortical volume, an enrichment of the repertoire of neural cell types and more elaborate patterns of neuronal connectivity. Retrotransposon-derived sequences have recently been implicated in multiple layers of gene regulation in the brain, from transcriptional and post-transcriptional control to both local and large-scale three-dimensional chromatin organization. Accordingly, an increasing variety of neurodevelopmental and neurodegenerative conditions are being recognized to be associated with retrotransposon dysregulation. We review here a large body of recent studies lending support to the idea that retrotransposon-dependent evolutionary novelties were crucial for the emergence of mammalian, primate and human peculiarities of brain morphology and function.

Investigating the Potential Roles of SINEs in the Human Genome

Short interspersed nuclear elements (SINEs) are nonautonomous retrotransposons that occupy approximately 13% of the human genome. They are transcribed by RNA polymerase III and can be retrotranscribed and inserted back into the genome with the help of other autonomous retroelements. Because they are preferentially located close to or within gene-rich regions, they can regulate gene expression by various mechanisms that act at both the DNA and the RNA levels. In this review, we summarize recent findings on the involvement of SINEs in different types of gene regulation and discuss the potential regulatory functions of SINEs that are in close proximity to genes, Pol III-transcribed SINE RNAs, and embedded SINE sequences within Pol II-transcribed genes in the human genome. These discoveries illustrate how the human genome has exapted some SINEs into functional regulatory elements.

New Ther1-derived SINE Squam3 in scaled reptiles

BACKGROUND

SINEs comprise a significant part of animal genomes and are used to study the evolution of diverse taxa. Despite significant advances in SINE studies in vertebrates and higher eukaryotes in general, their own evolution is poorly understood.

RESULTS

We have discovered and described in detail a new Squam3 SINE specific for scaled reptiles (Squamata). The subfamilies of this SINE demonstrate different distribution in the genomes of squamates, which together with the data on similar SINEs in the tuatara allowed us to propose a scenario of their evolution in the context of reptilian evolution.

CONCLUSIONS

Ancestral SINEs preserved in small numbers in most genomes can give rise to taxa-specific SINE families. Analysis of this aspect of SINEs can shed light on the history and mechanisms of SINE variation in reptilian genomes.

Restriction of an intron size en route to endothermy

Ca2+-insensitive and -sensitive E1 subunits of the 2-oxoglutarate dehydrogenase complex (OGDHC) regulate tissue-specific NADH and ATP supply by mutually exclusive OGDH exons 4a and 4b. Here we show that their splicing is enforced by distant lariat branch points (dBPs) located near the 5' splice site of the intervening intron. dBPs restrict the intron length and prevent transposon insertions, which can introduce or eliminate dBP competitors. The size restriction was imposed by a single dominant dBP in anamniotes that expanded into a conserved constellation of four dBP adenines in amniotes. The amniote clusters exhibit taxon-specific usage of individual dBPs, reflecting accessibility of their extended motifs within a stable RNA hairpin rather than U2 snRNA:dBP base-pairing. The dBP expansion took place in early terrestrial species and was followed by a uridine enrichment of large downstream polypyrimidine tracts in mammals. The dBP-protected megatracts permit reciprocal regulation of exon 4a and 4b by uridine-binding proteins, including TIA-1/TIAR and PUF60, which promote U1 and U2 snRNP recruitment to the 5' splice site and BP, respectively, but do not significantly alter the relative dBP usage. We further show that codons for residues critically contributing to protein binding sites for Ca2+ and other divalent metals confer the exon inclusion order that mirrors the Irving-Williams affinity series, linking the evolution of auxiliary splicing motifs in exons to metallome constraints. Finally, we hypothesize that the dBP-driven selection for Ca2+-dependent ATP provision by E1 facilitated evolution of endothermy by optimizing the aerobic scope in target tissues.

DNA methylation enables transposable element-driven genome expansion

Multicellular eukaryotic genomes show enormous differences in size. A substantial part of this variation is due to the presence of transposable elements (TEs). They contribute significantly to a cell's mass of DNA and have the potential to become involved in host gene control. We argue that the suppression of their activities by methylation of the C-phosphate-G (CpG) dinucleotide in DNA is essential for their long-term accommodation in the host genome and, therefore, to its expansion. An inevitable consequence of cytosine methylation is an increase in C-to-T transition mutations via deamination, which causes CpG loss. Cytosine deamination is often needed for TEs to take on regulatory functions in the host genome. Our study of the whole-genome sequences of 53 organisms showed a positive correlation between the size of a genome and the percentage of TEs it contains, as well as a negative correlation between size and the CpG observed/expected (O/E) ratio in both TEs and the host DNA. TEs are seldom found at promoters and transcription start sites, but they are found more at enhancers, particularly after they have accumulated C-to-T and other mutations. Therefore, the methylation of TE DNA allows for genome expansion and also leads to new opportunities for gene control by TE-based regulatory sites.

Our Conflict with Transposable Elements and Its Implications for Human Disease

Our genome is a historic record of successive invasions of mobile genetic elements. Like other eukaryotes, we have evolved mechanisms to limit their propagation and minimize the functional impact of new insertions. Although these mechanisms are vitally important, they are imperfect, and a handful of retroelement families remain active in modern humans. This review introduces the intrinsic functions of transposons, the tactics employed in their restraint, and the relevance of this conflict to human pathology. The most straightforward examples of disease-causing transposable elements are germline insertions that disrupt a gene and result in a monogenic disease allele. More enigmatic are the abnormal patterns of transposable element expression in disease states. Changes in transposon regulation and cellular responses to their expression have implicated these sequences in diseases as diverse as cancer, autoimmunity, and neurodegeneration. Distinguishing their epiphenomenal from their pathogenic effects may provide wholly new perspectives on our understanding of disease.

Evidence for a novel overlapping coding sequence in POLG initiated at a CUG start codon

BACKGROUND

POLG, located on nuclear chromosome 15, encodes the DNA polymerase γ(Pol γ). Pol γ is responsible for the replication and repair of mitochondrial DNA (mtDNA). Pol γ is the only DNA polymerase found in mitochondria for most animal cells. Mutations in POLG are the most common single-gene cause of diseases of mitochondria and have been mapped over the coding region of the POLG ORF.

RESULTS

Using PhyloCSF to survey alternative reading frames, we found a conserved coding signature in an alternative frame in exons 2 and 3 of POLG, herein referred to as ORF-Y that arose de novo in placental mammals. Using the synplot2 program, synonymous site conservation was found among mammals in the region of the POLG ORF that is overlapped by ORF-Y. Ribosome profiling data revealed that ORF-Y is translated and that initiation likely occurs at a CUG codon. Inspection of an alignment of mammalian sequences containing ORF-Y revealed that the CUG codon has a strong initiation context and that a well-conserved predicted RNA stem-loop begins 14 nucleotides downstream. Such features are associated with enhanced initiation at near-cognate non-AUG codons. Reanalysis of the Kim et al. (2014) draft human proteome dataset yielded two unique peptides that map unambiguously to ORF-Y. An additional conserved uORF, herein referred to as ORF-Z, was also found in exon 2 of POLG. Lastly, we surveyed Clinvar variants that are synonymous with respect to the POLG ORF and found that most of these variants cause amino acid changes in ORF-Y or ORF-Z.

CONCLUSIONS

We provide evidence for a novel coding sequence, ORF-Y, that overlaps the POLG ORF. Ribosome profiling and mass spectrometry data show that ORF-Y is expressed. PhyloCSF and synplot2 analysis show that ORF-Y is subject to strong purifying selection. An abundance of disease-correlated mutations that map to exons 2 and 3 of POLG but also affect ORF-Y provides potential clinical significance to this finding.

Short interspersed nuclear element (SINE)-mediated post-transcriptional effects on human and mouse gene expression: SINE-UP for active duty

Primate-specific Alu short interspersed nuclear elements (SINEs) and rodent-specific B and ID (B/ID) SINEs are non-autonomous and generally non-coding retrotransposons that have been copied and pasted into the respective genomes so as to constitute what is estimated to be a remarkable 13% and 8% of those genomes. In the context of messenger RNAs (mRNAs), those residing within 3'-untranslated regions (3'UTRs) can influence mRNA export from the nucleus to the cytoplasm, mRNA translation and/or mRNA decay via proteins with which they associate either individually or base-paired in cis or in trans with a partially complementary SINE. Each of these influences impinges on the primary function of mRNA, which is to serve as a template for protein synthesis. This review describes how human cells have used 3'UTR Alu elements to mediate post-transcriptional gene regulation and also describes examples of convergent evolution between human and mouse 3'UTR SINEs. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

Structural variation and its potential impact on genome instability: Novel discoveries in the EGFR landscape by long-read sequencing

Structural variation (SV) is typically defined as variation within the human genome that exceeds 50 base pairs (bp). SV may be copy number neutral or it may involve duplications, deletions, and complex rearrangements. Recent studies have shown SV to be associated with many human diseases. However, studies of SV have been challenging due to technological constraints. With the advent of third generation (long-read) sequencing technology, exploration of longer stretches of DNA not easily examined previously has been made possible. In the present study, we utilized third generation (long-read) sequencing techniques to examine SV in the EGFR landscape of four haplotypes derived from two human samples. We analyzed the EGFR gene and its landscape (+/- 500,000 base pairs) using this approach and were able to identify a region of non-coding DNA with over 90% similarity to the most common activating EGFR mutation in non-small cell lung cancer. Based on previously published Alu-element genome instability algorithms, we propose a molecular mechanism to explain how this non-coding region of DNA may be interacting with and impacting the stability of the EGFR gene and potentially generating this cancer-driver gene. By these techniques, we were also able to identify previously hidden structural variation in the four haplotypes and in the human reference genome (hg38). We applied previously published algorithms to compare the relative stabilities of these five different EGFR gene landscape haplotypes to estimate their relative potentials to generate the EGFR exon 19, 15 bp canonical deletion. To our knowledge, the present study is the first to use the differences in genomic architecture between targeted cancer-linked phased haplotypes to estimate their relative potentials to form a common cancer-linked driver mutation.

Cooperative evolution of two different TEs results in lineage-specific novel transcripts in the BLOC1S2 gene

BACKGROUND

The BLOC1S2 gene encodes the multifunctional protein BLOS2, a shared subunit of two lysosomal trafficking complexes: i) biogenesis of lysosome-related organelles complex-1 and i) BLOC-1-related complex. In our previous study, we identified an intriguing unreported transcript of the BLOC1S2 gene that has a novel exon derived from two transposable elements (TEs), MIR and AluSp. To investigate the evolutionary footprint and molecular mechanism of action of this transcript, we performed PCR and RT-PCR experiments and sequencing analyses using genomic DNA and RNA samples from humans and various non-human primates.

RESULTS

The results showed that the MIR element had integrated into the genome of our common ancestor, specifically in the BLOC1S2 gene region, before the radiation of all primate lineages and that the AluSp element had integrated into the genome of our common ancestor, fortunately in the middle of the MIR sequences, after the divergence of Old World monkeys and New World monkeys. The combined MIR and AluSp sequences provide a 3' splice site (AG) and 5' splice site (GT), respectively, and generate the Old World monkey-specific transcripts. Moreover, branch point sequences for the intron removal process are provided by the MIR and AluSp combination.

CONCLUSIONS

We show for the first time that sequential integration into the same location and sequence divergence events of two different TEs generated lineage-specific transcripts through sequence collaboration during primate evolution.

MIR sequences recruit zinc finger protein ZNF768 to expressed genes

Mammalian-wide interspersed repeats (MIRs) are retrotransposed elements of mammalian genomes. Here, we report the specific binding of zinc finger protein ZNF768 to the sequence motif GCTGTGTG (N20) CCTCTCTG in the core region of MIRs. ZNF768 binding is preferentially associated with euchromatin and promoter regions of genes. Binding was observed for genes expressed in a cell type-specific manner in human B cell line Raji and osteosarcoma U2OS cells. Mass spectrometric analysis revealed binding of ZNF768 to Elongator components Elp1, Elp2 and Elp3 and other nuclear factors. The N-terminus of ZNF768 contains a heptad repeat array structurally related to the C-terminal domain (CTD) of RNA polymerase II. This array evolved in placental animals but not marsupials and monotreme species, displays species-specific length variations, and possibly fulfills CTD related functions in gene regulation. We propose that the evolution of MIRs and ZNF768 has extended the repertoire of gene regulatory mechanisms in mammals and that ZNF768 binding is associated with cell type-specific gene expression.

Synteny search identifies carnivore Y chromosome for evolution of male specific genes

The explosive accumulation of mammalian genomes has provided a valuable resource to characterize the evolution of the Y chromosome. Unexpectedly, the Y-chromosome sequence has been characterized in only a small handful of species, with the majority being model organisms. Thus, identification of Y-linked scaffolds from unordered genome sequences is becoming more important. Here, we used a syntenic-based approach to generate the scaffolds of the male-specific region of the Y chromosome (MSY) from the genome sequence of 6 male carnivore species. Our results identified 14, 15, 9, 28, 14 and 11 Y-linked scaffolds in polar bears, pacific walruses, red pandas, cheetahs, ferrets and tigers, covering 1.55 Mbp, 2.62 Mbp, 964 Kb, 1.75 Mb, 2.17 Mbp and 1.84 Mb MSY, respectively. All the candidate Y-linked scaffolds in 3 selected species (red pandas, polar bears and tigers) were successfully verified using polymerase chain reaction. We re-annotated 8 carnivore MSYs including these 6 Y-linked scaffolds and domestic dog and cat MSY; a total of 11 orthologous genes conserved in at least 7 of the 8 carnivores were identified. These 11 Y-linked genes have significantly higher evolutionary rates compared with their X-linked counterparts, indicating less purifying selection for MSY genes. Taken together, our study shows that the approach of synteny search is a reliable and easily affordable strategy to identify Y-linked scaffolds from unordered carnivore genomes and provides a preliminary evolutionary study for carnivore MSY genes.

The Ribosome as a Missing Link in Prebiotic Evolution III: Over-Representation of tRNA- and rRNA-Like Sequences and Plieofunctionality of Ribosome-Related Molecules Argues for the Evolution of Primitive Genomes from Ribosomal RNA Modules

We propose that ribosomal RNA (rRNA) formed the basis of the first cellular genomes, and provide evidence from a review of relevant literature and proteonomic tests. We have proposed previously that the ribosome may represent the vestige of the first self-replicating entity in which rRNAs also functioned as genes that were transcribed into functional messenger RNAs (mRNAs) encoding ribosomal proteins. rRNAs also encoded polymerases to replicate itself and a full complement of the transfer RNAs (tRNAs) required to translate its genes. We explore here a further prediction of our “ribosome-first” theory: the ribosomal genome provided the basis for the first cellular genomes. Modern genomes should therefore contain an unexpectedly large percentage of tRNA- and rRNA-like modules derived from both sense and antisense reading frames, and these should encode non-ribosomal proteins, as well as ribosomal ones with key cell functions. Ribosomal proteins should also have been co-opted by cellular evolution to play extra-ribosomal functions. We review existing literature supporting these predictions. We provide additional, new data demonstrating that rRNA-like sequences occur at significantly higher frequencies than predicted on the basis of mRNA duplications or randomized RNA sequences. These data support our “ribosome-first” theory of cellular evolution.

Comparative Analysis of Immune Cells Reveals a Conserved Regulatory Lexicon

Most well-characterized enhancers are deeply conserved. In contrast, genome-wide comparative studies of steady-state systems showed that only a small fraction of active enhancers are conserved. To better understand conservation of enhancer activity, we used a comparative genomics approach that integrates temporal expression and epigenetic profiles in an innate immune system. We found that gene expression programs diverge among mildly induced genes, while being highly conserved for strongly induced genes. The fraction of conserved enhancers varies greatly across gene expression programs, with induced genes and early-response genes, in particular, being regulated by a higher fraction of conserved enhancers. Clustering of conserved accessible DNA sequences within enhancers resulted in over 60 sequence motifs including motifs for known factors, as well as many with unknown function. We further show that the number of instances of these motifs is a strong predictor of the responsiveness of a gene to pathogen detection.

Liver X Receptor-Binding DNA Motif Associated With Atherosclerosis-Specific DNA Methylation Profiles of Alu Elements and Neighboring CpG Islands

BACKGROUND

The signals that determine atherosclerosis-specific DNA methylation profiles are only partially known. We previously identified a 29-bp DNA motif (differential methylation motif [DMM]) proximal to CpG islands (CGIs) that undergo demethylation in advanced human atheromas. Those data hinted that the DMM docks modifiers of DNA methylation and transcription.

METHODS AND RESULTS

We sought to functionally characterize the DMM. We showed that the DMM overlaps with the RNA polymerase III-binding B box of Alu short interspersed nuclear elements and contains a DR2 nuclear receptor response element. Pointing to a possible functional role for an Alu DMM, CGIs proximal (<100 bp) to near-intact DMM-harboring Alu are significantly less methylated relative to CGIs proximal to degenerate DMM-harboring Alu or to DMM-devoid mammalian-wide interspersed repeat short interspersed nuclear elements in human arteries. As for DMM-binding factors, LXRB (liver X receptor β) binds the DMM in a DR2-dependent fashion, and LXR (liver X receptor) agonists induce significant hypermethylation of the bulk of Alu in THP-1 cells. Furthermore, we describe 3 intergenic long noncoding RNAs that harbor a DMM, are under transcriptional control by LXR agonists, and are differentially expressed between normal and atherosclerotic human aortas. Notably, CGIs adjacent to those long noncoding RNAs tend to be hypomethylated in symptomatic relative to stable human atheromas.

CONCLUSIONS

Collectively, the data suggest that a DMM is associated with 2 distinct methylation states: relatively low methylation of in cis CGIs and Alu element hypermethylation. Based on the known atheroprotective role of LXRs, we propose that LXR agonist-induced Alu hypermethylation, a landmark of atherosclerosis, is a compensatory rather than proatherogenic response.

Human transposable elements in Repbase: genomic footprints from fish to humans

Repbase is a comprehensive database of eukaryotic transposable elements (TEs) and repeat sequences, containing over 1300 human repeat sequences. Recent analyses of these repeat sequences have accumulated evidences for their contribution to human evolution through becoming functional elements, such as protein-coding regions or binding sites of transcriptional regulators. However, resolving the origins of repeat sequences is a challenge, due to their age, divergence, and degradation. Ancient repeats have been continuously classified as TEs by finding similar TEs from other organisms. Here, the most comprehensive picture of human repeat sequences is presented. The human genome contains traces of 10 clades (L1, CR1, L2, Crack, RTE, RTEX, R4, Vingi, Tx1 and Penelope) of non-long terminal repeat (non-LTR) retrotransposons (long interspersed elements, LINEs), 3 types (SINE1/7SL, SINE2/tRNA, and SINE3/5S) of short interspersed elements (SINEs), 1 composite retrotransposon (SVA) family, 5 classes (ERV1, ERV2, ERV3, Gypsy and DIRS) of LTR retrotransposons, and 12 superfamilies (Crypton, Ginger1, Harbinger, hAT, Helitron, Kolobok, Mariner, Merlin, MuDR, P, piggyBac and Transib) of DNA transposons. These TE footprints demonstrate an evolutionary continuum of the human genome.

A tissue-specific promoter derived from a SINE retrotransposon drives biallelic expression of PLAGL1 in human lymphocytes

The imprinted gene PLAGL1 is an important regulator of apoptosis and cell cycle arrest. Loss of its expression has been implicated in tumorigenesis in a range of different cancers, and overexpression during fetal development causes transient neonatal diabetes mellitus (TNDM). PLAGL1 lies within an imprinted region of chromosome 6q24, and monoallelic expression from the major, differentially methylated promoter (P1) occurs in most human tissues. However, in peripheral blood leukocytes, the active promoter (P2) is non-imprinted and drives biallelic transcription. We report here a novel PLAGL1 promoter (P5) derived from the insertion of a primate-specific, MIR3 SINE retrotransposon. P5 is highly utilized in lymphocytes, particularly in T cells, and like P2, directs biallelic transcription. Our results show that it is important to consider P5 in relation to PLAGL1 function in T cells when investigating the dysregulation of this gene.

Transposable elements in cancer

Transposable elements give rise to interspersed repeats, sequences that comprise most of our genomes. These mobile DNAs have been historically underappreciated - both because they have been presumed to be unimportant, and because their high copy number and variability pose unique technical challenges. Neither impediment now seems steadfast. Interest in the human mobilome has never been greater, and methods enabling its study are maturing at a fast pace. This Review describes the activity of transposable elements in human cancers, particularly long interspersed element-1 (LINE-1). LINE-1 sequences are self-propagating, protein-coding retrotransposons, and their activity results in somatically acquired insertions in cancer genomes. Altered expression of transposable elements and animation of genomic LINE-1 sequences appear to be hallmarks of cancer, and can be responsible for driving mutations in tumorigenesis.

Transposable elements in cancer

Transposable elements give rise to interspersed repeats, sequences that comprise most of our genomes. These mobile DNAs have been historically underappreciated - both because they have been presumed to be unimportant, and because their high copy number and variability pose unique technical challenges. Neither impediment now seems steadfast. Interest in the human mobilome has never been greater, and methods enabling its study are maturing at a fast pace. This Review describes the activity of transposable elements in human cancers, particularly long interspersed element-1 (LINE-1). LINE-1 sequences are self-propagating, protein-coding retrotransposons, and their activity results in somatically acquired insertions in cancer genomes. Altered expression of transposable elements and animation of genomic LINE-1 sequences appear to be hallmarks of cancer, and can be responsible for driving mutations in tumorigenesis.

Whole-genome expression analysis of mammalian-wide interspersed repeat elements in human cell lines

With more than 500,000 copies, mammalian-wide interspersed repeats (MIRs), a sub-group of SINEs, represent ∼2.5% of the human genome and one of the most numerous family of potential targets for the RNA polymerase (Pol) III transcription machinery. Since MIR elements ceased to amplify ∼130 myr ago, previous studies primarily focused on their genomic impact, while the issue of their expression has not been extensively addressed. We applied a dedicated bioinformatic pipeline to ENCODE RNA-Seq datasets of seven human cell lines and, for the first time, we were able to define the Pol III-driven MIR transcriptome at single-locus resolution. While the majority of Pol III-transcribed MIR elements are cell-specific, we discovered a small set of ubiquitously transcribed MIRs mapping within Pol II-transcribed genes in antisense orientation that could influence the expression of the overlapping gene. We also identified novel Pol III-transcribed ncRNAs, deriving from transcription of annotated MIR fragments flanked by unique MIR-unrelated sequences, and confirmed the role of Pol III-specific internal promoter elements in MIR transcription. Besides demonstrating widespread transcription at these retrotranspositionally inactive elements in human cells, the ability to profile MIR expression at single-locus resolution will facilitate their study in different cell types and states including pathological alterations.

Identification of RNA Polymerase III-Transcribed SINEs at Single-Locus Resolution from RNA Sequencing Data

Short Interspersed Element (SINE) retrotransposons are one of the most abundant DNA repeat elements in the human genome. They have been found to impact the expression of protein-coding genes, but the possible roles in cell physiology of their noncoding RNAs, generated by RNA polymerase (Pol) III, are just starting to be elucidated. For this reason, Short Interspersed Element (SINE) expression profiling is becoming mandatory to obtain a comprehensive picture of their regulatory roles. However, their repeated nature and frequent location within Pol II-transcribed genes represent a serious obstacle to the identification and quantification of genuine, Pol III-derived SINE transcripts at single-locus resolution on a genomic scale. Among the recent Next Generation Sequencing technologies, only RNA sequencing (RNA-Seq) holds the potential to solve these issues, even though both technical and biological matters need to be taken into account. A bioinformatic pipeline has been recently set up that, by exploiting RNA-seq features and knowledge of SINE transcription mechanisms, allows for easy identification and profiling of transcriptionally active genomic loci which are a source of genuine Pol III SINE transcripts.

Transcriptional interference by RNA polymerase III affects expression of the Polr3e gene

A mammalian interspersed repeat (MIR) nested in antisense orientation within the first intron of the Polr3e gene, encoding an RNA polymerase (Pol) III subunit, is conserved in mammals and highly occupied by Pol III. Here, Yeganeh et al. show that the MIR affects Polr3e expression through transcriptional interference.

Overlapping gene arrangements can potentially contribute to gene expression regulation. A mammalian interspersed repeat (MIR) nested in antisense orientation within the first intron of the Polr3e gene, encoding an RNA polymerase III (Pol III) subunit, is conserved in mammals and highly occupied by Pol III. Using a fluorescence assay, CRISPR/Cas9-mediated deletion of the MIR in mouse embryonic stem cells, and chromatin immunoprecipitation assays, we show that the MIR affects Polr3e expression through transcriptional interference. Our study reveals a mechanism by which a Pol II gene can be regulated at the transcription elongation level by transcription of an embedded antisense Pol III gene.

Short Interspersed Nuclear Element (SINE) Sequences in the Genome of the Human Pathogenic Fungus Aspergillus fumigatus Af293

Novel families of short interspersed nuclear element (SINE) sequences in the human pathogenic fungus Aspergillus fumigatus, clinical isolate Af293, were identified and categorised into tRNA-related and 5S rRNA-related SINEs. Eight predicted tRNA-related SINE families originating from different tRNAs, and nominated as AfuSINE2 sequences, contained target site duplications of short direct repeat sequences (4–14 bp) flanking the elements, an extended tRNA-unrelated region and typical features of RNA polymerase III promoter sequences. The elements ranged in size from 140–493 bp and were present in low copy number in the genome and five out of eight were actively transcribed. One putative tRNA^Arg-derived sequence, AfuSINE2-1a possessed a unique feature of repeated trinucleotide ACT residues at its 3’-terminus. This element was similar in sequence to the I-4_AO element found in A. oryzae and an I-1_AF long nuclear interspersed element-like sequence identified in A. fumigatus Af293. Families of 5S rRNA-related SINE sequences, nominated as AfuSINE3, were also identified and their 5'-5S rRNA-related regions show 50–65% and 60–75% similarity to respectively A. fumigatus 5S rRNAs and SINE3-1_AO found in A. oryzae. A. fumigatus Af293 contains five copies of AfuSINE3 sequences ranging in size from 259–343 bp and two out of five AfuSINE3 sequences were actively transcribed. Investigations on AfuSINE distribution in the fungal genome revealed that the elements are enriched in pericentromeric and subtelomeric regions and inserted within gene-rich regions. We also demonstrated that some, but not all, AfuSINE sequences are targeted by host RNA silencing mechanisms. Finally, we demonstrated that infection of the fungus with mycoviruses had no apparent effects on SINE activity.

Coevolution Theory of the Genetic Code at Age Forty: Pathway to Translation and Synthetic Life

The origins of the components of genetic coding are examined in the present study. Genetic information arose from replicator induction by metabolite in accordance with the metabolic expansion law. Messenger RNA and transfer RNA stemmed from a template for binding the aminoacyl-RNA synthetase ribozymes employed to synthesize peptide prosthetic groups on RNAs in the Peptidated RNA World. Coevolution of the genetic code with amino acid biosynthesis generated tRNA paralogs that identify a last universal common ancestor (LUCA) of extant life close to Methanopyrus, which in turn points to archaeal tRNA introns as the most primitive introns and the anticodon usage of Methanopyrus as an ancient mode of wobble. The prediction of the coevolution theory of the genetic code that the code should be a mutable code has led to the isolation of optional and mandatory synthetic life forms with altered protein alphabets.

MetaSINEs: Broad Distribution of a Novel SINE Superfamily in Animals

SINEs (short interspersed elements) are transposable elements that typically originate independently in each taxonomic clade (order/family). However, some SINE families share a highly similar central sequence and are thus categorized as a SINE superfamily. Although only four SINE superfamilies (CORE-SINEs, V-SINEs, DeuSINEs, and Ceph-SINEs) have been reported so far, it is expected that new SINE superfamilies would be discovered by deep exploration of new SINEs in metazoan genomes. Here we describe 15 SINEs, among which 13 are novel, that have a similar 66-bp central region and therefore constitute a new SINE superfamily, MetaSINEs. MetaSINEs are distributed from fish to cnidarians, suggesting their common evolutionary origin at least 640 Ma. Because the 3′ tails of MetaSINEs are variable, these SINEs most likely survived by changing their partner long interspersed elements for retrotransposition during evolution. Furthermore, we examined the presence of members of other SINE superfamilies in bivalve genomes and characterized eight new SINEs belonging to the CORE-SINEs, V-SINEs, and DeuSINEs, in addition to the MetaSINEs. The broad distribution of bivalve SINEs suggests that at least three SINEs originated in the common ancestor of Bivalvia. Our comparative analysis of the central domains of the SINEs revealed that, in each superfamily, only a restricted region is shared among all of its members. Because the functions of the central domains of the SINE superfamilies remain unknown, such structural information of SINE superfamilies will be useful for future experimental and comparative analyses to reveal why they have been retained in metazoan genomes during evolution.

The 3D folding of metazoan genomes correlates with the association of similar repetitive elements

The potential roles of the numerous repetitive elements found in the genomes of multi-cellular organisms remain speculative. Several studies have suggested a role in stabilizing specific 3D genomic contacts. To test this hypothesis, we exploited inter-chromosomal contacts frequencies obtained from Hi-C experiments and show that the folding of the human, mouse and Drosophila genomes is associated with a significant co-localization of several specific repetitive elements, notably many elements of the SINE family. These repeats tend to be the oldest ones and are enriched in transcription factor binding sites. We propose that the co-localization of these repetitive elements may explain the global conservation of genome folding observed between homologous regions of the human and mouse genome. Taken together, these results support a contribution of specific repetitive elements in maintaining and/or reshaping genome architecture over evolutionary times.

The Dfam database of repetitive DNA families

Repetitive DNA, especially that due to transposable elements (TEs), makes up a large fraction of many genomes. Dfam is an open access database of families of repetitive DNA elements, in which each family is represented by a multiple sequence alignment and a profile hidden Markov model (HMM). The initial release of Dfam, featured in the 2013 NAR Database Issue, contained 1143 families of repetitive elements found in humans, and was used to produce more than 100 Mb of additional annotation of TE-derived regions in the human genome, with improved speed. Here, we describe recent advances, most notably expansion to 4150 total families including a comprehensive set of known repeat families from four new organisms (mouse, zebrafish, fly and nematode). We describe improvements to coverage, and to our methods for identifying and reducing false annotation. We also describe updates to the website interface. The Dfam website has moved to http://dfam.org. Seed alignments, profile HMMs, hit lists and other underlying data are available for download.

MIR retrotransposon sequences provide insulators to the human genome

Insulators are regulatory elements that help to organize eukaryotic chromatin via enhancer-blocking and chromatin barrier activity. Although there are several examples of transposable element (TE)-derived insulators, the contribution of TEs to human insulators has not been systematically explored. Mammalian-wide interspersed repeats (MIRs) are a conserved family of TEs that have substantial regulatory capacity and share sequence characteristics with tRNA-related insulators. We sought to evaluate whether MIRs can serve as insulators in the human genome. We applied a bioinformatic screen using genome sequence and functional genomic data from CD4(+) T cells to identify a set of 1,178 predicted MIR insulators genome-wide. These predicted MIR insulators were computationally tested to serve as chromatin barriers and regulators of gene expression in CD4(+) T cells. The activity of predicted MIR insulators was experimentally validated using in vitro and in vivo enhancer-blocking assays. MIR insulators are enriched around genes of the T-cell receptor pathway and reside at T-cell-specific boundaries of repressive and active chromatin. A total of 58% of the MIR insulators predicted here show evidence of T-cell-specific chromatin barrier and gene regulatory activity. MIR insulators appear to be CCCTC-binding factor (CTCF) independent and show a distinct local chromatin environment with marked peaks for RNA Pol III and a number of histone modifications, suggesting that MIR insulators recruit transcriptional complexes and chromatin modifying enzymes in situ to help establish chromatin and regulatory domains in the human genome. The provisioning of insulators by MIRs across the human genome suggests a specific mechanism by which TE sequences can be used to modulate gene regulatory networks.

Genome-Wide Search Identifies 1.9 Mb from the Polar Bear Y Chromosome for Evolutionary Analyses

The male-inherited Y chromosome is the major haploid fraction of the mammalian genome, rendering Y-linked sequences an indispensable resource for evolutionary research. However, despite recent large-scale genome sequencing approaches, only a handful of Y chromosome sequences have been characterized to date, mainly in model organisms. Using polar bear (Ursus maritimus) genomes, we compare two different in silico approaches to identify Y-linked sequences: 1) Similarity to known Y-linked genes and 2) difference in the average read depth of autosomal versus sex chromosomal scaffolds. Specifically, we mapped available genomic sequencing short reads from a male and a female polar bear against the reference genome and identify 112 Y-chromosomal scaffolds with a combined length of 1.9 Mb. We verified the in silico findings for the longer polar bear scaffolds by male-specific in vitro amplification, demonstrating the reliability of the average read depth approach. The obtained Y chromosome sequences contain protein-coding sequences, single nucleotide polymorphisms, microsatellites, and transposable elements that are useful for evolutionary studies. A high-resolution phylogeny of the polar bear patriline shows two highly divergent Y chromosome lineages, obtained from analysis of the identified Y scaffolds in 12 previously published male polar bear genomes. Moreover, we find evidence of gene conversion among ZFX and ZFY sequences in the giant panda lineage and in the ancestor of ursine and tremarctine bears. Thus, the identification of Y-linked scaffold sequences from unordered genome sequences yields valuable data to infer phylogenomic and population-genomic patterns in bears.

Interspersed repeats in the horse (Equus caballus); spatial correlations highlight conserved chromosomal domains

The interspersed repeat content of mammalian genomes has been best characterized in human, mouse and cow. In this study, we carried out de novo identification of repeated elements in the equine genome and identified previously unknown elements present at low copy number. The equine genome contains typical eutherian mammal repeats, but also has a significant number of hybrid repeats in addition to clade-specific Long Interspersed Nuclear Elements (LINE). Equus caballus clade specific LINE 1 (L1) repeats can be classified into approximately five subfamilies, three of which have undergone significant expansion. There are 1115 full-length copies of these equine L1, but of the 103 presumptive active copies, 93 fall within a single subfamily, indicating a rapid recent expansion of this subfamily. We also analysed both interspersed and simple sequence repeats (SSR) genome-wide, finding that some repeat classes are spatially correlated with each other as well as with G+C content and gene density. Based on these spatial correlations, we have confirmed that recently-described ancestral vs. clade-specific genome territories can be defined by their repeat content. The clade-specific Short Interspersed Nuclear Element correlations were scattered over the genome and appear to have been extensively remodelled. In contrast, territories enriched for ancestral repeats tended to be contiguous domains. To determine if the latter territories were evolutionarily conserved, we compared these results with a similar analysis of the human genome, and observed similar ancestral repeat enriched domains. These results indicate that ancestral, evolutionarily conserved mammalian genome territories can be identified on the basis of repeat content alone. Interspersed repeats of different ages appear to be analogous to geologic strata, allowing identification of ancient vs. newly remodelled regions of mammalian genomes.

Ancient traces of tailless retropseudogenes in therian genomes

Transposable elements, once described by Barbara McClintock as controlling genetic units, not only occupy the largest part of our genome but are also a prominent moving force of genomic plasticity and innovation. They usually replicate and reintegrate into genomes silently, sometimes causing malfunctions or misregulations, but occasionally millions of years later, a few may evolve into new functional units. Retrotransposons make their way into the genome following reverse transcription of RNA molecules and chromosomal insertion. In therian mammals, long interspersed elements 1 (LINE1s) self-propagate but also coretropose many RNAs, including mRNAs and small RNAs that usually exhibit an oligo(A) tail. The revitalization of specific LINE1 elements in the mammalian lineage about 150 Ma parallels the rise of many other nonautonomous mobilized genomic elements. We previously identified and described hundreds of tRNA-derived retropseudogenes missing characteristic oligo(A) tails consequently termed tailless retropseudogenes. Additional analyses now revealed hundreds of thousands of tailless retropseudogenes derived from nearly all types of RNAs. We extracted 2,402 perfect tailless sequences (with discernible flanking target site duplications) originating from tRNAs, spliceosomal RNAs, 5S rRNAs, 7SK RNAs, mRNAs, and others. Interestingly, all are truncated at one or more defined positions that coincide with internal single-stranded regions. 5S ribosomal and U2 spliceosomal RNAs were analyzed in the context of mammalian phylogeny to discern the origin of the therian LINE1 retropositional system that evolved in our 150-Myr-old ancestor.

Identification of a novel large deletion in a patient with severe factor V deficiency using an in-house F5 MLPA assay

Factor V (FV) deficiency is a rare autosomal recessive bleeding disorder caused by mutations in the F5 gene. FV-deficient patients in whom no mutation or only one mutation is found may harbour large gene rearrangements, which are not detected by conventional mutation screening strategies. The aim of this study was to develop and validate a multiplex ligation-dependent probe amplification (MLPA) assay for the detection of large deletions and duplications in the F5 gene. Twenty-two MLPA probes targeting 19 of the 25 exons and the upstream and downstream regions of the F5 gene were designed and tested in 10 normal controls, a patient with a known heterozygous deletion of F5 exons 1-7 (positive control) and 14 genetically unexplained FV-deficient patients. MLPA results were confirmed by digital PCR on a QuantStudio(™) 3D Digital PCR System. The F5-specific probes yielded a reproducible peak profile in normal controls, correctly detected the known deletion in the positive control and suggested the presence of a novel deletion of exons 9-10 in a patient with undetectable FV levels and only one identified mutation. Follow-up by chip-based digital PCR, long-range PCR and direct sequencing confirmed that this patient carried a heterozygous F5 deletion of 1823 bp extending from intron 8 to intron 10. Bioinformatics sequence analysis pinpointed repetitive elements that might have originated the deletion. In conclusion, we have developed and validated an MLPA assay for the detection of gross F5 gene rearrangements. This assay may represent a valuable tool for the molecular diagnosis of FV deficiency.

Mammalian-wide interspersed repeat (MIR)-derived enhancers and the regulation of human gene expression

Background

Mammalian-wide interspersed repeats (MIRs) are the most ancient family of transposable elements (TEs) in the human genome. The deep conservation of MIRs initially suggested the possibility that they had been exapted to play functional roles for their host genomes. MIRs also happen to be the only TEs whose presence in-and-around human genes is positively correlated to tissue-specific gene expression. Similar associations of enhancer prevalence within genes and tissue-specific expression, along with MIRs’ previous implication as providing regulatory sequences, suggested a possible link between MIRs and enhancers.

Results

To test the possibility that MIRs contribute functional enhancers to the human genome, we evaluated the relationship between MIRs and human tissue-specific enhancers in terms of genomic location, chromatin environment, regulatory function, and mechanistic attributes. This analysis revealed MIRs to be highly concentrated in enhancers of the K562 and HeLa human cell-types. Significantly more enhancers were found to be linked to MIRs than would be expected by chance, and putative MIR-derived enhancers are characterized by a chromatin environment highly similar to that of canonical enhancers. MIR-derived enhancers show strong associations with gene expression levels, tissue-specific gene expression and tissue-specific cellular functions, including a number of biological processes related to erythropoiesis. MIR-derived enhancers were found to be a rich source of transcription factor binding sites, underscoring one possible mechanistic route for the element sequences co-option as enhancers. There is also tentative evidence to suggest that MIR-enhancer function is related to the transcriptional activity of non-coding RNAs.

Conclusions

Taken together, these data reveal enhancers to be an important cis-regulatory platform from which MIRs can exercise a regulatory function in the human genome and help to resolve a long-standing conundrum as to the reason for MIRs’ deep evolutionary conservation.

Evolution of Alu elements toward enhancers

The human genome contains approximately one million Alu repetitive elements comprising 10% of the genome, yet their functions are not well understood. Here, we show that Alu elements resemble enhancers. Alu elements are bound by two well-phased nucleosomes that contain histones bearing marks of active chromatin, and they show tissue-specific enrichment for the enhancer mark H3K4me1. A proportion of Alu elements were experimentally validated as bona fide active enhancers with an in vitro reporter assay. In addition, Hi-C data indicate that Alus show long-range interactions with gene promoters. We also find that Alus are generally more conserved when located in the proximal upstream region of genes. Their similarity to enhancers becomes more prominent with their age in the human genome, following a clear evolutionary continuum reminiscent of the evolutionary pattern of proto-genes. Therefore, we conclude that some Alu elements can function as enhancers and propose that many more may be proto-enhancers that serve as a repertoire for the de novo birth of enhancers.

A fast-evolving human NPAS3 enhancer gained reporter expression in the developing forebrain of transgenic mice

The developmental brain gene NPAS3 stands out as a hot spot in human evolution because it contains the largest number of human-specific, fast-evolving, conserved, non-coding elements. In this paper we studied 2xHAR142, one of these elements that is located in the fifth intron of NPAS3. Using transgenic mice, we show that the mouse and chimp 2xHAR142 orthologues behave as transcriptional enhancers driving expression of the reporter gene lacZ to a similar NPAS3 expression subdomain in the mouse central nervous system. Interestingly, the human 2xHAR142 orthologue drives lacZ expression to an extended expression pattern in the nervous system. Thus, molecular evolution of 2xHAR142 provides the first documented example of human-specific heterotopy in the forebrain promoted by a transcriptional enhancer and suggests that it may have contributed to assemble the unique properties of the human brain.

A novel application of ecological analyses to assess transposable element distributions in the genome of the domestic cow, Bos taurus

Transposable elements (TEs) are among the most abundant components of many eukaryotic genomes. Efforts to explain TE abundance, as well as TE diversity among genomes, have led some researchers to draw an analogy between genomic and ecological processes. Adopting this perspective, we conducted an analysis of the cow (Bos taurus) genome using techniques developed by community ecologists to determine whether environmental factors influence community composition. Specifically, each chromosome within the Bos taurus genome was treated as a "linear transect", and a multivariate redundancy analysis (RDA) was used to identify large-scale spatial patterns in TE communities associated with 10 TE families. The position of each TE community on the chromosome accounted for ∼50% of the variation along the chromosome "transect". Multivariate analysis further revealed an effect of gene density on TE communities that is influenced by several other factors in the (genomic) environment, including chromosome length and TE density. The results of this analysis demonstrate that ecological methods can be applied successfully to help answer genomic questions.

RNA-Mediated Gene Duplication and Retroposons: Retrogenes, LINEs, SINEs, and Sequence Specificity

A substantial number of “retrogenes” that are derived from the mRNA of various intron-containing genes have been reported. A class of mammalian retroposons, long interspersed element-1 (LINE1, L1), has been shown to be involved in the reverse transcription of retrogenes (or processed pseudogenes) and non-autonomous short interspersed elements (SINEs). The 3′-end sequences of various SINEs originated from a corresponding LINE. As the 3′-untranslated regions of several LINEs are essential for retroposition, these LINEs presumably require “stringent” recognition of the 3′-end sequence of the RNA template. However, the 3′-ends of mammalian L1s do not exhibit any similarity to SINEs, except for the presence of 3′-poly(A) repeats. Since the 3′-poly(A) repeats of L1 and Alu SINE are critical for their retroposition, L1 probably recognizes the poly(A) repeats, thereby mobilizing not only Alu SINE but also cytosolic mRNA. Many flowering plants only harbor L1-clade LINEs and a significant number of SINEs with poly(A) repeats, but no homology to the LINEs. Moreover, processed pseudogenes have also been found in flowering plants. I propose that the ancestral L1-clade LINE in the common ancestor of green plants may have recognized a specific RNA template, with stringent recognition then becoming relaxed during the course of plant evolution.

Control of myogenesis by rodent SINE-containing lncRNAs

Staufen1-mediated mRNA decay (SMD) degrades mRNAs that harbor a Staufen1-binding site (SBS) in their 3' untranslated regions (UTRs). Human SBSs can form by intermolecular base-pairing between a 3' UTR Alu element and an Alu element within a long noncoding RNA (lncRNA) called a &frac12;-sbsRNA. Since Alu elements are confined to primates, it was unclear how SMD occurs in rodents. Here we identify mouse mRNA 3' UTRs and lncRNAs that contain a B1, B2, B4, or identifier (ID) element. We show that SMD occurs in mouse cells via mRNA-lncRNA base-pairing of partially complementary elements and that mouse &frac12;-sbsRNA (m&frac12;-sbsRNA)-triggered SMD regulates C2C12 cell myogenesis. Our findings define new roles for lncRNAs as well as B and ID short interspersed elements (SINEs) in mice that undoubtedly influence many developmental and homeostatic pathways.

SINEBase: a database and tool for SINE analysis

SINEBase (http://sines.eimb.ru) integrates the revisited body of knowledge about short interspersed elements (SINEs). A set of formal definitions concerning SINEs was introduced. All available sequence data were screened through these definitions and the genetic elements misidentified as SINEs were discarded. As a result, 175 SINE families have been recognized in animals, flowering plants and green algae. These families were classified by the modular structure of their nucleotide sequences and the frequencies of different patterns were evaluated. These data formed the basis for the database of SINEs. The SINEBase website can be used in two ways: first, to explore the database of SINE families, and second, to analyse candidate SINE sequences using specifically developed tools. This article presents an overview of the database and the process of SINE identification and analysis.

Extrathymic generation of regulatory T cells in placental mammals mitigates maternal-fetal conflict

Regulatory T (Treg) cells, whose differentiation and function are controlled by X chromosome-encoded transcription factor Foxp3, are generated in the thymus (tTreg) and extrathymically (peripheral, pTreg), and their deficiency results in fatal autoimmunity. Here, we demonstrate that a Foxp3 enhancer, conserved noncoding sequence 1 (CNS1), essential for pTreg but dispensable for tTreg cell generation, is present only in placental mammals. CNS1 is largely composed of mammalian-wide interspersed repeats (MIR) that have undergone retrotransposition during early mammalian radiation. During pregnancy, pTreg cells specific to a model paternal alloantigen were generated in a CNS1-dependent manner and accumulated in the placenta. Furthermore, when mated with allogeneic, but not syngeneic, males, CNS1-deficient females showed increased fetal resorption accompanied by increased immune cell infiltration and defective remodeling of spiral arteries. Our results suggest that, during evolution, a CNS1-dependent mechanism of extrathymic differentiation of Treg cells emerged in placental animals to enforce maternal-fetal tolerance.