Genome-wide Distribution Research Articles

Distribution Pattern of N6-Methyladenine DNA Modification in the Seashore Paspalum (Paspalum vaginatum) Genome.

N6-methyladenine (6mA) DNA modification has been detected in several eukaryotic organisms, in some of them, it plays important role in the regulation process of stress-resistance response. However, the genome-wide distribution patterns and potential functions of 6mA DNA modification in halophyte Seashore paspalum (Paspalum vaginatum) remain largely unknown. Here, we examined the 6mA landscape in the P. vaginatum genome by adopting single molecule real-time sequencing technology and found that 6mA modification sites were broadly distributed across the P. vaginatum genome. We demonstrated distinct 6mA methylation levels and 6mA distribution patterns in different types of transcription genes, which hinted at different epigenetic rules. Furthermore, the moderate 6mA density genes in P. vaginatum functionally correlated with stress resistance, which also maintained a higher transcriptional level. On the other hand, a specific 6mA distribution pattern in the gene body and near TSS was observed in gene groups with higher RNA expression, which maybe implied some kind of regularity between 6mA site distribution and the protein coding genes transcription was possible. Our study provides new insights into the association between 6mA methylation and gene expression, which may also contribute to key agronomic traits in P. vaginatum.

Genetic dissection of the RNA polymerase II transcription cycle.

How DNA sequence affects the dynamics and position of RNA Polymerase II (Pol II) during transcription remains poorly understood. Here, we used naturally occurring genetic variation in F1 hybrid mice to explore how DNA sequence differences affect the genome-wide distribution of Pol II. We measured the position and orientation of Pol II in eight organs collected from heterozygous F1 hybrid mice using ChRO-seq. Our data revealed a strong genetic basis for the precise coordinates of transcription initiation and promoter proximal pause, allowing us to redefine molecular models of core transcriptional processes. Our results implicate DNA sequence, including both known and novel DNA sequence motifs, as key determinants of the position of Pol II initiation and pause. We report evidence that initiation site selection follows a stochastic process similar to Brownian motion along the DNA template. We found widespread differences in the position of transcription termination, which impact the primary structure and stability of mature mRNA. Finally, we report evidence that allelic changes in transcription often affect mRNA and ncRNA expression across broad genomic domains. Collectively, we reveal how DNA sequences shape core transcriptional processes at single nucleotide resolution in mammals.

Genomic analysis of Rad26 and Rad1-Rad10 reveals differences in their dependence on Mediator and RNA polymerase II.

Mediator is a conserved coregulator playing a key role in RNA polymerase (Pol) II transcription. Mediator also links transcription and nucleotide excision repair (NER) via a direct contact with Rad2/ERCC5(XPG) endonuclease. In this work, we analyzed the genome-wide distribution of Rad26/ERCC6(CSB) and Rad1-Rad10/ERCC4(XPF)-ERCC1, addressing the question of a potential link of these proteins with Mediator and Pol II in yeast Saccharomyces cerevisiae Our genomic analyses reveal that Rad1-Rad10 and Rad26 are present on the yeast genome in the absence of genotoxic stress, especially at highly transcribed regions, with Rad26 binding strongly correlating with that of Pol II. Moreover, we show that Rad1-Rad10 and Rad26 colocalize with Mediator at intergenic regions and physically interact with this complex. Using kin28 TFIIH mutant, we found that Mediator stabilization on core promoters leads to an increase in Rad1-Rad10 chromatin binding, whereas Rad26 occupancy follows mainly a decrease in Pol II transcription. Combined with multivariate analyses, our results show the relationships between Rad1-Rad10, Rad26, Mediator, and Pol II, modulated by the changes in binding dynamics of Mediator and Pol II transcription. In conclusion, we extend the Mediator link to Rad1-Rad10 and Rad26 NER proteins and reveal important differences in their dependence on Mediator and Pol II. Rad2 is the most dependent on Mediator, followed by Rad1-Rad10, whereas Rad26 is the most closely related to Pol II. Our work thus contributes to new concepts of the functional interplay between transcription and DNA repair machineries, which are relevant for human diseases including cancer and XP/CS syndromes.

Hierarchical Accumulation of Histone Variant H2A.Z Regulates Transcriptional States and Histone Modifications in Early Mammalian Embryos.

Early embryos undergo extensive epigenetic reprogramming to achieve gamete‐to‐embryo transition, which involves the loading and removal of histone variant H2A.Z on chromatin. However, how does H2A.Z regulate gene expression and histone modifications during preimplantation development remains unrevealed. Here, by using ultra‐low‐input native chromatin immunoprecipitation and sequencing, the genome‐wide distribution of H2A.Z is delineated in mouse oocytes and early embryos. These landscapes indicate that paternal H2A.Z is removed upon fertilization, followed by unbiased accumulation on parental genomes during zygotic genome activation (ZGA). Remarkably, H2A.Z exhibits hierarchical accumulation as different peak types at promoters: promoters with double H2A.Z peaks are colocalized with H3K4me3 and indicate transcriptional activation; promoters with a single H2A.Z peak are more likely to occupy bivalent marks (H3K4me3+H3K27me3) and indicate development gene suppression; promoters with no H2A.Z accumulation exhibit persisting gene silencing in early embryos. Moreover, H2A.Z depletion changes the enrichment of histone modifications and RNA polymerase II binding at promoters, resulting in abnormal gene expression and developmental arrest during lineage commitment. Furthermore, similar transcription and accumulation patterns between mouse and porcine embryos indicate that a dual role of H2A.Z in regulating the epigenome required for proper gene expression is conserved during mammalian preimplantation development.

Time makes histone H3 modifications drift in mouse liver.

To detect the epigenetic drift of time passing, we determined the genome-wide distributions of mono- and tri-methylated lysine 4 and acetylated and tri-methylated lysine 27 of histone H3 in the livers of healthy 3, 6 and 12 months old C57BL/6 mice. The comparison of different age profiles of histone H3 marks revealed global redistribution of histone H3 modifications with time, in particular in intergenic regions and near transcription start sites, as well as altered correlation between the profiles of different histone modifications. Moreover, feeding mice with caloric restriction diet, a treatment known to retard aging, reduced the extent of changes occurring during the first year of life in these genomic regions.

DIPG-06. Uncovering the FACTs in Diffuse Midline Glioma (DMG)

Effective treatments are urgently needed for the incurable paediatric brainstem tumour Diffuse Midline Glioma (DMG). Most DMGs contain Histone H3 mutations (H3K27M), which produce extensive epigenetic dysregulation by inhibiting EZH2-mediated trimethylation of H3K27 (H3K27me3). Global depletion of the repressive H3K27me3 modification results in aberrantly open chromatin and is central to DMG tumorigenesis. Thus, targeting the epigenome is a promising avenue of treatment for DMG. We found that targeting the histone chaperone complex Facilitates Chromatin Transcription (FACT) with the curaxin drug CBL0137 to have potent pre-clinical efficacy against DMG, leading to a paediatric Phase I/II clinical trial for CBL0137 which includes a DMG/DIPG cohort (NCT04870944). In this project we aim to elucidate CBL0137’s molecular mechanism in DMG. FACT is critical for maintaining chromatin homeostasis during DNA replication, transcription, and repair. We therefore hypothesised that FACT maintains the aberrant epigenetic landscape resulting from H3K27M. Consistently, we found CBL0137 to be more cytotoxic against H3K27M-mutant, compared to H3-WT or isogenic DMG cells lacking the mutation. Furthermore, FACT directly interacts with H3K27M, and FACT inhibition increases both EZH2 catalytic activity and global H3K27me3 levels. We are now using ChIP-seq to discover the genome-wide distribution of FACT and H3K27M, and will interrogate the consequence of CBL0137 treatment on the epigenome and transcriptome. Preliminary results suggest that FACT is located at certain genes co-occupied by H3K27M, the active histone mark H3K27ac, and the BET protein BRD4. This implies that FACT is involved in maintaining open chromatin and perhaps transcription at these regions. Combining CBL0137 with other epigenetic drugs, such as BET inhibitors, could therefore represent a rational therapeutic opportunity for DMG. This work will ultimately inform mechanism-based targeted therapies to combine with CBL0137 to improve its efficacy and uncover valuable new insights into DMG epigenetics and pathobiology.

310-OR: Histone Variant H2A.Z-Mediated Chromatin Organization Regulates Brown Adipocyte Thermogenic Function

Obesity accounts for a substantial and growing burden of global disease, leading to high morbidity and mortality in type 2 diabetes, cardiovascular disease, and certain cancers. Brown adipose tissue (BAT) dissipates energy as heat and has been proposed as a therapeutic target to combat obesity. The unique thermogenic capacity of BAT is primarily conferred by uncoupling protein 1 (Ucp1) , which dissipates energy by uncoupling the proton motive force from ATP generation. To identify novel transcriptional regulators of Ucp1, we performed an unbiased CRISPR/Cas9-based chromatin immunoprecipitation assay (enChIP) and identified H2A.Z as a novel trans-acting factor recruited to the promoter and enhancer of Ucp1 by β3-adrenergic receptor (β3-AR) stimulation. H2A.Z is an essential and evolutionary conserved histone variant of the canonical histone H2A. Knockdown of H2A.Z in brown adipocytes led to lower β3-AR-induced Ucp1 level and impeded the expression of other known thermogenic genes (e.g., Cidea, Pgc1a) , leading to impaired cellular thermogenesis. Mechanistically, by profiling the genome-wide distribution of H2A.Z via ChIP-seq, we observed H2A.Z was enriched at promoter and enhancer of Ucp1 upon β3-AR stimulation, which is consistent with the enChIP result. Intriguingly, we found that the average occupancy of H2A.Z in the promoter and enhancer regions of other β3-AR-upregulated genes (e.g., Zfp516, Vdr, Pck1, Gys2) was also higher than that in non-stimulated condition, indicating a global regulatory role of H2A.Z in the thermogenic program in response to β3-AR stimulation. Additionally, in the enhancer regions of β3-AR-induced genes, we also detected an increased level of H2A.Z acetylation, which is a mark of active chromatin and confers an open chromatin conformation. Taken together, these findings reveal the essential role of H2A.Z-mediated chromatin organization in BAT thermogenic function. Disclosure Y.Zhang: None. R.Zheng: None. C.Wang: None. J.Darcy: None. K.Chen: None. Y.Tseng: Consultant; Cellarity. Funding National Institutes of Health (R01DK102898)

USP16-mediated histone H2A lysine-119 deubiquitination during oocyte maturation is a prerequisite for zygotic genome activation.

Maternal-to-zygotic transition (MZT) is the first and key step in the control of animal development and intimately related to changes in chromatin structure and histone modifications. H2AK119ub1, an important epigenetic modification in regulating chromatin configuration and function, is primarily catalyzed by PRC1 and contributes to resistance to transcriptional reprogramming in mouse embryos. In this study, the genome-wide dynamic distribution of H2AK119ub1 during MZT in mice was investigated using chromosome immunoprecipitation and sequencing. The results indicated that H2AK119ub1 accumulated in fully grown oocytes and was enriched at the TSSs of maternal genes, but was promptly declined after meiotic resumption at genome-wide including the TSSs of early zygotic genes, by a previously unidentified mechanism. Genetic evidences indicated that ubiquitin-specific peptidase 16 (USP16) is the major deubiquitinase for H2AK119ub1 in mouse oocytes. Conditional knockout of Usp16 in oocytes did not impair their survival, growth, or meiotic maturation. However, oocytes lacking USP16 have defects when undergoing zygotic genome activation or gaining developmental competence after fertilization, potentially associated with high levels of maternal H2AK119ub1 deposition on the zygotic genomes. Taken together, H2AK119ub1 level is declined during oocyte maturation by an USP16-dependent mechanism, which ensures zygotic genome reprogramming and transcriptional activation of essential early zygotic genes.

Encyclopaedia of eukaryotic DNA methylation: from patterns to mechanisms and functions.

DNA methylation is an epigenetic modification with a very long evolutionary history. However, DNA methylation evolves surprisingly rapidly across eukaryotes. The genome-wide distribution of methylation diversifies rapidly in different lineages, and DNA methylation is lost altogether surprisingly frequently. The growing availability of genomic and epigenomic sequencing across organisms highlights this diversity but also illuminates potential factors that could explain why both the DNA methylation machinery and its genome-wide distribution evolve so rapidly. Key to this are new discoveries about the fitness costs associated with DNA methylation, and new theories about how the fundamental biochemical mechanisms of DNA methylation introduction and maintenance could explain how new genome-wide patterns of methylation evolve.

Comparative Transcriptome Analysis Reveals Regulatory Mechanism of Long Non-Coding RNAs during Abdominal Preadipocyte Adipogenic Differentiation in Chickens.

Simple SummaryExcessive fat deposition, particularly abdominal fat, negatively affects meat quality, feed utilization, and reproduction performance in chickens. Adipogenesis is a well-orchestrated process that involves a network of regulatory molecules. Emerging evidence has shown that long non-coding RNAs (lncRNAs) participate in mammalian adipogenesis. However, limited information is available about their role in adipogenesis in chickens. Therefore, we investigated the genome-wide identification, expression profiles, and regulatory networks of lncRNAs during adipogenic differentiation of chicken abdominal preadipocytes. We found that lncRNAs were expressed in chicken abdominal adipocytes, characterized by their differentiation stage-specific expression pattern and were involved in multiple lipid metabolism-related pathways. Several lncRNAs were identified as key lncRNAs responsible for adipogenic differentiation. These key lncRNAs could regulate lipid-related protein-encoding genes in cis- and trans-regulation manners, and/or competitively sponge miRNA to sequester it away from lipid-related genes, thereby functioning as regulators of chicken adipogenesis. These findings provide valuable information that improves our understanding of functional mechanisms of lncRNAs underlying adipogenesis and contributes to the genetic improvement of excessive abdominal fat in chickens.Long non-coding RNAs (lncRNAs) are implicated in mammalian adipogenesis and obesity. However, their genome-wide distribution, expression profiles, and regulatory mechanisms during chicken adipogenesis remain rarely understood. In the present study, lncRNAs associated with adipogenesis were identified from chicken abdominal adipocytes at multiple differentiation stages using Ribo-Zero RNA-seq. A total of 15,179 lncRNAs were identified and characterized by stage-specific expression patterns. Of these, 840 differentially expressed lncRNAs were detected, and their cis- and trans-target genes were significantly enriched in multiple lipid-related pathways. Through weighted gene co-expression network analysis (WGCNA) and time-series expression profile clustering analysis, 14 key lncRNAs were identified as candidate regulatory lncRNAs in chicken adipogenic differentiation. The cis- and trans-regulatory interactions of key lncRNAs were constructed based on their differentially expressed cis- and trans-target genes, respectively. We also constructed a competing endogenous RNA (ceRNA) network based on the key lncRNAs, differentially expressed miRNAs, and differentially expressed mRNAs. MSTRG.25116.1 was identified as a potential regulator of chicken abdominal preadipocyte adipogenic differentiation by acting as a transcriptional trans-regulator of fatty acid amide hydrolase (FAAH) gene expression and/or a ceRNA that post-transcriptionally mediates FAAH gene expression by sponging gga-miR-1635.

Genome-Wide TSS Distribution in Three Related Clostridia with Normalized Capp-Switch Sequencing.

ABSTRACTTranscription initiation is a tightly regulated process that is crucial for many aspects of prokaryotic physiology. High-throughput transcription start site (TSS) mapping can shed light on global and local regulation of transcription initiation, which in turn may help us understand and predict microbial behavior. In this study, we used Capp-Switch sequencing to determine the TSS positions in the genomes of three model solventogenic clostridia: Clostridium acetobutylicum ATCC 824, C. beijerinckii DSM 6423, and C. beijerinckii NCIMB 8052. We first refined the approach by implementing a normalization pipeline accounting for gene expression, yielding a total of 12,114 mapped TSSs across the species. We further compared the distributions of these sites in the three strains. Results indicated similar distribution patterns at the genome scale, but also some sharp differences, such as for the butyryl-CoA synthesis operon, particularly when comparing C. acetobutylicum to the C. beijerinckii strains. Lastly, we found that promoter structure is generally poorly conserved between C. acetobutylicum and C. beijerinckii. A few conserved promoters across species are discussed, showing interesting examples of how TSS determination and comparison can improve our understanding of gene expression regulation at the transcript level.IMPORTANCE Solventogenic clostridia have been employed in industry for more than a century, initially being used in the acetone-butanol-ethanol (ABE) fermentation process for acetone and butanol production. Interest in these bacteria has recently increased in the context of green chemistry and sustainable development. However, our current understanding of their genomes and physiology limits their optimal use as industrial solvent production platforms. The gene regulatory mechanisms of solventogenesis are still only partly understood, impeding efforts to increase rates and yields. Genome-wide mapping of transcription start sites (TSSs) for three model solventogenic Clostridium strains is an important step toward understanding mechanisms of gene regulation in these industrially important bacteria.

Transcription factors KLF15 and PPARδ cooperatively orchestrate genome-wide regulation of lipid metabolism in skeletal muscle

Skeletal muscle dynamically regulates systemic nutrient homeostasis through transcriptional adaptations to physiological cues. In response to changes in the metabolic environment (e.g., alterations in circulating glucose or lipid levels), networks of transcription factors and coregulators are recruited to specific genomic loci to fine-tune homeostatic gene regulation. Elucidating these mechanisms is of particular interest as these gene regulatory pathways can serve as potential targets to treat metabolic disease. The zinc-finger transcription factor Krüppel-like factor 15 (KLF15) is a critical regulator of metabolic homeostasis; however, its genome-wide distribution in skeletal muscle has not been previously identified. Here, we characterize the KLF15 cistrome in vivo in skeletal muscle and find that the majority of KLF15 binding is localized to distal intergenic regions and associated with genes related to circadian rhythmicity and lipid metabolism. We also identify critical interdependence between KLF15 and the nuclear receptor PPARδ in the regulation of lipid metabolic gene programs. We further demonstrate that KLF15 and PPARδ colocalize genome-wide, physically interact, and are dependent on one another to exert their transcriptional effects on target genes. These findings reveal that skeletal muscle KLF15 plays a critical role in metabolic adaptation through its direct actions on target genes and interactions with other nodal transcription factors such as PPARδ.

Th1 immune maturation effects of Nocardia rubra cell-wall skeleton via PI3K/Akt/PAX8 regulatory axis.

Nocardia rubra cell-wall skeleton (Nr-CWS) is reported as an external immunotherapeutic enhancer with the advantage of antitumor effect on human cancers. However, the immune regulatory role of Nr-CWS is not fully illustrated. We studied mouse CD4+ T lymphocytes isolated from mice spleen were induced by Nr-CWS and observed that the differentiation of Th1 CD4+ T cells and the cytokines of IL-2, TNF-α, IFN-γ were all enhanced by Nr-CWS. Furthermore, RNA sequencing was conducted to investigate the different mRNA profiling induced by Nr-CWS. We observed that paired box 8 (PAX8) was significantly up-regulated in Nr-CWS-treated Th1 cells compared to control. As a transcription factor, chromatin immunoprecipitation sequencing was carried out to study the genome-wide distribution of PAX8. Interestingly, we found that the binding domain of PAX8 was elevated by Nr-CWS, and the target genes associated with these binding sites showed a positive correlation between their transcription and PAX8 binding strength. Finally, we determined that Nr-CWS could enhance the activity of the PI3 K/Akt signaling pathway. Akt agonist could mimic the effect of Nr-CWS for PAX8 up-regulation, while Akt inhibitor compromised the expression of PAX8. Taken together, we determined a novel role of Nr-CWS in boosting the activity of Th1 maturation via the PI3 K/Akt/PAX8 axis.

Investigation of genome-wide InDel distribution and segregation in Populus with restriction-site associated DNA sequencing data

Insertion or deletion polymorphism (InDel) is one of the main genetic variations in plant genomes. However, there are few studies on InDels across the whole genome in Populus. In this study, we investigated genome-wide InDels in Populus deltoides and Populus simonii and InDel segregation in their F1 hybrid population with restriction-site associated DNA sequencing (RAD-seq) data. A total of 119,066 InDels were identified in P. deltoides and P. simonii according to the reference genome of Populus trichocarpa, including 58,532 unique InDels in P. deltoides, 54,469 unique InDels in P. simonii, and 6,065 common InDels in both. Meanwhile, the distribution of these InDels was analyzed along chromosomes, indicating that the distribution patterns for both species were largely similar, but the average InDel density was slightly higher in P. deltoides than in P. simonii. GO annotation and enrichment analysis of those genes harboring InDels showed the same patterns between the two poplar species. It is interesting to find that the ratio (~ 46%) of the common InDels within genes to all common InDels was higher than that of the InDels within genes to all InDels in P. deltoides (~ 35%) or in P. simonii (~ 34%), possibly indicating that those InDels are more conservative between poplar species. Moreover, investigation of the InDel segregation patterns demonstrated that a large number of Mendelian InDels could be selected for genetic mapping in the F1 hybrid population. RAD-seq provides genome-wide insights into the InDel distributions in P. deltoides and P. simonii and the segregation patterns in their progeny, providing valuable genomic variation information for genetic and evolutionary studies in Populus.

The dark proteome: translation from noncanonical open reading frames.

Omics-based technologies have revolutionized our understanding of the coding potential of the genome. In particular, these studies revealed widespread unannotated open reading frames (ORFs) throughout genomes and that these regions have the potential to encode novel functional (micro-)proteins and/or hold regulatory roles. However, despite their genomic prevalence, relatively few of these noncanonical ORFs have been functionally characterized, likely in part due to their under-recognition by the broader scientific community. The few that have been investigated in detail have demonstrated their essentiality in critical and divergent biological processes. As such, here we aim to discuss recent advances in understanding the diversity of noncanonical ORFs and their roles, as well as detail biologically important examples within the context of the mammalian genome.

Identification of the genome-wide distribution of cyanobacterial group-2 sigma factor SigE, accountable for its regulon.

SUMMARYSigma factors are the subunits of bacterial RNA polymerase that govern the expression of genes by recognizing the promoter sequence. Cyanobacteria, which are oxygenic phototrophic eubacteria, have multiple alternative sigma factors that respond to various environmental stresses. The subgroup highly homologous to the primary sigma factor (SigA) is called the group‐2 sigma factor. The model cyanobacteriumSynechocystis sp. PCC 6803 has four group‐2 sigma factors (SigB–SigE) conserved within the phylum Cyanobacteria. Among the group‐2 sigma factors in Synechocystis sp. PCC 6803, SigE is unique because it alters metabolism by inducing the expression of genes related to sugar catabolism and nitrogen metabolism. However, the features of the promoter sequences of the SigE regulon remain elusive. Here, we identified the direct targets of SigA and SigE by chromatin immunoprecipitation sequencing (ChIP‐seq). We then showed that the binding sites of SigE and SigA overlapped substantially, but that SigE localized exclusively to SigE‐dependent promoters. We also found consensus sequences from SigE‐dependent promoters and confirmed their importance. ChIP‐seq analysis showed both the redundancy and the specificity of SigE, compared with SigA, integrating information obtained from a previously adopted genetic approach and in vitro assays. The features of SigE elucidated in our study indicate its similarity with group‐2 sigma factors of other bacteria, even though they are evolutionarily irrelevant. Our approach is also applicable to other organisms and organelles, such as plant plastids, which have multiple group‐2 sigma factors.

Epigenetic reprogramming by TET enzymes impacts co-transcriptional R-loops.

DNA oxidation by ten-eleven translocation (TET) family enzymes is essential for epigenetic reprogramming. The conversion of 5-methylcytosine (5mC) into 5-hydroxymethylcytosine (5hmC) initiates developmental and cell-type-specific transcriptional programs through mechanisms that include changes in the chromatin structure. Here, we show that the presence of 5hmC in the transcribed gene promotes the annealing of the nascent RNA to the template DNA strand, leading to the formation of an R-loop. Depletion of TET enzymes reduced global R-loops in the absence of gene expression changes, whereas CRISPR-mediated tethering of TET to an active gene promoted the formation of R-loops. The genome-wide distribution of 5hmC and R-loops shows a positive correlation in mouse and human stem cells and overlap in half of all active genes. Moreover, R-loop resolution leads to differential expression of a subset of genes that are involved in crucial events during stem cell proliferation. Altogether, our data reveal that epigenetic reprogramming via TET activity promotes co-transcriptional R-loop formation, disclosing new mechanisms of gene expression regulation.

Histone H3 Lysine 4 and 27 Trimethylation Landscape of Human Alzheimer's Disease.

Background: Epigenetic remodeling is emerging as a critical process for both the onset and progression of Alzheimer’s disease (AD), the most common form of neurodegenerative dementia. However, it is not clear to what extent the distribution of histone modifications is involved in AD. Methods: To investigate histone H3 modifications in AD, we compared the genome-wide distributions of H3K4me3 and H3K27me3 in entorhinal cortices from severe sporadic AD patients and from age-matched healthy individuals of both sexes. Results: AD samples were characterized by typical average levels and distributions of the H3K4me3 and H3K27me3 signals. However, AD patients showed a lower H3K4me3 and higher H3K27me3 signal, particularly in males. Interestingly, the genomic sites found differentially trimethylated at the H3K4 between healthy and AD samples involve promoter regions of genes belonging to AD-related pathways such as glutamate receptor signaling. Conclusions: The signatures of H3K4me3 and H3K27me3 identified in AD patients validate the role of epigenetic chromatin remodeling in neurodegenerative disease and shed light on the genomic adaptive mechanisms involved in AD.

Transposon- and Genome Dynamics in the Fungal Genus Neurospora: Insights from Nearly Gapless Genome Assemblies

A large portion of nuclear DNA is composed of transposable element (TE) sequences, whose transposition is controlled by diverse host defense strategies in order to maintain genomic integrity. One such strategy is the fungal-specific Repeat-Induced Point mutation (RIP) that hyper-mutates repetitive DNA sequences. While RIP is found across Fungi, it has been shown to vary in efficiency. The filamentous ascomycete Neurospora crassa has been a pioneer in the study of RIP, but data on TEs and RIP from other species in the genus is limited. In this study, we investigated 18 nearly gapless genome assemblies of ten Neurospora species, which diverged from a common ancestor about 7 MYA, to determine and compare genome-wide TE distribution and their associated RIP patterns. Four of these assemblies, generated by PacBio technology, represent new genomic datasets. We showed that the TE contents between 8.7-18.9% covary with genome sizes that range between 37.8-43.9 Mb. Degraded copies of Long Terminal Repeat (LTR) retrotransposons were abundant among the identified TEs, and these are distributed across the genome at varying frequencies. In all investigated Neurospora genomes, TE sequences had signs of numerous C-to-T substitutions, suggesting that RIP occurred in all species, and accordingly, RIP signatures correlated with TE-dense regions in all genomes. In conclusion, essentially gapless genome assemblies allowed us to identify TEs in Neurospora genomes, and reveal that TEs contribute to genome size variation in this group. Our study suggests that TEs and RIP are highly correlated in each examined Neurospora species, and hence, the pattern of interaction is conserved over the investigated evolutionary timescale. Finally, with our results, we verify that RIP signatures can be used to facilitate the identification of TE-rich region in the genome. The comprehensive genomic dataset of Neurospora is a rich resource for further in-depth analyses of fungal genomes by the community.

A Chromosome-Level Genome Assembly of the European Beech (Fagus sylvatica) Reveals Anomalies for Organelle DNA Integration, Repeat Content and Distribution of SNPs.

The European Beech is the dominant climax tree in most regions of Central Europe and valued for its ecological versatility and hardwood timber. Even though a draft genome has been published recently, higher resolution is required for studying aspects of genome architecture and recombination. Here, we present a chromosome-level assembly of the more than 300 year-old reference individual, Bhaga, from the Kellerwald-Edersee National Park (Germany). Its nuclear genome of 541 Mb was resolved into 12 chromosomes varying in length between 28 and 73 Mb. Multiple nuclear insertions of parts of the chloroplast genome were observed, with one region on chromosome 11 spanning more than 2 Mb which fragments up to 54,784 bp long and covering the whole chloroplast genome were inserted randomly. Unlike in Arabidopsis thaliana, ribosomal cistrons are present in Fagus sylvatica only in four major regions, in line with FISH studies. On most assembled chromosomes, telomeric repeats were found at both ends, while centromeric repeats were found to be scattered throughout the genome apart from their main occurrence per chromosome. The genome-wide distribution of SNPs was evaluated using a second individual from Jamy Nature Reserve (Poland). SNPs, repeat elements and duplicated genes were unevenly distributed in the genomes, with one major anomaly on chromosome 4. The genome presented here adds to the available highly resolved plant genomes and we hope it will serve as a valuable basis for future research on genome architecture and for understanding the past and future of European Beech populations in a changing climate.