Survey of extrachromosomal circular DNA derived from plant satellite repeats

Centromere diversity: How different repeat-based holocentromeres may have evolved

In addition to monocentric eukaryotes, which have a single localized centromere on each chromosome, there are holocentric species, with extended repeat-based or repeat-less centromeres distributed over the entire chromosome length. At least two types of repeat-based holocentromeres exist, one composed of many small repeat-based centromere units (small unit-type), and another one characterized by a few large centromere units (large unit-type). We hypothesize that the transposable element-mediated dispersal of hundreds of short satellite arrays formed the small centromere unit-type holocentromere in Rhynchospora pubera. The large centromere unit-type of the plant Chionographis japonica is likely a product of simultaneous DNA double-strand breaks (DSBs), which initiated the de novo formation of repeat-based holocentromeres via insertion of satellite DNA, derived from extra-chromosomal circular DNAs (eccDNAs). The number of initial DSBs along the chromosomes must be higher than the number of centromere units since only a portion of the breaks will have incorporated eccDNA at an appropriate position to serve as future centromere unit sites. Subsequently, preferential incorporation of the centromeric histone H3 variant at these positions is assumed. The identification of repeat-based holocentromeres across lineages will unveil the centromere plasticity and elucidate the mechanisms underlying the diverse formation of holocentromeres.

Extrachromosomal circular DNA (eccDNA): from carcinogenesis to drug resistance

Extrachromosomal circular DNA (eccDNA) is a circular form of DNA that exists outside of the chromosome. Although it has only been a few decades since its discovery, in recent years, it has been found to have a close relationship with cancer, which has attracted widespread attention from researchers. Thus far, under the persistent research of researchers from all over the world, eccDNA has been found to play an important role in a variety of tumors, including breast cancer, lung cancer, ovarian cancer, etc. Herein, we review the sources of eccDNA, classifications, and the mechanisms responsible for their biogenesis. In addition, we introduce the relationship between eccDNA and various cancers and the role of eccDNA in the generation and evolution of cancer. Finally, we summarize the research significance and importance of eccDNA in cancer, and highlight new prospects for the application of eccDNA in the future detection and treatment of cancer.

Repeat turnover meets stable chromosomes: repetitive DNA sequences mark speciation and gene pool boundaries in sugar beet and wild beets

Sugar beet and its wild relatives share a base chromosome number of nine and similar chromosome morphologies. Yet, interspecific breeding is impeded by chromosome and sequence divergence that is still not fully understood. Since repetitive DNAs are among the fastest evolving parts of the genome, we investigated, if repeatome innovations and losses are linked to chromosomal differentiation and speciation. We traced genome and chromosome-wide evolution across 13 beet species comprising all sections of the genera Beta and Patellifolia. For this, we combined short and long read sequencing, flow cytometry, and cytogenetics to build a comprehensive framework that spans the complete scale from DNA to chromosome to genome. Genome sizes and repeat profiles reflect the separation into three gene pools with contrasting evolutionary patterns. Among all repeats, satellite DNAs harbor most genomic variability, leading to fundamentally different centromere architectures, ranging from chromosomal uniformity in Beta and Patellifolia to the formation of patchwork chromosomes in Corollinae/Nanae. We show that repetitive DNAs are causal for the genome expansions and contractions across the beet genera, providing insights into the genomic underpinnings of beet speciation. Satellite DNAs in particular vary considerably between beet genomes, leading to the evolution of distinct chromosomal setups in the three gene pools, likely contributing to the barriers in beet breeding. Thus, with their isokaryotypic chromosome sets, beet genomes present an ideal system for studying the link between repeats, genomic variability, and chromosomal differentiation and provide a theoretical fundament for understanding barriers in any crop breeding effort.

The evolution and formation of centromeric repeats analysis in Vitis vinifera

MAIN CONCLUSION

Six grape centromere-specific markers for cytogenetics were mined by combining genetic and immunological assays, and the possible evolution mechanism of centromeric repeats was analyzed. Centromeric histone proteins are functionally conserved; however, centromeric repetitive DNA sequences may represent considerable diversity in related species. Therefore, studying the characteristics and structure of grape centromere repeat sequences contributes to a deeper understanding of the evolutionary process of grape plants, including their origin and mechanisms of polyploidization. Plant centromeric regions are mainly composed of repetitive sequences, including SatDNA and transposable elements (TE). In this research, the characterization of centromere sequences in the whole genome of grapevine (Vitis vinifera L.) has been conducted. Five centromeric tandem repeat sequences (Vv1, Vv2, Vv5, Vv6, and Vv8) and one long terminal repeat (LTR) sequence Vv24 were isolated. These sequences had different centromeric distributions, which indicates that grape centromeric sequences may undergo rapid evolution. The existence of extrachromosomal circular DNA (eccDNA) and gene expression in CenH3 subdomain region may provide various potential mechanisms for the generation of new centromeric regions.

Dynamics of extrachromosomal circular DNA in rice

The genome's dynamic nature, exemplified by elements like extrachromosomal circular DNA (eccDNA), is crucial for biodiversity and adaptation. Yet, the role of eccDNA in plants, particularly rice, remains underexplored. Here, we identify 25,598 eccDNAs, unveiling the widespread presence of eccDNA across six rice tissues and revealing its formation as a universal and random process. Interestingly, we discover that direct repeats play a pivotal role in eccDNA formation, pointing to a unique origin mechanism. Despite eccDNA's prevalence in coding sequences, its impact on gene expression is minimal, implying its roles beyond gene regulation. We also observe the association between eccDNA's formation and minor chromosomal deletions, providing insights of its possible function in regulating genome stability. Further, we discover eccDNA specifically accumulated in rice leaves, which may be associated with DNA damage caused by environmental stressors like intense light. In summary, our research advances understanding of eccDNA's role in the genomic architecture and offers valuable insights for rice cultivation and breeding.

Holocentromeres can consist of merely a few megabase-sized satellite arrays

The centromere is the chromosome region where microtubules attach during cell division. In contrast to monocentric chromosomes with one centromere, holocentric species usually distribute hundreds of centromere units along the entire chromatid. We assembled the chromosome-scale reference genome and analyzed the holocentromere and (epi)genome organization of the lilioid Chionographis japonica. Remarkably, each of its holocentric chromatids consists of only 7 to 11 evenly spaced megabase-sized centromere-specific histone H3-positive units. These units contain satellite arrays of 23 and 28 bp-long monomers capable of forming palindromic structures. Like monocentric species, C. japonica forms clustered centromeres in chromocenters at interphase. In addition, the large-scale eu- and heterochromatin arrangement differs between C. japonica and other known holocentric species. Finally, using polymer simulations, we model the formation of prometaphase line-like holocentromeres from interphase centromere clusters. Our findings broaden the knowledge about centromere diversity, showing that holocentricity is not restricted to species with numerous and small centromere units.

The Chlamydomonas Genome Project, version 6: Reference assemblies for mating-type plus and minus strains reveal extensive structural mutation in the laboratory

Five versions of the Chlamydomonas reinhardtii reference genome have been produced over the last two decades. Here we present version 6, bringing significant advances in assembly quality and structural annotations. PacBio-based chromosome-level assemblies for two laboratory strains, CC-503 and CC-4532, provide resources for the plus and minus mating-type alleles. We corrected major misassemblies in previous versions and validated our assemblies via linkage analyses. Contiguity increased over ten-fold and >80% of filled gaps are within genes. We used Iso-Seq and deep RNA-seq datasets to improve structural annotations, and updated gene symbols and textual annotation of functionally characterized genes via extensive manual curation. We discovered that the cell wall-less classical reference strain CC-503 exhibits genomic instability potentially caused by deletion of the helicase RECQ3, with major structural mutations identified that affect >100 genes. We therefore present the CC-4532 assembly as the primary reference, although this strain also carries unique structural mutations and is experiencing rapid proliferation of a Gypsy retrotransposon. We expect all laboratory strains to harbor gene-disrupting mutations, which should be considered when interpreting and comparing experimental results. Collectively, the resources presented here herald a new era of Chlamydomonas genomics and will provide the foundation for continued research in this important reference organism.

Sex chromosome differentiation via changes in the Y chromosome repeat landscape in African annual killifishes Nothobranchius furzeri and N. kadleci

Homomorphic sex chromosomes and their turnover are common in teleosts. We investigated the evolution of nascent sex chromosomes in several populations of two sister species of African annual killifishes, Nothobranchius furzeri and N. kadleci, focusing on their under-studied repetitive landscape. We combined bioinformatic analyses of the repeatome with molecular cytogenetic techniques, including comparative genomic hybridization, fluorescence in situ hybridization with satellite sequences, ribosomal RNA genes (rDNA) and bacterial artificial chromosomes (BACs), and immunostaining of SYCP3 and MLH1 proteins to mark lateral elements of synaptonemal complexes and recombination sites, respectively. Both species share the same heteromorphic XY sex chromosome system, which thus evolved prior to their divergence. This was corroborated by sequence analysis of a putative master sex determining (MSD) gene gdf6Y in both species. Based on their divergence, differentiation of the XY sex chromosome pair started approximately 2 million years ago. In all populations, the gdf6Y gene mapped within a region rich in satellite DNA on the Y chromosome long arms. Despite their heteromorphism, X and Y chromosomes mostly pair regularly in meiosis, implying synaptic adjustment. In N. kadleci, Y-linked paracentric inversions like those previously reported in N. furzeri were detected. An inversion involving the MSD gene may suppress occasional recombination in the region, which we otherwise evidenced in the N. furzeri population MZCS-121 of the Limpopo clade lacking this inversion. Y chromosome centromeric repeats were reduced compared with the X chromosome and autosomes, which points to a role of relaxed meiotic drive in shaping the Y chromosome repeat landscape. We speculate that the recombination rate between sex chromosomes was reduced due to heterochiasmy. The observed differences between the repeat accumulations on the X and Y chromosomes probably result from high repeat turnover and may not relate closely to the divergence inferred from earlier SNP analyses.

The wtf meiotic driver gene family has unexpectedly persisted for over 100 million years

Meiotic drivers are selfish elements that bias their own transmission into more than half of the viable progeny produced by a driver+/driver- heterozygote. Meiotic drivers are thought to exist for relatively short evolutionary timespans because a driver gene or gene family is often found in a single species or in a group of very closely related species. Additionally, drivers are generally considered doomed to extinction when they spread to fixation or when suppressors arise. In this study, we examine the evolutionary history of the wtf meiotic drivers first discovered in the fission yeast Schizosaccharomyces pombe. We identify homologous genes in three other fission yeast species, S. octosporus, S. osmophilus, and S. cryophilus, which are estimated to have diverged over 100 million years ago from the S. pombe lineage. Synteny evidence supports that wtf genes were present in the common ancestor of these four species. Moreover, the ancestral genes were likely drivers as wtf genes in S. octosporus cause meiotic drive. Our findings indicate that meiotic drive systems can be maintained for long evolutionary timespans.

Extrachromosomal circular DNA: A neglected nucleic acid molecule in plants

Throughout the years, most plant genomic studies were focused on nuclear chromosomes. Extrachromosomal circular DNA (eccDNA) has largely been neglected for decades since its discovery in 1965. While initial research showed that eccDNAs can originate from highly repetitive sequences, recent findings show that many regions of the genome can contribute to the eccDNA pool. Currently, the biological functions of eccDNAs, if any, are a mystery but recent studies have indicated that they can be regulated by different genomic loci and contribute to stress response and adaptation. In this review, we outline current relevant technological developments facilitating eccDNA identification and the latest discoveries about eccDNAs in plants. Finally, we explore the probable functions and future research directions that could be undertaken with respect to different eccDNA sources.

Sequence characterization of eccDNA content in glyphosate sensitive and resistant Palmer amaranth from geographically distant populations

The discovery of non-chromosomal circular DNA offers new directions in linking genome structure with function in plant biology. Glyphosate resistance through EPSPS gene copy amplification in Palmer amaranth was due to an autonomously replicating extra-chromosomal circular DNA mechanism (eccDNA). CIDER-Seq analysis of geographically distant glyphosate sensitive (GS) and resistant (GR) Palmer Amaranth (Amaranthus palmeri) revealed the presence of numerous small extra-chromosomal circular DNAs varying in size and with degrees of repetitive content, coding sequence, and motifs associated with autonomous replication. In GS biotypes, only a small portion of these aligned to the 399 kb eccDNA replicon, the vehicle underlying gene amplification and genetic resistance to the herbicide glyphosate. The aligned eccDNAs from GS were separated from one another by large gaps in sequence. In GR biotypes, the eccDNAs were present in both abundance and diversity to assemble into a nearly complete eccDNA replicon. Mean sizes of eccDNAs were similar in both biotypes and were around 5kb with larger eccDNAs near 25kb. Gene content for eccDNAs ranged from 0 to 3 with functions that include ribosomal proteins, transport, metabolism, and general stress response genetic elements. Repeat content among smaller eccDNAs indicate a potential for recombination into larger structures. Genomic hotspots were also identified in the Palmer amaranth genome with a disposition for gene focal amplifications as eccDNA. The presence of eccDNA may serve as a reservoir of genetic heterogeneity in this species and may be functionally important for survival.

Diverse and mobile: eccDNA-based identification of carrot low-copy-number LTR retrotransposons active in callus cultures

Long terminal repeat retrotransposons (LTR-RTs) are mobilized via an RNA intermediate using a 'copy and paste' mechanism, and account for the majority of repetitive DNA in plant genomes. As a side effect of mobilization, the formation of LTR-RT-derived extrachromosomal circular DNAs (eccDNAs) occurs. Thus, high-throughput sequencing of eccDNA can be used to identify active LTR-RTs in plant genomes. Despite the release of a reference genome assembly, carrot LTR-RTs have not yet been thoroughly characterized. LTR-RTs are abundant and diverse in the carrot genome. We identified 5976 carrot LTR-RTs, 2053 and 1660 of which were attributed to Copia and Gypsy superfamilies, respectively. They were further classified into lineages, families and subfamilies. More diverse LTR-RT lineages, i.e. lineages comprising many low-copy-number subfamilies, were more frequently associated with genic regions. Certain LTR-RT lineages have been recently active in Daucus carota. In particular, low-copy-number LTR-RT subfamilies, e.g. those belonging to the DcAle lineage, have significantly contributed to carrot genome diversity as a result of continuing activity. We utilized eccDNA sequencing to identify and characterize two DcAle subfamilies, Alex1 and Alex3, active in carrot callus. We documented 14 and 32 de novo insertions of Alex1 and Alex3, respectively, which were positioned in non-repetitive regions.

eccDNAdb: a database of extrachromosomal circular DNA profiles in human cancers

Extrachromosomal circular DNA (eccDNA) elements are circular DNA molecules that are derived from but are independent of chromosomal DNA. EccDNA is emerging as a rising star because of its ubiquitous existence in cancers and its crucial role in oncogene amplification and tumor progression. In the present study, whole-genome sequencing (WGS) data of cancer samples were downloaded from public repositories. Afterwards, eccDNAs were identified from WGS data via bioinformatic analyses. To leverage database coverage, eccDNAs were also collected by manual curation of literatures. Gene expression and clinical data were downloaded from TCGA and CCLE and then used to investigate the roles of eccDNAs in cancers. Finally, the first integrated database of eccDNAs, eccDNAdb, was developed. eccDNAdb currently includes 1270 eccDNAs, which were identified in 480 samples (of 42 cancers) after analyzing a total number of 3395 tumor samples (of 57 cancers) including patient tissues, patient-derived xenografts, and cancer cell lines. A total number of 54,901 eccDNA genes were annotated and included in the database as well. With the integration of gene expression, clinical information and chromatin accessibility data, eccDNAdb enables users to easily determine the biological function and clinical relevance of eccDNAs in human cancers. In conclusion, eccDNAdb is freely accessible at http://www.eccdnadb.org . To our knowledge, eccDNAdb is the first database in the eccDNA research field. It is expected to provide insight for novel cancer therapies.

ecc_finder: A Robust and Accurate Tool for Detecting Extrachromosomal Circular DNA From Sequencing Data

Extrachromosomal circular DNA (eccDNA) has been observed in different species for decades, and more and more evidence shows that this specific type of DNA molecules may play an important role in rapid adaptation. Therefore, characterizing the full landscape of eccDNA has become critical, and there are several protocols for enriching eccDNAs and performing short-read or long-read sequencing. However, there is currently no available bioinformatic tool to identify eccDNAs from Nanopore reads. More importantly, the current tools based on Illumina short reads lack an efficient standardized pipeline notably to identify eccDNA originating from repeated loci and cannot be applied to very large genomes. Here, we introduce a comprehensive tool to solve both of these two issues. Applying ecc_finder to eccDNA-seq data (either mobilome-seq, Circle-Seq and CIDER-seq) from Arabidopsis, human, and wheat (with genome sizes ranging from 120Mb to 17 Gb), we document the improvement of computational time, sensitivity, and accuracy and demonstrate ecc_finder wide applicability and functionality.

Quantitative assessment reveals the dominance of duplicated sequences in germline-derived extrachromosomal circular DNA

Extrachromosomal circular DNA (eccDNA) plays a role in human diseases such as cancer, but little is known about the impact of eccDNA in healthy human biology. Since eccDNA is a tiny fraction of nuclear DNA, artificial amplification has been employed to increase eccDNA amounts, resulting in the loss of native compositions. We developed an approach to enrich eccDNA populations at the native state (naïve small circular DNA, nscDNA) and investigated their origins in the human genome. We found that, in human sperm, the vast majority of nscDNA came from high-copy genomic regions, including the most variable regions between individuals. Because eccDNA can be incorporated back into chromosomes, eccDNA may promote human genetic variation.

Extrachromosomal circular DNA (eccDNA) originates from linear chromosomal DNA in various human tissues under physiological and disease conditions. The genomic origins of eccDNA have largely been investigated using in vitro–amplified DNA. However, in vitro amplification obscures quantitative information by skewing the total population stoichiometry. In addition, the analyses have focused on eccDNA stemming from single-copy genomic regions, leaving eccDNA from multicopy regions unexamined. To address these issues, we isolated eccDNA without in vitro amplification (naïve small circular DNA, nscDNA) and assessed the populations quantitatively by integrated genomic, molecular, and cytogenetic approaches. nscDNA of up to tens of kilobases were successfully enriched by our approach and were predominantly derived from multicopy genomic regions including segmental duplications (SDs). SDs, which account for 5% of the human genome and are hotspots for copy number variations, were significantly overrepresented in sperm nscDNA, with three times more sequencing reads derived from SDs than from the entire single-copy regions. SDs were also overrepresented in mouse sperm nscDNA, which we estimated to comprise 0.2% of nuclear DNA. Considering that eccDNA can be integrated into chromosomes, germline-derived nscDNA may be a mediator of genome diversity.

Interstitial Telomeric-like Repeats (ITR) in Seed Plants as Assessed by Molecular Cytogenetic Techniques: A Review

The discovery of telomeric repeats in interstitial regions of plant chromosomes (ITRs) through molecular cytogenetic techniques was achieved several decades ago. However, the information is scattered and has not been critically evaluated from an evolutionary perspective. Based on the analysis of currently available data, it is shown that ITRs are widespread in major evolutionary lineages sampled. However, their presence has been detected in only 45.6% of the analysed families, 26.7% of the sampled genera, and in 23.8% of the studied species. The number of ITR sites greatly varies among congeneric species and higher taxonomic units, and range from one to 72 signals. ITR signals mostly occurs as homozygous loci in most species, however, odd numbers of ITR sites reflecting a hemizygous state have been reported in both gymnosperm and angiosperm groups. Overall, the presence of ITRs appears to be poor predictors of phylogenetic and taxonomic relatedness at most hierarchical levels. The presence of ITRs and the number of sites are not significantly associated to the number of chromosomes. The longitudinal distribution of ITR sites along the chromosome arms indicates that more than half of the ITR presences are between proximal and terminal locations (49.5%), followed by proximal (29.0%) and centromeric (21.5%) arm regions. Intraspecific variation concerning ITR site number, chromosomal locations, and the differential presence on homologous chromosome pairs has been reported in unrelated groups, even at the population level. This hypervariability and dynamism may have likely been overlooked in many lineages due to the very low sample sizes often used in cytogenetic studies.

The Genomics of Plant Satellite DNA

The twenty-first century began with a certain indifference to the research of satellite DNA (satDNA). Neither genome sequencing projects were able to accurately encompass the study of satDNA nor classic methodologies were able to go further in undertaking a better comprehensive study of the whole set of satDNA sequences of a genome. Nonetheless, knowledge of satDNA has progressively advanced during this century with the advent of new analytical techniques. The enormous advantages that genome-wide approaches have brought to its analysis have now stimulated a renewed interest in the study of satDNA. At this point, we can look back and try to assess more accurately many of the key questions that were left unsolved in the past about this enigmatic and important component of the genome. I review here the understanding gathered on plant satDNAs over the last few decades with an eye on the near future.

Evolutionary History of Alpha Satellite DNA Repeats Dispersed within Human Genome Euchromatin

Major human alpha satellite DNA repeats are preferentially assembled within (peri)centromeric regions but are also dispersed within euchromatin in the form of clustered or short single repeat arrays. To study the evolutionary history of single euchromatic human alpha satellite repeats (ARs), we analyzed their orthologous loci across the primate genomes. The continuous insertion of euchromatic ARs throughout the evolutionary history of primates starting with the ancestors of Simiformes (45-60 Ma) and continuing up to the ancestors of Homo is revealed. Once inserted, the euchromatic ARs were stably transmitted to the descendant species, some exhibiting copy number variation, whereas their sequence divergence followed the species phylogeny. Many euchromatic ARs have sequence characteristics of (peri)centromeric alpha repeats suggesting heterochromatin as a source of dispersed euchromatic ARs. The majority of euchromatic ARs are inserted in the vicinity of other repetitive elements such as L1, Alu, and ERV or are embedded within them. Irrespective of the insertion context, each AR insertion seems to be unique and once inserted, ARs do not seem to be subsequently spread to new genomic locations. In spite of association with (retro)transposable elements, there is no indication that such elements play a role in ARs proliferation. The presence of short duplications at most of ARs insertion sites suggests site-directed recombination between homologous motifs in ARs and in the target genomic sequence, probably mediated by extrachromosomal circular DNA, as a mechanism of spreading within euchromatin.

Small extrachromosomal circular DNA (eccDNA): major functions in evolution and cancer

Extrachromosomal circular DNA (eccDNA) refers to a type of circular DNA that originate from but are likely independent of chromosomes. Due to technological advancements, eccDNAs have recently emerged as multifunctional molecules with numerous characteristics. The unique topological structure and genetic characteristics of eccDNAs shed new light on the monitoring, early diagnosis, treatment, and prediction of cancer. EccDNAs are commonly observed in both normal and cancer cells and function via different mechanisms in the stress response to exogenous and endogenous stimuli, aging, and carcinogenesis and in drug resistance during cancer treatment. The structural diversity of eccDNAs contributes to the function and numerical diversity of eccDNAs and thereby endows eccDNAs with powerful roles in evolution and in cancer initiation and progression by driving genetic plasticity and heterogeneity from extrachromosomal sites, which has been an ignored function in evolution in recent decades. EccDNAs show great potential in cancer, and we summarize the features, biogenesis, evaluated functions, functional mechanisms, related methods, and clinical utility of eccDNAs with a focus on their role in evolution and cancer.

The decreased exclusion of nuclear eccDNA: From molecular and subcellular levels to human aging and age-related diseases

Extrachromosomal circular DNA (eccDNA) accumulates within the nucleus of eukaryotic cells during physiological aging and in age-related diseases (ARDs) and the accumulation could be caused by the declined exclusion of nuclear eccDNA in these states. This review focuses on the formation of eccDNA and the roles of some main factors, such as nuclear pore complexes (NPCs), nucleoplasmic reticulum (NR), and nuclear actin, in eccDNA exclusion. eccDNAs are mostly formed from non-coding DNA during DNA damage repair. They move to NPCs along nuclear actin and are excluded out of the nucleus through functional NPCs in young and healthy cells. However, it has been demonstrated that defective NPCs, abnormal NPC components and nuclear actin rods are increased in aged cells, various cancers and certain other ARDs such as cardiovascular diseases, premature aging, neurodegenerative diseases and myopathies. Therefore, mainly resulting from the increase of dysfunctional NPCs, the exclusion of nuclear eccDNAs may be reduced and eccDNAs thus accumulate within the nucleus in aging and the aforementioned ARDs. In addition, the protective function of non-coding DNA in tumorigenesis is further discussed.

The Role of Satellite DNAs in Genome Architecture and Sex Chromosome Evolution in Crambidae Moths

Tandem repeats are important parts of eukaryotic genomes being crucial e.g., for centromere and telomere function and chromatin modulation. In Lepidoptera, knowledge of tandem repeats is very limited despite the growing number of sequenced genomes. Here we introduce seven new satellite DNAs (satDNAs), which more than doubles the number of currently known lepidopteran satDNAs. The satDNAs were identified in genomes of three species of Crambidae moths, namely Ostrinia nubilalis, Cydalima perspectalis, and Diatraea postlineella, using graph-based computational pipeline RepeatExplorer. These repeats varied in their abundance and showed high variability within and between species, although some degree of conservation was noted. The satDNAs showed a scattered distribution, often on both autosomes and sex chromosomes, with the exception of both satellites in D. postlineella, in which the satDNAs were located at a single autosomal locus. Three satDNAs were abundant on the W chromosomes of O. nubilalis and C. perspectalis, thus contributing to their differentiation from the Z chromosomes. To provide background for the in situ localization of the satDNAs, we performed a detailed cytogenetic analysis of the karyotypes of all three species. This comparative analysis revealed differences in chromosome number, number and location of rDNA clusters, and molecular differentiation of sex chromosomes.

Satellite DNA-Mediated Gene Expression Regulation: Physiological and Evolutionary Implication

Satellite DNAs are tandemly repeated sequences organized in large clusters within (peri)centromeric and/or subtelomeric heterochromatin. However, in many species, satellite DNAs are not restricted to heterochromatin but are also dispersed as short arrays within euchromatin. Such genomic organization together with transcriptional activity seems to be a prerequisite for the gene-modulatory effect of satellite DNAs which was first demonstrated in the beetle Tribolium castaneum upon heat stress. Namely, enrichment of a silent histone mark at euchromatic repeats of a major beetle satellite DNA results in epigenetic silencing of neighboring genes. In addition, human satellite III transcripts induced by heat shock contribute to genome-wide gene silencing, providing protection against stress-induced cell death. Gene silencing mediated by satellite RNA was also shown to be fundamental for the early embryonic development of the mosquito Aedes aegypti. Apart from a physiological role during embryogenesis and heat stress response, activation of satellite DNAs in terms of transcription and proliferation can have an evolutionary impact. Spreading of satellite repeats throughout euchromatin promotes the variation of epigenetic landscapes and gene expression diversity, contributing to the evolution of gene regulatory networks and to genome adaptation in fluctuating environmental conditions.

Dynamic Evolution of Euchromatic Satellites on the X Chromosome in Drosophila melanogaster and the simulans Clade

Satellite DNAs (satDNAs) are among the most dynamically evolving components of eukaryotic genomes and play important roles in genome regulation, genome evolution, and speciation. Despite their abundance and functional impact, we know little about the evolutionary dynamics and molecular mechanisms that shape satDNA distributions in genomes. Here, we use high-quality genome assemblies to study the evolutionary dynamics of two complex satDNAs, Rsp-like and 1.688 g/cm3, in Drosophila melanogaster and its three nearest relatives in the simulans clade. We show that large blocks of these repeats are highly dynamic in the heterochromatin, where their genomic location varies across species. We discovered that small blocks of satDNA that are abundant in X chromosome euchromatin are similarly dynamic, with repeats changing in abundance, location, and composition among species. We detail the proliferation of a rare satellite (Rsp-like) across the X chromosome in D. simulans and D. mauritiana. Rsp-like spread by inserting into existing clusters of the older, more abundant 1.688 satellite, in events likely facilitated by microhomology-mediated repair pathways. We show that Rsp-like is abundant on extrachromosomal circular DNA in D. simulans, which may have contributed to its dynamic evolution. Intralocus satDNA expansions via unequal exchange and the movement of higher order repeats also contribute to the fluidity of the repeat landscape. We find evidence that euchromatic satDNA repeats experience cycles of proliferation and diversification somewhat analogous to bursts of transposable element proliferation. Our study lays a foundation for mechanistic studies of satDNA proliferation and the functional and evolutionary consequences of satDNA movement.

Near-Random Distribution of Chromosome-Derived Circular DNA in the Condensed Genome of Pigeons and the Larger, More Repeat-Rich Human Genome

Extrachromosomal circular DNA (eccDNA) elements of chromosomal origin are known to be common in a number of eukaryotic species. However, it remains to be addressed whether genomic features such as genome size, the load of repetitive elements within a genome, and/or animal physiology affect the number of eccDNAs. Here, we investigate the distribution and numbers of eccDNAs in a condensed and less repeat-rich genome compared with the human genome, using Columba livia domestica (domestic rock pigeon) as a model organism. By sequencing eccDNA in blood and breast muscle from three pigeon breeds at various ages and with different flight behavior, we characterize 30,000 unique eccDNAs. We identify genomic regions that are likely hotspots for DNA circularization in breast muscle, including genes involved in muscle development. We find that although eccDNA counts do not correlate with the biological age in pigeons, the number of unique eccDNAs in a nonflying breed (king pigeons) is significantly higher (9-fold) than homing pigeons. Furthermore, a comparison between eccDNA from skeletal muscle in pigeons and humans reveals ∼9-10 times more unique eccDNAs per human nucleus. The fraction of eccDNA sequences, derived from repetitive elements, exist in proportions to genome content, that is, human 72.4% (expected 52.5%) and pigeon 8.7% (expected 5.5%).

Overall, our results support that eccDNAs are common in pigeons, that the amount of unique eccDNA types per nucleus can differ between species as well as subspecies, and suggest that eccDNAs from repeats are found in proportions relative to the content of repetitive elements in a genome.

Heteromorphic Sex Chromosomes and Their DNA Content in Fish: An Insight through Satellite DNA Accumulation in Megaleporinus elongatus

The repetitive DNA content of fish sex chromosomes provides valuable insights into specificities and patterns of their genetic sex determination systems. In this study, we revealed the genomic satellite DNA (satDNA) content of Megaleporinuselongatus, a Neotropical fish species with Z1Z1Z2Z2/Z1W1Z2W2 multiple sex chromosomes, through high-throughput analysis and graph-based clustering, isolating 68 satDNA families. By physically mapping these sequences in female metaphases, we discovered 15 of the most abundant satDNAs clustered in its chromosomes, 9 of which were found exclusively in the highly heterochromatic W1. This heteromorphic sex chromosome showed the highest amount of satDNA accumulations in this species. The second most abundant family, MelSat02-26, shared FISH signals with the NOR-bearing pair in similar patterns and is linked to the multiple sex chromosome system. Our results demonstrate the diverse satDNA content in M. elongatus, especially in its heteromorphic sex chromosome. Additionally, we highlighted the different accumulation patterns and distribution of these sequences across species by physically mapping these satDNAs in other Anostomidae, Megaleporinusmacrocephalus and Leporinusfriderici (a species without differentiated sex chromosomes).

What makes a centromere?

Centromeres are the eukaryotic chromosomal sites at which the kinetochore forms and attaches to spindle microtubules to orchestrate chromosomal segregation in mitosis and meiosis. Although centromeres are essential for cell division, their sequences are not conserved and evolve rapidly. Centromeres vary dramatically in size and organization. Here we categorize their diversity and explore the evolutionary forces shaping them. Nearly all centromeres favor AT-rich DNA that is gene-free and transcribed at a very low level. Repair of frequent centromere-proximal breaks probably contributes to their rapid sequence evolution. Point centromeres are only ~125 bp and are specified by common protein-binding motifs, whereas short regional centromeres are 1-5 kb, typically have unique sequences, and may have pericentromeric repeats adapted to facilitate centromere clustering. Transposon-rich centromeres are often ~100-300 kb and are favored by RNAi machinery that silences transposons, by suppression of meiotic crossovers at centromeres, and by the ability of some transposons to target centromeres. Megabase-length satellite centromeres arise in plants and animals with asymmetric female meiosis that creates centromere competition, and favors satellite monomers one or two nucleosomes in length that position and stabilize centromeric nucleosomes. Holocentromeres encompass the length of a chromosome and may differ dramatically between mitosis and meiosis. We propose a model in which low level transcription of centromeres facilitates the formation of non-B DNA that specifies centromeres and promotes loading of centromeric nucleosomes.

Characterization of repeat arrays in ultra-long nanopore reads reveals frequent origin of satellite DNA from retrotransposon-derived tandem repeats

Amplification of monomer sequences into long contiguous arrays is the main feature distinguishing satellite DNA from other tandem repeats, yet it is also the main obstacle in its investigation because these arrays are in principle difficult to assemble. Here we explore an alternative, assembly-free approach that utilizes ultra-long Oxford Nanopore reads to infer the length distribution of satellite repeat arrays, their association with other repeats and the prevailing sequence periodicities. Using the satellite DNA-rich legume plant Lathyrus sativus as a model, we demonstrated this approach by analyzing 11 major satellite repeats using a set of nanopore reads ranging from 30 to over 200 kb in length and representing 0.73× genome coverage. We found surprising differences between the analyzed repeats because only two of them were predominantly organized in long arrays typical for satellite DNA. The remaining nine satellites were found to be derived from short tandem arrays located within LTR-retrotransposons that occasionally expanded in length. While the corresponding LTR-retrotransposons were dispersed across the genome, this array expansion occurred mainly in the primary constrictions of the L. sativus chromosomes, which suggests that these genome regions are favourable for satellite DNA accumulation.

Programmed genome rearrangements in Oxytricha produce transcriptionally active extrachromosomal circular DNA

Extrachromosomal circular DNA (eccDNA) is both a driver of eukaryotic genome instability and a product of programmed genome rearrangements, but its extent had not been surveyed in Oxytricha, a ciliate with elaborate DNA elimination and translocation during development. Here, we captured rearrangement-specific circular DNA molecules across the genome to gain insight into its processes of programmed genome rearrangement. We recovered thousands of circularly excised Tc1/mariner-type transposable elements and high confidence non-repetitive germline-limited loci. We verified their bona fide circular topology using circular DNA deep-sequencing, 2D gel electrophoresis and inverse polymerase chain reaction. In contrast to the precise circular excision of transposable elements, we report widespread heterogeneity in the circular excision of non-repetitive germline-limited loci. We also demonstrate that circular DNAs are transcribed in Oxytricha, producing rearrangement-specific long non-coding RNAs. The programmed formation of thousands of eccDNA molecules makes Oxytricha a model system for studying nucleic acid topology. It also suggests involvement of eccDNA in programmed genome rearrangement.

Adding color to a century-old enigma: multi-color chromosome identification unravels the autotriploid nature of saffron (Crocus sativus) as a hybrid of wild Crocus cartwrightianus cytotypes

Saffron crocus (Crocus sativus) is the source of the most expensive spice of the world, produced from manually harvested stigmas, thus serving as a cash crop for rural communities. However, despite its economic importance, its genome and chromosomes are poorly studied. C. sativus is a sterile triploid species harboring eight chromosome triplets, and propagated only as a clonal lineage by corms. Saffron's evolutionary origin, parental species and allo- or autotriploidy has been a matter of discussion for almost a century. We performed a survey sequencing of the saffron genome and selected cytogenetic landmark sequences consisting of major tandem repeats, which we used as probes in comparative multicolor fluorescent in situ hybridization (FISH). We tagged 92 chromosomal positions and resolved the chromosomal composition of saffron triplets. By comparative FISH of six Crocus species from 11 accessions, we demonstrate that C. sativus is an autotriploid hybrid derived from heterogeneous Crocus cartwrightianus cytotypes. The FISH reference karyotype of saffron is crucial for integrating genome sequencing data with chromosomes and for investigating the relationship among Crocus species. We provide an evolutionary model of the saffron emergence; the knowledge of the parental origin offers a route towards the resynthesis of C. sativus from C. cartwrightianus to broaden saffron's gene pool.

Characterization of repeated DNA sequences in genomes of blue-flowered flax

BACKGROUND

Members of different sections of the genus Linum are characterized by wide variability in size, morphology and number of chromosomes in karyotypes. Since such variability is determined mainly by the amount and composition of repeated sequences, we conducted a comparative study of the repeatomes of species from four sections forming a clade of blue-flowered flax. Based on the results of high-throughput genome sequencing performed in this study as well as available WGS data, bioinformatic analyses of repeated sequences from 12 flax samples were carried out using a graph-based clustering method.

RESULTS

It was found that the genomes of closely related species, which have a similar karyotype structure, are also similar in the repeatome composition. In contrast, the repeatomes of karyologically distinct species differed significantly, and no similar tandem-organized repeats have been identified in their genomes. At the same time, many common mobile element families have been identified in genomes of all species, among them, Athila Ty3/gypsy LTR retrotransposon was the most abundant. The 30-chromosome members of the sect. Linum (including the cultivated species L. usitatissimum) differed significantly from other studied species by a great number of satellite DNA families as well as their relative content in genomes.

CONCLUSIONS

The evolution of studied flax species was accompanied by waves of amplification of satellite DNAs and LTR retrotransposons. The observed inverse correlation between the total contents of dispersed repeats and satellite DNAs allowed to suggest a relationship between both classes of repeating sequences. Significant interspecific differences in satellite DNA sets indicated a high rate of evolution of this genomic fraction. The phylogenetic relationships between the investigated flax species, obtained by comparison of the repeatomes, agreed with the results of previous molecular phylogenetic studies.

Satellite DNA in Vicia faba is characterized by remarkable diversity in its sequence composition, association with centromeres, and replication timing

Satellite DNA, a class of repetitive sequences forming long arrays of tandemly repeated units, represents substantial portions of many plant genomes yet remains poorly characterized due to various methodological obstacles. Here we show that the genome of the field bean (Vicia faba, 2n = 12), a long-established model for cytogenetic studies in plants, contains a diverse set of satellite repeats, most of which remained concealed until their present investigation. Using next-generation sequencing combined with novel bioinformatics tools, we reconstructed consensus sequences of 23 novel satellite repeats representing 0.008–2.700% of the genome and mapped their distribution on chromosomes. We found that in addition to typical satellites with monomers hundreds of nucleotides long, V. faba contains a large number of satellite repeats with unusually long monomers (687–2033 bp), which are predominantly localized in pericentromeric regions. Using chromatin immunoprecipitation with CenH3 antibody, we revealed an extraordinary diversity of centromeric satellites, consisting of seven repeats with chromosome-specific distribution. We also found that in spite of their different nucleotide sequences, all centromeric repeats are replicated during mid-S phase, while most other satellites are replicated in the first part of late S phase, followed by a single family of FokI repeats representing the latest replicating chromatin.

Discoveries of Extrachromosomal Circles of DNA in Normal and Tumor Cells

While the vast majority of cellular DNA in eukaryotes is contained in long linear strands in chromosomes, we have long recognized some exceptions like mitochondrial DNA, plasmids in yeasts, and double minutes (DMs) in cancer cells where the DNA is present in extrachromosomal circles. In addition, specialized extrachromosomal circles of DNA (eccDNA) have been noted to arise from repetitive genomic sequences like telomeric DNA or rDNA. Recently eccDNA arising from unique (nonrepetitive) DNA have been discovered in normal and malignant cells, raising interesting questions about their biogenesis, function and clinical utility. Here, we review recent results and future directions of inquiry on these new forms of eccDNA.

High-throughput analysis of satellite DNA in the grasshopper Pyrgomorpha conica reveals abundance of homologous and heterologous higher-order repeats

Satellite DNA (satDNA) constitutes an important fraction of repetitive DNA in eukaryotic genomes, but it is barely known in most species. The high-throughput analysis of satDNA in the grasshopper Pyrgomorpha conica revealed 87 satDNA variants grouped into 76 different families, representing 9.4% of the genome. Fluorescent in situ hybridization (FISH) analysis of the 38 most abundant satDNA families revealed four different patterns of chromosome distribution. Homology search between the 76 satDNA families showed the existence of 15 superfamilies, each including two or more families, with the most abundant superfamily representing more than 80% of all satDNA found in this species. This also revealed the presence of two types of higher-order repeats (HORs), one showing internal homologous subrepeats, as conventional HORs, and an additional type showing non-homologous internal subrepeats, the latter arising by the combination of a given satDNA family with a non-annotated sequence, or with telomeric DNA. Interestingly, the heterologous subrepeats included in these HORs showed higher divergence within the HOR than outside it, suggesting that heterologous HORs show poor homogenization, in high contrast with conventional (homologous) HORs. Finally, heterologous HORs can show high differences in divergence between their constituent subrepeats, suggesting the possibility of regional homogenization.

Circular DNA elements of chromosomal origin are common in healthy human somatic tissue

The human genome is generally organized into stable chromosomes, and only tumor cells are known to accumulate kilobase (kb)-sized extrachromosomal circular DNA elements (eccDNAs). However, it must be expected that kb eccDNAs exist in normal cells as a result of mutations. Here, we purify and sequence eccDNAs from muscle and blood samples from 16 healthy men, detecting ~100,000 unique eccDNA types from 16 million nuclei. Half of these structures carry genes or gene fragments and the majority are smaller than 25 kb. Transcription from eccDNAs suggests that eccDNAs reside in nuclei and recurrence of certain eccDNAs in several individuals implies DNA circularization hotspots. Gene-rich chromosomes contribute to more eccDNAs per megabase and the most transcribed protein-coding gene in muscle, TTN (titin), provides the most eccDNAs per gene. Thus, somatic genomes are rich in chromosome-derived eccDNAs that may influence phenotypes through altered gene copy numbers and transcription of full-length or truncated genes.

Somatic cells can accumulate structural variations such as deletions. Here, Møller et al. show that normal human cells generate large extrachromosomal circular DNAs (eccDNAs), most likely the products of excised DNA, that can be transcriptionally active and, thus, may have phenotypic consequences.

Dispersion Profiles and Gene Associations of Repetitive DNAs in the Euchromatin of the Beetle Tribolium castaneum

Satellite DNAs are tandemly repeated sequences clustered within heterochromatin. However, in some cases, such as the major TCAST1 satellite DNA from the beetle Tribolium castaneum, they are found partially dispersed within euchromatin. Such organization together with transcriptional activity enables TCAST1 to modulate the activity of neighboring genes. In order to explore if other T. castaneum repetitive families have features that could provide them with a possible gene-modulatory role, we compare here the structure, organization, dispersion profiles, and transcription activity of 10 distinct TCAST repetitive families including TCAST1. The genome organization of TCAST families exhibit either satellite-like or transposon-like characteristics. In addition to heterochromatin localization, bioinformatic searches of the assembled genome have revealed dispersion of all families within euchromatin, preferentially in the form of single repeats. Dispersed TCAST repeats are mutually correlated in distribution and are grouped in distinct regions of euchromatin. The repeats are associated with genes, are enriched in introns relative to intergenic regions, and very rarely overlap exons. In spite of the different mechanisms of repeat proliferation, such as transposition and homologous recombination, all TCAST families share a similar frequency of spreading as well as dispersion and gene association profiles. Additionally, TCAST families are transcribed and their transcription is significantly activated by heat stress. A possibility that such common features of TCAST families might be related to their potential gene-modulatory role is discussed.

Fine organization of genomic regions tagged to the 5S rDNA locus of the bread wheat 5B chromosome

BACKGROUND

The multigene family encoding the 5S rRNA, one of the most important structurally-functional part of the large ribosomal subunit, is an obligate component of all eukaryotic genomes. 5S rDNA has long been a favored target for cytological and phylogenetic studies due to the inherent peculiarities of its structural organization, such as the tandem arrays of repetitive units and their high interspecific divergence. The complex polyploid nature of the genome of bread wheat, Triticum aestivum, and the technically difficult task of sequencing clusters of tandem repeats mean that the detailed organization of extended genomic regions containing 5S rRNA genes remains unclear. This is despite the recent progress made in wheat genomic sequencing. Using pyrosequencing of BAC clones, in this work we studied the organization of two distinct 5S rDNA-tagged regions of the 5BS chromosome of bread wheat.

RESULTS

Three BAC-clones containing 5S rDNA were identified in the 5BS chromosome-specific BAC-library of Triticum aestivum. Using the results of pyrosequencing and assembling, we obtained six 5S rDNA- containing contigs with a total length of 140,417 bp, and two sets (pools) of individual 5S rDNA sequences belonging to separate, but closely located genomic regions on the 5BS chromosome. Both regions are characterized by the presence of approximately 70-80 copies of 5S rDNA, however, they are completely different in their structural organization. The first region contained highly diverged short-type 5S rDNA units that were disrupted by multiple insertions of transposable elements. The second region contained the more conserved long-type 5S rDNA, organized as a single tandem array. FISH using probes specific to both 5S rDNA unit types showed differences in the distribution and intensity of signals on the chromosomes of polyploid wheat species and their diploid progenitors.

CONCLUSION

A detailed structural organization of two closely located 5S rDNA-tagged genomic regions on the 5BS chromosome of bread wheat has been established. These two regions differ in the organization of both 5S rDNA and the neighboring sequences comprised of transposable elements, implying different modes of evolution for these regions.

Molecular characterization of cell-free eccDNAs in human plasma

Extrachromosomal circular DNAs (eccDNAs) have been reported in most eukaryotes. However, little is known about the cell-free eccDNA profiles in circulating system such as blood. To characterize plasma cell-free eccDNAs, we performed sequencing analysis in 26 libraries from three blood donors and negative controls. We identified thousands of unique plasma eccDNAs in the three subjects. We observed proportional eccDNA increase with initial DNA input. The detected eccDNAs were also associated with circular DNA enrichment efficiency. Increasing the sequencing depth in an additional sample identified many more eccDNAs with highly heterogenous molecular structure. Size distribution of eccDNAs varied significantly from 31 bp to 19,989 bp. We found significantly higher GC content in smaller eccDNAs (<500 bp) than the larger ones (>500 bp) (p < 0.01). We also found an enrichment of eccDNAs at exons and 3′UTR (enrichment folds from 1.36 to 3.1) as well as the DNase hypersensitive sites (1.58–2.42 fold), H3K4Me1 (1.23–1.42 fold) and H3K27Ac (1.33–1.62 fold) marks. Junction sequence analysis suggested fundamental role of nonhomologous end joining mechanism during eccDNA formation. Further characterization of the extracellular eccDNAs in peripheral blood will facilitate understanding of their molecular mechanisms and potential clinical utilities.

Characterization of the microDNA through the response to chemotherapeutics in lymphoblastoid cell lines

Recently, a new class of extrachromosomal circular DNA, called microDNA, was identified. They are on average 100 to 400 bp long and are derived from unique non-repetitive genomic regions with high gene density. MicroDNAs are thought to arise from DNA breaks associated with RNA metabolism or replication slippage. Given the paucity of information on this entirely novel phenomenon, we aimed to get an additional insight into microDNA features by performing the microDNA analysis in 20 independent human lymphoblastoid cell lines (LCLs) prior and after treatment with chemotherapeutic drugs. The results showed non-random genesis of microDNA clusters from the active regions of the genome. The size periodicity of 190 bp was observed, which matches DNA fragmentation typical for apoptotic cells. The chemotherapeutic drug-induced apoptosis of LCLs increased both number and size of clusters further suggesting that part of microDNAs could result from the programmed cell death. Interestingly, proportion of identified microDNA sequences has common loci of origin when compared between cell line experiments. While compatible with the original observation that microDNAs originate from a normal physiological process, obtained results imply complementary source of its production. Furthermore, non-random genesis of microDNAs depicted by redundancy between samples makes these entities possible candidates for new biomarker generation.

The satellite DNA AflaSAT-1 in the A and B chromosomes of the grasshopper Abracris flavolineata

BACKGROUND

Satellite DNAs (satDNAs) are organized in repetitions directly contiguous to one another, forming long arrays and composing a large portion of eukaryote genomes. These sequences evolve according to the concerted evolution model, and homogenization of repeats is observed at the intragenomic level. Satellite DNAs are the primary component of heterochromatin, located primarily in centromeres and telomeres. Moreover, satDNA enrichment in specific chromosomes has been observed, such as in B chromosomes, that can provide clues about composition, origin and evolution of this chromosome. In this study, we isolated and characterized a satDNA in A and B chromosomes of Abracris flavolineata by integrating cytogenetic, molecular and genomics approaches at intra- and inter-population levels, with the aim to understand the evolution of satDNA and composition of B chromosomes.

RESULTS

AflaSAT-1 satDNA was shared with other species and in A. flavolineata, was associated with another satDNA, AflaSAT-2. Chromosomal mapping revealed centromeric blocks variable in size in almost all chromosomes (except pair 11) of A complement for both satDNAs, whereas for B chromosome, only a small centromeric signal occurred. In distinct populations, variable number of AflaSAT-1 chromosomal sites correlated with variability in copy number. Instead of such variability, low sequence diversity was observed in A complement, but monomers from B chromosome were more variable, presenting also exclusive mutations. AflaSAT-1 was transcribed in five tissues of adults in distinct life cycle phases.

CONCLUSIONS

The sharing of AflaSAT-1 with other species is consistent with the library hypothesis and indicates common origin in a common ancestor; however, AflaSAT-1 was highly amplified in the genome of A. flavolineata. At the population level, homogenization of repeats in distinct populations was documented, but dynamic expansion or elimination of repeats was also observed. Concerning the B chromosome, our data provided new information on the composition in A. flavolineata. Together with previous results, the sequences of heterochromatic nature were not likely highly amplified in the entire B chromosome. Finally, the constitutive transcriptional activity suggests a possible unknown functional role, which should be further investigated.

Higher-order organisation of extremely amplified, potentially functional and massively methylated 5S rDNA in European pikes (Esox sp.)

Background

Pikes represent an important genus (Esox) harbouring a pre-duplication karyotype (2n = 2x = 50) of economically important salmonid pseudopolyploids. Here, we have characterized the 5S ribosomal RNA genes (rDNA) in Esox lucius and its closely related E. cisalpinus using cytogenetic, molecular and genomic approaches. Intragenomic homogeneity and copy number estimation was carried out using Illumina reads. The higher-order structure of rDNA arrays was investigated by the analysis of long PacBio reads. Position of loci on chromosomes was determined by FISH. DNA methylation was analysed by methylation-sensitive restriction enzymes.

Results

The 5S rDNA loci occupy exclusively (peri)centromeric regions on 30–38 acrocentric chromosomes in both E. lucius and E. cisalpinus. The large number of loci is accompanied by extreme amplification of genes (>20,000 copies), which is to the best of our knowledge one of the highest copy number of rRNA genes in animals ever reported. Conserved secondary structures of predicted 5S rRNAs indicate that most of the amplified genes are potentially functional. Only few SNPs were found in genic regions indicating their high homogeneity while intergenic spacers were more heterogeneous and several families were identified. Analysis of 10–30 kb-long molecules sequenced by the PacBio technology (containing about 40% of total 5S rDNA) revealed that the vast majority (96%) of genes are organised in large several kilobase-long blocks. Dispersed genes or short tandems were less common (4%). The adjacent 5S blocks were directly linked, separated by intervening DNA and even inverted. The 5S units differing in the intergenic spacers formed both homogeneous and heterogeneous (mixed) blocks indicating variable degree of homogenisation between the loci. Both E. lucius and E. cisalpinus 5S rDNA was heavily methylated at CG dinucleotides.

Conclusions

Extreme amplification of 5S rRNA genes in the Esox genome occurred in the absence of significant pseudogenisation suggesting its recent origin and/or intensive homogenisation processes. The dense methylation of units indicates that powerful epigenetic mechanisms have evolved in this group of fish to silence amplified genes. We discuss how the higher-order repeat structures impact on homogenisation of 5S rDNA in the genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3774-7) contains supplementary material, which is available to authorized users.

Gene conversion events and variable degree of homogenization of rDNA loci in cultivars of Brassica napus

BACKGROUND AND AIMS

Brassica napus (AACC, 2n = 38, oilseed rape) is a relatively recent allotetraploid species derived from the putative progenitor diploid species Brassica rapa (AA, 2n = 20) and Brassica oleracea (CC, 2n = 18). To determine the influence of intensive breeding conditions on the evolution of its genome, we analysed structure and copy number of rDNA in 21 cultivars of B. napus, representative of genetic diversity.

METHODS

We used next-generation sequencing genomic approaches, Southern blot hybridization, expression analysis and fluorescence in situ hybridization (FISH). Subgenome-specific sequences derived from rDNA intergenic spacers (IGS) were used as probes for identification of loci composition on chromosomes.

KEY RESULTS

Most B. napus cultivars (18/21, 86 %) had more A-genome than C-genome rDNA copies. Three cultivars analysed by FISH ('Darmor', 'Yudal' and 'Asparagus kale') harboured the same number (12 per diploid set) of loci. In B. napus 'Darmor', the A-genome-specific rDNA probe hybridized to all 12 rDNA loci (eight on the A-genome and four on the C-genome) while the C-genome-specific probe showed weak signals on the C-genome loci only. Deep sequencing revealed high homogeneity of arrays suggesting that the C-genome genes were largely overwritten by the A-genome variants in B. napus 'Darmor'. In contrast, B. napus 'Yudal' showed a lack of gene conversion evidenced by additive inheritance of progenitor rDNA variants and highly localized hybridization signals of subgenome-specific probes on chromosomes. Brassica napus 'Asparagus kale' showed an intermediate pattern to 'Darmor' and 'Yudal'. At the expression level, most cultivars (95 %) exhibited stable A-genome nucleolar dominance while one cultivar ('Norin 9') showed co-dominance.

CONCLUSIONS

The B. napus cultivars differ in the degree and direction of rDNA homogenization. The prevalent direction of gene conversion (towards the A-genome) correlates with the direction of expression dominance indicating that gene activity may be needed for interlocus gene conversion.

Structural and functional liaisons between transposable elements and satellite DNAs

Transposable elements (TEs) and satellite DNAs (satDNAs) are typically identified as major repetitive DNA components in eukaryotic genomes. TEs are DNA segments able to move throughout a genome while satDNAs are tandemly repeated sequences organized in long arrays. Both classes of repetitive sequences are extremely diverse, and many TEs and satDNAs exist within a genome. Although they differ in structure, genomic organization, mechanisms of spread, and evolutionary dynamics, TEs and satDNAs can share sequence similarity and organizational patterns, thus indicating that complex mutual relationships can determine their evolution, and ultimately define roles they might have on genome architecture and function. Motivated by accumulating data about sequence elements that incorporate features of both TEs and satDNAs, here we present an overview of their structural and functional liaisons.

High-throughput analysis of the satellitome illuminates satellite DNA evolution

Satellite DNA (satDNA) is a major component yet the great unknown of eukaryote genomes and clearly underrepresented in genome sequencing projects. Here we show the high-throughput analysis of satellite DNA content in the migratory locust by means of the bioinformatic analysis of Illumina reads with the RepeatExplorer and RepeatMasker programs. This unveiled 62 satDNA families and we propose the term "satellitome" for the whole collection of different satDNA families in a genome. The finding that satDNAs were present in many contigs of the migratory locust draft genome indicates that they show many genomic locations invisible by fluorescent in situ hybridization (FISH). The cytological pattern of five satellites showing common descent (belonging to the SF3 superfamily) suggests that non-clustered satDNAs can become into clustered through local amplification at any of the many genomic loci resulting from previous dissemination of short satDNA arrays. The fact that all kinds of satDNA (micro- mini- and satellites) can show the non-clustered and clustered states suggests that all these elements are mostly similar, except for repeat length. Finally, the presence of VNTRs in bacteria, showing similar properties to non-clustered satDNAs in eukaryotes, suggests that this kind of tandem repeats show common properties in all living beings.

Identification of Extrachromosomal Circular DNA in Hop via Rolling Circle Amplification

During a survey for new viruses affecting hop plants, a circular DNA molecule was identified via rolling circle amplification (RCA) and later characterized. A small region of the 5.7-kb long molecule aligned with a microsatellite region in the Humulus lupulus genome, and no coding sequence was identified. Sequence analysis and literature review suggest that the small DNA molecule is an extranuclear DNA element, specifically, an extrachromosomal circular DNA (eccDNA), and its presence was confirmed by electron microscopy. This work is the first report of eccDNAs in the family Cannabaceae. Additionally, this work highlights the advantages of using RCA to study extrachromosomal DNA in higher plants.

Structural polymorphisms and distinct genomic composition suggest recurrent origin and ongoing evolution of B chromosomes in the Prospero autumnale complex (Hyacinthaceae)

Supernumerary B chromosomes (Bs) are genomic parasitic components, originating from the A complement via chromosomal rearrangements, which follow their own evolutionary trajectories. They often contain repetitive DNAs, some shared with regular chromosomes and some newly evolved. Genomic composition, origin and evolution of Bs have been analysed in the chromosomally variable Prospero autumnale complex. Two rDNAs and a satellite DNA (PaB6) from regular chromosomes were mapped to Bs of 26 plants from three diploid cytotypes, their hybrids and polyploid derivatives. In homoploid diploid hybrids, genomic in situ hybridization (GISH) allowed B painting with the parental DNAs. Bs were structurally variable and highly enriched in 5S rDNA and satDNA PaB6, and rarely in 35S rDNA. Eleven combinations of rDNA and PaB6 localization were observed. The quantities of PaB6 in Bs and regular chromosomes were not correlated, suggesting amplification mechanisms other than recombination. PaB6 and 5S rDNA amounts increased with increasing ploidy level. GISH revealed two independent origins of Bs. The structural variation, repeat content, repeat-type fluctuations and differing genomic affinities of Bs in different cytotypes suggest that they represent young proto-B chromosomes. Bs in P. autumnale probably form recurrently as by-products of the extensive genome restructuring within this chromosomally variable species complex.

Gene stacking by recombinases

Efficient methods of stacking genes into plant genomes are needed to expedite transfer of multigenic traits to crop varieties of diverse ecosystems. Over two decades of research has identified several DNA recombinases that carryout efficient cis and trans recombination between the recombination sites artificially introduced into the plant chromosome. The specificity and efficiency of recombinases make them extremely attractive for genome engineering. In plant biotechnology, recombinases have mostly been used for removing selectable marker genes and have rarely been extended to more complex applications. The reversibility of recombination, a property of the tyrosine family of recombinases, does not lend itself to gene stacking approaches that involve rounds of transformation for integrating genes into the engineered sites. However, recent developments in the field of recombinases have overcome these challenges and paved the way for gene stacking. Some of the key advancements include the application of unidirectional recombination systems, modification of recombination sites and transgene site modifications to allow repeated site-specific integrations into the selected site. Gene stacking is relevant to agriculturally important crops, many of which are difficult to transform; therefore, development of high-efficiency gene stacking systems will be important for its application on agronomically important crops, and their elite varieties. Recombinases, by virtue of their specificity and efficiency in plant cells, emerge as powerful tools for a variety of applications including gene stacking.

Genome-wide analysis of tandem repeats in Tribolium castaneum genome reveals abundant and highly dynamic tandem repeat families with satellite DNA features in euchromatic chromosomal arms

Although satellite DNAs are well-explored components of heterochromatin and centromeres, little is known about emergence, dispersal and possible impact of comparably structured tandem repeats (TRs) on the genome-wide scale. Our bioinformatics analysis of assembled Tribolium castaneum genome disclosed significant contribution of TRs in euchromatic chromosomal arms and clear predominance of satellite DNA-typical 170 bp monomers in arrays of ≥5 repeats. By applying different experimental approaches, we revealed that the nine most prominent TR families Cast1-Cast9 extracted from the assembly comprise ∼4.3% of the entire genome and reside almost exclusively in euchromatic regions. Among them, seven families that build ∼3.9% of the genome are based on ∼170 and ∼340 bp long monomers. Results of phylogenetic analyses of 2500 monomers originating from these families show high-sequence dynamics, evident by extensive exchanges between arrays on non-homologous chromosomes. In addition, our analysis shows that concerted evolution acts more efficiently on longer than on shorter arrays. Efficient genome-wide distribution of nine TR families implies the role of transposition only in expansion of the most dispersed family, and involvement of other mechanisms is anticipated. Despite similarities in sequence features, FISH experiments indicate high-level compartmentalization of centromeric and euchromatic tandem repeats.

High evolutionary turnover of satellite families in Caenorhabditis

Background

The high density of tandem repeat sequences (satellites) in nematode genomes and the availability of genome sequences from several species in the group offer a unique opportunity to better understand the evolutionary dynamics and the functional role of these sequences. We take advantage of the previously developed SATFIND program to study the satellites in four Caenorhabditis species and investigate these questions.

Methods

The identification and comparison of satellites is carried out in three steps. First we find all the satellites present in each species with the SATFIND program. Each satellite is defined by its length, number of repeats, and repeat sequence. Only satellites with at least ten repeats are considered. In the second step we build satellite families with a newly developed alignment program. Satellite families are defined by a consensus sequence and the number of satellites in the family. Finally we compare the consensus sequence of satellite families in different species.

Results

We give a catalog of individual satellites in each species. We have also identified satellite families with a related sequence and compare them in different species. We analyze the turnover of satellites: they increased in size through duplications of fragments of 100-300 bases. It appears that in many cases they have undergone an explosive expansion. In C. elegans we have identified a subset of large satellites that have strong affinity for the centromere protein CENP-A. We have also compared our results with those obtained from other species, including one nematode and three mammals.

Conclusions

Most satellite families found in Caenorhabditis are species-specific; in particular those with long repeats. A subset of these satellites may facilitate the formation of kinetochores in mitosis. Other satellite families in C. elegans are either related to Helitron transposons or to meiotic pairing centers.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0495-x) contains supplementary material, which is available to authorized users.

Satellite DNA in Plants: More than Just Rubbish

For decades, satellite DNAs have been the hidden part of genomes. Initially considered as junk DNA, there is currently an increasing appreciation of the functional significance of satellite DNA repeats and of their sequences. Satellite DNA families accumulate in the heterochromatin in different parts of the eukaryotic chromosomes, mainly in pericentromeric and subtelomeric regions, but they also span the functional centromere. Tandem repeat sequences may spread from subtelomeric to interstitial loci, leading to the formation of chromosome-specific loci or to the accumulation in equilocal sites in different chromosomes. They also appear as the main components of the heterochromatin in the sex-specific region of sex chromosomes. Satellite DNA, required for chromosome organization, also plays a role in pairing and segregation. Some satellite repeats are transcribed and can participate in the formation and maintenance of heterochromatin structure and in the modulation of gene expression. In addition to the identification of the different satellite DNA families, their characteristics and location, we are interested in determining their impact on the genomes, by identifying the mechanisms leading to their appearance and amplification as well as in understanding how they change over time, the factors affecting these changes, and the influence exerted by the evolutionary history of the organisms. On the other hand, satellite DNA sequences are rapidly evolving sequences that may cause reproductive barriers between organisms and promote speciation. The accumulation of experimental data collected in recent years and the emergence of new approaches based on next-generation sequencing and high-throughput genome analysis are opening new perspectives that are changing our understanding of satellite DNA. This review examines recent data to provide a timely update on the overall information gathered about this part of the genome, focusing on the advances in the knowledge of its origin, its evolution, and its potential functional roles.