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 wide distribution and horizontal transfers of beta satellite DNA in eukaryotes.

Beta satellite DNA (satDNA), also known as Sau3A sequences, are repetitive DNA sequences reported in human and primate genomes. It is previously thought that beta satDNAs originated in old world monkeys and bursted in great apes. In this study, we searched 7821 genome assemblies of 3767 eukaryotic species and found that beta satDNAs are widely distributed across eukaryotes. The four major branches of eukaryotes, animals, fungi, plants and Harosa/SAR, all have multiple clades containing beta satDNAs. These results were also confirmed by searching whole genome sequencing data (SRA) and PCR assay. Beta satDNA sequences were found in all the primate clades, as well as in Dermoptera and Scandentia, indicating that the beta satDNAs in primates might originate in the common ancestor of Primatomorpha or Euarchonta. In contrast, the widely patchy distribution of beta satDNAs across eukaryotes presents a typical scenario of multiple horizontal transfers.

Characterization of centromeric satellite DNAs (MALREP) in the Asian swamp eel (Monopterus albus) suggests the possible origin of repeats from transposable elements.

Centromeric satellite DNA (cen-satDNA) sequences of the Asian swamp eel (Monopterus albus) were characterized. Three GC-rich cen-satDNA sequences were detected as a 233 bp MALREP-A and a 293 bp MALREP-B localized to all chromosomes, and a 293 bp MALREP-C distributed on eight chromosome pairs. Sequence lengths of MALREP-B and MALREP-C were 60 bp larger than that of MALREP-A, showing partial homology with core sequences (233 bp). Size differences between MALREP-A and MALREP-B/C suggest the possible occurrence of two satDNA families. The presence of an additional 60 bp in MALREP-B/C resulted from an ancient dimer of 233 bp monomers and subsequent mutation and homogenization between the two monomers. All MALREPs showed partial homology with transposable elements (TEs), suggesting that the MALREPs originated from the TEs. The MALREPs might have been acquired in the Asian swamp eel, thereby promoting fixation in the species.

Extraordinary Sequence Diversity and Promiscuity of Centromeric Satellites in the Legume Tribe Fabeae.

Satellite repeats are major sequence constituents of centromeres in many plant and animal species. Within a species, a single family of satellite sequences typically occupies centromeres of all chromosomes and is absent from other parts of the genome. Due to their common origin, sequence similarities exist among the centromere-specific satellites in related species. Here, we report a remarkably different pattern of centromere evolution in the plant tribe Fabeae, which includes genera Pisum, Lathyrus, Vicia, and Lens. By immunoprecipitation of centromeric chromatin with CENH3 antibodies, we identified and characterized a large and diverse set of 64 families of centromeric satellites in 14 species. These families differed in their nucleotide sequence, monomer length (33-2,979 bp), and abundance in individual species. Most families were species-specific, and most species possessed multiple (2-12) satellites in their centromeres. Some of the repeats that were shared by several species exhibited promiscuous patterns of centromere association, being located within CENH3 chromatin in some species, but apart from the centromeres in others. Moreover, FISH experiments revealed that the same family could assume centromeric and noncentromeric positions even within a single species. Taken together, these findings suggest that Fabeae centromeres are not shaped by the coevolution of a single centromeric satellite with its interacting CENH3 proteins, as proposed by the centromere drive model. This conclusion is also supported by the absence of pervasive adaptive evolution of CENH3 sequences retrieved from Fabeae species.

Eight Million Years of Satellite DNA Evolution in Grasshoppers of the Genus Schistocerca Illuminate the Ins and Outs of the Library Hypothesis.

Satellite DNA (satDNA) is an abundant class of tandemly repeated noncoding sequences, showing high rate of change in sequence, abundance, and physical location. However, the mechanisms promoting these changes are still controversial. The library model was put forward to explain the conservation of some satDNAs for long periods, predicting that related species share a common collection of satDNAs, which mostly experience quantitative changes. Here, we tested the library model by analyzing three satDNAs in ten species of Schistocerca grasshoppers. This group represents a valuable material because it diversified during the last 7.9 Myr across the American continent from the African desert locust (Schistocerca gregaria), and this thus illuminates the direction of evolutionary changes. By combining bioinformatic and cytogenetic, we tested whether these three satDNA families found in S. gregaria are also present in nine American species, and whether differential gains and/or losses have occurred in the lineages. We found that the three satDNAs are present in all species but display remarkable interspecies differences in their abundance and sequences while being highly consistent with genus phylogeny. The number of chromosomal loci where satDNA is present was also consistent with phylogeny for two satDNA families but not for the other. Our results suggest eminently chance events for satDNA evolution. Several evolutionary trends clearly imply either massive amplifications or contractions, thus closely fitting the library model prediction that changes are mostly quantitative. Finally, we found that satDNA amplifications or contractions may influence the evolution of monomer consensus sequences and by chance playing a major role in driftlike dynamics.

The CENP-B box, a nucleotide motif involved in centromere formation, has multiple origins in New World monkeys.

Centromere protein B (CENP-B), a protein participating in centromere formation, binds to centromere satellite DNA by recognizing a 17-bp motif called the CENP-B box. This motif is found in hominids (humans and great apes) at an identical location in repeat units of their centromere satellite DNA. We have recently reported that the CENP-B box exists at diverse locations in three New World monkey species (marmoset, squirrel monkey and tamarin). However, the evolutionary origin of the CENP-B box in these species was not determined. It could have been present in a common ancestor, or emerged multiple times in different lineages. Here we present results of a phylogenetic analysis of centromere satellite DNA that support the multiple emergence hypothesis. Repeat units almost invariably formed monophyletic groups in each species and the CENP-B box location was unique for each species. The CENP-B box is not essential for the immediate survival of its host organism. On the other hand, it is known to be required for de novo centromere assembly. Our results suggest that the CENP-B box confers a long-term selective advantage. For example, it may play a pivotal role when a centromere is accidentally lost or impaired.

Expression dynamics of repetitive DNA in early human embryonic development.

BACKGROUND: The last decade witnessed a number of genome-wide studies on human pre-implantation, which mostly focused on genes and provided only limited information on repeats, excluding the satellites. Considering the fact that repeats constitute a large portion of our genome with reported links to human physiology and disease, a thorough understanding of their spatiotemporal regulation during human embryogenesis will give invaluable clues on chromatin dynamics across time and space. Therefore, we performed a detailed expression analysis of all repetitive DNA elements including the satellites across stages of human pre-implantation and embryonic stem cells. RESULTS: We uncovered stage-specific expressions of more than a thousand repeat elements whose expressions fluctuated with a mild global decrease at the blastocyst stage. Most satellites were highly expressed at the 4-cell level and expressions of ACRO1 and D20S16 specifically peaked at this point. Whereas all members of the SVA elements were highly upregulated at 8-cell and morula stages, other transposons and small RNA repeats exhibited a high level of variation among their specific subtypes. Our repeat enrichment analysis in gene promoters coupled with expression correlations highlighted potential links between repeat expressions and nearby genes, emphasising mostly 8-cell and morula specific genes together with SVA_D, LTR5_Hs and LTR70 transposons. The DNA methylation analysis further complemented the understanding on the mechanistic aspects of the repeatome's regulation per se and revealed critical stages where DNA methylation levels are negatively correlating with repeat expression. CONCLUSIONS: Taken together, our study shows that specific expression patterns are not exclusive to genes and long non-coding RNAs but the repeatome also exhibits an intriguingly dynamic pattern at the global scale. Repeats identified in this study; particularly satellites, which were historically associated with heterochromatin, and those with potential links to nearby gene expression provide valuable insights into the understanding of key events in genomic regulation and warrant further research in epigenetics, genomics and developmental biology.

Genome invasion by a hypomethylated satellite repeat in Australian crucifer Ballantinia antipoda.

Repetitive sequences are ubiquitous components of all eukaryotic genomes. They contribute to genome evolution and the regulation of gene transcription. However, the uncontrolled activity of repetitive sequences can negatively affect genome functions and stability. Therefore, repetitive DNAs are embedded in a highly repressive heterochromatic environment in plant cell nuclei. Here, we analyzed the sequence, composition and the epigenetic makeup of peculiar non-pericentromeric heterochromatic segments in the genome of the Australian crucifer Ballantinia antipoda. By the combination of high throughput sequencing, graph-based clustering and cytogenetics, we found that the heterochromatic segments consist of a mixture of unique sequences and an A-T-rich 174 bp satellite repeat (BaSAT1). BaSAT1 occupies about 10% of the B. antipoda nuclear genome in >250 000 copies. Unlike many other highly repetitive sequences, BaSAT1 repeats are hypomethylated; this contrasts with the normal patterns of DNA methylation in the B. antipoda genome. Detailed analysis of several copies revealed that these non-methylated BaSAT1 repeats were also devoid of heterochromatic histone H3K9me2 methylation. However, the factors decisive for the methylation status of BaSAT1 repeats remain currently unknown. In summary, we show that even highly repetitive sequences can exist as hypomethylated in the plant nuclear genome.

Unexpected conformational variations of the human centromeric chromatin complex.

We combined classical salt fractionation with chromatin immunoprecipitation to recover human centromeric chromatin under native conditions. We found that >85% of the total centromeric chromatin is insoluble under conditions typically used for native chromatin extraction. To map both soluble and insoluble chromatin in situ, we combined CUT&RUN (cleavage under targets and release using nuclease), a targeted nuclease method, with salt fractionation. Using this approach, we observed unexpected structural and conformational variations of centromere protein A (CENP-A)-containing complexes on different α-satellite dimeric units within highly homogenous arrays. Our results suggest that slight α-satellite sequence differences control the structure and occupancy of the associated centromeric chromatin complex.

DNA-Assembled Core-Satellite Upconverting-Metal-Organic Framework Nanoparticle Superstructures for Efficient Photodynamic Therapy.

DNA-mediated assembly of core-satellite structures composed of Zr(IV)-based porphyrinic metal-organic framework (MOF) and NaYF4 ,Yb,Er upconverting nanoparticles (UCNPs) for photodynamic therapy (PDT) is reported. MOF NPs generate singlet oxygen (1 O2 ) upon photoirradiation with visible light without the need for additional small molecule, diffusional photosensitizers such as porphyrins. Using DNA as a templating agent, well-defined MOF-UCNP clusters are produced where UCNPs are spatially organized around a centrally located MOF NP. Under NIR irradiation, visible light emitted from the UCNPs is absorbed by the core MOF NP to produce 1 O2 at significantly greater amounts than what can be produced from simply mixing UCNPs and MOF NPs. The MOF-UCNP core-satellite superstructures also induce strong cell cytotoxicity against cancer cells, which are further enhanced by attaching epidermal growth factor receptor targeting affibodies to the PDT clusters, highlighting their promise as theranostic photodynamic agents.

Single-molecule sequencing resolves the detailed structure of complex satellite DNA loci in Drosophila melanogaster.

Highly repetitive satellite DNA (satDNA) repeats are found in most eukaryotic genomes. SatDNAs are rapidly evolving and have roles in genome stability and chromosome segregation. Their repetitive nature poses a challenge for genome assembly and makes progress on the detailed study of satDNA structure difficult. Here, we use single-molecule sequencing long reads from Pacific Biosciences (PacBio) to determine the detailed structure of all major autosomal complex satDNA loci in Drosophila melanogaster, with a particular focus on the 260-bp and Responder satellites. We determine the optimal de novo assembly methods and parameter combinations required to produce a high-quality assembly of these previously unassembled satDNA loci and validate this assembly using molecular and computational approaches. We determined that the computationally intensive PBcR-BLASR assembly pipeline yielded better assemblies than the faster and more efficient pipelines based on the MHAP hashing algorithm, and it is essential to validate assemblies of repetitive loci. The assemblies reveal that satDNA repeats are organized into large arrays interrupted by transposable elements. The repeats in the center of the array tend to be homogenized in sequence, suggesting that gene conversion and unequal crossovers lead to repeat homogenization through concerted evolution, although the degree of unequal crossing over may differ among complex satellite loci. We find evidence for higher-order structure within satDNA arrays that suggest recent structural rearrangements. These assemblies provide a platform for the evolutionary and functional genomics of satDNAs in pericentric heterochromatin.

Reorganization of chromosome architecture in replicative cellular senescence.

Replicative cellular senescence is a fundamental biological process characterized by an irreversible arrest of proliferation. Senescent cells accumulate a variety of epigenetic changes, but the three-dimensional (3D) organization of their chromatin is not known. We applied a combination of whole-genome chromosome conformation capture (Hi-C), fluorescence in situ hybridization, and in silico modeling methods to characterize the 3D architecture of interphase chromosomes in proliferating, quiescent, and senescent cells. Although the overall organization of the chromatin into active (A) and repressive (B) compartments and topologically associated domains (TADs) is conserved between the three conditions, a subset of TADs switches between compartments. On a global level, the Hi-C interaction matrices of senescent cells are characterized by a relative loss of long-range and gain of short-range interactions within chromosomes. Direct measurements of distances between genetic loci, chromosome volumes, and chromatin accessibility suggest that the Hi-C interaction changes are caused by a significant reduction of the volumes occupied by individual chromosome arms. In contrast, centromeres oppose this overall compaction trend and increase in volume. The structural model arising from our study provides a unique high-resolution view of the complex chromosomal architecture in senescent cells.

Influence of DNA methylation on positioning and DNA flexibility of nucleosomes with pericentric satellite DNA.

DNA methylation occurs on CpG sites and is important to form pericentric heterochromatin domains. The satellite 2 sequence, containing seven CpG sites, is located in the pericentric region of human chromosome 1 and is highly methylated in normal cells. In contrast, the satellite 2 region is reportedly hypomethylated in cancer cells, suggesting that the methylation status may affect the chromatin structure around the pericentric regions in tumours. In this study, we mapped the nucleosome positioning on the satellite 2 sequence in vitro and found that DNA methylation modestly affects the distribution of the nucleosome positioning. The micrococcal nuclease assay revealed that the DNA end flexibility of the nucleosomes changes, depending on the DNA methylation status. However, the structures and thermal stabilities of the nucleosomes are unaffected by DNA methylation. These findings provide new information to understand how DNA methylation functions in regulating pericentric heterochromatin formation and maintenance in normal and malignant cells.

Molecular cloning and characterization of satellite DNA sequences from constitutive heterochromatin of the habu snake (Protobothrops flavoviridis, Viperidae) and the Burmese python (Python bivittatus, Pythonidae).

Highly repetitive DNA sequences of the centromeric heterochromatin provide valuable molecular cytogenetic markers for the investigation of genomic compartmentalization in the macrochromosomes and microchromosomes of sauropsids. Here, the relationship between centromeric heterochromatin and karyotype evolution was examined using cloned repetitive DNA sequences from two snake species, the habu snake (Protobothrops flavoviridis, Crotalinae, Viperidae) and Burmese python (Python bivittatus, Pythonidae). Three satellite DNA (stDNA) families were isolated from the heterochromatin of these snakes: 168-bp PFL-MspI from P. flavoviridis and 196-bp PBI-DdeI and 174-bp PBI-MspI from P. bivittatus. The PFL-MspI and PBI-DdeI sequences were localized to the centromeric regions of most chromosomes in the respective species, suggesting that the two sequences were the major components of the centromeric heterochromatin in these organisms. The PBI-MspI sequence was localized to the pericentromeric region of four chromosome pairs. The PFL-MspI and the PBI-DdeI sequences were conserved only in the genome of closely related species, Gloydius blomhoffii (Crotalinae) and Python molurus, respectively, although their locations on the chromosomes were slightly different. In contrast, the PBI-MspI sequence was also in the genomes of P. molurus and Boa constrictor (Boidae), and additionally localized to the centromeric regions of eight chromosome pairs in B. constrictor, suggesting that this sequence originated in the genome of a common ancestor of Pythonidae and Boidae, approximately 86 million years ago. The three stDNA sequences showed no genomic compartmentalization between the macrochromosomes and microchromosomes, suggesting that homogenization of the centromeric and/or pericentromeric stDNA sequences occurred in the macrochromosomes and microchromosomes of these snakes.

Molecular characterization of Chilli leaf curl viruses infecting new host plant Petunia hybrida in India.

Petunia hybrida is an important ornamental plant grown in many countries including India. It is a good model plant for the study of genetics and molecular biology. During a survey in 2013-2014, severe leaf curling was observed on most of the P. hybrida grown in the Sikar district, Rajasthan. The infected plants were analyzed for begomovirus infection by rolling circular amplification (RCA) and sequenced. Full length sequences confirmed the association of monopartite begomovirus with betasatellites. Phylogenetic analysis showed the highest percentage of identity with Chilli leaf curl virus (ChLCuV) and therefore considered to be an isolate of ChLCuV. Recombination analysis showed that ChLCuV has broadened its host range by recombination process. To the best our knowledge, this is the first report of natural occurrence of ChLCuV on P. hybrida in India.

Altered imprinted gene expression and methylation patterns in mid-gestation aborted cloned porcine fetuses and placentas.

PURPOSE: To determine the expression patterns of imprinted genes and their methylation status in aborted cloned porcine fetuses and placentas. METHODS: RNA and DNA were prepared from fetuses and placentas that were produced by SCNT and controls from artificial insemination. The expression of 18 imprinted genes was determined by quantitative real-time PCR (q-PCR). Bisulfite sequencing PCR (BSP) was conducted to determine the methylation status of PRE-1 short interspersed repetitive element (SINE), satellite DNA and H19 differentially methylated region 3 (DMR3). RESULTS: The weight, imprinted gene expression and genome-wide DNA methylation patterns were compared between the mid-gestation aborted and normal control samples. The results showed hypermethylation of PRE-1 and satellite sequences, the aberrant expression of imprinted genes, and the hypomethylation of H19 DMR3 occurred in mid-gestation aborted fetuses and placentas. CONCLUSIONS: Cloned pigs generated by somatic cell nuclear transfer (SCNT) showed a greater ratio of early abortion during mid-gestation than did normal controls because of the incomplete epigenetic reprogramming of the donor cells. Altered expression of imprinted genes and the hypermethylation profile of the repetitive regions (PRE-1 and satellite DNA) may be associated with defective development and early abortion of cloned pigs, emphasizing the importance of epigenetics during pregnancy and implications thereof for patient-specific embryonic stem cells for human therapeutic cloning and improvement of human assisted reproduction.

Estimation of the RNU2 macrosatellite mutation rate by BRCA1 mutation tracing.

Large tandem repeat sequences have been poorly investigated as severe technical limitations and their frequent absence from the genome reference hinder their analysis. Extensive allelotyping of this class of variation has not been possible until now and their mutational dynamics are still poorly known. In order to estimate the mutation rate of a macrosatellite, we analysed in detail the RNU2 locus, which displays at least 50 different alleles containing 5-82 copies of a 6.1 kb repeat unit. Mining data from the 1000 Genomes Project allowed us to precisely estimate copy numbers of the RNU2 repeat unit using read depth of coverage. This further revealed significantly different mean values in various recent modern human populations, favoring a scenario of fast evolution of this locus. Its proximity to a disease gene with numerous founder mutations, BRCA1, within the same linkage disequilibrium block, offered the unique opportunity to trace RNU2 arrays over a large timescale. Analysis of the transmission of RNU2 arrays associated with one 'private' mutation in an extended kindred and four founder mutations in multiple kindreds gave an estimation by maximum likelihood of 5 × 10(-3) mutations per generation, which is close to that of microsatellites.

Circomics of Cuban geminiviruses reveals the first alpha-satellite DNA in the Caribbean.

Circomics (circular DNA genomics), the combination of rolling circle amplification (RCA), restriction fragment length polymorphism (RFLP) analysis and pyro-sequencing, has been used recently to identify geminiviruses with high efficiency and low costs. Circular DNAs associated with Cuban geminiviruses were characterised by RCA/RFLP analysis and 454 sequencing of two batches of DNA amplified from selected plant samples as well as individual cloning and Sanger sequencing of DNA components and compared to other geminiviral DNAs by phylogenetic analysis. Cuban geminiviruses that were closely related to each other challenged the circomics approach. Ten geminiviral components and one alpha-satellite DNA were determined and compared to three geminiviral components obtained by conventional cloning. New strains of Sida yellow mottle virus (SiYMoV), tomato yellow distortion leaf virus (ToYDLV), Sida golden mosaic Florida virus (SiGMFV) and Sida golden mosaic Liguanea virus (SiGMLV) are described with host plant species being classified by molecular PCR-based bar coding. A new virus species is named Peristrophe mosaic virus. The first alpha-satellite found in Middle America establishes the New World branch of these elements which are related to nanoviruses and were previously thought to be restricted to the Old World. In conclusion, circomics is efficient for complex infections and closely related viruses to detected unexpected viral DNAs, but may need some scrutinisation by direct sequencing and cloning of individual components for certain cases.

Regional changes in the sequence of cotton leaf curl multan betasatellite.

Cotton leaf curl disease (CLCuD) in Pakistan and northwestern India is caused by monopartite begomoviruses in association with an essential, disease-specific satellite, Cotton leaf curl Multan betasatellite (CLCuMB). Following a recent upsurge in CLCuD problems in Sindh province (southern Pakistan), sequences of clones of CLCuMB were obtained from Sindh and Punjab province (central Pakistan), where CLCuD has been a problem since the mid-1980s. The sequences were compared to all sequences of CLCuMB available in the databases. Analysis of the sequences shows extensive sequence variation in CLCuMB, most likely resulting from recombination. The range of sequence variants differ between Sindh, the Punjab and northwestern India. The possible significance of the findings with respect to movement of the CLCuD between the three regions is discussed. Additionally, the lack of sequence variation within the only coding sequence of CLCuMB suggests that the betasatellite is not involved in resistance breaking which became a problem after 2001 in the Punjab and subsequently also in northwestern India.

Viral metagenomics: analysis of begomoviruses by illumina high-throughput sequencing.

Traditional DNA sequencing methods are inefficient, lack the ability to discern the least abundant viral sequences, and ineffective for determining the extent of variability in viral populations. Here, populations of single-stranded DNA plant begomoviral genomes and their associated beta- and alpha-satellite molecules (virus-satellite complexes) (genus, Begomovirus; family, Geminiviridae) were enriched from total nucleic acids isolated from symptomatic, field-infected plants, using rolling circle amplification (RCA). Enriched virus-satellite complexes were subjected to Illumina-Next Generation Sequencing (NGS). CASAVA and SeqMan NGen programs were implemented, respectively, for quality control and for de novo and reference-guided contig assembly of viral-satellite sequences. The authenticity of the begomoviral sequences, and the reproducibility of the Illumina-NGS approach for begomoviral deep sequencing projects, were validated by comparing NGS results with those obtained using traditional molecular cloning and Sanger sequencing of viral components and satellite DNAs, also enriched by RCA or amplified by polymerase chain reaction. As the use of NGS approaches, together with advances in software development, make possible deep sequence coverage at a lower cost; the approach described herein will streamline the exploration of begomovirus diversity and population structure from naturally infected plants, irrespective of viral abundance. This is the first report of the implementation of Illumina-NGS to explore the diversity and identify begomoviral-satellite SNPs directly from plants naturally-infected with begomoviruses under field conditions.