Ribosomal DNA (rDNA) repeats are situated in the nucleolus organizer regions (NOR) of chromosomes and transcribed into rRNA for ribosome biogenesis. Thus, they are an essential component of eukaryotic genomes. rDNA repeat units consist of rRNA gene clusters that are transcribed into single pre-rRNA molecules, each separated by intergenic spacers (IGS) that contain regulatory elements for rRNA gene cluster transcription. Because of their high repeat content, rDNA sequences are usually absent from genome assemblies. In this work, we used the long-read sequencing technology to describe the chicken IGS and fill the knowledge gap on rDNA sequences of one of the key domesticated animals.

Factors affecting the organization and spacing of functionally unrelated genes in metazoan genomes are not well understood. Because of the vast size of a typical metazoan genome compared to known regulatory and protein-coding regions, functional DNA is generally considered to have a negligible impact on gene spacing and genome organization. In particular, it has been impossible to estimate the global impact, if any, of regulatory elements on genome architecture.

The geminivirus and nanovirus families of DNA plant viruses have proved to be a fertile source of viral genomic sequences, clearly demonstrated by the large number of sequence entries within public DNA sequence databases. Due to considerable conservation in genome organization, these viruses contain easily identifiable intergenic regions that have been found to contain multiple DNA sequence elements important to viral replication and gene regulation. As a first step in a broad screen of geminivirus and nanovirus intergenic sequences for DNA segments important in controlling viral gene expression, we have 'mined' a large set of viral intergenic regions for transcriptional enhancers. Viral sequences that are found to act as enhancers of transcription in plants are likely to contribute to viral gene activity during infection.

Meiotic recombination is initiated by DNA double-strand breaks (DSBs) made by Spo11 (Rec12 in fission yeast), which becomes covalently linked to the DSB ends. Like recombination events, DSBs occur at hotspots in the genome, but the genetic factors responsible for most hotspots have remained elusive. Here we describe in fission yeast the genome-wide distribution of meiosis-specific Rec12-DNA linkages, which closely parallel DSBs measured by conventional Southern blot hybridization. Prominent DSB hotspots are located approximately 65 kb apart, separated by intervals with little or no detectable breakage. Most hotspots lie within exceptionally large intergenic regions. Thus, the chromosomal architecture responsible for hotspots in fission yeast is markedly different from that of budding yeast, in which DSB hotspots are much more closely spaced and, in many regions of the genome, occur at each promoter. Our analysis in fission yeast reveals a clearly identifiable chromosomal feature that can predict the majority of recombination hotspots across a whole genome and provides a basis for searching for the chromosomal features that dictate hotspots of meiotic recombination in other organisms, including humans.

Carboxy-terminal binding protein (CtBP) is a well-known corepressor of several DNA binding transcription factors in Drosophila as well as in mammals. CtBP is implicated in Polycomb Group (PcG) complex-mediated transcriptional repression because it can bind to some PcG proteins, and mutation of the ctbp gene in flies results in lost PcG protein recruitment to Polycomb Response Elements (PREs) and lost PcG repression. However, the mechanism of reduced PcG DNA binding in CtBP mutant backgrounds is unknown. We show here that in a Drosophila CtBP mutant background, intergenic transcripts are induced across several PRE sequences and this corresponds to reduced DNA binding by PcG proteins Pleiohomeotic (PHO) and Polycomb (Pc), and reduced trimethylation of histone H3 on lysine 27, a hallmark of PcG repression. Restoration of CtBP levels by expression of a CtBP transgene results in repression of intergenic transcripts, restored PcG binding, and elevated trimethylation of H3 on lysine 27. Our results support a model in which CtBP regulates expression of intergenic transcripts that controls DNA binding by PcG proteins and subsequent histone modifications and transcriptional activity.

Deep sequencing of mammalian DNA methylomes has uncovered a previously unpredicted number of discrete hypomethylated regions in intergenic space (iHMRs). Here, we combined whole-genome bisulfite sequencing data with extensive gene expression and chromatin-state data to define functional classes of iHMRs, and to reconstruct the dynamics of their establishment in a developmental setting. Comparing HMR profiles in embryonic stem and primary blood cells, we show that iHMRs mark an exclusive subset of active DNase hypersensitive sites (DHS), and that both developmentally constitutive and cell-type-specific iHMRs display chromatin states typical of distinct regulatory elements. We also observe that iHMR changes are more predictive of nearby gene activity than the promoter HMR itself, and that expression of noncoding RNAs within the iHMR accompanies full activation and complete demethylation of mature B cell enhancers. Conserved sequence features corresponding to iHMR transcript start sites, including a discernible TATA motif, suggest a conserved, functional role for transcription in these regions. Similarly, we explored both primate-specific and human population variation at iHMRs, finding that while enhancer iHMRs are more variable in sequence and methylation status than any other functional class, conservation of the TATA box is highly predictive of iHMR maintenance, reflecting the impact of sequence plasticity and transcriptional signals on iHMR establishment. Overall, our analysis allowed us to construct a three-step timeline in which (1) intergenic DHS are pre-established in the stem cell, (2) partial demethylation of blood-specific intergenic DHSs occurs in blood progenitors, and (3) complete iHMR formation and transcription coincide with enhancer activation in lymphoid-specified cells.

Enzymes that catalyse CpG methylation in DNA, including the DNA methyltransferases 1 (DNMT1), 3A (DNMT3A) and 3B (DNMT3B), are indispensable for mammalian tissue development and homeostasis1-4. They are also implicated in human developmental disorders and cancers5-8, supporting the critical role of DNA methylation in the specification and maintenance of cell fate. Previous studies have suggested that post-translational modifications of histones are involved in specifying patterns of DNA methyltransferase localization and DNA methylation at promoters and actively transcribed gene bodies9-11. However, the mechanisms that control the establishment and maintenance of intergenic DNA methylation remain poorly understood. Tatton-Brown-Rahman syndrome (TBRS) is a childhood overgrowth disorder that is defined by germline mutations in DNMT3A. TBRS shares clinical features with Sotos syndrome (which is caused by haploinsufficiency of NSD1, a histone methyltransferase that catalyses the dimethylation of histone H3 at K36 (H3K36me2)8,12,13), which suggests that there is a mechanistic link between these two diseases. Here we report that NSD1-mediated H3K36me2 is required for the recruitment of DNMT3A and maintenance of DNA methylation at intergenic regions. Genome-wide analysis shows that the binding and activity of DNMT3A colocalize with H3K36me2 at non-coding regions of euchromatin. Genetic ablation of Nsd1 and its paralogue Nsd2 in mouse cells results in a redistribution of DNMT3A to H3K36me3-modified gene bodies and a reduction in the methylation of intergenic DNA. Blood samples from patients with Sotos syndrome and NSD1-mutant tumours also exhibit hypomethylation of intergenic DNA. The PWWP domain of DNMT3A shows dual recognition of H3K36me2 and H3K36me3 in vitro, with a higher binding affinity towards H3K36me2 that is abrogated by TBRS-derived missense mutations. Together, our study reveals a trans-chromatin regulatory pathway that connects aberrant intergenic CpG methylation to human neoplastic and developmental overgrowth.

Dunaliella salina Teodoresco, a unicellular, halophilic green alga belonging to the Chlorophyceae, is among the most industrially important microalgae. This is because D. salina can produce massive amounts of beta-carotene, which can be collected for commercial purposes, and because of its potential as a feedstock for biofuels production. Although the biochemistry and physiology of D. salina have been studied in great detail, virtually nothing is known about the genomes it carries, especially those within its mitochondrion and plastid. This study presents the complete mitochondrial and plastid genome sequences of D. salina and compares them with those of the model green algae Chlamydomonas reinhardtii and Volvox carteri.

Long intergenic noncoding RNAs (lincRNAs) are one of the major unexplored components of genomes. Here we re-analyzed a published methylated DNA immunoprecipitation sequencing (MeDIP-seq) dataset to characterize the DNA methylation pattern of pig lincRNA genes in adipose and muscle tissues. Our study showed that the methylation level of lincRNA genes was higher than that of mRNA genes, with similar trends observed in comparisons of the promoter, exon or intron regions. Different methylation pattern were observed across the transcription start sites (TSS) of lincRNA and protein-coding genes. Furthermore, an overlap was observed between many lincRNA genes and differentially methylated regions (DMRs) identified among different breeds of pigs, which show different fat contents, sexes and anatomic locations of tissues. We identify a lincRNA gene, linc-sscg3623, that displayed differential methylation levels in backfat between Min and Large White pigs at 60 and 120 days of age. We found that a demethylation process occurred between days 150 and 180 in the Min and Large White pigs, which was followed by remethylation between days 180 and 210. These results contribute to our understanding of the domestication of domestic animals and identify lincRNA genes involved in adipogenesis and muscle development.

Hypomethylation of LINE-1 repeats in cancer has been proposed as the main mechanism behind their activation; this assumption, however, was based on findings from early studies that were biased toward young and transpositionally active elements. Here, we investigate the relationship between methylation of 2 intergenic, transpositionally inactive LINE-1 elements and expression of the LINE-1 chimeric transcript (LCT) 13 and LCT14 driven by their antisense promoters (L1-ASP). Our data from DNA modification, expression, and 5'RACE analyses suggest that colorectal cancer methylation in the regions analyzed is not always associated with LCT repression. Consistent with this, in HCT116 colorectal cancer cells lacking DNA methyltransferases DNMT1 or DNMT3B, LCT13 expression decreases, while cells lacking both DNMTs or treated with the DNMT inhibitor 5-azacytidine (5-aza) show no change in LCT13 expression. Interestingly, levels of the H4K20me3 histone modification are inversely associated with LCT13 and LCT14 expression. Moreover, at these LINE-1s, H4K20me3 levels rather than DNA methylation seem to be good predictor of their sensitivity to 5-aza treatment. Therefore, by studying individual LINE-1 promoters we have shown that in some cases these promoters can be active without losing methylation; in addition, we provide evidence that other factors (e.g., H4K20me3 levels) play prominent roles in their regulation.

The miR-200 family members have been implicated in stress responses and ovarian tumorigenesis. Here, we find that miR-200c/141 transcription is intimately linked to the transcription of the proximal upstream gene PTPN6 (SHP1) in all physiological conditions tested. PTPN6 and miR-200c/141 are transcriptionally co-regulated by two complementary mechanisms. First, a bypass of the regular PTPN6 polyadenylation signal allows the transcription of the downstream miR-200c/141. Second, the promoters of the PTPN6 and miR-200c/141 transcription units physically interact through a 3-dimensional DNA loop and exhibit similar epigenetic regulation. Our findings highlight that transcription of intergenic miRNAs is a novel outcome of transcriptional read-through and reveal a yet unexplored type of DNA loop associating two closely located promoters. These mechanisms have significant relevance in ovarian cancers and stress response, pathophysiological conditions in which miR-200c/141 exert key functions.

DNA from a single bacterial artificial chromosome clone was used to sequence the mouse ribosomal DNA intergenic spacer from the 3' end of the 45S pre-RNA to the spacer promoter (Accession No. AF441733). This made possible the assembly of a complete mouse ribosomal DNA repeat unit (45309 bp long, TPA Accession No. BK000964). Analysis of the intergenic spacer (IGS) showed a high density of simple sequence repeats and transposable elements. The IGS contains two long sequence blocks, which are repeated tandemly. Some of the sequences in these blocks are also present in other parts of the IGS. A difference in the mutation rate along the mouse IGS was observed. The significance of sequence motifs in the IGS for transcription enhancement, transcription termination, origin of replication, and nucleolar organization is discussed.

Mammalian genomes comprise largely intergenic noncoding DNA with numerous cis-regulatory elements. Whether and how the size of intergenic DNA affects gene expression in a tissue-specific manner remain unknown. Here we show that genes with extended intergenic regions are preferentially expressed in neural tissues but repressed in other tissues in mice and humans. Extended intergenic regions contain twice as many active enhancers in neural tissues compared to other tissues. Neural genes with extended intergenic regions are globally co-expressed with neighboring neural genes controlled by distinct enhancers in the shared intergenic regions. Moreover, generic neural genes expressed in multiple tissues have significantly longer intergenic regions than neural genes expressed in fewer tissues. The intergenic regions of the generic neural genes have many tissue-specific active enhancers containing distinct transcription factor binding sites specific to each neural tissue. We also show that genes with extended intergenic regions are enriched for neural genes only in vertebrates. The expansion of intergenic regions may reflect the regulatory complexity of tissue-type-specific gene expression in the nervous system.

The intergenic spacer (IGS) of rDNA is frequently built of long blocks of tandem repeats. To estimate the intragenomic variability of such knotty regions, we employed PacBio sequencing of the Cucurbita moschata genome, in which thousands of rDNA copies are distributed across a number of loci. The rRNA coding regions are highly conserved, indicating intensive interlocus homogenization and/or high selection pressure. However, the IGS exhibits high intragenomic structural diversity. Two repeated blocks, R1 (300-1250 bp) and R2 (290-643 bp), account for most of the IGS variation. They exhibit minisatellite-like features built of multiple periodically spaced short GC-rich sequence motifs with the potential to adopt non-canonical DNA conformations, G-quadruplex-folded and left-handed Z-DNA. The mutual arrangement of these motifs can be used to classify IGS variants into five structural families. Subtle polymorphisms exist within each family due to a variable number of repeats, suggesting the coexistence of an enormous number of IGS variants. The substantial length and structural heterogeneity of IGS minisatellites suggests that the tempo of their divergence exceeds the tempo of the homogenization of rDNA arrays. As frequently occurring among plants, we hypothesize that their instability may influence transcription regulation and/or destabilize rDNA units, possibly spreading them across the genome.

Concerted evolution refers to the pattern in which copies of multigene families show high intraspecific sequence homogeneity but high interspecific sequence diversity. Sequence homogeneity of these copies depends on relative rates of mutation and recombination, including gene conversion and unequal crossing over, between misaligned copies. The internally repetitive intergenic spacer (IGS) is located between the genes for the 28S and 18S ribosomal RNAs. To identify patterns of recombination and/or homogenization within IGS repeat arrays, and to identify regions of the IGS that are under functional constraint, we analyzed 13 complete IGS sequences from 10 individuals representing four species in the Daphnia pulex complex.

Two DNA-based methods were compared for the ability to assign serotype to 139 isolates of Salmonella enterica ssp. I. Intergenic sequence ribotyping (ISR) evaluated single nucleotide polymorphisms occurring in a 5S ribosomal gene region and flanking sequences bordering the gene dkgB. A DNA microarray hybridization method that assessed the presence and the absence of sets of genes was the second method. Serotype was assigned for 128 (92.1%) of submissions by the two DNA methods. ISR detected mixtures of serotypes within single colonies and it cost substantially less than Kauffmann-White serotyping and DNA microarray hybridization. Decreasing the cost of serotyping S. enterica while maintaining reliability may encourage routine testing and research.

Mitochondrial DNA variations of Peruvian honey bee populations were surveyed by using the tRNAleu-cox2 intergenic region. Only two studies have characterized these populations, indicating the presence of Africanized honey bee colonies in different regions of Peru and varied levels of Africanization, but the current status of its genetic diversity is unknown. A total of 512 honey bee colonies were sampled from three regions to characterize them. Our results revealed the presence of European and African haplotypes: the African haplotypes identified belong to sub-lineage AI (13) and sub-lineage AIII (03), and the European haplotypes to lineages C (06) and M (02). Of 24 haplotypes identified, 15 new sequences are reported here (11 sub-lineage AI, 2 sub-lineage AIII, and 2 lineage M). Peruvian honey bee populations presented a higher proportion from African than European haplotypes. High proportions of African haplotype were reported for Piura and Junín, unlike Lima, which showed more European haplotypes from lineage C. Few colonies belonging to lineage M would represent accidental purchase or traces of the introduction into Peru in the 19th century.

DNA barcoding technology, which uses a short piece of DNA sequence to identify species, has wide ranges of applications. Until today, a universal DNA barcode marker for plants remains elusive. The rbcL and matK regions have been proposed as the "core barcode" for plants and the ITS2 and psbA-trnH intergenic spacer (PTIGS) regions were later added as supplemental barcodes. The use of PTIGS region as a supplemental barcode has been limited by the lack of computational tools that can handle significant insertions and deletions in the PTIGS sequences. Here, we compared the most commonly used alignment-based and alignment-free methods and developed a web server to allow the biologists to carry out PTIGS-based DNA barcoding analyses.

Tissue-specific enhancers are critical for gene regulation. In this study, we help elucidate the contribution of muscle-associated differential DNA methylation to the enhancer activity of highly muscle-specific genes. By bioinformatic analysis of 44 muscle-associated genes, we show that preferential gene expression in skeletal muscle (SkM) correlates with SkM-specific intragenic and intergenic enhancer chromatin and overlapping foci of DNA hypomethylation. Some genes, e.g., CASQ1 and FBXO32, displayed broad regions of both SkM- and heart-specific enhancer chromatin but exhibited focal SkM-specific DNA hypomethylation. Half of the genes had SkM-specific super-enhancers. In contrast to simple enhancer/gene-expression correlations, a super-enhancer was associated with the myogenic MYOD1 gene in both SkM and myoblasts even though SkM has < 1 percent as much MYOD1 expression. Local chromatin differences in this super-enhancer probably contribute to the SkM/myoblast differential expression. Transfection assays confirmed the tissue-specificity of the 0.3-kb core enhancer within MYOD1's super-enhancer and demonstrated its repression by methylation of its three CG dinucleotides. Our study suggests that DNA hypomethylation increases enhancer tissue-specificity and that SkM super-enhancers sometimes are poised for physiologically important, rapid up-regulation.

Despite growing consensus that long intergenic non-coding ribonucleic acids (lincRNAs) are modulators of cancer, the knowledge about the deoxyribonucleic acid (DNA) methylation patterns of lincRNAs in cancers remains limited. In this study, we constructed DNA methylation profiles for 4629 tumors and 705 normal tissue samples from 20 different types of human cancer by reannotating data of DNA methylation arrays. We found that lincRNAs had different promoter methylation patterns in cancers. We classified 2461 lincRNAs into two categories and three subcategories, according to their promoter methylation patterns in tumors. LincRNAs with resistant methylation patterns in tumors had conserved transcriptional regulation regions and were ubiquitously expressed across normal tissues. By integrating cancer subtype data and patient clinical information, we identified lincRNAs with promoter methylation patterns that were associated with cancer status, subtype or prognosis for several cancers. Network analysis of aberrantly methylated lincRNAs in cancers showed that lincRNAs with aberrant methylation patterns might be involved in cancer development and progression. The methylated and demethylated lincRNAs identified in this study provide novel insights for developing cancer biomarkers and potential therapeutic targets.