Regular Higher Order Repeat Structures in Beetle Tribolium castaneum Genome

Precise identification of cascading alpha satellite higher order repeats in T2T-CHM13 assembly of human chromosome 3

AIM

To precisely identify and analyze alpha-satellite higher-order repeats (HORs) in T2T-CHM13 assembly of human chromosome 3.

METHODS

From the recently sequenced complete T2T-CHM13 assembly of human chromosome 3, the precise alpha satellite HOR structure was computed by using the novel high-precision GRM2023 algorithm with global repeat map (GRM) and monomer distance (MD) diagrams.

RESULTS

The major alpha satellite HOR array in chromosome 3 revealed a novel cascading HOR, housing 17mer HOR copies with subfragments of periods 15 and 2. Within each row in the cascading HOR, the monomers were of different types, but different rows within the same cascading 17mer HOR contained more than one monomer of the same type. Each canonical 17mer HOR copy comprised 17 monomers belonging to 16 different monomer types. Another pronounced 10mer HOR array was of the regular Willard's type.

CONCLUSION

Our findings emphasize the complexity within the chromosome 3 centromere as well as deviations from expected highly regular patterns.

Tandem NBPF 3mer HORs (Olduvai triplets) in Neanderthal and two novel HOR tandem arrays in human chromosome 1 T2T-CHM13 assembly

It is known that the ~ 1.6 kb Neuroblastoma BreakPoint Family (NBPF) repeats are human specific and contributing to cognitive capabilities, with increasing frequency in higher order repeat 3mer HORs (Olduvai triplets). From chimpanzee to modern human there is a discontinuous jump from 0 to ~ 50 tandemly organized 3mer HORs. Here we investigate the structure of NBPF 3mer HORs in the Neanderthal genome assembly of Pääbo et al., comparing it to the results obtained for human hg38.p14 chromosome 1. Our findings reveal corresponding NBPF 3mer HOR arrays in Neanderthals with slightly different monomer structures and numbers of HOR copies compared to humans. Additionally, we compute the NBPF 3mer HOR pattern for the complete telomere-to-telomere human genome assembly (T2T-CHM13) by Miga et al., identifying two novel tandem arrays of NBPF 3mer HOR repeats with 5 and 9 NBPF 3mer HOR copies. We hypothesize that these arrays correspond to novel NBPF genes (here referred to as NBPFA1 and NBPFA2). Further improving the quality of the Neanderthal genome using T2T-CHM13 as a reference would be of great interest in determining the presence of such distant novel NBPF genes in the Neanderthal genome and enhancing our understanding of human evolution.

Tandemly repeated NBPF HOR copies (Olduvai triplets): Possible impact on human brain evolution

Previously it was found that the neuroblastoma breakpoint family (NBPF) gene repeat units of ∼1.6 kb have an important role in human brain evolution and function. The higher order organization of these repeat units has been discovered by both methods, the higher order repeat (HOR)-searching method and the HLS searching method. Using the HOR searching method with global repeat map algorithm, here we identified the tandemly organized NBPF HORs in the human and nonhuman primate NCBI reference genomes. We identified 50 tandemly organized canonical 3mer NBPF HOR copies (Olduvai triplets), but none in nonhuman primates chimpanzee, gorilla, orangutan, and Rhesus macaque. This discontinuous jump in tandemly organized HOR copy number is in sharp contrast to the known gradual increase in the number of Olduvai domains (NBPF monomers) from nonhuman primates to human, especially from ∼138 in chimpanzee to ∼300 in human genome. Using the same global repeat map algorithm method we have also determined the 3mer tandems of canonical 3mer HOR copies in 20 randomly chosen human genomes (10 male and 10 female). In all cases, we found the same 3mer HOR copy numbers as in the case of the reference human genome, with no mutation. On the other hand, some point mutations with respect to reference genome are found for some NBPF monomers which are not tandemly organized in canonical HORs.

Comparative Analysis of Transposable Elements in Genus Calliptamus Grasshoppers Revealed That Satellite DNA Contributes to Genome Size Variation

Simple Summary

Calliptamus is a genus of grasshoppers belonging to the family Acrididae. The genus Calliptamus includes approximately 17 recognized species. Calliptamus abbreviatus, Calliptamus italicus, and Calliptamus barbarus are three species that are widely found in northern China. These species are polyphagous, feeding on a variety of wild plants as well as crops, particularly legumes. The genome sizes, phylogenetic position, and transcriptome analysis of the genus Calliptamus were already known previous to this research. The repeatome analysis of these species was missing, which is directly linked to the larger genome sizes of the grasshoppers. Here, we classified repetitive DNA sequences at the level of superfamilies and sub-families, and found that LINE, TcMar-Tc1 and Ty3-gypsy LTR retrotransposons dominated the repeatomes of all genomes, accounting for 16–34% of the total genomes of these species. Satellite DNA dynamic evolutionary changes in all three genomes played a role in genome size evolution. This study would be a valuable source for future genome assemblies.

Abstract

Transposable elements (TEs) play a significant role in both eukaryotes and prokaryotes genome size evolution, structural changes, duplication, and functional variabilities. However, the large number of different repetitive DNA has hindered the process of assembling reference genomes, and the genus level TEs diversification of the grasshopper massive genomes is still under investigation. The genus Calliptamus diverged from Peripolus around 17 mya and its species divergence dated back about 8.5 mya, but their genome size shows rather large differences. Here, we used low-coverage Illumina unassembled short reads to investigate the effects of evolutionary dynamics of satDNAs and TEs on genome size variations. The Repeatexplorer2 analysis with 0.5X data resulted in 52%, 56%, and 55% as repetitive elements in the genomes of Calliptamus barbarus, Calliptamus italicus, and Calliptamus abbreviatus, respectively. The LINE and Ty3-gypsy LTR retrotransposons and TcMar-Tc1 dominated the repeatomes of all genomes, accounting for 16–35% of the total genomes of these species. Comparative analysis unveiled that most of the transposable elements (TEs) except satDNAs were highly conserved across three genomes in the genus Calliptamus grasshoppers. Out of a total of 20 satDNA families, 17 satDNA families were commonly shared with minor variations in abundance and divergence between three genomes, and 3 were Calliptamus barbarus specific. Our findings suggest that there is a significant amplification or contraction of satDNAs at genus phylogeny which is the main cause that made genome size different.

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.

First Description of a Satellite DNA in Manatees' Centromeric Regions

Trichechus manatus and Trichechus inunguis are the two Sirenia species that occur in the Americas. Despite their increasing extinction risk, many aspects of their biology remain understudied, including the repetitive DNA fraction of their genomes. Here we used the sequenced genome of T. manatus and TAREAN to identify satellite DNAs (satDNAs) in this species. We report the first description of TMAsat, a satDNA comprising ~0.87% of the genome, with ~684bp monomers and centromeric localization. In T. inunguis, TMAsat showed similar monomer length, chromosome localization and conserved CENP-B box-like motifs as in T. manatus. We also detected this satDNA in the Dugong dugon and in the now extinct Hydrodamalis gigas genomes. The neighbor-joining tree shows that TMAsat sequences from T. manatus, T. inunguis, D. dugon, and H. gigas lack species-specific clusters, which disagrees with the predictions of concerted evolution. We detected a divergent TMAsat-like homologous sequence in elephants and hyraxes, but not in other mammals, suggesting this sequence was already present in the common ancestor of Paenungulata, and later became a satDNA in the Sirenians. This is the first description of a centromeric satDNA in manatees and will facilitate the inclusion of Sirenia in future studies of centromeres and satDNA biology.

Genomic Tackling of Human Satellite DNA: Breaking Barriers through Time

(Peri)centromeric repetitive sequences and, more specifically, satellite DNA (satDNA) sequences, constitute a major human genomic component. SatDNA sequences can vary on a large number of features, including nucleotide composition, complexity, and abundance. Several satDNA families have been identified and characterized in the human genome through time, albeit at different speeds. Human satDNA families present a high degree of sub-variability, leading to the definition of various subfamilies with different organization and clustered localization. Evolution of satDNA analysis has enabled the progressive characterization of satDNA features. Despite recent advances in the sequencing of centromeric arrays, comprehensive genomic studies to assess their variability are still required to provide accurate and proportional representation of satDNA (peri)centromeric/acrocentric short arm sequences. Approaches combining multiple techniques have been successfully applied and seem to be the path to follow for generating integrated knowledge in the promising field of human satDNA biology.

The Pontastacus leptodactylus (Astacidae) Repeatome Provides Insight Into Genome Evolution and Reveals Remarkable Diversity of Satellite DNA

Pontastacus leptodactylus is a native European crayfish species found in both freshwater and brackish environments. It has commercial importance for fisheries and aquaculture industries. Up till now, most studies concerning P. leptodactylus have focused onto gaining knowledge about its phylogeny and population genetics. However, little is known about the chromosomal evolution and genome organization of this species. Therefore, we performed clustering analysis of a low coverage genomic dataset to identify and characterize repetitive DNA in the P. leptodactylus genome. In addition, the karyogram of P. leptodactylus (2n = 180) is presented here for the first time consisting of 75 metacentric, 14 submetacentric, and a submetacentric/metacentric heteromorphic chromosome pair. We determined the genome size to be at ~18.7 gigabase pairs. Repetitive DNA represents about 54.85% of the genome. Satellite DNA repeats are the most abundant type of repetitive DNA, making up to ~28% of the total amount of repetitive elements, followed by the Ty3/Gypsy retroelements (~15%). Our study established a surprisingly high diversity of satellite repeats in P. leptodactylus. The genome of P. leptodactylus is by far the most satellite-rich genome discovered to date with 258 satellite families described. Of the five mapped satellite DNA families on chromosomes, PlSAT3-411 co-localizes with the AT-rich DAPI positive probable (peri)centromeric heterochromatin on all chromosomes, while PlSAT14-79 co-localizes with the AT-rich DAPI positive (peri)centromeric heterochromatin on one chromosome and is also located subterminally and intercalary on some chromosomes. PlSAT1-21 is located intercalary in the vicinity of the (peri)centromeric heterochromatin on some chromosomes, while PlSAT6-70 and PlSAT7-134 are located intercalary on some P. leptodactylus chromosomes. The FISH results reveal amplification of interstitial telomeric repeats (ITRs) in P. leptodactylus. The prevalence of repetitive elements, especially the satellite DNA repeats, may have provided a driving force for the evolution of the P. leptodactylus genome.

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.

CenH3 distribution reveals extended centromeres in the model beetle Tribolium castaneum

Centromeres are chromosomal domains essential for kinetochore assembly and correct chromosome segregation. Inconsistent in their underlying DNA sequences, centromeres are defined epigenetically by the presence of the centromere-specific histone H3 variant CenH3. Most of the analyzed eukaryotes have monocentric chromosomes in which CenH3 proteins deposit into a single, primary constriction visible at metaphase chromosomes. Contrary to monocentrics, evolutionary sporadic holocentric chromosomes lack a primary constriction and have kinetochore activity distributed along the entire chromosome length. In this work, we identified cCENH3 protein, the centromeric H3 histone of the coleopteran model beetle Tribolium castaneum. By ChIP-seq analysis we disclosed that cCENH3 chromatin assembles upon a repertoire of repetitive DNAs. cCENH3 in situ mapping revealed unusually elongated T. castaneum centromeres that comprise approximately 40% of the chromosome length. Being the longest insect regional centromeres evidenced so far, T. castaneum centromeres are characterized by metapolycentric structure composed of several individual cCENH3-containing domains. We suggest that the model beetle T. castaneum with its metapolycentromeres could represent an excellent model for further studies of non-canonical centromeres in insects.

Decoding the Role of Satellite DNA in Genome Architecture and Plasticity-An Evolutionary and Clinical Affair

Repetitive DNA is a major organizational component of eukaryotic genomes, being intrinsically related with their architecture and evolution. Tandemly repeated satellite DNAs (satDNAs) can be found clustered in specific heterochromatin-rich chromosomal regions, building vital structures like functional centromeres and also dispersed within euchromatin. Interestingly, despite their association to critical chromosomal structures, satDNAs are widely variable among species due to their high turnover rates. This dynamic behavior has been associated with genome plasticity and chromosome rearrangements, leading to the reshaping of genomes. Here we present the current knowledge regarding satDNAs in the light of new genomic technologies, and the challenges in the study of these sequences. Furthermore, we discuss how these sequences, together with other repeats, influence genome architecture, impacting its evolution and association with disease.

Discovery of 33mer in chromosome 21 - the largest alpha satellite higher order repeat unit among all human somatic chromosomes

The centromere is important for segregation of chromosomes during cell division in eukaryotes. Its destabilization results in chromosomal missegregation, aneuploidy, hallmarks of cancers and birth defects. In primate genomes centromeres contain tandem repeats of ~171 bp alpha satellite DNA, commonly organized into higher order repeats (HORs). In spite of crucial importance, satellites have been understudied because of gaps in sequencing - genomic “black holes”. Bioinformatical studies of genomic sequences open possibilities to revolutionize understanding of repetitive DNA datasets. Here, using robust (Global Repeat Map) algorithm we identified in hg38 sequence of human chromosome 21 complete ensemble of alpha satellite HORs with six long repeat units (≥20 mers), five of them novel. Novel 33mer HOR has the longest HOR unit identified so far among all somatic chromosomes and novel 23mer reverse HOR is distant far from the centromere. Also, we discovered that for hg38 assembly the 33mer sequences in chromosomes 21, 13, 14, and 22 are 100% identical but nearby gaps are present; that seems to require an additional more precise sequencing. Chromosome 21 is of significant interest for deciphering the molecular base of Down syndrome and of aneuploidies in general. Since the chromosome identifier probes are largely based on the detection of higher order alpha satellite repeats, distinctions between alpha satellite HORs in chromosomes 21 and 13 here identified might lead to a unique chromosome 21 probe in molecular cytogenetics, which would find utility in diagnostics. It is expected that its complete sequence analysis will have profound implications for understanding pathogenesis of diseases and development of new therapeutic approaches.

Satellite DNA evolution: old ideas, new approaches

A substantial portion of the genomes of most multicellular eukaryotes consists of large arrays of tandemly repeated sequence, collectively called satellite DNA. The processes generating and maintaining different satellite DNA abundances across lineages are important to understand as satellites have been linked to chromosome mis-segregation, disease phenotypes, and reproductive isolation between species. While much theory has been developed to describe satellite evolution, empirical tests of these models have fallen short because of the challenges in assessing satellite repeat regions of the genome. Advances in computational tools and sequencing technologies now enable identification and quantification of satellite sequences genome-wide. Here, we describe some of these tools and how their applications are furthering our knowledge of satellite evolution and function.

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.

Satellite DNA: An Evolving Topic

Satellite DNA represents one of the most fascinating parts of the repetitive fraction of the eukaryotic genome. Since the discovery of highly repetitive tandem DNA in the 1960s, a lot of literature has extensively covered various topics related to the structure, organization, function, and evolution of such sequences. Today, with the advent of genomic tools, the study of satellite DNA has regained a great interest. Thus, Next-Generation Sequencing (NGS), together with high-throughput in silico analysis of the information contained in NGS reads, has revolutionized the analysis of the repetitive fraction of the eukaryotic genomes. The whole of the historical and current approaches to the topic gives us a broad view of the function and evolution of satellite DNA and its role in chromosomal evolution. Currently, we have extensive information on the molecular, chromosomal, biological, and population factors that affect the evolutionary fate of satellite DNA, knowledge that gives rise to a series of hypotheses that get on well with each other about the origin, spreading, and evolution of satellite DNA. In this paper, I review these hypotheses from a methodological, conceptual, and historical perspective and frame them in the context of chromosomal organization and evolution.

Light Sheet-based Fluorescence Microscopy of Living or Fixed and Stained Tribolium castaneum Embryos

The red flour beetle Tribolium castaneum has become an important insect model organism in developmental genetics and evolutionary developmental biology. The observation of Tribolium embryos with light sheet-based fluorescence microscopy has multiple advantages over conventional widefield and confocal fluorescence microscopy. Due to the unique properties of a light sheet-based microscope, three dimensional images of living specimens can be recorded with high signal-to-noise ratios and significantly reduced photo-bleaching as well as photo-toxicity along multiple directions over periods that last several days. With more than four years of methodological development and a continuous increase of data, the time seems appropriate to establish standard operating procedures for the usage of light sheet technology in the Tribolium community as well as in the insect community at large. This protocol describes three mounting techniques suitable for different purposes, presents two novel custom-made transgenic Tribolium lines appropriate for long-term live imaging, suggests five fluorescent dyes to label intracellular structures of fixed embryos and provides information on data post-processing for the timely evaluation of the recorded data. Representative results concentrate on long-term live imaging, optical sectioning and the observation of the same embryo along multiple directions. The respective datasets are provided as a downloadable resource. Finally, the protocol discusses quality controls for live imaging assays, current limitations and the applicability of the outlined procedures to other insect species. This protocol is primarily intended for developmental biologists who seek imaging solutions that outperform standard laboratory equipment. It promotes the continuous attempt to close the gap between the technically orientated laboratories/communities, which develop and refine microscopy methodologically, and the life science laboratories/communities, which require 'plug-and-play' solutions to technical challenges. Furthermore, it supports an axiomatic approach that moves the biological questions into the center of attention.