DNA Sequence Composition Research Articles

Epigenetics in rare neurological diseases.

Rare neurological diseases include a vast group of heterogenous syndromes with primary impairment(s) in the peripheral and/or central nervous systems. Such rare disorders may have overlapping phenotypes, despite their distinct genetic etiology. One unique aspect of rare neurological diseases is their potential common association with altered epigenetic mechanisms. Epigenetic mechanisms include regulatory processes that control gene expression and cellular phenotype without changing the composition of the corresponding DNA sequences. Epigenetic factors include three types of proteins, the "readers, writers, and erasers" of DNA and DNA-bound proteins. Thus, epigenetic impairments of many neurological diseases may contribute to their pathology and manifested phenotypes. Here, we aim to provide a comprehensive review on the general etiology of selected rare neurological diseases, that include Rett Syndrome, Prader-Willi Syndrome, Rubinstein-Taybi Syndrome, Huntington's disease, and Angelman syndrome, with respect to their associated aberrant epigenetic mechanisms.

Discerning the Role of DNA Sequence, Shape, and Flexibility in Recognition by Drosophila Transcription Factors.

The precise spatial and temporal orchestration of gene expression is crucial for the ontogeny of an organism and is mainly governed by transcription factors (TFs). The mechanism of recognition of cognate sites amid millions of base pairs in the genome by TFs is still incompletely understood. In this study, we focus on DNA sequence composition, shape, and flexibility preferences of 28 quintessential TFs from Drosophila melanogaster that are critical to development and body patterning mechanisms. Our study finds that TFs exhibit distinct predilections for DNA shape, flexibility, and sequence compositions in the proximity of transcription factor binding sites (TFBSs). Notably, certain zinc finger proteins prefer GC-rich areas with less negative propeller twist, while homeodomains mainly seek AT-rich regions with a more negative propeller twist at their sites. Intriguingly, while numerous cofactors share similar binding site preferences and bind closer to each other in the genome, some cofactors that have different preferences bind farther apart. These findings shed light on TF DNA recognition and provide novel insights into possible cofactor binding and transcriptional regulation mechanisms.

Unveiling the distinctive traits of functional rye centromeres: minisatellites, retrotransposons, and R-loop formation.

Centromeres play a vital role in cellular division by facilitating kinetochore assembly and spindle attachments. Despite their conserved functionality, centromeric DNA sequences exhibit rapid evolution, presenting diverse sizes and compositions across species. The functional significance of rye centromeric DNA sequences, particularly in centromere identity, remains unclear. In this study, we comprehensively characterized the sequence composition and organization of rye centromeres. Our findings revealed that these centromeres are primarily composed of long terminal repeat retrotransposons (LTR-RTs) and interspersed minisatellites. We systematically classified LTR-RTs into five categories, highlighting the prevalence of younger CRS1, CRS2, and CRS3 of CRSs (centromeric retrotransposons of Secale cereale) were primarily located in the core centromeres and exhibited a higher association with CENH3 nucleosomes. The minisatellites, mainly derived from retrotransposons, along with CRSs, played a pivotal role in establishing functional centromeres in rye. Additionally, we observed the formation of R-loops at specific regions of CRS1, CRS2, and CRS3, with both rye pericentromeres and centromeres exhibiting enrichment in R-loops. Notably, these R-loops selectively formed at binding regions of the CENH3 nucleosome in rye centromeres, suggesting a potential role in mediating the precise loading of CENH3 to centromeres and contributing to centromere specification. Our work provides insights into the DNA sequence composition, distribution, and potential function of R-loops in rye centromeres. This knowledge contributes valuable information to understanding the genetics and epigenetics of rye centromeres, offering implications for the development of synthetic centromeres in future plant modifications and beyond.

Developing an eDNA approach for wetland biomonitoring: Insights on technical and conventional approaches

AbstractWetlands are ecologically and culturally significant ecosystems that are experiencing biodiversity declines globally. Biomonitoring techniques that use environmental DNA (eDNA) to detect and monitor biodiversity are well established in lake, riverine, and marine ecosystems. However, their use in wetlands requires further development due to the presence of sediments that block eDNA filters to limit water filtration, alongside a lack of standardized methodology. In this study, we examined eDNA dynamics to understand spatiotemporal biodiversity patterns in an Aotearoa New Zealand wetland and to optimize their application to wetland‐specific challenges. We sampled four sites across Opuatia Wetland at three time points during an austral spring. We conducted conventional taxonomic surveys, tested three different filter sizes (1.2 μm, 5 μm, and semi‐quantitative dacron filters), and assessed our ability to detect foreign DNA (from kea; Nestor notabilis) at different time points and distances post‐release. We found significant differences in DNA sequence composition across time and space, and when using different sized filters. eDNA data generally complemented (versus replaced) conventional survey and identification methods, with certain species only detected by one method or the other. Taxonomic resolution of conventional sampling and identification methods often exceeded that of eDNA. Foreign DNA was detectable 10 m from its release point for up to 1 week post‐release. Our results provide new considerations for future eDNA research in wetland environments, where rapid biomonitoring techniques are needed to support conservation and preservation.

Stable bulged G-quadruplexes in the human genome: identification, experimental validation and functionalization.

DNA sequence composition determines the topology and stability of G-quadruplexes (G4s). Bulged G-quadruplex structures (G4-Bs) are a subset of G4s characterized by 3D conformations with bulges. Current search algorithms fail to capture stable G4-B, making their genome-wide study infeasible. Here, we introduced a large family of computationally defined and experimentally verified potential G4-B forming sequences (pG4-BS). We found 478263 pG4-BS regions that do not overlap 'canonical' G4-forming sequences in the human genome and are preferentially localized in transcription regulatory regions including R-loops and open chromatin. Over 90% of protein-coding genes contain pG4-BS in their promoter or gene body. We observed generally higher pG4-BS content inR-loops and their flanks,longer genes that are associated with brain tissue, immune and developmental processes. Also, the presence of pG4-BS on both template and non-template strands in promoters is associated with oncogenesis, cardiovascular disease and stemness. Our G4-BS models predicted G4-forming ability in vitro with 91.5% accuracy. Analysis of G4-seq and CUT&Tag data strongly supports the existence of G4-BS conformations genome-wide. We reconstructed a novel G4-B 3D structure located in the E2F8 promoter. This study defines a large family of G4-like sequences, offering new insights into the essential biological functions and potential future therapeutic uses of G4-B.

Transcriptional activators in the early Drosophila embryo perform different kinetic roles.

Combinatorial regulation of gene expression by transcription factors (TFs) may in part arise from kinetic synergy-wherein TFs regulate different steps in the transcription cycle. Kinetic synergy requires that TFs play distinguishable kinetic roles. Here, we used live imaging to determine the kinetic roles of three TFs that activate transcription in the Drosophila embryo-Zelda, Bicoid, and Stat92E-by introducing their binding sites into the even-skipped stripe 2 enhancer. These TFs influence different sets of kinetic parameters, and their influence can change over time. All three TFs increased the fraction of transcriptionally active nuclei; Zelda also shortened the first-passage time into transcription and regulated the interval between transcription events. Stat92E also increased the lifetimes of active transcription. Different TFs can therefore play distinct kinetic roles in activating the transcription. This has consequences for understanding the composition and flexibility of regulatory DNA sequences and the biochemical function of TFs. A record of this paper's transparent peer review process is included in the supplemental information.

The Sequence Dependent Nanoscale Structure of CENP-A Nucleosomes.

CENP-A is a histone variant found in high abundance at the centromere in humans. At the centromere, this histone variant replaces the histone H3 found throughout the bulk chromatin. Additionally, the centromere comprises tandem repeats of α-satellite DNA, which CENP-A nucleosomes assemble upon. However, the effect of the DNA sequence on the nucleosome assembly and centromere formation remains poorly understood. Here, we investigated the structure of nucleosomes assembled with the CENP-A variant using Atomic Force Microscopy. We assembled both CENP-A nucleosomes and H3 nucleosomes on a DNA substrate containing an α-satellite motif and characterized their positioning and wrapping efficiency. We also studied CENP-A nucleosomes on the 601-positioning motif and non-specific DNA to compare their relative positioning and stability. CENP-A nucleosomes assembled on α-satellite DNA did not show any positional preference along the substrate, which is similar to both H3 nucleosomes and CENP-A nucleosomes on non-specific DNA. The range of nucleosome wrapping efficiency was narrower on α-satellite DNA compared with non-specific DNA, suggesting a more stable complex. These findings indicate that DNA sequence and histone composition may be two of many factors required for accurate centromere assembly.

TBIO-05. R-loops in pediatric brain tumors

Abstract R-loops are structures containing DNA-RNA hybrids and displaced single-stranded DNA, which are widely distributed across genomic regions. The generation and removal of R-loops is dynamically regulated by several factors including helicases and topoisomerases. Previously, we have identified high levels of R-loops associated with genomic instability in Embryonal Tumors with Multilayered Rosettes (ETMR) (Lambo et al. 2019), which is associated with a sensitivity to drugs targeting topoisomerases and DNA repair. However, it is unknown whether there are other pediatric cancers with high levels of R-loops, and whether R-loops in these tumors are also associated with genomic instability and sensitivity to specific treatments. Here, we have used real R-loops sequencing data sets, including ETMR DRIP-seq data, and the bioinformatic tool QmRLFS-finder (Wongsurawat et al. 2012), which can predict R-loops forming sequences (RLFSs) based on DNA sequence composition. By performing bioinformatic analyses of DRIP-seq data sets and integrating these with RLFSs from human reference genome simulations, we aim to analyze how reliable the different methods are to measure and predict R-loops and what role they have in cancer genomes. In ETMR, R-loops signals from DRIP-seq data were enriched around transcriptional start sites (TSS), comparable to R-loops distribution and enrichment in other published cell models. However, the overlap with RLFSs was limited, indicating that more validations with real DRIP-seq data from more ETMRs and other pediatric cancers are needed to validate the different methods. Taken together, until now, we developed methods to (semi-)quantify R-loops by combining real sequencing data and computational simulation approaches to investigate R-loops distribution in cancer genomes and to investigate whether R-loops are associated with breakpoints or other genetic aberrations in pediatric cancer entities.

Y chromosome-related genetic diseases

As a male-specific chromosome, the structure of Y chromosome is complex and lacks of recombination, with numerous repeating, amplifying and palindromic sequences. The research of Y chromosome is difficult and slow since there are few protein coding genes and a large amount of heterochromatin which has caused extreme difficulty for sequencing. In recent years, an increasing number of studies have been focused on the Y chromosome. With the completion of the sequencing of human Y chromosome, the rapid development of sequencing technology, and the composition of DNA sequences in human Y chromosomes and the determination of gene content. This paper has summarized the structural composition and genes function of human Y chromosome, as well as the related hereditary diseases, with an aim to provide reference for Y chromosome-related genetic research.

Enhancer-silencer transitions in the human genome.

Dual-function regulatory elements (REs), acting as enhancers in some cellular contexts and as silencers in others, have been reported to facilitate the precise gene regulatory response to developmental signals in Drosophila melanogaster. However, with few isolated examples detected, dual-function REs in mammals have yet to be systematically studied. We herein investigated this class of REs in the human genome and profiled their activity across multiple cell types. Focusing on enhancer–silencer transitions specific to the development of T cells, we built an accurate deep learning classifier of REs and identified about 12,000 silencers active in primary peripheral blood T cells that act as enhancers in embryonic stem cells. Compared with regular silencers, these dual-function REs are evolving under stronger purifying selection and are enriched for mutations associated with disease phenotypes and altered gene expression. In addition, they are enriched in the loci of transcriptional regulators, such as transcription factors (TFs) and chromatin remodeling genes. Dual-function REs consist of two intertwined but largely distinct sets of binding sites bound by either activating or repressing TFs, depending on the type of RE function in a given cell line. This indicates the recruitment of different TFs for different regulatory modes and a complex DNA sequence composition of these REs with dual activating and repressive encoding. With an estimated >6% of cell type–specific human silencers acting as dual-function REs, this overlooked class of REs requires a specific investigation on how their inherent functional plasticity might be a contributing factor to human diseases.

Formation of artificial chromosomes in Caenorhabditis elegans and analyses of their segregation in mitosis, DNA sequence composition and holocentromere organization.

To investigate how exogenous DNA concatemerizes to form episomal artificial chromosomes (ACs), acquire equal segregation ability and maintain stable holocentromeres, we injected DNA sequences with different features, including sequences that are repetitive or complex, and sequences with different AT-contents, into the gonad of Caenorhabditis elegans to form ACs in embryos, and monitored AC mitotic segregation. We demonstrated that AT-poor sequences (26% AT-content) delayed the acquisition of segregation competency of newly formed ACs. We also co-injected fragmented Saccharomyces cerevisiae genomic DNA, differentially expressed fluorescent markers and ubiquitously expressed selectable marker to construct a less repetitive, more complex AC. We sequenced the whole genome of a strain which propagates this AC through multiple generations, and de novo assembled the AC sequences. We discovered CENP-AHCP-3 domains/peaks are distributed along the AC, as in endogenous chromosomes, suggesting a holocentric architecture. We found that CENP-AHCP-3 binds to the unexpressed marker genes and many fragmented yeast sequences, but is excluded in the yeast extremely high-AT-content centromeric and mitochondrial DNA (> 83% AT-content) on the AC. We identified A-rich motifs in CENP-AHCP-3 domains/peaks on the AC and on endogenous chromosomes, which have some similarity with each other and similarity to some non-germline transcription factor binding sites.

Genome-wide profiling reveals functional interplay of DNA sequence composition, transcriptional activity, and nucleosome positioning in driving DNA supercoiling and helix destabilization in C. elegans.

DNA topology and alternative DNA structures are implicated in regulating diverse biological processes. Although biomechanical properties of these structures have been studied extensively in vitro, characterization in vivo, particularly in multicellular organisms, is limited. We devised new methods to map DNA supercoiling and single-stranded DNA in Caenorhabditis elegans embryos and diapause larvae. To map supercoiling, we quantified the incorporation of biotinylated psoralen into DNA using high-throughput sequencing. To map single-stranded DNA, we combined permanganate treatment with genome-wide sequencing of induced double-stranded breaks. We found high levels of negative supercoiling at transcription start sites (TSSs) in embryos. GC-rich regions flanked by a sharp GC-to-AT transition delineate boundaries of supercoil propagation. In contrast to TSSs in embryos, TSSs in diapause larvae showed dramatic reductions in negative supercoiling without concomitant attenuation of transcription, suggesting developmental-stage-specific regulation. To assess whether alternative DNA structures control chromosome architecture and gene expression, we examined DNA supercoiling in the context of X-Chromosome dosage compensation. We showed that the condensin dosage compensation complex creates negative supercoils locally at its highest-occupancy binding sites but found no evidence for large-scale supercoiling domains along X Chromosomes. In contrast to transcription-coupled negative supercoiling, single-strandedness, which is most pronounced at transcript end sites, is dependent on high AT content and symmetrically positioned nucleosomes. We propose that sharp transitions in sequence composition at functional genomic elements constitute a common regulatory code and that DNA structure and propagation of torsional stress at regulatory elements are critical parameters in shaping important developmental events.

Human nucleoli comprise multiple constrained territories, tethered to individual chromosomes.

It is unknown how ribosomal gene (rDNA) arrays from multiple chromosomal nucleolar organizers (NORs) partition within human nucleoli. Exploration of this paradigm for chromosomal organization is complicated by the shared DNA sequence composition of five NOR-bearing acrocentric chromosome p-arms. Here, we devise a methodology for genetic manipulation of individual NORs. Efficient "scarless" genome editing of rDNA repeats is achieved on "poised" human NORs held within monochromosomal cell hybrids. Subsequent transfer to human cells introduces "active" NORs yielding readily discernible functional customized ribosomes. We reveal that ribosome biogenesis occurs entirely within constrained territories, tethered to individual NORs inside a larger nucleolus.

Impact of parasitic lifestyle and different types of centromere organization on chromosome and genome evolution in the plant genus Cuscuta.

The parasitic genus Cuscuta (Convolvulaceae) is exceptional among plants with respect to centromere organization, including both monocentric and holocentric chromosomes, and substantial variation in genome size and chromosome number. We investigated 12 species representing the diversity of the genus in a phylogenetic context to reveal the molecular and evolutionary processes leading to diversification of their genomes. We measured genome sizes and investigated karyotypes and centromere organization using molecular cytogenetic techniques. We also performed low-pass whole genome sequencing and comparative analysis of repetitive DNA composition. A remarkable 102-fold variation in genome sizes (342-34734 Mbp/1C) was detected for monocentric Cuscuta species, while genomes of holocentric species were of moderate sizes (533-1545 Mbp/1C). The genome size variation was primarily driven by the differential accumulation of LTR-retrotransposons and satellite DNA. The transition to holocentric chromosomes in the subgenus Cuscuta was associated with loss of histone H2A phosphorylation and elimination of centromeric retrotransposons. In addition, basic chromosome number of holocentric species (x=7) was smaller than in monocentrics (x=15 or 16). We demonstrated that the transition to holocentricity in Cuscuta was accompanied by significant changes in epigenetic marks, chromosome number and the repetitive DNA sequence composition.

Sequence Composition Underlying Centromeric and Heterochromatic Genome Compartments of the Pacific Oyster Crassostrea gigas.

Segments of the genome enriched in repetitive sequences still present a challenge and are omitted in genome assemblies. For that reason, the exact composition of DNA sequences underlying the heterochromatic regions and the active centromeres are still unexplored for many organisms. The centromere is a crucial region of eukaryotic chromosomes responsible for the accurate segregation of genetic material. The typical landmark of centromere chromatin is the rapidly-evolving variant of the histone H3, CenH3, while DNA sequences packed in constitutive heterochromatin are associated with H3K9me3-modified histones. In the Pacific oyster Crassostrea gigas we identified its centromere histone variant, Cg-CenH3, that shows stage-specific distribution in gonadal cells. In order to investigate the DNA composition of genomic regions associated with the two specific chromatin types, we employed chromatin immunoprecipitation followed by high-throughput next-generation sequencing of the Cg-CenH3- and H3K9me3-associated sequences. CenH3-associated sequences were assigned to six groups of repetitive elements, while H3K9me3-associated-ones were assigned only to three. Those associated with CenH3 indicate the lack of uniformity in the chromosomal distribution of sequences building the centromeres, being also in the same time dispersed throughout the genome. The heterochromatin of C. gigas exhibited general paucity and limited chromosomal localization as predicted, with H3K9me3-associated sequences being predominantly constituted of DNA transposons.

Circulating cell-free DNA is predominantly composed of retrotransposable elements and non-telomeric satellite DNA

Circulating cell-free DNAs (cfDNAs) are DNA fragments which can be isolated from mammalian blood serum or plasma. In order to gain deeper insight into their origin(s), we have characterized the composition of human and cattle cfDNA via large-scale analyses of high-throughput sequencing data. We observed significant differences between the composition of cfDNA in serum/plasma and the corresponding DNA sequence composition of the human genome. Retrotransposable elements and non-telomeric satellite DNA were particularly overrepresented in the cfDNA population, while telomeric satellite DNA was underrepresented. This was consistently observed for human plasma, bovine serum and for the supernatant of human cancer cell cultures. Our results suggest that reverse transcription of retrotransposable elements and secondary-structure formation during the replication of satellite DNA are contributing to the composition of the cfDNA molecules in the mammalian blood stream. We believe that our work is an important step towards the understanding of the biogenesis of cfDNAs and thus may also facilitate the future exploitation of their diagnostic potential.

DNA sequence-dependent chromatin architecture and nuclear hubs formation

In this study, by exploring chromatin conformation capture data, we show that the nuclear segregation of Topologically Associated Domains (TADs) is contributed by DNA sequence composition. GC-peaks and valleys of TADs strongly influence interchromosomal interactions and chromatin 3D structure. To gain insight on the compositional and functional constraints associated with chromatin interactions and TADs formation, we analysed intra-TAD and intra-loop GC variations. This led to the identification of clear GC-gradients, along which, the density of genes, super-enhancers, transcriptional activity, and CTCF binding sites occupancy co-vary non-randomly. Further, the analysis of DNA base composition of nucleolar aggregates and nuclear speckles showed strong sequence-dependant effects. We conjecture that dynamic DNA binding affinity and flexibility underlay the emergence of chromatin condensates, their growth is likely promoted in mechanically soft regions (GC-rich) of the lowest chromatin and nucleosome densities. As a practical perspective, the strong linear association between sequence composition and interchromosomal contacts can help define consensus chromatin interactions, which in turn may be used to study alternative states of chromatin architecture.

A common genomic code for chromatin architecture and recombination landscape.

Recent findings established a link between DNA sequence composition and interphase chromatin architecture and explained the evolutionary conservation of TADs (Topologically Associated Domains) and LADs (Lamina Associated Domains) in mammals. This prompted us to analyse conformation capture and recombination rate data to study the relationship between chromatin architecture and recombination landscape of human and mouse genomes. The results reveal that: (1) low recombination domains and blocks of elevated linkage disequilibrium tend to coincide with TADs and isochores, indicating co-evolving regulatory elements and genes in insulated neighbourhoods; (2) double strand break (DSB) and recombination frequencies increase in the short loops of GC-rich TADs, whereas recombination cold spots are typical of LADs and (3) the binding and loading of proteins, which are critical for DSB and meiotic recombination (SPO11, DMC1, H3K4me3 and PRMD9) are higher in GC-rich TADs. One explanation for these observations is that the occurrence of DSB and recombination in meiotic cells are associated with compositional and epigenetic features (genomic code) that influence DNA stiffness/flexibility and appear to be similar to those guiding the chromatin architecture in the interphase nucleus of pre-leptotene cells.

Molecular cytogenetic analysis of genome-specific repetitive elements in Citrus clementina Hort. Ex Tan. and its taxonomic implications

BackgroundClementine mandarin (Citrus clementina Hort. ex Tan.) is one of the most famous and widely grown citrus cultivars worldwide. Variations in relation to the composition and distribution of repetitive DNA sequences that dominate greatly in eukaryote genomes are considered to be species-, genome-, or even chromosome-specific. Repetitive DNA-based fluorescence in situ hybridization (FISH) is a powerful tool for molecular cytogenetic study. However, to date few studies have involved in the repetitive elements and cytogenetic karyotype of Clementine.ResultsA graph-based similarity sequence read clustering methodology was performed to analyze the repetitive DNA families in the Clementine genome. The bioinformatics analysis showed that repetitive DNAs constitute 41.95% of the Clementine genome, and the majority of repetitive elements are retrotransposons and satellite DNAs. Sequential multicolor FISH using a probe mix that contained CL17, four satellite DNAs, two rDNAs and an oligonucleotide of (TTTAGGG)3 was performed with Clementine somatic metaphase chromosomes. An integrated karyotype of Clementine was established based on unequivocal and reproducible chromosome discriminations. The distribution patterns of these probes in several Citrus, Poncirus and Fortunella species were summarized through extensive FISH analyses. Polymorphism and heterozygosity were commonly observed in the three genera. Some asymmetrical FISH loci in Clementine were in agreement with its hybrid origin.ConclusionsThe composition and abundance of repetitive elements in the Clementine genome were reanalyzed. Multicolor FISH-based karyotyping provided direct visual proof of the heterozygous nature of Clementine chromosomes with conspicuous asymmetrical FISH hybridization signals. We detected some similar and variable distribution patterns of repetitive DNAs in Citrus, Poncirus, and Fortunella, which revealed notable conservation among these genera, as well as obvious polymorphism and heterozygosity, indicating the potential utility of these repetitive element markers for the study of taxonomic, phylogenetic and evolutionary relationships in the future.

DNA Sequence and Histone Core Composition Control the Unwrapping of DNA from Nucleosome Core Particles

As controlling access to DNA is a major component of transcriptional regulation, understanding the effects of interactions between core proteins and the DNA is crucial to gaining an overall understanding of transcription. We report on the use of small angle X-ray scattering with contrast variation (CV-SAXS) to analyze DNA unwrapping from nucleosome core particles. Interactions between the histone core and the DNA are destabilized using increasing concentrations of neutral salt, lowering the strength of electrostatic interactions that bind the components together. As the complex is destabilized many unwrapping states are accessible and the DNA explores possible unwrapping pathways. The CV-SAXS data reports on full range of structures present in the sample. They are fit using sequence-dependent models of unwrapped structures with thermal motions, selected by an ensemble optimization method [1]. Different combinations of DNA sequence and histone composition also alter their interactions, and the unwrapped structures reflect biophysical differences in the nucleosome compositions. 1. Mauney A., Tokuda J.M., Gloss L.M., Gonzalez O., Pollack L. Local DNA sequence controls the cooperativity and asymmetry of DNA unwrapping from nucleosome core particles. Biophys J. 2018 Sep 4;115(5):773-781