Elimination of subtelomeric repeat sequences exerts little effect on telomere essential functions in Saccharomyces cerevisiae

Telomeres, which are chromosomal end structures, play a crucial role in maintaining genome stability and integrity in eukaryotes. In the baker's yeast Saccharomyces cerevisiae, the X- and Y'-elements are subtelomeric repetitive sequences found in all 32 and 17 telomeres, respectively. While the Y'-elements serve as a backup for telomere functions in cells lacking telomerase, the function of the X-elements remains unclear. This study utilized the S. cerevisiae strain SY12, which has three chromosomes and six telomeres, to investigate the role of X-elements (as well as Y'-elements) in telomere maintenance. Deletion of Y'-elements (SY12YΔ), X-elements (SY12XYΔ+Y), or both X- and Y'-elements (SY12XYΔ) did not impact the length of the terminal TG1-3 tracks or telomere silencing. However, inactivation of telomerase in SY12YΔ, SY12XYΔ+Y, and SY12XYΔ cells resulted in cellular senescence and the generation of survivors. These survivors either maintained their telomeres through homologous recombination-dependent TG1-3 track elongation or underwent microhomology-mediated intra-chromosomal end-to-end joining. Our findings indicate the non-essential role of subtelomeric X- and Y'-elements in telomere regulation in both telomerase-proficient and telomerase-null cells and suggest that these elements may represent remnants of S. cerevisiae genome evolution. Furthermore, strains with fewer or no subtelomeric elements exhibit more concise telomere structures and offer potential models for future studies in telomere biology.

Mathematical model for the role of multiple pericentromeric repeats on heterochromatin assembly

Although the length and constituting sequences for pericentromeric repeats are highly variable across eukaryotes, the presence of multiple pericentromeric repeats is one of the conserved features of the eukaryotic chromosomes. Pericentromeric heterochromatin is often misregulated in human diseases, with the expansion of pericentromeric repeats in human solid cancers. In this article, we have developed a mathematical model of the RNAi-dependent methylation of H3K9 in the pericentromeric region of fission yeast. Our model, which takes copy number as an explicit parameter, predicts that the pericentromere is silenced only if there are many copies of repeats. It becomes bistable or desilenced if the copy number of repeats is reduced. This suggests that the copy number of pericentromeric repeats alone can determine the fate of heterochromatin silencing in fission yeast. Through sensitivity analysis, we identified parameters that favor bistability and desilencing. Stochastic simulation shows that faster cell division and noise favor the desilenced state. These results show the unexpected role of pericentromeric repeat copy number in gene silencing and provide a quantitative basis for how the copy number allows or protects repetitive and unique parts of the genome from heterochromatin silencing, respectively.

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.

Selfish conflict underlies RNA-mediated parent-of-origin effects

Genomic imprinting-the non-equivalence of maternal and paternal genomes-is a critical process that has evolved independently in many plant and mammalian species1,2. According to kinship theory, imprinting is the inevitable consequence of conflictive selective forces acting on differentially expressed parental alleles3,4. Yet, how these epigenetic differences evolve in the first place is poorly understood3,5,6. Here we report the identification and molecular dissection of a parent-of-origin effect on gene expression that might help to clarify this fundamental question. Toxin-antidote elements (TAs) are selfish elements that spread in populations by poisoning non-carrier individuals7-9. In reciprocal crosses between two Caenorhabditis tropicalis wild isolates, we found that the slow-1/grow-1 TA is specifically inactive when paternally inherited. This parent-of-origin effect stems from transcriptional repression of the slow-1 toxin by the PIWI-interacting RNA (piRNA) host defence pathway. The repression requires PIWI Argonaute and SET-32 histone methyltransferase activities and is transgenerationally inherited via small RNAs. Remarkably, when slow-1/grow-1 is maternally inherited, slow-1 repression is halted by a translation-independent role of its maternal mRNA. That is, slow-1 transcripts loaded into eggs-but not SLOW-1 protein-are necessary and sufficient to counteract piRNA-mediated repression. Our findings show that parent-of-origin effects can evolve by co-option of the piRNA pathway and hinder the spread of selfish genes that require sex for their propagation.

Enzymatic Assembly of Small Synthetic Genes with Repetitive Elements

Gene synthesis efficiency has greatly improved in recent years but is limited when it comes to repetitive sequences, which results in synthesis failure or delays by DNA synthesis vendors. This represents a major obstacle for the development of synthetic biology since repetitive elements are increasingly being used in the design of genetic circuits and design of biomolecular nanostructures. Here, we describe a method for the assembly of small synthetic genes with repetitive elements: First, a gene of interest is split in silico into small synthons of up to 80 base pairs flanked by Golden-Gate-compatible overhangs. Then, synthons are made by oligo extension and finally assembled into a synthetic gene by Golden Gate Assembly. We demonstrate the method by constructing eight challenging genes with repetitive elements, e.g., multiple repeats of RNA aptamers and RNA origami scaffolds with multiple identical aptamers. The genes range in size from 133 to 456 base pairs and are assembled with fidelities of up to 87.5%. The method was developed to facilitate our own specific research but may be of general use for constructing challenging and repetitive genes and, thus, a valuable addition to the molecular cloning toolbox.

Dimeric structures of DNA ATTTC repeats promoted by divalent cations

Structural studies of repetitive DNA sequences may provide insights why and how certain repeat instabilities in their number and nucleotide sequence are managed or even required for normal cell physiology, while genomic variability associated with repeat expansions may also be disease-causing. The pentanucleotide ATTTC repeats occur in hundreds of genes important for various cellular processes, while their insertion and expansion in noncoding regions are associated with neurodegeneration, particularly with subtypes of spinocerebellar ataxia and familial adult myoclonic epilepsy. We describe a new striking domain-swapped DNA-DNA interaction triggered by the addition of divalent cations, including Mg2+ and Ca2+. The results of NMR characterization of d(ATTTC)3 in solution show that the oligonucleotide folds into a novel 3D architecture with two central C:C+ base pairs sandwiched between a couple of T:T base pairs. This structural element, referred to here as the TCCTzip, is characterized by intercalative hydrogen-bonding, while the nucleobase moieties are poorly stacked. The 5'- and 3'-ends of TCCTzip motif are connected by stem-loop segments characterized by A:T base pairs and stacking interactions. Insights embodied in the non-canonical DNA structure are expected to advance our understanding of why only certain pyrimidine-rich DNA repeats appear to be pathogenic, while others can occur in the human genome without any harmful consequences.

Tigerfish designs oligonucleotide-based in situ hybridization probes targeting intervals of highly repetitive DNA at the scale of genomes

Fluorescent in situ hybridization (FISH) is a powerful method for the targeted visualization of nucleic acids in their native contexts. Recent technological advances have leveraged computationally designed oligonucleotide (oligo) probes to interrogate > 100 distinct targets in the same sample, pushing the boundaries of FISH-based assays. However, even in the most highly multiplexed experiments, repetitive DNA regions are typically not included as targets, as the computational design of specific probes against such regions presents significant technical challenges. Consequently, many open questions remain about the organization and function of highly repetitive sequences. Here, we introduce Tigerfish, a software tool for the genome-scale design of oligo probes against repetitive DNA intervals. We showcase Tigerfish by designing a panel of 24 interval-specific repeat probes specific to each of the 24 human chromosomes and imaging this panel on metaphase spreads and in interphase nuclei. Tigerfish extends the powerful toolkit of oligo-based FISH to highly repetitive DNA.

Higher evolutionary dynamics of gene copy number for Drosophila glue genes located near short repeat sequences

BACKGROUND

During evolution, genes can experience duplications, losses, inversions and gene conversions. Why certain genes are more dynamic than others is poorly understood. Here we examine how several Sgs genes encoding glue proteins, which make up a bioadhesive that sticks the animal during metamorphosis, have evolved in Drosophila species.

RESULTS

We examined high-quality genome assemblies of 24 Drosophila species to study the evolutionary dynamics of four glue genes that are present in D. melanogaster and are part of the same gene family - Sgs1, Sgs3, Sgs7 and Sgs8 - across approximately 30 millions of years. We annotated a total of 102 Sgs genes and grouped them into 4 subfamilies. We present here a new nomenclature for these Sgs genes based on protein sequence conservation, genomic location and presence/absence of internal repeats. Two types of glue genes were uncovered. The first category (Sgs1, Sgs3x, Sgs3e) showed a few gene losses but no duplication, no local inversion and no gene conversion. The second group (Sgs3b, Sgs7, Sgs8) exhibited multiple events of gene losses, gene duplications, local inversions and gene conversions. Our data suggest that the presence of short "new glue" genes near the genes of the latter group may have accelerated their dynamics.

CONCLUSIONS

Our comparative analysis suggests that the evolutionary dynamics of glue genes is influenced by genomic context. Our molecular, phylogenetic and comparative analysis of the four glue genes Sgs1, Sgs3, Sgs7 and Sgs8 provides the foundation for investigating the role of the various glue genes during Drosophila life.

X Chromosome-Specific Repeats in Non-Domestic Bovidae

Repetitive sequences form a substantial and still enigmatic part of the mammalian genome. We isolated repetitive DNA blocks of the X chromosomes of three species of the family Bovidae: Kobus defassa (KDEXr sequence), Bos taurus (BTAXr sequence) and Antilope cervicapra (ACEXr sequence). The copy numbers of the isolated sequences were assessed using qPCR, and their chromosomal localisations were analysed using FISH in ten bovid tribes and in outgroup species. Besides their localisation on the X chromosome, their presence was also revealed on the Y chromosome and autosomes in several species. The KDEXr sequence abundant in most Bovidae species also occurs in distant taxa (Perissodactyla and Carnivora) and seems to be evolutionarily older than BTAXr and ACEXr. The ACEXr sequence, visible only in several Antilopini species using FISH, is probably the youngest, and arised in an ancestor common to Bovidae and Cervidae. All three repetitive sequences analysed in this study are interspersed among gene-rich regions on the X chromosomes, apparently preventing the crossing-over in their close vicinity. This study demonstrates that repetitive sequences on the X chromosomes have undergone a fast evolution, and their variation among related species can be beneficial for evolutionary studies.

Structure, Topology, and Stability of Multiple G-quadruplexes in Long Telomeric Overhangs

Telomeres and their single stranded overhangs gradually shorten with successive cell divisions, as part of the natural aging process, but can be elongated by telomerase, a nucleoprotein complex which is activated in the majority of cancers. This prominent implication in cancer and aging has made the repetitive telomeric sequences (TTAGGG repeats) and the G-quadruplex structures that form in their overhangs the focus of intense research in the past several decades. However, until recently most in vitro efforts to understand the structure, stability, dynamics, and interactions of telomeric overhangs had been focused on short sequences that are not representative of longer sequences encountered in a physiological setting. In this review, we will provide a broad perspective about telomeres and associated factors, and introduce the agents and structural characteristics involved in organizing, maintaining, and protecting telomeric DNA. We will also present a summary of recent research performed on long telomeric sequences, nominally defined as those that can form two or more tandem G-quadruplexes, i.e., which contain eight or more TTAGGG repeats. Results of experimental studies using a broad array of experimental tools, in addition to recent computational efforts will be discussed, particularly in terms of their implications for the stability, folding topology, and compactness of the tandem G-quadruplexes that form in long telomeric overhangs.

Evolution of the spinach sex-linked region within a rarely recombining pericentromeric region

Sex chromosomes have evolved independently in many different plant lineages. Here, we describe reference genomes for spinach (Spinacia oleracea) X and Y haplotypes by sequencing homozygous XX females and YY males. The long arm of 185-Mb chromosome 4 carries a 13-Mb X-linked region (XLR) and 24.1-Mb Y-linked region (YLR), of which 10 Mb is Y specific. We describe evidence that this reflects insertions of autosomal sequences creating a "Y duplication region" or "YDR" whose presence probably directly reduces genetic recombination in the immediately flanking regions, although both the X and Y sex-linked regions are within a large pericentromeric region of chromosome 4 that recombines rarely in meiosis of both sexes. Sequence divergence estimates using synonymous sites indicate that YDR genes started diverging from their likely autosomal progenitors about 3 MYA, around the time when the flanking YLR stopped recombining with the XLR. These flanking regions have a higher density of repetitive sequences in the YY than the XX assembly and include slightly more pseudogenes compared with the XLR, and the YLR has lost about 11% of the ancestral genes, suggesting some degeneration. Insertion of a male-determining factor would have caused Y linkage across the entire pericentromeric region, creating physically small, highly recombining, terminal pseudoautosomal regions. These findings provide a broader understanding of the origin of sex chromosomes in spinach.

exo-FISH: Protocol for detecting DNA breaks in repetitive regions of mammalian genomes

Detecting DNA breaks in defined regions of the genome is critical to advancing our understanding of genome stability maintenance. Here, we present exo-FISH, a protocol to label exposed single-stranded DNA in defined repetitive regions of mammalian genomes by combining in vitro restriction enzyme digestion on fixed cells with fluorescence in situ hybridization (FISH). We describe steps for cell harvesting and fixation, slide treatments, and FISH probe hybridization. We then detail procedures for imaging and analysis. For complete details on the use and execution of this protocol, please refer to Saayman et al. (2023).1.

Diversity and Evolution of Highly Repetitive DNA Sequences Constituting Chromosome Site-Specific Heterochromatin in Two Gerbillinae Species

Constitutive heterochromatin, consisting of repetitive sequences, diverges very rapidly; therefore, its nucleotide sequences and chromosomal distributions are often largely different, even between closely related species. The chromosome C-banding patterns of two Gerbillinae species, Meriones unguiculatus and Gerbillus perpallidus, vary greatly, even though they belong to the same subfamily. To understand the evolution of C-positive heterochromatin in these species, we isolated highly repetitive sequences, determined their nucleotide sequences, and characterized them using chromosomal and filter hybridization. We obtained a centromeric repeat (MUN-HaeIII) and a chromosome 13-specific repeat (MUN-EcoRI) from M. unguiculatus. We also isolated a centromeric/pericentromeric repeat (GPE-MBD) and an interspersed-type repeat that was predominantly amplified in the X and Y chromosomes (GPE-EcoRI) from G. perpallidus. GPE-MBD was found to contain a 17-bp motif that is essential for binding to the centromere-associated protein CENP-B. This indicates that it may play a role in the formation of a specified structure and/or function of centromeres. The nucleotide sequences of the three sequence families, except GPE-EcoRI, were conserved only in Gerbillinae. GPE-EcoRI was derived from the long interspersed nuclear elements 1 retrotransposon and showed sequence homology throughout Muridae and Cricetidae species, indicating that the repeat sequence occurred at least in the common ancestor of Muridae and Cricetidae. Due to a lack of assembly data of highly repetitive sequences constituting heterochromatin in whole-genome sequences of vertebrate species published to date, the knowledge obtained in this study provides useful information for a deep understanding of the evolution of repetitive sequences in not only rodents but also in mammals.

Continuous synthesis of E. coli genome sections and Mb-scale human DNA assembly

Whole-genome synthesis provides a powerful approach for understanding and expanding organism function1-3. To build large genomes rapidly, scalably and in parallel, we need (1) methods for assembling megabases of DNA from shorter precursors and (2) strategies for rapidly and scalably replacing the genomic DNA of organisms with synthetic DNA. Here we develop bacterial artificial chromosome (BAC) stepwise insertion synthesis (BASIS)-a method for megabase-scale assembly of DNA in Escherichia coli episomes. We used BASIS to assemble 1.1 Mb of human DNA containing numerous exons, introns, repetitive sequences, G-quadruplexes, and long and short interspersed nuclear elements (LINEs and SINEs). BASIS provides a powerful platform for building synthetic genomes for diverse organisms. We also developed continuous genome synthesis (CGS)-a method for continuously replacing sequential 100 kb stretches of the E. coli genome with synthetic DNA; CGS minimizes crossovers1,4 between the synthetic DNA and the genome such that the output for each 100 kb replacement provides, without sequencing, the input for the next 100 kb replacement. Using CGS, we synthesized a 0.5 Mb section of the E. coli genome-a key intermediate in its total synthesis1-from five episomes in 10  days. By parallelizing CGS and combining it with rapid oligonucleotide synthesis and episome assembly5,6, along with rapid methods for compiling a single genome from strains bearing distinct synthetic genome sections1,7,8, we anticipate that it will be possible to synthesize entire E. coli genomes from functional designs in less than 2 months.

A new type of tandem repeats in Myotis petax (Chiroptera, Vespertilionidae) mitochondrial control region

BACKGROUND

Tandem repeats in mitochondrial DNA control region are known to different animal taxa, including bat species of the family Vespertilionidae. The long R1-repeats in the bat ETAS-domain are often presented in a variable copy number and may exhibit both inter-individual and intra-individual sequence diversity. The function of repeats in the control region is still unclear, but it has been shown that repetitive sequences in some animal groups (shrews, cats and sheep) may include parts of ETAS1 and ETAS2 conservative blocks of mitochondrial DNA.

METHODS AND RESULTS

Analysis of the control region sequences for 31 Myotis petax specimens allowed the identification of the inter-individual variability and clarification of the composition of the R1-repeats. The copy number of the R1-repeats varies from 4 to 7 in individuals. The specimens examined do not exhibit a size heteroplasmy previously described for Myotis species. The unusual short 30 bp R1-repeats have been detected in M. petax for the first time. The ten specimens from Amur Region and Primorsky Territory have one or two copies of these additional repeats.

CONCLUSIONS

It was determined that the R1-repeats in M. petax control region consist of parts of the ETAS1 and ETAS2 blocks. The origin of the additional repeats seems to be related to the 51 bp deletion in the central part of the R1-repeat unit and subsequent duplication. Comparison of repetitive sequences in the control region of closely-related Myotis species identified the occurrence of incomplete repeats also resulting from the short deletions, but distinct from additional repeats of M. petax.

Molecular Cytogenetic Characterization of C-Band-Positive Heterochromatin of the Greater Long-Tailed Hamster (<b><i>Tscherskia triton</i></b>, Cricetinae)

The greater long-tailed hamster (<i>Tscherskia triton</i>, Cricetinae) has a unique karyotype (2n = 28), containing 11 pairs of acrocentric chromosomes with large C-band-positive centromeric heterochromatin blocks. To understand the origin and evolutionary process of heterochromatin in this species, we isolated novel families of chromosome site-specific highly repetitive DNA sequences from <i>Taq</i>I-digested genomic DNA and then characterized them by chromosome in situ and filter hybridization. The <i>Taq</i>I-families of repetitive sequences were classified into 2 types by their genome organization and chromosomal distribution: the 110-bp repeated sequence organized in large tandem arrays (as satellite DNA), localized to centromeric C-positive heterochromatin of acrocentric autosomes (chromosomes 1–11) and submetacentric X chromosome, and the 405-bp repeated sequence that was composed of 30–32-bp internal repeats, distributed in the pericentromeric region on the short arms of X and Y chromosomes. The repetitive sequences did not cross-hybridize with genomic DNA of any genera of Cricetinae (<i>Mesocricetus</i>, <i>Cricetulus</i>, and <i>Phodopus</i>). These results suggest that the 110-bp and 405-bp repeats rapidly diverged in the lineage of <i>T. triton</i>, evolving in a concerted manner among autosomes and X chromosome and within X and Y chromosomes, respectively. The 110-bp centromeric repeat contained a 17-bp motif in which 9 bases are essential for binding with the centromere-associated protein CENP-B, suggesting the possibility that the 110-bp major satellite DNA carrying the 17-bp motif may have a role in the formation of specified structure and/or function of centromeres in <i>T. triton</i>.

A Spontaneous Inversion of the X Chromosome Heterochromatin Provides a Tool for Studying the Structure and Activity of the Nucleolus in Drosophila melanogaster

The pericentromeric heterochromatin is largely composed of repetitive sequences, making it difficult to analyze with standard molecular biological methods. At the same time, it carries many functional elements with poorly understood mechanisms of action. The search for new experimental models for the analysis of heterochromatin is an urgent task. In this work, we used the Rif1 mutation, which suppresses the underreplication of all types of repeated sequences, to analyze heterochromatin regions in polytene chromosomes of Drosophila melanogaster. In the Rif1 background, we discovered and described in detail a new inversion, In(1)19EHet, which arose on a chromosome already carrying the In(1)sc⁸ inversion and transferred a large part of X chromosome heterochromatin, including the nucleolar organizer to a new euchromatic environment. Using nanopore sequencing and FISH, we have identified the eu- and heterochromatin breakpoints of In(1)19EHet. The combination of the new inversion and the Rif1 mutation provides a promising tool for studies of X chromosome heterochromatin structure, nucleolar organization, and the nucleolar dominance phenomenon. In particular, we found that, with the complete polytenization of rDNA repeats, the nucleolus consists of a cloud-like structure corresponding to the classical nucleolus of polytene chromosomes, as well as an unusual intrachromosomal structure containing alternating transcriptionally active and inactive regions.

Cloning and expansion of repetitive DNA sequences

Repeat and structure-prone DNA sequences comprise a large proportion of the human genome. The instability of these sequences has been implicated in a range of diseases, including cancers and neurodegenerative disorders. However, the mechanism of pathogenicity is poorly understood. As such, further studies on repetitive DNA are required. Cloning and maintaining repeat-containing substrates is challenging due to their inherent ability to form non-B DNA secondary structures which are refractory to DNA polymerases and prone to undergo rearrangements. Here, we describe an approach to clone and expand tandem-repeat DNA without interruptions, thereby allowing for its manipulation and subsequent investigation.

The Complete Mitochondrial Genome of Ophioglossum vulgatum L. Is with Highly Repetitive Sequences: Intergenomic Fragment Transfer and Phylogenetic Analysis

Many plant mitochondrial (mt) genomes have been sequenced but few in ferns. Ophioglossum vulgatum represents a typical species of fern genus Ophioglossum with medicinal and scientific value. However, its mt genome structure remains to be characterized. This study assembled and annotated the complete O. vulgatum mt genome and presented its structural characters and repeat sequences firstly. Its mt and chloroplast (cp) transfer sequences were explored, and the phylogenetic significance of both mt and cp genomes was also evaluated at the family level. Our results showed that the complete mt genome of O. vulgatum is a single circular genome of 369,673 bp in length, containing 5000 dispersed repetitive sequences. Phylogenetic trees reconstructed from cp and mt genomes displayed similar topologies, but also showed subtle differences at certain nodes. There exist 4818 bp common gene fragments between cp and mt genomes, of which more than 70% are located in tRNA intergenic regions (in mt). In conclusion, we assembled the complete mt genome of O. vulgatum, identified its remarkable structural characters, and provided new insights on ferns. The complementary results derived from mt and cp phylogeny highlighted that some higher taxonomic-level phylogenetic relationships among ferns remain to be resolved.

3D genome, on repeat: Higher-order folding principles of the heterochromatinized repetitive genome

Nearly half of the human genome is comprised of diverse repetitive sequences ranging from satellite repeats to retrotransposable elements. Such sequences are susceptible to stepwise expansions, duplications, inversions, and recombination events which can compromise genome function. In this review, we discuss the higher-order folding mechanisms of compartmentalization and loop extrusion and how they shape, and are shaped by, heterochromatin. Using primarily mammalian model systems, we contrast mechanisms governing H3K9me3-mediated heterochromatinization of the repetitive genome and highlight emerging links between repetitive elements and chromatin folding.

Long-read mapping to repetitive reference sequences using Winnowmap2

Approximately 5-10% of the human genome remains inaccessible due to the presence of repetitive sequences such as segmental duplications and tandem repeat arrays. We show that existing long-read mappers often yield incorrect alignments and variant calls within long, near-identical repeats, as they remain vulnerable to allelic bias. In the presence of a nonreference allele within a repeat, a read sampled from that region could be mapped to an incorrect repeat copy. To address this limitation, we developed a new long-read mapping method, Winnowmap2, by using minimal confidently alignable substrings. Winnowmap2 computes each read mapping through a collection of confident subalignments. This approach is more tolerant of structural variation and more sensitive to paralog-specific variants within repeats. Our experiments highlight that Winnowmap2 successfully addresses the issue of allelic bias, enabling more accurate downstream variant calls in repetitive sequences.

Interaction of Proteins with Inverted Repeats and Cruciform Structures in Nucleic Acids

Cruciforms occur when inverted repeat sequences in double-stranded DNA adopt intra-strand hairpins on opposing strands. Biophysical and molecular studies of these structures confirm their characterization as four-way junctions and have demonstrated that several factors influence their stability, including overall chromatin structure and DNA supercoiling. Here, we review our understanding of processes that influence the formation and stability of cruciforms in genomes, covering the range of sequences shown to have biological significance. It is challenging to accurately sequence repetitive DNA sequences, but recent advances in sequencing methods have deepened understanding about the amounts of inverted repeats in genomes from all forms of life. We highlight that, in the majority of genomes, inverted repeats are present in higher numbers than is expected from a random occurrence. It is, therefore, becoming clear that inverted repeats play important roles in regulating many aspects of DNA metabolism, including replication, gene expression, and recombination. Cruciforms are targets for many architectural and regulatory proteins, including topoisomerases, p53, Rif1, and others. Notably, some of these proteins can induce the formation of cruciform structures when they bind to DNA. Inverted repeat sequences also influence the evolution of genomes, and growing evidence highlights their significance in several human diseases, suggesting that the inverted repeat sequences and/or DNA cruciforms could be useful therapeutic targets in some cases.

PlantRep: a database of plant repetitive elements

We re-annotated repeats of 459 plant genomes and released a new database: PlantRep ( http://www.plantrep.cn/ ). PlantRep sheds lights of repeat evolution and provides fundamental data for deep exploration of genome.

Polypolish: Short-read polishing of long-read bacterial genome assemblies

Long-read-only bacterial genome assemblies usually contain residual errors, most commonly homopolymer-length errors. Short-read polishing tools can use short reads to fix these errors, but most rely on short-read alignment which is unreliable in repeat regions. Errors in such regions are therefore challenging to fix and often remain after short-read polishing. Here we introduce Polypolish, a new short-read polisher which uses all-per-read alignments to repair errors in repeat sequences that other polishers cannot. Polypolish performed well in benchmarking tests using both simulated and real reads, and it almost never introduced errors during polishing. The best results were achieved by using Polypolish in combination with other short-read polishers.

Recent improvements in Oxford Nanopore Technologies sequencing platforms and assembly algorithms have made it easier than ever to generate complete bacterial genome sequences. However, Oxford Nanopore genome sequences suffer from errors that limit their utility in downstream analyses. To fix these errors, one can ‘polish’ the genome with Illumina sequencing, exploiting the fact that Oxford Nanopore and Illumina sequencing have different error profiles. There are several polishing tools which can fix most errors in an Oxford Nanopore genome, but they struggle with errors in repetitive regions of the genome. With this in mind, we have developed a polisher, Polypolish, which uses a novel approach that allows it to fix more errors in genomic repeats. Our results show that Polypolish is both effective at repairing sequence errors and very unlikely to introduce new errors. Polypolish can often fix errors that other polishers cannot and vice versa, so the best results come from using a combination of tools. Polypolish therefore has an important role in bacterial genome assembly methods that aim for the highest possible sequence accuracy.

Bioinformatics Approaches for Determining the Functional Impact of Repetitive Elements on Non-coding RNAs

With a large number of annotated non-coding RNAs (ncRNAs), repetitive sequences are found to constitute functional components (termed as repetitive elements) in ncRNAs that perform specific biological functions. Bioinformatics analysis is a powerful tool for improving our understanding of the role of repetitive elements in ncRNAs. This chapter summarizes recent findings that reveal the role of repetitive elements in ncRNAs. Furthermore, relevant bioinformatics approaches are systematically reviewed, which promises to provide valuable resources for studying the functional impact of repetitive elements on ncRNAs.

Genome Sequence of Chrysotila roscoffensis, a Coccolithphore Contributed to Global Biogeochemical Cycles

Chrysotila is a genus of coccolithophores. Together with Emiliania, it is one of the representative genera in the Haptophyta which have been extensively studied. They are photosynthetic unicellular marine algae sharing the common characteristic of the production of CaCO₃ platelets (coccoliths) on the surface of their cells and are crucial contributors to global biogeochemical cycles. Here, we report the genome assembly of Chrysotila roscoffensis. The assembled genome size was ~636 Mb distributed across 769 scaffolds with N50 of 1.63 Mb, and maximum contig length of ~2.6 Mb. Repetitive elements accounted for approximately 59% of the genome. A total of 23,341 genes were predicted from C. roscoffensis genome. The divergence time between C. roscoffensis and Emiliania huxleyi was estimated to be around 537.6 Mya. Gene families related to cytoskeleton, cellular motility and morphology, and ion transport were expanded. The genome of C. roscoffensis will provide a foundation for understanding the genetic and phenotypic diversification and calcification mechanisms of coccolithophores.

Assembly of the Complete Mitochondrial Genome of Chinese Plum (Prunus salicina): Characterization of Genome Recombination and RNA Editing Sites

Despite the significant progress that has been made in the genome sequencing of Prunus, this area of research has been lacking a systematic description of the mitochondrial genome of this genus for a long time. In this study, we assembled the mitochondrial genome of the Chinese plum (Prunus salicina) using Illumina and Oxford Nanopore sequencing data. The mitochondrial genome size of P. salicina was found to be 508,035 base pair (bp), which is the largest reported in the Rosaceae family to date, and P. salicina was shown to be 63,453 bp longer than sweet cherry (P. avium). The P. salicina mitochondrial genome contained 37 protein-coding genes (PCGs), 3 ribosomal RNA (rRNA) genes, and 16 transfer RNA (tRNA) genes. Two plastid-derived tRNA were identified. We also found two short repeats that captured the nad3 and nad6 genes and resulted in two copies. In addition, nine pairs of repeat sequences were identified as being involved in the mediation of genome recombination. This is crucial for the formation of subgenomic configurations. To characterize RNA editing sites, transcriptome data were used, and we identified 480 RNA editing sites in protein-coding sequences. Among them, the initiation codon of the nad1 gene confirmed that an RNA editing event occurred, and the genomic encoded ACG was edited as AUG in the transcript. Combined with previous reports on the chloroplast genome, our data complemented our understanding of the last part of the organelle genome of plum, which will facilitate our understanding of the evolution of organelle genomes.

Efficient Compression and Indexing for Highly Repetitive DNA Sequence Collections

In this paper, we focus upon the important problem of indexing and searching highly repetitive DNA sequence collections. Given a collection G of t sequences S_i of length n each, we can represent G succinctly in 2nH_k(T) + O(n^' loglogn) + o(q n^') + o(tn) bits using O(t n² + q n^') time, where H_k(T) is the kth-order empirical entropy of the sequence T ∈ G that is used as the reference sequence, n^' is the total number of variations between T and the sequences in G, and q is a small fixed constant. We can restore any length len substring S[ sp, ..., sp + len-1] of S ∈ G in O(n_s^' + len(logn)² / loglogn) time and report all positions where P occurs in G in O(m ·t + occ ·t ·(logn)²/loglogn ) time. In addition, we propose a dynamic programming method to find the variations between T and the sequences in G in a space-efficient way, with which we can build succinct structures to enable efficient search. For highly repetitive sequences, experimental results on the tested data demonstrate that the proposed method has significant advantages in space usage and retrieval time over the current state-of-the-art methods. The source code is available online.

The Organization of Pericentromeric Heterochromatin in Polytene Chromosome 3 of the Drosophilamelanogaster Line with the Rif11; SuURES Su(var)3-906 Mutations Suppressing Underreplication

Although heterochromatin makes up 40% of the Drosophila melanogaster genome, its organization remains little explored, especially in polytene chromosomes, as it is virtually not represented in them due to underreplication. Two all-new approaches were used in this work: (i) with the use of a newly synthesized Drosophila line that carries three mutations, Rif11, SuURES and Su(var)3-906, suppressing the underreplication of heterochromatic regions, we obtained their fullest representation in polytene chromosomes and described their structure; (ii) 20 DNA fragments with known positions on the physical map as well as molecular genetic features of the genome (gene density, histone marks, heterochromatin proteins, origin recognition complex proteins, replication timing sites and satellite DNAs) were mapped in the newly polytenized heterochromatin using FISH and bioinformatics data. The borders of the heterochromatic regions and variations in their positions on arm 3L have been determined for the first time. The newly polytenized heterochromatic material exhibits two main types of morphology: a banding pattern (locations of genes and short satellites) and reticular chromatin (locations of large blocks of satellite DNA). The locations of the banding and reticular polytene heterochromatin was determined on the physical map.

DNA Satellites Are Transcribed as Part of the Non-Coding Genome in Eukaryotes and Bacteria

It has been shown in recent years that many repeated sequences in the genome are expressed as RNA transcripts, although the role of such RNAs is poorly understood. Some isolated and tandem repeats (satellites) have been found to be transcribed, such as mammalian Alu sequences and telomeric/centromeric satellites in different species. However, there is no detailed study on the eventual transcription of the interspersed satellites found in many species. Therefore, we decided to study for the first time the transcription of the abundant DNA satellites in the bacterium Bacillus coagulans and in the nematode Caenorhabditis elegans. We have updated the data for C. elegans satellites using the latest version of the genome. We analyzed the transcription of satellites in both species in available RNA-seq results and found that they are widely transcribed. Our demonstration that satellite RNAs are transcribed adds a new family of non-coding RNAs. This is a field that requires further investigation and will provide a deeper understanding of gene expression and control.

A sensitive repeat identification framework based on short and long reads

Numerous studies have shown that repetitive regions in genomes play indispensable roles in the evolution, inheritance and variation of living organisms. However, most existing methods cannot achieve satisfactory performance on identifying repeats in terms of both accuracy and size, since NGS reads are too short to identify long repeats whereas SMS (Single Molecule Sequencing) long reads are with high error rates. In this study, we present a novel identification framework, LongRepMarker, based on the global de novo assembly and k-mer based multiple sequence alignment for precisely marking long repeats in genomes. The major characteristics of LongRepMarker are as follows: (i) by introducing barcode linked reads and SMS long reads to assist the assembly of all short paired-end reads, it can identify the repeats to a greater extent; (ii) by finding the overlap sequences between assemblies or chomosomes, it locates the repeats faster and more accurately; (iii) by using the multi-alignment unique k-mers rather than the high frequency k-mers to identify repeats in overlap sequences, it can obtain the repeats more comprehensively and stably; (iv) by applying the parallel alignment model based on the multi-alignment unique k-mers, the efficiency of data processing can be greatly optimized and (v) by taking the corresponding identification strategies, structural variations that occur between repeats can be identified. Comprehensive experimental results show that LongRepMarker can achieve more satisfactory results than the existing de novo detection methods (https://github.com/BioinformaticsCSU/LongRepMarker).

Extensive genomic rearrangements mediated by repetitive sequences in plastomes of Medicago and its relatives

BACKGROUND

Although plastomes are highly conserved with respect to gene content and order in most photosynthetic angiosperms, extensive genomic rearrangements have been reported in Fabaceae, particularly within the inverted repeat lacking clade (IRLC) of Papilionoideae. Two hypotheses, i.e., the absence of the IR and the increased repeat content, have been proposed to affect the stability of plastomes. However, this is still unclear for the IRLC species. Here, we aimed to investigate the relationships between repeat content and the degree of genomic rearrangements in plastomes of Medicago and its relatives Trigonella and Melilotus, which are nested firmly within the IRLC.

RESULTS

We detected abundant repetitive elements and extensive genomic rearrangements in the 75 newly assembled plastomes of 20 species, including gene loss, intron loss and gain, pseudogenization, tRNA duplication, inversion, and a second independent IR gain (IR ~ 15 kb in Melilotus dentata) in addition to the previous first reported cases in Medicago minima. We also conducted comparative genomic analysis to evaluate plastome evolution. Our results indicated that the overall repeat content is positively correlated with the degree of genomic rearrangements. Some of the genomic rearrangements were found to be directly linked with repetitive sequences. Tandem repeated sequences have been detected in the three genes with accelerated substitution rates (i.e., accD, clpP, and ycf1) and their length variation could be explained by the insertions of tandem repeats. The repeat contents of the three localized hypermutation regions around these three genes with accelerated substitution rates are also significantly higher than that of the remaining plastome sequences.

CONCLUSIONS

Our results suggest that IR reemergence in the IRLC species does not ensure their plastome stability. Instead, repeat-mediated illegitimate recombination is the major mechanism leading to genome instability, a pattern in agreement with recent findings in other angiosperm lineages. The plastome data generated herein provide valuable genomic resources for further investigating the plastome evolution in legumes.

Repeat proliferation and partial endoreplication jointly shape the patterns of genome size evolution in orchids

Although the evolutionary drivers of genome size change are known, the general patterns and mechanisms of plant genome size evolution are yet to be established. Here we aim to assess the relative importance of proliferation of repetitive DNA, chromosomal variation (including polyploidy), and the type of endoreplication for genome size evolution of the Pleurothallidinae, the most species-rich orchid lineage. Phylogenetic relationships between 341 Pleurothallidinae representatives were refined using a target enrichment hybrid capture combined with high-throughput sequencing approach. Genome size and the type of endoreplication were assessed using flow cytometry supplemented with karyological analysis and low-coverage Illumina sequencing for repeatome analysis on a subset of samples. Data were analyzed using phylogeny-based models. Genome size diversity (0.2-5.1 Gbp) was mostly independent of profound chromosome count variation (2n = 12-90) but tightly linked with the overall content of repetitive DNA elements. Species with partial endoreplication (PE) had significantly greater genome sizes, and genomic repeat content was tightly correlated with the size of the non-endoreplicated part of the genome. In PE species, repetitive DNA is preferentially accumulated in the non-endoreplicated parts of their genomes. Our results demonstrate that proliferation of repetitive DNA elements and PE together shape the patterns of genome size diversity in orchids.

EPGA-SC : A Framework for de novo Assembly of Single-Cell Sequencing Reads

Assembling genomes from single-cell sequencing data is essential for single-cell studies. However, single-cell assemblies are challenging due to (i) the highly non-uniform read coverage and (ii) the elevated levels of sequencing errors and chimeric reads. Although several assemblers for single-cell data have been proposed in recent years, most of them fail to construct correct long contigs. In this study, we present a new framework called EPGA-SC for de novo assembly of single-cell sequencing reads. The EPGA assembler has designed strategies to solve the problems caused by sequencing errors, sequencing biases, and repetitive regions. However, the extremely unbalanced and richer error types prevent EPGA to achieve high performance in single-cell sequencing data. In this study, we designed EPGA-SC based on EPGA. The main innovations of EPGA-SC are as follows: (i) classifying reads to reduce the proportion of false reads; (ii) using multiple sets of high precision paired-end reads generated from the high precision assemblies produced by other assembler such as SPAdes to overcome the impact of sequencing biases and repetitive regions; and (iii) developing novel algorithms for removing chimeric errors and extending contigs. We test EPGA-SC with seven datasets. The experimental results show that EPGA-SC can generate better assemblies than most current tools in most time in term of MAX contig, N50, NG50, NA50, and NGA50.

Single Copy Oligonucleotide Fluorescence In Situ Hybridization Probe Design Platforms: Development, Application and Evaluation

Oligonucleotides fluorescence in situ hybridization (Oligo-FISH) is an emerging technology and is an important tool in research areas such as detection of chromosome variation, identification of allopolyploid, and deciphering of three-dimensional (3D) genome structures. Based on the demand for highly efficient oligo probes for oligo-FISH experiments, increasing numbers of tools have been developed for probe design in recent years. Obsolete oligonucleotide design tools have been adapted for oligo-FISH probe design because of their similar considerations. With the development of DNA sequencing and large-scale synthesis, novel tools have been designed to increase the specificity of designed oligo probes and enable genome-scale oligo probe design, which has greatly improved the application of single copy oligo-FISH. Despite this, few studies have introduced the development of the oligo-FISH probe design tools and their application in FISH experiments systematically. Besides, a comprehensive comparison and evaluation is lacking for the available tools. In this review, we provide an overview of the oligo-FISH probe design process, summarize the development and application of the available tools, evaluate several state-of-art tools, and eventually provide guidance for single copy oligo-FISH probe design.

Eusociality in snapping shrimps is associated with larger genomes and an accumulation of transposable elements

Despite progress uncovering the genomic underpinnings of sociality, much less is known about how social living affects the genome. In different insect lineages, for example, eusocial species show both positive and negative associations between genome size and structure, highlighting the dynamic nature of the genome. Here, we explore the relationship between sociality and genome architecture in Synalpheus snapping shrimps that exhibit multiple origins of eusociality and extreme interspecific variation in genome size. Our goal is to determine whether eusociality leads to an accumulation of repetitive elements and an increase in genome size, presumably due to reduced effective population sizes resulting from a reproductive division of labor, or whether an initial accumulation of repetitive elements leads to larger genomes and independently promotes the evolution of eusociality through adaptive evolution. Using phylogenetically informed analyses, we find that eusocial species have larger genomes with more transposable elements (TEs) and microsatellite repeats than noneusocial species. Interestingly, different TE subclasses contribute to the accumulation in different species. Phylogenetic path analysis testing alternative causal relationships between sociality and genome architecture is most consistent with the hypothesis that TEs modulate the relationship between sociality and genome architecture. Although eusociality appears to influence TE accumulation, ancestral state reconstruction suggests moderate TE abundances in ancestral species could have fueled the initial transitions to eusociality. Ultimately, we highlight a complex and dynamic relationship between genome and social evolution, demonstrating that sociality can influence the evolution of the genome, likely through changes in demography related to patterns of reproductive skew.

Crystal structure of parallel G-quadruplex formed by the two-repeat ALS- and FTD-related GGGGCC sequence

The hexanucleotide repeat expansion, GGGGCC (G4C2), within the first intron of the C9orf72 gene is known to be the most common genetic cause of both amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). The G4C2 repeat expansions, either DNA or RNA, are able to form G-quadruplexes which induce toxicity leading to ALS/FTD. Herein, we report a novel crystal structure of d(G4C2)2 that self-associates to form an eight-layer parallel tetrameric G-quadruplex. Two d(G4C2)2 associate together as a parallel dimeric G-quadruplex which folds into a tetramer via 5'-to-5' arrangements. Each dimer consists of four G-tetrads connected by two CC propeller loops. Especially, the 3'-end cytosines protrude out and form C·C+•C·C+/ C·C•C·C+ quadruple base pair or C•C·C+ triple base pair stacking on the dimeric block. Our work sheds light on the G-quadruplexes adopted by d(G4C2) and yields the invaluable structural details for the development of small molecules to tackle neurodegenerative diseases, ALS and FTD.

Characterization of structural variation in Tibetans reveals new evidence of high-altitude adaptation and introgression

BACKGROUND

Structural variation (SV) acts as an essential mutational force shaping the evolution and function of the human genome. However, few studies have examined the role of SVs in high-altitude adaptation and little is known of adaptive introgressed SVs in Tibetans so far.

RESULTS

Here, we generate a comprehensive catalog of SVs in a Chinese Tibetan (n = 15) and Han (n = 10) population using nanopore sequencing technology. Among a total of 38,216 unique SVs in the catalog, 27% are sequence-resolved for the first time. We systematically assess the distribution of these SVs across repeat sequences and functional genomic regions. Through genotyping in additional 276 genomes, we identify 69 Tibetan-Han stratified SVs and 80 candidate adaptive genes. We also discover a few adaptive introgressed SV candidates and provide evidence for a deletion of 335 base pairs at 1p36.32.

CONCLUSIONS

Overall, our results highlight the important role of SVs in the evolutionary processes of Tibetans' adaptation to the Qinghai-Tibet Plateau and provide a valuable resource for future high-altitude adaptation studies.

Toxic Y chromosome: Increased repeat expression and age-associated heterochromatin loss in male Drosophila with a young Y chromosome

Sex-specific differences in lifespan are prevalent across the tree of life and influenced by heteromorphic sex chromosomes. In species with XY sex chromosomes, females often outlive males. Males and females can differ in their overall repeat content due to the repetitive Y chromosome, and repeats on the Y might lower survival of the heterogametic sex (toxic Y effect). Here, we take advantage of the well-assembled young Y chromosome of Drosophila miranda to study the sex-specific dynamics of chromatin structure and repeat expression during aging in male and female flies. Male D. miranda have about twice as much repetitive DNA compared to females, and live shorter than females. Heterochromatin is crucial for silencing of repetitive elements, yet old D. miranda flies lose H3K9me3 modifications in their pericentromere, with heterochromatin loss being more severe during aging in males than females. Satellite DNA becomes de-repressed more rapidly in old vs. young male flies relative to females. In contrast to what is observed in D. melanogaster, we find that transposable elements (TEs) are expressed at higher levels in male D. miranda throughout their life. We show that epigenetic silencing via heterochromatin formation is ineffective on the TE-rich neo-Y chromosome, presumably due to active transcription of a large number of neo-Y linked genes, resulting in up-regulation of Y-linked TEs already in young males. This is consistent with an interaction between the evolutionary age of the Y chromosome and the genomic effects of aging. Our data support growing evidence that “toxic Y chromosomes” can diminish male fitness and a reduction in heterochromatin can contribute to sex-specific aging.

Y chromosomes can be toxic. The Y chromosome of many species contains a large number of transposable elements (TEs), which are transcriptionally constrained by repressive chromatin marks. When relieved of these epigenetic modifications, many TEs can readily move from one genomic location to another. We show that TEs located on the Y chromosome are less effectively silenced in male Drosophila, and the toxic Y effect appears more pronounced in a species that contains a larger Y chromosome with more repeats and more actively transcribed genes. Our data demonstrate that repeat-rich Y chromosomes are a genomic liability for males.

BLM has Contrary Effects on Repeat-Mediated Deletions, based on the Distance of DNA DSBs to a Repeat and Repeat Divergence

Repeat-mediated deletions (RMDs) often involve repetitive elements (e.g., short interspersed elements) with sequence divergence that is separated by several kilobase pairs (kbps). We have examined RMDs induced by DNA double-strand breaks (DSBs) under varying conditions of repeat sequence divergence (identical versus 1% and 3% divergent) and DSB/repeat distance (16 bp–28.4 kbp). We find that the BLM helicase promotes RMDs with long DSB/repeat distances (e.g., 28.4 kbp), which is consistent with a role in extensive DSB end resection, because the resection nucleases EXO1 and DNA2 affect RMDs similarly to BLM. In contrast, BLM suppresses RMDs with sequence divergence and intermediate (e.g., 3.3 kbp) DSB/repeat distances, which supports a role in heteroduplex rejection. The role of BLM in heteroduplex rejection is not epistatic with MSH2 and is independent of the annealing factor RAD52. Accordingly, the role of BLM on RMDs is substantially affected by DSB/repeat distance and repeat sequence divergence.

Mendez-Dorantes et al. identify the BLM helicase as a key regulator of repeat-mediated deletions (RMDs). BLM, EXO1, and DNA2 mediate RMDs with remarkably long DNA break/repeat distances. BLM suppresses RMDs with sequence divergence that is optimal with a long non-homologous tail and is independent of MSH2 and RAD52.

The MTNR1A mRNA is stabilized by the cytoplasmic hnRNPL in renal tubular cells

The downregulation of melatonin receptor 1A (MTNR1A) is associated with a range of pathological conditions, including membranous nephropathy. Knowledge of the mechanism underlying MTNR1A expression has been limited to the transcriptional regulation level. Here, RNA interference screening in human kidney cells revealed that heterogeneous nuclear ribonucleoprotein L (hnRNPL) upregulated MTNR1A RNA post-transcriptionally. hnRNPL knockdown or overexpression led to increased or decreased levels of cyclic adenosine monophosphate-responsive element-binding protein phosphorylation, respectively. Molecular studies showed that cytoplasmic hnRNPL exerts a stabilizing effect on the MTNR1A transcript through CA-repeat elements in its coding region. Further studies revealed that the interaction between hnRNPL and MTNR1A serves to protect MNTR1A RNA degradation by the exosome component 10 protein. MTNR1A, but not hnRNPL, displays a diurnal rhythm in mouse kidneys. Enhanced levels of MTNR1A recorded at midnight correlated with robust binding activity between cytoplasmic hnRNPL and the MTNR1A transcript. Both hnRNPL and MTNR1A were decreased in the cytoplasm of tubular epithelial cells from experimental membranous nephropathy kidneys, supporting their clinical relevance. Collectively, our data identified cytoplasmic hnRNPL as a novel player in the upregulation of MTNR1A expression in renal tubular epithelial cells, and as a potential therapeutic target.

Universal fluorescence in situ hybridization (FISH) protocol for mapping repetitive DNAs in insects and other arthropods

Repetitive DNAs comprise large portion of eukaryote genomes. In genome projects, the assembly of repetitive DNAs is challenging due to the similarity between repeats, which generate ambiguities for alignment. Fluorescence in situ hybridization (FISH) is a powerful technique for the physical mapping of various sequences on chromosomes. This technique is thus very helpful in chromosome-based genome assemblies, providing information on the fine architecture of genomes and their evolution. However, various protocols are currently used for FISH mapping, most of which are relatively laborious and expensive, or work properly only with a specific type of probes or sequences, and there is a need for a universal and affordable FISH protocol. Here we tested a FISH protocol for mapping of different DNA repeats, such as multigene families (rDNAs, U snDNAs, histone genes), satellite DNAs, microsatellites, transposable elements, DOP-PCR products, and telomeric motif (TTAGG)_n, on the chromosomes of various insects and other arthropods. Different cell types and stages obtained from diverse tissues were used. The FISH procedure proved high quality and reliable results in all experiments performed. We obtained data on the chromosomal distribution of DNA repeats in representatives of insects and other arthropods. Thus, our results allow us to conclude that the protocol is universal and requires only time adjustment for chromosome/DNA denaturation. The use of this FISH protocol will facilitate studies focused on understanding the evolution and role of repetitive DNA in arthropod genomes.

Physical mapping of repetitive oligonucleotides facilitates the establishment of a genome map-based karyotype to identify chromosomal variations in peanut

BACKGROUND

Chromosomal variants play important roles in crop breeding and genetic research. The development of single-stranded oligonucleotide (oligo) probes simplifies the process of fluorescence in situ hybridization (FISH) and facilitates chromosomal identification in many species. Genome sequencing provides rich resources for the development of oligo probes. However, little progress has been made in peanut due to the lack of efficient chromosomal markers. Until now, the identification of chromosomal variants in peanut has remained a challenge.

RESULTS

A total of 114 new oligo probes were developed based on the genome-wide tandem repeats (TRs) identified from the reference sequences of the peanut variety Tifrunner (AABB, 2n = 4x = 40) and the diploid species Arachis ipaensis (BB, 2n = 2x = 20). These oligo probes were classified into 28 types based on their positions and overlapping signals in chromosomes. For each type, a representative oligo was selected and modified with green fluorescein 6-carboxyfluorescein (FAM) or red fluorescein 6-carboxytetramethylrhodamine (TAMRA). Two cocktails, Multiplex #3 and Multiplex #4, were developed by pooling the fluorophore conjugated probes. Multiplex #3 included FAM-modified oligo TIF-439, oligo TIF-185-1, oligo TIF-134-3 and oligo TIF-165. Multiplex #4 included TAMRA-modified oligo Ipa-1162, oligo Ipa-1137, oligo DP-1 and oligo DP-5. Each cocktail enabled the establishment of a genome map-based karyotype after sequential FISH/genomic in situ hybridization (GISH) and in silico mapping. Furthermore, we identified 14 chromosomal variants of the peanut induced by radiation exposure. A total of 28 representative probes were further chromosomally mapped onto the new karyotype. Among the probes, eight were mapped in the secondary constrictions, intercalary and terminal regions; four were B genome-specific; one was chromosome-specific; and the remaining 15 were extensively mapped in the pericentric regions of the chromosomes.

CONCLUSIONS

The development of new oligo probes provides an effective set of tools which can be used to distinguish the various chromosomes of the peanut. Physical mapping by FISH reveals the genomic organization of repetitive oligos in peanut chromosomes. A genome map-based karyotype was established and used for the identification of chromosome variations in peanut following comparisons with their reference sequence positions.

Repetitive Elements in Humans

Repetitive DNA in humans is still widely considered to be meaningless, and variations within this part of the genome are generally considered to be harmless to the carrier. In contrast, for euchromatic variation, one becomes more careful in classifying inter-individual differences as meaningless and rather tends to see them as possible influencers of the so-called 'genetic background', being able to at least potentially influence disease susceptibilities. Here, the known 'bad boys' among repetitive DNAs are reviewed. Variable numbers of tandem repeats (VNTRs = micro- and minisatellites), small-scale repetitive elements (SSREs) and even chromosomal heteromorphisms (CHs) may therefore have direct or indirect influences on human diseases and susceptibilities. Summarizing this specific aspect here for the first time should contribute to stimulating more research on human repetitive DNA. It should also become clear that these kinds of studies must be done at all available levels of resolution, i.e., from the base pair to chromosomal level and, importantly, the epigenetic level, as well.

IUPACpal: efficient identification of inverted repeats in IUPAC-encoded DNA sequences

BACKGROUND

An inverted repeat is a DNA sequence followed downstream by its reverse complement, potentially with a gap in the centre. Inverted repeats are found in both prokaryotic and eukaryotic genomes and they have been linked with countless possible functions. Many international consortia provide a comprehensive description of common genetic variation making alternative sequence representations, such as IUPAC encoding, necessary for leveraging the full potential of such broad variation datasets.

RESULTS

We present IUPACPAL, an exact tool for efficient identification of inverted repeats in IUPAC-encoded DNA sequences allowing also for potential mismatches and gaps in the inverted repeats.

CONCLUSION

Within the parameters that were tested, our experimental results show that IUPACPAL compares favourably to a similar application packaged with EMBOSS. We show that IUPACPAL identifies many previously unidentified inverted repeats when compared with EMBOSS, and that this is also performed with orders of magnitude improved speed.

Customized optical mapping by CRISPR-Cas9 mediated DNA labeling with multiple sgRNAs

Whole-genome mapping technologies have been developed as a complementary tool to provide scaffolds for genome assembly and structural variation analysis (1,2). We recently introduced a novel DNA labeling strategy based on a CRISPR-Cas9 genome editing system, which can target any 20bp sequences. The labeling strategy is specifically useful in targeting repetitive sequences, and sequences not accessible to other labeling methods. In this report, we present customized mapping strategies that extend the applications of CRISPR-Cas9 DNA labeling. We first design a CRISPR-Cas9 labeling strategy to interrogate and differentiate the single allele differences in NGG protospacer adjacent motifs (PAM sequence). Combined with sequence motif labeling, we can pinpoint the single-base differences in highly conserved sequences. In the second strategy, we design mapping patterns across a genome by selecting sets of specific single-guide RNAs (sgRNAs) for labeling multiple loci of a genomic region or a whole genome. By developing and optimizing a single tube synthesis of multiple sgRNAs, we demonstrate the utility of CRISPR-Cas9 mapping with 162 sgRNAs targeting the 2Mb Haemophilus influenzae chromosome. These CRISPR-Cas9 mapping approaches could be particularly useful for applications in defining long-distance haplotypes and pinpointing the breakpoints in large structural variants in complex genomes and microbial mixtures.

Visualization of lncRNA and mRNA Structure Models Within the Integrative Genomics Viewer

Every class of RNA forms base-paired structures that impact biological functions. Chemical probing of RNA structure, especially with the advent of strategies such as SHAPE-MaP, vastly expands the scale and quantitative accuracy over which RNA structure can be examined. These methods have enabled large-scale structural studies of mRNAs and lncRNAs, but the length and complexity of these RNAs makes interpretation of the data challenging. We have created modules available through the open-source Integrative Genomics Viewer (IGV) for straightforward visualization of RNA structures along with complementary experimental data. Here we present detailed and stepwise strategies for exploring and visualizing complex RNA structures in IGV. Individuals can use these instructions and supplied sample data to become adept at using IGV to visualize RNA structure models in conjunction with useful allied information.

Developmental remodelling of non-CG methylation at satellite DNA repeats

In vertebrates, DNA methylation predominantly occurs at CG dinucleotides however, widespread non-CG methylation (mCH) has been reported in mammalian embryonic stem cells and in the brain. In mammals, mCH is found at CAC trinucleotides in the nervous system, where it is associated with transcriptional repression, and at CAG trinucleotides in embryonic stem cells, where it positively correlates with transcription. Moreover, CAC methylation appears to be a conserved feature of adult vertebrate brains. Unlike any of those methylation signatures, here we describe a novel form of mCH that occurs in the TGCT context within zebrafish mosaic satellite repeats. TGCT methylation is inherited from both male and female gametes, remodelled during mid-blastula transition, and re-established during gastrulation in all embryonic layers. Moreover, we identify DNA methyltransferase 3ba (Dnmt3ba) as the primary enzyme responsible for the deposition of this mCH mark. Finally, we observe that TGCT-methylated repeats are specifically associated with H3K9me3-marked heterochromatin suggestive of a functional interplay between these two gene-regulatory marks. Altogether, this work provides insight into a novel form of vertebrate mCH and highlights the substrate diversity of vertebrate DNA methyltransferases.

The genomic structure of the highly-conserved dmrt1 gene in Solea senegalensis (Kaup, 1868) shows an unexpected intragenic duplication

Knowing the factors responsible for sex determination in a species has significant theoretical and practical implications; the dmrt1 gene (Doublesex and Mab-3 (DM)-related Transcription factor 1) plays this role in diverse animal species. Solea senegalensis is a commercially important flat fish in which females grow 30% faster than males. It has 2n = 42 chromosomes and an XX / XY chromosome system for sex determination, without heteromorph chromosomes but with sex proto-chromosome. In the present study, we are providing the genomic structure and nucleotide sequence of dmrt1 gene obtained from cDNA from male and female adult gonads. A cDNA of 2027 containing an open-reading frame (ORF) of 1206 bp and encoding a 402 aa protein it is described for dmrt1 gene of S. senegalensis. Multiple mRNA isoforms indicating a high variable system of alternative splicing in the expression of dmrt1 of the sole in gonads were studied. None isoforms could be related to sex of individuals. The genomic structure of the dmrt1 of S. senegalensis showed a gene of 31400 bp composed of 7 exons and 6 introns. It contains an unexpected duplication of more than 10399 bp, involving part of the exon I, exons II and III and a SINE element found in the sequence that it is proposed as responsible for the duplication. A mature miRNA of 21 bp in length was localized at 336 bp from exon V. Protein-protein interacting networks of the dmrt1 gene showed matches with dmrt1 protein from Cynoglossus semilaevis and a protein interaction network with 11 nodes (dmrt1 plus 10 other proteins). The phylogenetic relationship of the dmrt1 gene in S. senegalensis is consistent with the evolutionary position of its species. The molecular characterization of this gene will enhance its functional analysis and the understanding of sex differentiation in Solea senegalensis and other flatfish.

Repeat-sequence turnover shifts fundamentally in species with large genomes

Given the 2,400-fold range of genome sizes (0.06-148.9 Gbp (gigabase pair)) of seed plants (angiosperms and gymnosperms) with a broadly similar gene content (amounting to approximately 0.03 Gbp), the repeat-sequence content of the genome might be expected to increase with genome size, resulting in the largest genomes consisting almost entirely of repetitive sequences. Here we test this prediction, using the same bioinformatic approach for 101 species to ensure consistency in what constitutes a repeat. We reveal a fundamental change in repeat turnover in genomes above around 10 Gbp, such that species with the largest genomes are only about 55% repetitive. Given that genome size influences many plant traits, habits and life strategies, this fundamental shift in repeat dynamics is likely to affect the evolutionary trajectory of species lineages.