Search

Alu are high copy number interspersed repeats that have accumulated near genes during primate and human evolution. They are a pervasive source of structural variation in modern humans. Impacts that Alu insertions may have on gene expression are not well understood, although some have been associated with expression quantitative trait loci (eQTLs). Here, we directly test regulatory effects of polymorphic Alu insertions in isolation of other variants on the same haplotype. To screen insertion variants for those with such effects, we used ectopic luciferase reporter assays and evaluated 110 Alu insertion variants, including more than 40 with a potential role in disease risk. We observed a continuum of effects with significant outliers that up- or down-regulate luciferase activity. Using a series of reporter constructs, which included genomic context surrounding the Alu, we can distinguish between instances in which the Alu disrupts another regulator and those in which the Alu introduces new regulatory sequence. We next focused on three polymorphic Alu loci associated with breast cancer that display significant effects in the reporter assay. We used CRISPR to modify the endogenous sequences, establishing cell lines varying in the Alu genotype. Our findings indicate that Alu genotype can alter expression of genes implicated in cancer risk, including PTHLH, RANBP9, and MYC. These data show that commonly occurring polymorphic Alu elements can alter transcript levels and potentially contribute to disease risk.

Background Plasmodium interspersed repeat (pir) is the largest multigene family in the genomes of most Plasmodium species. A variety of functions for the PIR proteins which they encode have been proposed, including antigenic variation, immune evasion, sequestration and rosetting. However, direct evidence for these is lacking. The repetitive nature of the family has made it difficult to determine function experimentally. However, there has been some success in using gene expression studies to suggest roles for some members in virulence and chronic infection. Methods Here pir gene expression was examined across the life cycle of Plasmodium berghei using publicly available RNAseq data-sets, and at high resolution in the intraerythrocytic development cycle using new data from Plasmodium chabaudi. Results Expression of pir genes is greatest in stages of the parasite which invade and reside in red blood cells. The marked exception is that liver merozoites and male gametocytes produce a very large number of pir gene transcripts, notably compared to female gametocytes, which produce relatively few. Within the asexual blood stages different subfamilies peak at different times, suggesting further functional distinctions. Representing a subfamily of its own, the highly conserved ancestral pir gene warrants further investigation due to its potential tractability for functional investigation. It is highly transcribed in multiple life cycle stages and across most studied Plasmodium species and thus is likely to play an important role in parasite biology. Conclusions The identification of distinct expression patterns for different pir genes and subfamilies is likely to provide a basis for the design of future experiments to uncover their function.

Background Transposable element (TE) sequences are classified into families based on the reconstructed history of replication, and into subfamilies based on more fine-grained features that are often intended to capture family history. We evaluate the reliability of annotation with common subfamilies by assessing the extent to which subfamily annotation is reproducible in replicate copies created by segmental duplications in the human genome, and in homologous copies shared by human and chimpanzee. Results We find that standard methods annotate over 10% of replicates as belonging to different subfamilies, despite the fact that they are expected to be annotated as belonging to the same subfamily. Point mutations and homologous recombination appear to be responsible for some of this discordant annotation (particularly in the young Alu family), but are unlikely to fully explain the annotation unreliability. Conclusions The surprisingly high level of disagreement in subfamily annotation of homologous sequences highlights a need for further research into definition of TE subfamilies, methods for representing subfamily annotation confidence of TE instances, and approaches to better utilizing such nuanced annotation data in downstream analysis.

While near-cognate codons are frequently used for translation initiation in eukaryotes, their efficiencies are usually low (

DNA hypomethylation of long interspersed repetitive DNA retrotransposon (LINE-1) and Alu repeats elements of short interspersed elements family (SINEs) is an early event in carcinogenesis that causes transcriptional activation and leads to chromosomal instability. In the current study, DNA methylation levels of LINE-1 and Alu repeats were analyzed in tumoral tissues of invasive breast cancer in a Tunisian cohort and its association with the clinicopathological features of patients was defined. DNA methylation of LINE-1 and Alu repeats were analyzed using pyrosequencing in 61 invasive breast cancers. Median values observed for DNA methylation of LINE-1 and Alu repeats were considered as the cut-off (59.81 and 18.49%, respectively). The results of the current study demonstrated a positive correlation between DNA methylation levels of LINE-1 and Alu repeats (rho=0.284; P

Despite a vast expansion in the availability of epigenomic data, our knowledge of the chromatin landscape at interspersed repeats remains highly limited by difficulties in mapping short-read sequencing data to these regions. In particular, little is known about the locus-specific regulation of evolutionarily young transposable elements (TEs), which have been implicated in genome stability, gene regulation and innate immunity in a variety of developmental and disease contexts. Here we propose an approach for generating locus-specific protein–DNA binding profiles at interspersed repeats, which leverages information on the spatial proximity between repetitive and non-repetitive genomic regions. We demonstrate that the combination of HiChIP and a newly developed mapping tool (PAtChER) yields accurate protein enrichment profiles at individual repetitive loci. Using this approach, we reveal previously unappreciated variation in the epigenetic profiles of young TE loci in mouse and human cells. Insights gained using our method will be invaluable for dissecting the molecular determinants of TE regulation and their impact on the genome.

Endemic species are especially vulnerable to biodiversity loss caused by isolation or habitat specificity, small population size, and anthropogenic factors. Endemic species biodiversity analysis has a critically important global value for the development of conservation strategies. The rare onion Allium ledebourianum is a narrow-lined endemic species, with natural populations located in the extreme climatic conditions of the Kazakh Altai. A. ledebourianum populations are decreasing everywhere due to anthropogenic impact, and therefore, this species requires preservation and protection. Conservation of this rare species is associated with monitoring studies to investigate the genetic diversity of natural populations. Fundamental components of eukaryote genome include multiple classes of interspersed repeats. Various PCR-based DNA fingerprinting methods are used to detect chromosomal changes related to recombination processes of these interspersed elements. These methods are based on interspersed repeat sequences and are an effective approach for assessing the biological diversity of plants and their variability. We applied DNA profiling approaches based on conservative sequences of interspersed repeats to assess the genetic diversity of natural A. ledebourianum populations located in the territory of Kazakhstan Altai. The analysis of natural A. ledebourianum populations, carried out using the DNA profiling approach, allowed the effective differentiation of the populations and assessment of their genetic diversity. We used conservative sequences of tRNA primer binding sites (PBS) of the long-terminal repeat (LTR) retrotransposons as PCR primers. Amplification using the three most effective PBS primers generated 628 PCR amplicons, with an average of 209 amplicons. The average polymorphism level varied from 34% to 40% for all studied samples. Resolution analysis of the PBS primers showed all of them to have high or medium polymorphism levels, which varied from 0.763 to 0.965. Results of the molecular analysis of variance showed that the general biodiversity of A. ledebourianum populations is due to interpopulation (67%) and intrapopulation (33%) differences. The revealed genetic diversity was higher in the most distant population of A. ledebourianum LD64, located on the Sarymsakty ridge of Southern Altai. This is the first genetic diversity study of the endemic species A. ledebourianum using DNA profiling approaches. This work allowed us to collect new genetic data on the structure of A. ledebourianum populations in the Altai for subsequent development of preservation strategies to enhance the reproduction of this relict species. The results will be useful for the conservation and exploitation of this species, serving as the basis for further studies of its evolution and ecology.

Despite major progress in defining the functional roles of genes, a complete understanding of their influences is far from being realized, even in relatively simple organisms. A major milestone in this direction arose via the completion of the yeast Saccharomyces cerevisiae gene-knockout collection (YKOC), which has enabled high-throughput reverse genetics, phenotypic screenings and analyses of synthetic-genetic interactions1–3. Ensuing experimental work has also highlighted some inconsistencies and mistakes in the YKOC, or genome instability events that rebalance the effects of specific knockouts4–6, but a complete overview of these is lacking. The identification and analysis of genes that are required for maintaining genomic stability have traditionally relied on reporter assays and on the study of deletions of individual genes, but whole-genome-sequencing technologies now enable—in principle—the direct observation of genome instability globally and at scale. To exploit this opportunity, we sequenced the whole genomes of nearly all of the 4,732 strains comprising the homozygous diploid YKOC. Here, by extracting information on copy-number variation of tandem and interspersed repetitive DNA elements, we describe—for almost every single non-essential gene—the genomic alterations that are induced by its loss. Analysis of this dataset reveals genes that affect the maintenance of various genomic elements, highlights cross-talks between nuclear and mitochondrial genome stability, and shows how strains have genetically adapted to life in the absence of individual non-essential genes. Whole-genome sequencing of the strains of the Saccharomyces cerevisiae gene-knockout collection reveals the effects of the deletion of non-essential genes on genome stability.

Comprehensive and accurate annotation of the repeatome, including transposons, is critical for deepening our understanding of repeat origins, biogenesis, regulatory mechanisms, and roles. Here, we developed Generic Repeat Finder (GRF), a tool for genome-wide repeat detection based on fast, exhaustive numerical calculation algorithms integrated with optimized dynamic programming strategies. GRF sensitively identifies terminal inverted repeats (TIRs), terminal direct repeats (TDRs), and interspersed repeats that bear both inverted and direct repeats. GRF also detects DNA or RNA transposable elements characterized by these repeats in plant and animal genomes. For TIRs and TDRs, GRF identifies spacers in the middle and mismatches/insertions or deletions in terminal repeats, showing their alignment or base-pairing information. GRF helps improve the annotation for various DNA transposons and retrotransposons, such as miniature inverted-repeat transposable elements (MITEs), long terminal repeat (LTR) retrotransposons, and non-LTR retrotransposons, including long interspersed nuclear elements and short interspersed nuclear elements in plants. We used GRF to perform TIR/TDR, interspersed-repeat, and MITE detection in several species, including Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), and mouse (Mus musculus). As a generic bioinformatics tool in repeat finding implemented as a parallelized C++ program, GRF was faster and more sensitive than the existing inverted repeat/MITE detection tools based on numerical approaches (i.e. detectIR and detectMITE) in Arabidopsis and mouse. GRF is more sensitive than Inverted Repeat Finder in TIR detection, LTR_FINDER in short TDR detection (≤1,000 nt), and phRAIDER in interspersed repeat detection in Arabidopsis and rice. GRF is an open source available from Github.

Background Draft quality genomes for a multitude of organisms have become common due to the advancement of genome assemblers using long-read technologies with high error rates. Although current assemblies are substantially more contiguous than assemblies based on short reads, complete chromosomal assemblies are still challenging. Interspersed repeat families with multiple copy versions dominate the contig and scaffold ends of current long-read assemblies for complex genomes. These repeat families generally remain unresolved, as existing algorithmic solutions either do not scale to large copy numbers or can not handle the current high read error rates. Results We propose novel repeat resolution methods for large interspersed repeat families and assess their accuracy on simulated data sets with various distinct repeat structures and on drosophila melanogaster transposons. Additionally, we compare our methods to an existing long read repeat resolution tool and show the improved accuracy of our method. Conclusions Our results demonstrate the applicability of our methods for the improvement of the contiguity of genome assemblies. Electronic supplementary material The online version of this article (10.1186/s12859-019-2807-4) contains supplementary material, which is available to authorized users.

Atlantic salmon (Salmo salar) is a major species produced in world aquaculture and an important vertebrate model organism for studying the process of rediploidization following whole genome duplication events (Ss4R, 80 mya). The current Salmo salar transcriptome is largely generated from genome sequence based in silico predictions supported by ESTs and short-read sequencing data. However, recent progress in long-read sequencing technologies now allows for full-length transcript sequencing from single RNA-molecules. This study provides a de novo full-length mRNA transcriptome from liver, head-kidney and gill materials. A pipeline was developed based on Iso-seq sequencing of long-reads on the PacBio platform (HQ reads) followed by error-correction of the HQ reads by short-reads from the Illumina platform. The pipeline successfully processed more than 1.5 million long-reads and more than 900 million short-reads into error-corrected HQ reads. A surprisingly high percentage (32%) represented expressed interspersed repeats, while the remaining were processed into 71 461 full-length mRNAs from 23 071 loci. Each transcript was supported by several single-molecule long-read sequences and at least three short-reads, assuring a high sequence accuracy. On average, each gene was represented by three isoforms. Comparisons to the current Atlantic salmon transcripts in the RefSeq database showed that the long-read transcriptome validated 25% of all known transcripts, while the remaining full-length transcripts were novel isoforms, but few were transcripts from novel genes. A comparison to the current genome assembly indicates that the long-read transcriptome may aid in improving transcript annotation as well as provide long-read linkage information useful for improving the genome assembly. More than 80% of transcripts were assigned GO terms and thousands of transcripts were from genes or splice-variants expressed in an organ-specific manner demonstrating that hybrid error-corrected long-read transcriptomes may be applied to study genes and splice-variants expressed in certain organs or conditions (e.g., challenge materials). In conclusion, this is the single largest contribution of full-length mRNAs in Atlantic salmon. The results will be of great value to salmon genomics research, and the pipeline outlined may be applied to generate additional de novo transcriptomes in Atlantic Salmon or applied for similar projects in other species.

We present SquiggleNet, the first deep-learning model that can classify nanopore reads directly from their electrical signals. SquiggleNet operates faster than DNA passes through the pore, allowing real-time classification and read ejection. Using 1 s of sequencing data, the classifier achieves significantly higher accuracy than base calling followed by sequence alignment. Our approach is also faster and requires an order of magnitude less memory than alignment-based approaches. SquiggleNet distinguished human from bacterial DNA with over 90% accuracy, generalized to unseen bacterial species in a human respiratory meta genome sample, and accurately classified sequences containing human long interspersed repeat elements. Supplementary Information The online version contains supplementary material available at (10.1186/s13059-021-02511-y).

Analysis of the genetic diversity of natural populations of threatened and endangered species of plants is a main aspect of conservation strategy. The endangered species Allium altaicum is a relict plant of the Ice Age and natural populations are located in extreme climatic conditions of Kazakstan’s Altai Mountains. Mobile genetic elements and other interspersed repeats are basic components of a eukaryote genome, which can activate under stress conditions and indirectly promote the survival of an organism against environmental stresses. Detections of chromosomal changes related to recombination processes of mobile genetic elements are performed by various PCR methods. These methods are based on interspersed repeat sequences and are an effective tool for research of biological diversity of plants and their variability. In our research, we used conservative sequences of tRNA primer binding sites (PBS) when initializing the retrotransposon replication as PCR primers to research the genetic diversity of 12 natural populations of A. altaicum found in various ecogeographic conditions of the Kazakhstani Altai. High efficiency of the PBS amplification method used was observed already at the intrapopulation level. Unique amplicons representative of a certain population were found at the intrapopulation level. Analysis of molecular dispersion revealed that the biodiversity of populations of mountainous and lowland A. altaicum is due to intrapopulation differences for climatic zones of habitation. This is likely conditional upon predominance of vegetative reproduction over seed reproduction in some populations. In the case of vegetative reproduction, somatic recombination related to the activity of mobile genetic elements are preserved in subsequent generations. This leads to an increase of intrapopulation genetic diversity. Thus, high genetic diversity was observed in populations such as A. altaicum located in the territory of the Kalbinskii Altai, whereas the minimum diversity was observed in the populations of the Leninororsk ecogeographic group. Distinctions between these populations were also identified depending on the areas of their distribution. Low-land and mid-mountain living environments are characterized by a great variety of shapes and plasticity. This work allowed us to obtain new genetic data on the structure of A. altaicum populations on the territory of the Kazakhstan Altai for the subsequent development of preservation and reproduction strategies for this relict species.

Repetitive DNA sequences are the homologous DNA fragments with multiple copies and are the main components in higher plant genomes. Reliable chromosome-scaled genome assembly provided abundant genome information regarding genome structure research in sesame. In this chapter, we present the genome structure character and distribution of interspersed repeat sequences and the tandem repeats in the assembled genome of sesame var. Yuzhi 11 (version 3.0). For S. indicum, the content of interspersed repeat sequences is 46.11%. Of the five groups of interspersed repetitive sequences, i.e., SINEs, LINEs, long terminal repeat (LTR) element, DNA element, and unknown dispersed repeat elements, LTR group is the biggest known repeat type and occupies 11.41% of the assembled genome (335 Mb). Ty1-Copia is the major superfamily in LTR-RTs. In sesame genome, distribution of the tandem repeats, such as rDNA, telomere, and centromere repeats, reflects the character of species. Cytogenetic analyses and genome structure showed that the individual 26 chromosomes in S. indicum have telomeric repeats at the terminal position. The telomeres of sesame are conserved, and the repeat sequences are same with the common repeat unit of TTTAGGG. A total of 1235 telomeric repeat of (TTTAGGG)3 were found in the assembled genome. Meanwhile, about three thousand 45S and 5S rDNA repeats are detected in the assemble genome, which accord with the rDNA content range in other plants. Distribution characters and copy number of telomeres and rDNAs reflect the relative completeness of the assembled sesame genome. In addition, distribution of simple sequence repeats (SSRs) in sesame genome also has been precisely analyzed. A 153 bp centromeric repeat sequence is identified in sesame genome for the first time, which supplies more information for exploring the genome evolution of Sesamum.

Transposable elements (TEs), also known as mobile elements (MEs), are interspersed repeats that constitute a major fraction of the genomes of higher organisms. As one of their important functional impacts on gene function and genome evolution, TEs participate in regulating the expression of genes nearby and even far away at transcriptional and post-transcriptional levels. There are two principal ways by which TEs regulate expression of genes in the human genome. First, TEs provide cis-regulatory sequences in the genome. TEs’ intrinsic regulatory properties for their own expression make them potential factors for regulating the expression of the host genes. TE-derived cis-regulatory sites are found in promoter and enhancer elements, providing binding sites for a wide range of trans-acting factors. Second, TEs encode for regulatory RNAs. TEs sequences have been revealed to be present in a substantial fraction of miRNAs and long non-coding RNAs (lncRNAs), indicating their TE origin. Furthermore, TEs sequences were found to be critical for regulatory functions of these RNAs including binding to the target mRNA. TEs thus provide crucial regulatory roles by being part of cis-regulatory and regulatory RNA sequences. Moreover, both TE-derived cis-regulatory sequences and TE-derived regulatory RNAs, have been implicated to provide evolutionary novelty to gene regulation. These TE-derived regulatory mechanisms also tend to function in tissue-specific fashion. In this review, we aim to comprehensively cover the studies regarding these two aspects of TE-mediated gene regulation, mainly focusing on the mechanisms, contribution of different types of TEs, differential roles among tissue types, and lineage specificity, based on data mostly in humans.

The human CGGBP1 is implicated in a variety of cellular functions. It regulates genomic integrity, cell cycle, gene expression and cellular response to growth signals. Evidence suggests that these functions of CGGBP1 manifest through binding to GC-rich regions in the genome and regulation of interspersed repeats. Recent works show that CGGBP1 is needed for cytosine methylation homeostasis and genome-wide occupancy patterns of the epigenome regulator protein CTCF. It has remained unknown if cytosine methylation regulation and CTCF occupancy regulation by CGGBP1 are independent or interdependent processes. By sequencing immunoprecipitated methylated DNA, we have found that some transcription factor-binding sites resist stochastic changes in cytosine methylation. Of these, we have analyzed the CTCF-binding sites thoroughly and show that cytosine methylation regulation at CTCF-binding DNA sequence motifs by CGGBP1 is deterministic. These CTCF-binding sites are positioned at locations where the spread of cytosine methylation in cis depends on the levels of CGGBP1. Our findings suggest that CTCF occupancy and functions are determined by CGGBP1-regulated cytosine methylation patterns.

A great number of interspersed repetitive sequences in chromosomes make it difficult to identify chromosomal material via fluorescence in situ hybridization (FISH). The traditional approach to solve this problem is chromosome in situ suppression hybridization (CISS-hybridization). Unfortunately, it is impossible to be performed or fails with chromosomes of many eukaryote species. The aim of this study was to consider the image enhance procedure [1] and the in silico method of chromosome specific signal visualization (method VISSIS) [2] as alternatives to CISS-hybridization. The effectiveness of approaches for identification of specific signals was estimated by signal-to-background ratio (SNR). The computer methods were applied to images of human chromosomes, obtained with FISH of the whole chromosome painting DNA probes. Results showed that effectiveness of image processing methods depends on ratio of short and line interspersed elements (SINEs/LINEs) in DNA probes. The closer chromosomes in ratio of SINEs/LINEs, the higher specific signal intensities and signal-to-background ratios could be achieved. This suggests that computer methods can be efficient only with application of DNA probes derived from chromosomes characterized with similar ratio of SINE and LINE contents.

The common pheasant Phasianus colchicus, belonging to the order Galliformes and family Phasianidae, is the most widespread species. Despite a long history of captivity, the domestication of this bird is still at a preliminary stage. Recently, the demand for accelerating its transformation to poultry for meat and egg production has been increasing. In this study, we assembled high quality, chromosome scale genome of the common pheasant by using PacBio long reads, next-generation short reads, and Hi-C technology. The primary assembly has contig N50 size of 1.33 Mb and scaffold N50 size of 59.46 Mb, with a total size of 0.99 Gb, resolving most macrochromosomes into single scaffolds. A total of 23,058 genes and 10.71 Mb interspersed repeats were identified, constituting 30.31% and 10.71% of the common pheasant genome, respectively. Our phylogenetic analysis revealed that the common pheasant shared common ancestors with turkey about 24.7-34.5 million years ago (Ma). Rapidly evolved gene families, as well as branch-specific positively selected genes, indicate that calcium-related genes are potentially related to the adaptive and evolutionary change of the common pheasant. Interestingly, we found that the common pheasant has a unique major histocompatibility complex B locus (MHC-B) structure: three major inversions occurred in the sequence compared with chicken MHC-B. Furthermore, we detected signals of selection in five breeds of domestic common pheasant, several of which are production-oriented.

Long interspersed nuclear elements-1 (LINE-1) are mobile DNA elements that comprise the majority of interspersed repeats in the mammalian genome. During the last decade, these transposable sequences have been described as controlling elements involved in transcriptional regulation and genome plasticity. Recently, LINE-1 have been implicated in neurogenesis, but to date little is known about their nuclear organization in neurons. The olfactory epithelium is a site of continuous neurogenesis, and loci of olfactory receptor genes are enriched in LINE-1 copies. Olfactory neurons have a unique inverted nuclear architecture and constitutive heterochromatin forms a block in the center of the nuclei. Our DNA FISH images show that, even though LINE-1 copies are dispersed throughout the mice genome, they are clustered forming a cap around the central heterochromatin block and frequently occupy the same position as facultative heterochromatin in olfactory neurons nuclei. This specific LINE-1 organization could not be observed in other olfactory epithelium cell types. Analyses of H3K27me3 and H3K9me3 ChIP-seq data from olfactory epithelium revealed that LINE-1 copies located at OR gene loci show different enrichment for these heterochromatin marks. We also found that LINE-1 are transcribed in mouse olfactory epithelium. These results suggest that LINE-1 play a role in the olfactory neurons' nuclear architecture. SIGNIFICANCE STATEMENT: LINE-1 are mobile DNA elements and comprise almost 20% of mice and human genomes. These retrotransposons have been implicated in neurogenesis. We show for the first time that LINE-1 retrotransposons have a specific nuclear organization in olfactory neurons, forming aggregates concentric to the heterochromatin block and frequently occupying the same region as facultative heterochromatin. We found that LINE-1 at olfactory receptor gene loci are differently enriched for H3K9me3 and H3K27me3, but LINE-1 transcripts could be detected in the olfactory epithelium. We speculate that these retrotransposons play an active role in olfactory neurons' nuclear architecture.

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.

Diuraphis noxia Kurdjumov, 1913 (Hemiptera: Aphididae), commonly known as the Russian wheat aphid, is a devastating pest of wheat and barley. Although fourteen sources of resistance (Dn genes) have...

BACKGROUND Approximately 1% of normal tension glaucoma (NTG) cases are caused by TANK-binding kinase 1 (TBK1) gene duplications and triplications. However, the precise borders and orientation of these TBK1 gene copy number variations (CNVs) on chromosome 12 are unknown. METHODS We determined the exact borders of TBK1 CNVs and the orientation of duplicated or triplicated DNA segments in 5 NTG patients with different TBK1 mutations using whole-genome sequencing. RESULTS Tandemly duplicated chromosome segments spanning the TBK1 gene were detected in 4 NTG patients, each with unique borders. Four of 5 CNVs had borders located within interspersed repetitive DNA sequences (Alu and long interspersed nuclear element-L1 elements), suggesting that mismatched homologous recombinations likely generated these CNVs. A fifth NTG patient had a complex rearrangement including triplication of a chromosome segment spanning the TBK1 gene. CONCLUSIONS No specific mutation hotspots for TBK1 CNVs were detected, however, interspersed repetitive sequences (ie, Alu elements) were identified at the borders of TBK1 CNVs, which suggest that mismatch of these elements during meiosis may be the mechanism that generated TBK1 gene dosage mutations.

The factors guiding retrotransposon insertion site preference are not well understood. Different types of retrotransposons share common replication machinery and yet occupy distinct genomic domains. Autonomous long interspersed elements accumulate in gene-poor domains and their non-autonomous short interspersed elements accumulate in gene-rich domains. To determine genomic factors that contribute to this discrepancy we analysed the distribution of retrotransposons within the framework of chromosomal domains and regulatory elements. Using comparative genomics, we identified large-scale conserved patterns of retrotransposon accumulation across several mammalian genomes. Importantly, retrotransposons that were active after our sample-species diverged accumulated in orthologous regions. This suggested a similar evolutionary interaction between retrotransposon activity and conserved genome architecture across our species. In addition, we found that retrotransposons accumulated at regulatory element boundaries in open chromatin, where accumulation of particular retrotransposon types depended on insertion size and local regulatory element density. From our results, we propose a model where density and distribution of genes and regulatory elements canalise retrotransposon accumulation. Through conservation of synteny, gene regulation and nuclear organisation, mammalian genomes with dissimilar retrotransposons follow similar evolutionary trajectories.

Background DNA transposons helitrons are mobile genetic elements responsible for major movements of the genetic material within and across different genomes. This ability makes helitrons suitable candidate elements for the development of new approaches of multilocus genotyping of live-stock animals, along with the well-known microsatellite loci. Objective We aimed to estimate the informativeness of helitron and microsatellite markers in assessing the consolidation and the "gene pool" standards of two commercial dairy cattle breeds (Ayrshire breed and holsteinized Black-and-White cattle) and one local breed of Kalmyk cattle, and to reveal any inter-breed difference in the organization of genomic regions flanked by helitrons in the studied cattle breeds. Method We used the combination of two highly-polymorphic genomic elements - helitrons and trinu-cleotide microsatellites (AGC)6G and (GAG)6C, respectively - for genome scanning of the sampled groups of cattle. Also, we pyrosequenced the genomic regions flanked by the inverted repeats of 3'-end of Heligloria family of helitron fragments. Results Generally, the both combinations of markers generated polymorphic spectra, based on which certain interbreed differentiation could be observed. The analysis of the identified interspersed repeats suggests that in factory and local cattle the genomic regions flanked by helitron fragments are shaped differently and contain different superfamilies of transposable elements, especially retrotransposons. Conclusion Despite the well-known fact of retrotransposon-dependent microsatellite expansion, our data suggest that, in the cattle genome, the DNA transposons and microsatellites can also be found in close neighbourhood, and that helitrons and retrotransposons may form domains of increased variability - targets for factors of artificial selection.

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.

Plasmodium simium, a malaria parasite of non-human primates in the Atlantic forest region of Brazil was recently shown to cause zoonotic infection in humans in the region. Phylogenetic analyses based on the whole genome sequences of six P. simium isolates infecting humans and two isolates from brown howler monkeys revealed that P. simium is monophyletic within the broader diversity of South American Plasmodium vivax, consistent with the hypothesis that P. simium first infected non-human primates as a result of a host-switch from humans carrying P. vivax. We provide molecular evidence that the current zoonotic infections of people have likely resulted from multiple independent host switches, each seeded from a different monkey infection. Very low levels of genetic diversity within P. simium genomes and the absence of P. simium-P. vivax hybrids suggest that the P. simium population emerged recently and has subsequently experienced a period of independent evolution in Platyrrhini monkeys. We further find that Plasmodium Interspersed Repeat (PIR) genes, Plasmodium Helical Interspersed Subtelomeric (PHIST) genes and Tryptophan-Rich Antigens (TRAg) genes in P. siumium are genetically divergent from P. vivax and are enriched for non-synonymous single nucleotide polymorphisms, consistent with the rapid evolution of these genes. Analysis of genes involved in erythrocyte invasion revealed several notable differences between P. vivax and P. simium, including large deletions within the coding region of the Duffy Binding Protein 1 (DBP1) and Reticulocyte Binding Protein 2a (RBP2a) genes in P. simium. Genotyping of P. simium isolates from non-human primates (NHPs) and zoonotic human infections showed that a precise deletion of 38 amino acids in DBP1 is exclusively present in all human infecting isolates, whereas non-human primate infecting isolates were polymorphic for the deletion. We speculate that these deletions in the parasite-encoded key erythrocyte invasion ligands and the additional rapid genetic changes have facilitated zoonotic transfer to humans. Non-human primate malaria parasites can be considered a reservoir of potential infectious human parasites that must be considered in any attempt of malaria elimination. The genome of P. simium will thus form an important basis for future functional characterizations on the mechanisms underlying malaria zoonosis.

Background CGGBP1 is a repeat-binding protein with diverse functions in the regulation of gene expression, cytosine methylation, repeat silencing and genomic integrity. CGGBP1 has also been identified as a cooperator of histone-modifying enzymes and as a component of CTCF-containing complexes that regulate the enhancer–promoter looping. CGGBP1–CTCF cross talk in chromatin regulation has been hitherto unknown. Results Here, we report that the occupancy of CTCF at repeats depends on CGGBP1. Using ChIP-sequencing for CTCF, we describe its occupancy at repetitive DNA. Our results show that endogenous level of CGGBP1 ensures CTCF occupancy preferentially on repeats over canonical CTCF motifs. By combining CTCF ChIP-sequencing results with ChIP sequencing for three different kinds of histone modifications (H3K4me3, H3K9me3 and H3K27me3), we show that the CGGBP1-dependent repeat-rich CTCF-binding sites regulate histone marks in flanking regions. Conclusion CGGBP1 affects the pattern of CTCF occupancy. Our results posit CGGBP1 as a regulator of CTCF and its binding sites in interspersed repeats.

The complete mitochondrial genome of Morchella importuna, the famous edible and medicinal mushroom, was assembled as a 272,238 bp single circular dsDNA. As the largest mitogenome among fungi, it exhibits several distinct characteristics. The mitogenome of M. importuna encoded 14 core conserved mitochondrial protein-coding genes and 151 mitochondrial non-conserved open reading frames (ncORFs) were predicted, of which 61 were annotated as homing endonuclease genes, and 108 were confirmed to be expressed during the vegetative growth stages of M. importuna. In addition, 34 introns were identified in seven core genes (cob, cox1, cox2, cox3, nad1, nad4 and nad5) and two rRNA genes (rrnS and rrnL) with a length from 383 bp to 7453 bp, and eight large introns with a length range of 2340 bp to 7453 bp contained multiple intronic mtORFs. Moreover, 34 group I (IA, IB, IC1, IC2, ID and derived group I introns) and four group II intron domains were identified for the 34 introns, including five hybrid ones. Furthermore, the M. importuna mitogenome showed the presence of about 18.7% mitogenomic interspersed repeats. These and the aforementioned ncORFs and introns, contributed to the enlarged size of the mitogenome.

Mitochondria is a powerhouse of all eukaryotic cells that have its own circular DNA (mtDNA) encoding various RNAs and proteins. Somatic perturbations of mtDNA are accumulating with age thus it is of great importance to uncover the main sources of mtDNA instability. Recent analyses demonstrated that somatic mtDNA deletions depend on imperfect repeats of various nature between distant mtDNA segments. However, till now there are no comprehensive databases annotating all types of imperfect repeats in numerous species with sequenced complete mitochondrial genome as well as there are no algorithms capable to call all types of imperfect repeats in circular mtDNA. We implemented naive algorithm of pattern recognition by analogy to standard dot-plot construction procedures allowing us to find both perfect and imperfect repeats of four main types: direct, inverted, mirror and complementary. Our algorithm is adapted to specific characteristics of mtDNA such as circularity and an excess of short repeats - it calls imperfect repeats starting from the length of 10 b.p. We constructed interactive web available database ImtRDB depositing perfect and imperfect repeats positions in mtDNAs of more than 3500 Vertebrate species. Additional tools, such as visualization of repeats within a genome, comparison of repeat densities among different genomes and a possibility to download all results make this database useful for many biologists. Our first analyses of the database demonstrated that mtDNA imperfect repeats (i) are usually short; (ii) associated with unfolded DNA structures; (iii) four types of repeats positively correlate with each other forming two equivalent pairs: direct and mirror versus inverted and complementary, with identical nucleotide content and similar distribution between species; (iv) abundance of repeats is negatively associated with GC content; (v) dinucleotides GC versus CG are overrepresented on light chain of mtDNA covered by repeats. ImtRDB is available at http://bioinfodbs.kantiana.ru/ImtRDB/ . It is accompanied by the software calling all types of interspersed repeats with different level of degeneracy in circular DNA. This database and software can become a very useful tool in various areas of mitochondrial and chloroplast DNA research.

RNA editing is a post-transcription maturation process that diversifies genomically encoded information and can lead to transcriptome diversity. Thanks to next-generation sequencing technologies, a large number of editing sites have been identified in different species. Although this mechanism is well described in mammals, only a few studies have been performed in chicken. Here, candidate or potential RNA editing sites were identified in eight different tissues of chicken (brain, spleen, colon, lung, kidney, heart, testes and liver). We identified 68 A-to-G editing sites in 46 genes. Only two of these were previously reported in chicken. We found no C-to-T sites, attesting to the lack of this type of editing mechanism in chicken. Similar to mammals, the editing sites were enriched in non-coding regions, rarely resulted in a change in amino acids, showed a critical role in the nervous system and had a low guanosine level upstream of the editing site and some enrichment downstream from the site. Moreover, in contrast to mammals, editing sites were weakly enriched in interspersed repeats and the number and editing ratio of non-synonymous sites were higher than for those of synonymous sites. Interestingly, we found several tissue-specific edited genes, including GABRA3, SORL1 and HTR1D in brain and RYR2 and FHOD3 in heart, that were associated with functional processes relevant to the corresponding tissue. This finding highlights the importance of RNA editing in several chicken tissues, especially the brain, and establishes a foundation for further exploration of this process.

F-box proteins are a core component of Skp1-Cul1-F-box (SCF) ubiquitin/ligase complexes and are involved in a lot of cellular processes in yeasts. However, the current knowledge of the molecular evolution of the F-box gene family in yeasts remains unclear. In this study, 136 F-box genes were identified in 10 yeast species of the Saccharomycetaceae. In addition to the F-box domain, the other six domains were identified in these F-box proteins. The evolutionary history of F-box gene numbers in 10 Saccharomycetaceae yeasts was reconstructed. Whole-genome duplication, interspersed repeats, and gene loss events were inferred. These events contributed to F-box gene number variation in the 10 yeast species. Eighty-seven and 33 positively selected sites were detected in program Selecton and Datamonkey web-server, respectively. Three of them were considered the significant positively selected sites, and 23 of them had changed radically in amino acid properties by using TreeSAAP. We investigated F-box gene number variation and underlying mechanisms, and selection patterns, all of which were beneficial to deeply understand genome evolution and figure out the function of the F-box proteins.

Transposable elements (TEs) are interspersed repeat sequences that make up much of the human genome. Their expression has been implicated in development and disease. However, TE-derived RNA-seq reads are difficult to quantify. Past approaches have excluded these reads or aggregated RNA expression to subfamilies shared by similar TE copies, sacrificing quantitative accuracy or the genomic context necessary to understand the basis of TE transcription. As a result, the effects of TEs on gene expression and associated phenotypes are not well understood. Here, we present Software for Quantifying Interspersed Repeat Expression (SQuIRE), the first RNA-seq analysis pipeline that provides a quantitative and locus-specific picture of TE expression (https://github.com/wyang17/SQuIRE). SQuIRE is an accurate and user-friendly tool that can be used for a variety of species. We applied SQuIRE to RNA-seq from normal mouse tissues and a Drosophila model of amyotrophic lateral sclerosis. In both model organisms, we recapitulated previously reported TE subfamily expression levels and revealed locus-specific TE expression. We also identified differences in TE transcription patterns relating to transcript type, gene expression and RNA splicing that would be lost with other approaches using subfamily-level analyses. Altogether, our findings illustrate the importance of studying TE transcription with locus-level resolution.

BACKGROUND AND AIMS: The genus Solanum includes important vegetable crops and their wild relatives. Introgression of their useful traits into elite cultivars requires effective recombination between hom(e)ologues, which is partially determined by genome sequence differentiation. In this study we compared the repetitive genome fractions of wild and cultivated species of the potato and tomato clades in a phylogenetic context. METHODS: Genome skimming followed by a clustering approach was used as implemented in the RepeatExplorer pipeline. Repeat classes were annotated and the sequences of their main domains were compared. KEY RESULTS: Repeat abundance and genome size were correlated and the larger genomes of species in the tomato clade were found to contain a higher proportion of unclassified elements. Families and lineages of repetitive elements were largely conserved between the clades, but their relative proportions differed. The most abundant repeats were Ty3/Gypsy elements. Striking differences in abundance were found in the highly dynamic Ty3/Gypsy Chromoviruses and Ty1/Copia Tork elements. Within the potato clade, early branching Solanum cardiophyllum showed a divergent repeat profile. There were also contrasts between cultivated and wild potatoes, mostly due to satellite amplification in the cultivated species. Interspersed repeat profiles were very similar among potatoes. The repeat profile of Solanum etuberosum was more similar to that of the potato clade. CONCLUSIONS: The repeat profiles in Solanum seem to be very similar despite genome differentiation at the level of collinearity. Removal of transposable elements by unequal recombination may have been responsible for structural rearrangements across the tomato clade. Sequence variability in the tomato clade is congruent with clade-specific amplification of repeats after its divergence from S. etuberosum and potatoes. The low differentiation among potato and its wild relatives at the level of interspersed repeats may explain the difficulty in discriminating their genomes by genomic in situ hybridization techniques.

Genes that are commonly deregulated in cancer are clinically attractive as candidate pan-diagnostic markers and therapeutic targets. To globally identify such targets, we compared Cap Analysis of Gene Expression profiles from 225 different cancer cell lines and 339 corresponding primary cell samples to identify transcripts that are deregulated recurrently in a broad range of cancer types. Comparing RNA-seq data from 4,055 tumors and 563 normal tissues profiled in the The Cancer Genome Atlas and FANTOM5 datasets, we identified a core transcript set with theranostic potential. Our analyses also revealed enhancer RNAs, which are upregulated in cancer, defining promoters that overlap with repetitive elements (especially SINE/Alu and LTR/ERV1 elements) that are often upregulated in cancer. Lastly, we documented for the first time upregulation of multiple copies of the REP522 interspersed repeat in cancer. Overall, our genome-wide expression profiling approach identified a comprehensive set of candidate biomarkers with pan-cancer potential, and extended the perspective and pathogenic significance of repetitive elements that are frequently activated during cancer progression. Cancer Res; 76(2); 216–26. ©2015 AACR.

Alport syndrome-diffuse leiomyomatosis (AS-DL, OMIM: 308940) is a rare variant of the X-linked Alport syndrome that shows overgrowth of visceral smooth muscles in the gastrointestinal, respiratory and female reproductive tracts in addition to renal symptoms. AS-DL results from deletions that encompass the 5′ ends of the COL4A5 and COL4A6 genes, but deletion breakpoints between COL4A5 and COL4A6 have been determined in only four cases. Here, we characterize deletion breakpoints in five AS-DL patients and show a contiguous COL4A6/COL4A5 deletion in each case. We also demonstrate that eight out of nine deletion alleles involved sequences homologous between COL4A5 and COL4A6. Most breakpoints took place in recognizable transposed elements, including long and short interspersed repeats, DNA transposons and long-terminal repeat retrotransposons. Because deletions involved the bidirectional promoter region in each case, we suggest that the occurrence of leiomyomatosis in AS-DL requires inactivation of both genes. Altogether, our study highlights the importance of homologous recombination involving multiple transposed elements for the development of this continuous gene syndrome and other atypical loss-of-function phenotypes.

Summary Circular RNAs (circRNAs) are stable non-coding RNAs with a closed circular structure. One of the best studied circRNAs is ciRS-7 (CDR1as), which acts as a regulator of the microRNA miR-7; however, its biosynthetic pathway has remained an enigma. Here we delineate the biosynthetic pathway of ciRS-7. The back-splicing events that form circRNAs are often facilitated by flanking inverted repeats of the primate-specific Alu elements. The ciRS-7 gene lacks these elements, but, instead, we identified a set of flanking inverted elements belonging to the mammalian-wide interspersed repeat (MIR) family. Splicing reporter assays in HEK293 cells demonstrated that these inverted MIRs are required to generate ciRS-7 through back-splicing, and CRISPR/Cas9-mediated deletions confirmed the requirement of the endogenous MIR elements in SH-SY5Y cells. Using bioinformatic searches, we identified several other MIR-dependent circRNAs and confirmed them experimentally. We propose that MIR-mediated RNA circularization is used to generate a subset of mammalian circRNAs., Graphical Abstract, Highlights • The circular RNA, ciRS-7 (CDR1as), functions as a regulator of miR-7 • ciRS-7 is generated by back-splicing, not via intra-lariat splicing • Back-splicing of ciRS-7 is promoted by the flanking inverted MIR elements • The biosynthesis of a subset of mammalian circRNAs could be mediated by MIRs, Biological Sciences; Molecular Biology

A genomic differential display method was developed that analyzes many restriction fragment length polymorphisms simultaneously. Interspersed repeat sequences were used to reduce DNA sample complexity and to target genomic subsets of interest. This work focused on trinucleotide repeats because of their importance in human inherited diseases. Immobilized repeat-containing oligonucleotides were used to capture genomic DNA fragments containing sequences complementary to the oligonucleotide. Captured fragments were amplified by PCR and fluorescently labeled using primers complementary to the repeat sequence and/or to the known sequences ligated to the ends of the restriction fragments. The labeled PCR fragments were displayed by size on a high-resolution automated fluorescent DNA sequencing instrument. Although there was a conservation in the overall pattern of displayed genome subsets, many clear and reproducible differences were detected when genomes from different individuals were compared. Fewer differences were detected within, than between, monozygotic twin pair genomes. In control experiments, the method distinguished between Huntington disease alleles with normal and expanded CAG repeat lengths.

SUMMARYCircular RNAs (circRNAs) are stable noncoding RNAs with a closed circular structure. One of the first and best studied circRNAs is ciRS-7 (CDR1as) that acts as a regulator of the microRNA miR-7, however, the biosynthesis pathway has remained an enigma. Here we delineate the biosynthesis pathway of ciRS-7. The back-splicing events that form circRNAs are often facilitated by flanking inverted repeats of the primate-specific Alu elements. ciRS-7 gene lacks these elements but, instead, we identified a set of flanking inverted elements belonging to the mammalian-wide interspersed repeat (MIR) family. Splicing reporter assays in HEK293 cells demonstrated that these inverted MIRs are required to generate ciRS-7 through a back-splicing and CRISPR/Cas9-mediated deletions confirmed the requirement of the endogenous MIR elements in SH-SY5Y cells. Using bioinformatics searches, we identified several other MIR-dependent circRNAs that we confirmed experimentally. We propose that MIR-mediated RNA circularization constitutes a new widespread biosynthesis principle for mammalian circRNAs.

Assembling a plant genome is challenging due to the abundance of repetitive sequences, yet no standard is available to evaluate the assembly of repeat space. LTR retrotransposons (LTR-RTs) are the predominant interspersed repeat that is poorly assembled in draft genomes. Here, we propose a reference-free genome metric called LTR Assembly Index (LAI) that evaluates assembly continuity using LTR-RTs. After correcting for LTR-RT amplification dynamics, we show that LAI is independent of genome size, genomic LTR-RT content, and gene space evaluation metrics (i.e., BUSCO and CEGMA). By comparing genomic sequences produced by various sequencing techniques, we reveal the significant gain of assembly continuity by using long-read-based techniques over short-read-based methods. Moreover, LAI can facilitate iterative assembly improvement with assembler selection and identify low-quality genomic regions. To apply LAI, intact LTR-RTs and total LTR-RTs should contribute at least 0.1% and 5% to the genome size, respectively. The LAI program is freely available on GitHub: https://github.com/oushujun/LTR_retriever.

Determining cell lineage and function is critical to understanding human physiology and pathology. Although advances in lineage tracing methods provide new insight into cell fate, defining cellular diversity at the mammalian level remains a challenge. Here, we develop a genome editing strategy using a cytidine deaminase fused with nickase Cas9 (nCas9) to specifically target endogenous interspersed repeat regions in mammalian cells. The resulting mutation patterns serve as a genetic barcode, which is induced by targeted mutagenesis with single-guide RNA (sgRNA), leveraging substitution events, and subsequent read out by a single primer pair. By analyzing interspersed mutation signatures, we show the accurate reconstruction of cell lineage using both bulk cell and single-cell data. We envision that our genetic barcode system will enable fine-resolution mapping of organismal development in healthy and diseased mammalian states., Lineage tracing has provided new insights into cell fate but defining cellular diversity remains a challenge. Here the authors target endogenous repeat regions in mammalian cells with cytidine deaminase fused to nCas9 to create genetic barcodes for fine-resolution mapping.

RNA editing is a post-transcription maturation process that diversifies genomically encoded information and can lead to diversity and complexity of transcriptome, especially in the brain. Thanks to next-generation sequencing technologies, a large number of editing sites have been identified in different species, especially in human, mouse and rat. While this mechanism is well described in mammals, only a few studies have been performed in the chicken. Here, we developed a rigorous computational strategy to identify RNA editing sites in eight different tissues of the chicken (brain, spleen, colon, lung, kidney, heart, testes and liver), based on RNA sequencing data alone. We identified 68 A-to-G editing sites in 46 genes. Only two of these were previously reported in chicken. We found no C-to-U sites, attesting the lack of this type of editing mechanism in the chicken. Similar to mammals, the editing sites were enriched in non-coding regions, rarely resulted in change of amino acids, showed a critical role in nervous system and had a low guanosine level upstream of the editing site and some enrichment downstream from the site. Moreover, in contrast to mammals, editing sites were weakly enriched in interspersed repeats and the frequency and editing ratio of non-synonymous sites were higher than those of synonymous sites.Interestingly, we found several tissue-specific edited genes including GABRA3, SORL1 and HTR1D in brain and RYR2 and FHOD3 in heart that were associated with functional processes relevant to the corresponding tissue. This finding highlighted the importance of the RNA editing in several chicken tissues, especially the brain. This study extends our understanding of RNA editing in chicken tissues and establish a foundation for further exploration of this process.

Transposable elements are interspersed repeat sequences that make up much of the human genome. Conventional approaches to RNA-seq analysis often exclude these sequences, fail to optimally adjudicate read alignments, or align reads to interspersed repeat consensus sequences without considering these transcripts in their genomic contexts. As a result, repetitive sequence contributions to transcriptomes are not well understood. Here, we present Software for Quantifying Interspersed Repeat Expression (SQuIRE), an RNA-seq analysis pipeline that integrates repeat and genome annotation (RepeatMasker), read alignment (STAR), gene expression (StringTie) and differential expression (DESeq2). SQuIRE uniquely provides a locus-specific picture of interspersed repeat-encoded RNA expression. SQuIRE can be downloaded at (github.com/wyang17/SQuIRE).

Transpositions transfer DNA segments between different loci within a genome; in particular, when a transposition is found in a sample but not in a reference genome, it is called a non-reference transposition. They are important structural variations that have clinical impact. Transpositions can be called by analyzing second generation high-throughput sequencing datasets. Current methods follow either a database-based or a database-free approach. Database-based methods require a database of transposable elements. Some of them have good specificity; however this approach cannot detect novel transpositions, and it requires a good database of transposable elements, which is not yet available for many species. Database-free methods perform de novo calling of transpositions, but their accuracy is low. We observe that this is due to the misalignment of the reads; since reads are short and the human genome has many repeats, false alignments create false positive predictions while missing alignments reduce the true positive rate. This paper proposes new techniques to improve database-free non-reference transposition calling: first, we propose a realignment strategy called one-end remapping that corrects the alignments of reads in interspersed repeats; second, we propose a SNV-aware filter that removes some incorrectly aligned reads. By combining these two techniques and other techniques like clustering and positive-to-negative ratio filter, our proposed transposition caller TranSurVeyor shows at least 3.1-fold improvement in terms of F1-score over existing database-free methods. More importantly, even though TranSurVeyor does not use databases of prior information, its performance is at least as good as existing database-based methods such as MELT, Mobster and Retroseq. We also illustrate that TranSurVeyor can discover transpositions that are not known in the current database.

Assessing correctness of an assembled chromosome architecture is a central challenge. We create a geometric analysis method (called GenomeLandscaper) to conduct landscape analysis of genome-fingerprints maps (GFM), trace large-scale repetitive regions, and assess their impacts on the global architectures of assembled chromosomes. We develop an alignment-free method for phylogenetics analysis. The human Y chromosomes (GRCh.chrY, HuRef.chrY and YH.chrY) are analysed as a proof-of-concept study. We construct a galaxy of genome-fingerprints maps (GGFM) for them, and a landscape compatibility among relatives is observed. But a long sharp straight line on the GGFM breaks such a landscape compatibility, distinguishing GRCh38p1.chrY (and throughout GRCh38p7.chrY) from GRCh37p13.chrY, HuRef.chrY and YH.chrY. We delete a 1.30-Mbp target segment to rescue the landscape compatibility, matching the antecedent GRCh37p13.chrY. We re-locate it into the modelled centromeric and pericentromeric region of GRCh38p10.chrY, matching a gap placeholder of GRCh37p13.chrY. We decompose it into sub-constituents (such as BACs, interspersed repeats, and tandem repeats) and trace their homologues by phylogenetics analysis. We elucidate that most examined tandem repeats are of reasonable quality, but the BAC-sized repeats, 173U1020C (176.46 Kbp) and 5U41068C (205.34 Kbp), are likely over-repeated. These results offer unique insights into the centromeric and pericentromeric regions of the human Y chromosomes.

Repetitive DNA--sequence motifs repeated hundreds or thousands of times in the genome--makes up the major proportion of all the nuclear DNA in most eukaryotic genomes. However, the significance of repetitive DNA in the genome is not completely understood, and it has been considered to have both structural and functional roles, or perhaps even no essential role. High-throughput DNA sequencing reveals huge numbers of repetitive sequences. Most bioinformatic studies focus on low-copy DNA including genes, and hence, the analyses collapse repeats in assemblies presenting only one or a few copies, often masking out and ignoring them in both DNA and RNA read data. Chromosomal studies are proving vital to examine the distribution and evolution of sequences because of the challenges of analysis of sequence data. Many questions are open about the origin, evolutionary mode and functions that repetitive sequences might have in the genome. Some, the satellite DNAs, are present in long arrays of similar motifs at a small number of sites, while others, particularly the transposable elements (DNA transposons and retrotranposons), are dispersed over regions of the genome; in both cases, sequence motifs may be located at relatively specific chromosome domains such as centromeres or subtelomeric regions. Here, we overview a range of works involving detailed characterization of the nature of all types of repetitive sequences, in particular their organization, abundance, chromosome localization, variation in sequence within and between chromosomes, and, importantly, the investigation of their transcription or expression activity. Comparison of the nature and locations of sequences between more, and less, related species is providing extensive information about their evolution and amplification. Some repetitive sequences are extremely well conserved between species, while others are among the most variable, defining differences between even closely relative species. These data suggest contrasting modes of evolution of repetitive DNA of different types, including selfish sequences that propagate themselves and may even be transferred horizontally between species rather than by descent, through to sequences that have a tendency to amplification because of their sequence motifs, to those that have structural significance because of their bulk rather than precise sequence. Functional consequences of repeats include generation of variability by movement and insertion in the genome (giving useful genetic markers), the definition of centromeres, expression under stress conditions and regulation of gene expression via RNA moieties. Molecular cytogenetics and bioinformatic studies in a comparative context are now enabling understanding of the nature and behaviour of this major genomic component.

The main classes of repetitive DNA sequences, including coding (rRNA genes) and noncoding (tandem and interspersed repeats) sequences are reviewed. Emphasis is placed on their special role in the formation of the structural and functional organization of the genomes of higher plants and in the support of their higher genetic variation, compared to animal genomes, at the levels of individual sequences and of the whole genome.