Higher-Order Inter-chromosomal Hubs Shape 3D Genome Organization in the Nucleus.

Eukaryotic genomes are packaged into a 3-dimensional structure in the nucleus. Current methods for studying genome-wide structure are based on proximity ligation. However, this approach can fail to detect known structures, such as interactions with nuclear bodies, because these DNA regions can be too far apart to directly ligate. Accordingly, our overall understanding of genome organization remains incomplete. Here, we develop split-pool recognition of interactions by tag extension (SPRITE), a method that enables genome-wide detection of higher-order interactions within the nucleus. Using SPRITE, we recapitulate known structures identified by proximity ligation and identify additional interactions occurring across larger distances, including two hubs of inter-chromosomal interactions that are arranged around the nucleolus and nuclear speckles. We show that a substantial fraction of the genome exhibits preferential organization relative to these nuclear bodies. Our results generate a global model whereby nuclear bodies act as inter-chromosomal hubs that shape the overall packaging of DNA in the nucleus.

Enhancer-promoter interactions are reconfigured through the formation of long-range multiway hubs as mouse ES cells exit pluripotency.

Enhancers bind transcription factors, chromatin regulators, and non-coding transcripts to modulate the expression of target genes. Here, we report 3D genome structures of single mouse ES cells as they are induced to exit pluripotency and transition through a formative stage prior to undergoing neuroectodermal differentiation. We find that there is a remarkable reorganization of 3D genome structure where inter-chromosomal intermingling increases dramatically in the formative state. This intermingling is associated with the formation of a large number of multiway hubs that bring together enhancers and promoters with similar chromatin states from typically 5-8 distant chromosomal sites that are often separated by many Mb from each other. In the formative state, genes important for pluripotency exit establish contacts with emerging enhancers within these multiway hubs, suggesting that the structural changes we have observed may play an important role in modulating transcription and establishing new cell identities.

TFIIIC Binding to Alu Elements Controls Gene Expression via Chromatin Looping and Histone Acetylation.

How repetitive elements, epigenetic modifications, and architectural proteins interact ensuring proper genome expression remains poorly understood. Here, we report regulatory mechanisms unveiling a central role of Alu elements (AEs) and RNA polymerase III transcription factor C (TFIIIC) in structurally and functionally modulating the genome via chromatin looping and histone acetylation. Upon serum deprivation, a subset of AEs pre-marked by the activity-dependent neuroprotector homeobox Protein (ADNP) and located near cell-cycle genes recruits TFIIIC, which alters their chromatin accessibility by direct acetylation of histone H3 lysine-18 (H3K18). This facilitates the contacts of AEs with distant CTCF sites near promoter of other cell-cycle genes, which also become hyperacetylated at H3K18. These changes ensure basal transcription of cell-cycle genes and are critical for their re-activation upon serum re-exposure. Our study reveals how direct manipulation of the epigenetic state of AEs by a general transcription factor regulates 3D genome folding and expression.

A three-dimensional genomics view for speciation research.

Genomic information is folded in a three-dimensional (3D) structure, a rarely explored evolutionary driver of speciation. Technological advances now enable the study of 3D genome structures (3DGSs) across the Tree of Life. At the onset of 3D speciation genomics, we discuss the putative roles of 3DGSs in speciation.

Reversible acetylation of HDAC8 regulates cell cycle.

HDAC8, a member of class I HDACs, plays a pivotal role in cell cycle regulation by deacetylating the cohesin subunit SMC3. While cyclins and CDKs are well-established cell cycle regulators, our knowledge of other regulators remains limited. Here we reveal the acetylation of K202 in HDAC8 as a key cell cycle regulator responsive to stress. K202 acetylation in HDAC8, primarily catalyzed by Tip60, restricts HDAC8 activity, leading to increased SMC3 acetylation and cell cycle arrest. Furthermore, cells expressing the mutant form of HDAC8 mimicking K202 acetylation display significant alterations in gene expression, potentially linked to changes in 3D genome structure, including enhanced chromatid loop interactions. K202 acetylation impairs cell cycle progression by disrupting the expression of cell cycle-related genes and sister chromatid cohesion, resulting in G2/M phase arrest. These findings indicate the reversible acetylation of HDAC8 as a cell cycle regulator, expanding our understanding of stress-responsive cell cycle dynamics.

Assembly, comparative analysis, and utilization of a single haplotype reference genome for soybean.

Cultivar Williams 82 has served as the reference genome for the soybean research community since 2008, but is known to have areas of genomic heterogeneity among different sub-lines. This work provides an updated assembly (version Wm82.a6) derived from a specific sub-line known as Wm82-ISU-01 (seeds available under USDA accession PI 704477). The genome was assembled using Pacific BioSciences HiFi reads and integrated into chromosomes using HiC. The 20 soybean chromosomes assembled into a genome of 1.01Gb, consisting of 36 contigs. The genome annotation identified 48 387 gene models, named in accordance with previous assembly versions Wm82.a2 and Wm82.a4. Comparisons of Wm82.a6 with other near-gapless assemblies of Williams 82 reveal large regions of genomic heterogeneity, including regions of differential introgression from the cultivar Kingwa within approximately 30 Mb and 25 Mb segments on chromosomes 03 and 07, respectively. Additionally, our analysis revealed a previously unknown large (>20 Mb) heterogeneous region in the pericentromeric region of chromosome 12, where Wm82.a6 matches the 'Williams' haplotype while the other two near-gapless assemblies do not match the haplotype of either parent of Williams 82. In addition to the Wm82.a6 assembly, we also assembled the genome of 'Fiskeby III,' a rich resource for abiotic stress resistance genes. A genome comparison of Wm82.a6 with Fiskeby III revealed the nucleotide and structural polymorphisms between the two genomes within a QTL region for iron deficiency chlorosis resistance. The Wm82.a6 and Fiskeby III genomes described here will enhance comparative and functional genomics capacities and applications in the soybean community.

Evolution of sea-surfing plant propagule as revealed by the genomes of Heritiera mangroves.

Coastal forests, such as mangroves, protect much of the tropical and subtropical coasts. Long-distance dispersal via sea-surfing propagules is essential for coastal plants, but the genomic and molecular basis of sea-surfing plant propagule evolution remains unclear. Heritiera fomes and Heritiera littoralis are two coastal plants with typical buoyant fruits. We de novo sequenced and assembled their high-quality genomes. Our phylogenomic analysis indicates H. littoralis and H. fomes originated (at ~6.08 Mya) just before the start of Quaternary sea-level fluctuations. Whole-genome duplication occurred earlier, permitting gene copy gains in the two species. Many of the expanded gene families are involved in lignin and flavonoid biosynthesis, likely contributing to buoyant fruit emergence. It is repeatedly revealed that one duplicated copy to be under positive selection while the other is not. By examining H. littoralis fruits at three different developmental stages, we found that gene expression levels remain stable from young to intermediate. However, ~1000 genes are up-regulated and ~ 3000 genes are down-regulated as moving to mature. Particularly in fruit epicarps, the upregulation of WRKY12 and E2Fc likely constrains the production of p-Coumaroyl-CoA, the key internal substrate for lignin biosynthesis. Hence, to increase fruit impermeability, methylated lignin biosynthesis is shut down by down-regulating the genes CCoAOMT, F5H, COMT, and CSE, while unmethylated lignins are preferentially produced by upregulating CAD and CCR. Similarly, cutin polymers and cuticular waxes accumulate with high levels before maturation in epicarps. Overall, our genome assemblies and analyses uncovered the genomic evolution and temporal transcriptional regulation of sea-surfing propagule.

InferLoop: leveraging single-cell chromatin accessibility for the signal of chromatin loop.

Deciphering cell-type-specific 3D structures of chromatin is challenging. Here, we present InferLoop, a novel method for inferring the strength of chromatin interaction using single-cell chromatin accessibility data. The workflow of InferLoop is, first, to conduct signal enhancement by grouping nearby cells into bins, and then, for each bin, leverage accessibility signals for loop signals using a newly constructed metric that is similar to the perturbation of the Pearson correlation coefficient. In this study, we have described three application scenarios of InferLoop, including the inference of cell-type-specific loop signals, the prediction of gene expression levels and the interpretation of intergenic loci. The effectiveness and superiority of InferLoop over other methods in those three scenarios are rigorously validated by using the single-cell 3D genome structure data of human brain cortex and human blood, the single-cell multi-omics data of human blood and mouse brain cortex, and the intergenic loci in the GWAS Catalog database as well as the GTEx database, respectively. In addition, InferLoop can be applied to predict loop signals of individual spots using the spatial chromatin accessibility data of mouse embryo. InferLoop is available at https://github.com/jumphone/inferloop.

A Kuhnian revolution in molecular biology: Most genes in complex organisms express regulatory RNAs.

Thomas Kuhn described the progress of science as comprising occasional paradigm shifts separated by interludes of 'normal science'. The paradigm that has held sway since the inception of molecular biology is that genes (mainly) encode proteins. In parallel, theoreticians posited that mutation is random, inferred that most of the genome in complex organisms is non-functional, and asserted that somatic information is not communicated to the germline. However, many anomalies appeared, particularly in plants and animals: the strange genetic phenomena of paramutation and transvection; introns; repetitive sequences; a complex epigenome; lack of scaling of (protein-coding) genes and increase in 'noncoding' sequences with developmental complexity; genetic loci termed 'enhancers' that control spatiotemporal gene expression patterns during development; and a plethora of 'intergenic', overlapping, antisense and intronic transcripts. These observations suggest that the original conception of genetic information was deficient and that most genes in complex organisms specify regulatory RNAs, some of which convey intergenerational information. Also see the video abstract here: https://youtu.be/qxeGwahBANw.

Transposable elements drive intron gain in diverse eukaryotes.

There is massive variation in intron numbers across eukaryotic genomes, yet the major drivers of intron content during evolution remain elusive. Rapid intron loss and gain in some lineages contrast with long-term evolutionary stasis in others. Episodic intron gain could be explained by recently discovered specialized transposons called Introners, but so far Introners are only known from a handful of species. Here, we performed a systematic search across 3,325 eukaryotic genomes and identified 27,563 Introner-derived introns in 175 genomes (5.2%). Species with Introners span remarkable phylogenetic diversity, from animals to basal protists, representing lineages whose last common ancestor dates to over 1.7 billion years ago. Aquatic organisms were 6.5 times more likely to contain Introners than terrestrial organisms. Introners exhibit mechanistic diversity but most are consistent with DNA transposition, indicating that Introners have evolved convergently hundreds of times from nonautonomous transposable elements. Transposable elements and aquatic taxa are associated with high rates of horizontal gene transfer, suggesting that this combination of factors may explain the punctuated and biased diversity of species containing Introners. More generally, our data suggest that Introners may explain the episodic nature of intron gain across the eukaryotic tree of life. These results illuminate the major source of ongoing intron creation in eukaryotic genomes.

Host-specific co-evolution likely driven by diet in Buchnera aphidicola.

BACKGROUND: Russian wheat aphid (Diuraphis noxia Kurd.) is a severe pest to wheat, and even though resistance varieties are available to curb this pest, they are becoming obsolete with the development of new virulent aphid populations. Unlike many other aphids, D noxia only harbours a single endosymbiont, Buchnera aphidicola. Considering the importance of Buchnera, this study aimed to elucidate commonalities and dissimilarities between various hosts, to better understand its distinctiveness within its symbiotic relationship with D. noxia. To do so, the genome of the D. noxia's Buchnera was assembled and compared to those of other aphid species that feed on diverse host species. RESULTS: The overall importance of several features such as gene length and percentage GC content was found to be critical for the maintenance of Buchnera genes when compared to their closest free-living relative, Escherichia coli. Buchnera protein coding genes were found to have percentage GC contents that tended towards a mean of ~ 26% which had strong correlation to their identity to their E. coli homologs. Several SNPs were identified between different aphid populations and multiple isolates of Buchnera were confirmed in single aphids. CONCLUSIONS: Establishing the strong correlation of percentage GC content of protein coding genes and gene identity will allow for identifying which genes will be lost in the continually shrinking Buchnera genome. This is also the first report of a parthenogenically reproducing aphid that hosts multiple Buchnera strains in a single aphid, raising questions regarding the benefits of maintaining multiple strains. We also found preliminary evidence for post-transcriptional regulation of Buchnera genes in the form of polyadenylation.

Alteration of chromosome structure impacts gene expressions implicated in pancreatic ductal adenocarcinoma cells.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) is a lethal malignancy with a five-year survival rate of approximately 10%. Genetic mutations are pivotal drivers in PDAC pathogenesis, but recent investigations also revealed the involvement of non-genetic alterations in the disease development. In this study, we undertook a multi-omics approach, encompassing ATAC-seq, RNA-seq, ChIP-seq, and Hi-C methodologies, to dissect gene expression alterations arising from changes in chromosome accessibility and chromatin three-dimensional interactions in PDAC. RESULTS: Our findings indicate that chromosomal structural alterations can lead to abnormal expressions on key genes during PDAC development. Notably, overexpression of oncogenes FGFR2, FOXA2, CYP2R1, and CPOX can be attributed to the augmentation of promoter accessibility, coupled with long-range interactions with distal elements. Additionally, our findings indicate that chromosomal structural alterations caused by genomic instability can lead to abnormal expressions in PDACs. As an example, by analyzing chromosomal changes, we identified a putative oncogenic gene, LPAR1, which shows upregulated expression in both PDAC cell lines and clinical samples. The overexpression is correlated with alterations in LPAR1-associated 3D genome structure and chromatin state. We further demonstrated that high LPAR1 activity is required for enhanced PDAC cell migration in vitro. CONCLUSIONS: Collectively, our findings reveal that the chromosomal conformational alterations, in addition to the well-known genetic mutations, are critical for PDAC tumorigenesis.

Mitochondrial genome structure and composition in 70 fishes: a key resource for fisheries management in the South Atlantic.

BACKGROUND: Phylogenetic gaps of public databases of reference sequences are a major obstacle for comparative genomics and management of marine resources, particularly in the Global South, where economically important fisheries and conservation flagship species often lack closely-related references. We applied target-enrichment to obtain complete mitochondrial genomes of marine ichthyofauna from the Brazilian coast selected based on economic significance, conservation status and lack of phylogenetically-close references. These included sardines (Dorosomatidae, Alosidae), mackerels (Scombridae) croakers (Sciaenidae), groupers (Epinephelidae) and snappers (Lutjanidae). RESULTS: Custom baits were designed to enrich mitochondrial DNA across a broad phylogenetic range of fishes. Sequencing generated approximately 100k reads per sample, which were assembled in a total of 70 complete mitochondrial genomes and include fifty-two new additions to GenBank, including five species with no previous mitochondrial data. Departures from the typical gene content and order occurred in only three taxa and mostly involved tRNA gene duplications. Start-codons for all genes, except Cytochrome C Oxidase subunit I (COI), were consistently ATG, whilst a wide range of stop-codons deviated from the prevailing TAA. Phylogenetic analysis confirmed assembly accuracy and revealed signs of cryptic diversification within the Mullus genus. Lineage delimitation methods using Sardinella aurita and S. brasiliensis mitochondrial genomes support a single Operational Taxonomic Unit. CONCLUSIONS: Target enrichment was highly efficient, providing complete novel mitochondrial genomes with little sequencing effort. These sequences are deposited in public databases to enable subsequent studies in population genetics and adaptation of Latin American fish species and serve as a vital resource for conservation and management programs that rely on molecular data for species and genus-level identification.

Molecular Characterization and Phylogenetic Analysis of Centipedegrass [Eremochloa ophiuroides (Munro) Hack.] Based on the Complete Chloroplast Genome Sequence.

Centipedegrass (Eremochloa ophiuroides) is an important warm-season grass plant used as a turfgrass as well as pasture grass in tropical and subtropical regions, with wide application in land surface greening and soil conservation in South China and southern United States. In this study, the complete cp genome of E. ophiuroides was assembled using high-throughput Illumina sequencing technology. The circle pseudomolecule for E. ophiuroides cp genome is 139,107 bp in length, with a quadripartite structure consisting of a large single copyregion of 82,081 bp and a small single copy region of 12,566 bp separated by a pair of inverted repeat regions of 22,230 bp each. The overall A + T content of the whole genome is 61.60%, showing an asymmetric nucleotide composition. The genome encodes a total of 131 gene species, composed of 20 duplicated genes within the IR regions and 111 unique genes comprising 77 protein-coding genes, 30 transfer RNA genes, and 4 ribosome RNA genes. The complete cp genome sequence contains 51 long repeats and 197 simple sequence repeats, and a high degree of collinearity among E. ophiuroide and other Gramineae plants was disclosed. Phylogenetic analysis showed E. ophiuroides, together with the other two Eremochloa species, is closely related to Mnesithea helferi within the subtribe Rottboelliinae. These findings will be beneficial for the classification and identification of the Eremochloa taxa, phylogenetic resolution, novel gene discovery, and functional genomic studies for the genus Eremochloa.

Plastid DNA is a major source of nuclear genome complexity and of RNA genes in the orphan crop moringa.

BACKGROUND: Unlike Transposable Elements (TEs) and gene/genome duplication, the role of the so-called nuclear plastid DNA sequences (NUPTs) in shaping the evolution of genome architecture and function remains poorly studied. We investigate here the functional and evolutionary fate of NUPTs in the orphan crop Moringa oleifera (moringa), featured by the highest fraction of plastid DNA found so far in any plant genome, focusing on (i) any potential biases in their distribution in relation to specific nuclear genomic features, (ii) their contribution to the emergence of new genes and gene regions, and (iii) their impact on the expression of target nuclear genes. RESULTS: In agreement with their potential mutagenic effect, NUPTs are underrepresented among structural genes, although their overall transcription levels and broadness were only lower when involved exonic regions; the occurrence of plastid DNA generally did not result in a broader expression, except among those affected in introns by older NUPTs. In contrast, we found a strong enrichment of NUPTs among specific superfamilies of retrotransposons and several classes of RNA genes, including those participating in the protein biosynthetic machinery (i.e., rRNA and tRNA genes) and a specific class of regulatory RNAs. A significant fraction of NUPT RNA genes was found to be functionally expressed, thus potentially contributing to the nuclear pool. CONCLUSIONS: Our results complete our view of the molecular factors driving the evolution of nuclear genome architecture and function, and support plastid DNA in moringa as a major source of (i) genome complexity and (ii) the nuclear pool of RNA genes.

Modulated illumination microscopy: Application perspectives in nuclear nanostructure analysis.

The structure of the cell nucleus of higher organisms has become a major topic of advanced light microscopy. So far, a variety of methods have been applied, including confocal laser scanning fluorescence microscopy, 4Pi, STED and localisation microscopy approaches, as well as different types of patterned illumination microscopy, modulated either laterally (in the object plane) or axially (along the optical axis). Based on our experience, we discuss here some application perspectives of Modulated Illumination Microscopy (MIM) and its combination with single-molecule localisation microscopy (SMLM). For example, spatially modulated illumination microscopy/SMI (illumination modulation along the optical axis) has been used to determine the axial extension (size) of small, optically isolated fluorescent objects between ≤ 200 nm and ≥ 40 nm diameter with a precision down to the few nm range; it also allows the axial positioning of such structures down to the 1 nm scale; combined with laterally structured illumination/SIM, a 3D localisation precision of ≤1 nm is expected using fluorescence yields typical for SMLM applications. Together with the nanosizing capability of SMI, this can be used to analyse macromolecular nuclear complexes with a resolution approaching that of cryoelectron microscopy.

Is RNA the working genome in eukaryotes ? The 60 year evolution of a conceptual challenge.

About 80 years ago, in 1943, after a century of biochemical and genetic research, DNA was established as the carrier of genetic information. At the onset of Molecular Biology around 1960, the genome of living organisms embodied 3 basic, still unknown paradigms: its composition, organisation and expression. Between 1980 and 1990, its replication was understood, and ideas about its 3D-organisation were suggested and finally confirmed by 2010. The basic mechanisms of gene expression in higher organisms, the synthesis of precursor RNAs and their processing into functional RNAs, were also discovered about 60 years ago in 1961/62. However, some aspects were then, and are still now debated, although the latest results in post-genomic research have confirmed the basic principles. When my history-essay was published in 2003, describing the discovery of RNA processing 40 years earlier, the main facts were not yet generally confirmed or acknowledged. The processing of pre-rRNA to 28 S and 18 S rRNA was clearly demonstrated, confirmed by others and generally accepted as a fact. However, the "giant" size of pre-mRNA 10-100 kb-long and pervasive DNA transcription were still to be confirmed by post-genomic methods. It was found, surprisingly, that up to 90% of DNA is transcribed in the life cycle of eukaryotic organisms thus showing that pervasive transcription was the general rule. In this essay, we shall take a journey through the 60-year history of evolving paradigms of gene expression which followed the emergence of Molecular Biology, and we will also evoke some of the "folklore" in research throughout this period. Most important was the growing recognition that although the genome is encoded in DNA, the Working Genome in eukaryotic organisms is RNA.

Spatial genome organization, TGFß, and biomolecular condensates: Do they talk during development?

Cis-regulatory elements govern gene expression programs to determine cell identity during development. Recently, the possibility that multiple enhancers are orchestrated in clusters of enhancers has been suggested. How these elements are arranged in the 3D space to control the activation of a specific promoter remains unclear. Our recent work revealed that the TGFß pathway drives the assembly of enhancer clusters and precise gene activation during neurogenesis. We discovered that the TGFß pathway coactivator JMJD3 was essential in maintaining these structures in the 3D space. To do that, JMJD3 required an intrinsically disordered region involved in forming phase-separated biomolecular condensates found in the enhancer clusters. Our data support the existence of a relationship between 3D-conformation of the chromatin, biomolecular condensates, and TGFß-driven response during mammalian neurogenesis. In this review, we discuss how signaling (TGFß), epigenetics (JMJD3), and biochemical properties (biomolecular condensates nucleation) are coordinated to modulate the genome structure to guarantee proper neural development. Moreover, we comment on the potential underlying mechanisms and implications of the enhancer-mediated regulation. Finally, we point out the knowledge gaps that still need to be addressed.

Plastome structure, phylogeny and evolution of plastid genes in Reevesia (Helicteroideae, Malvaceae).

Reevesia is an eastern Asian-eastern North American disjunction genus in the family Malvaceae s.l. and comprises approximately 25 species. The relationships within the genus are not well understood. Here, 15 plastomes representing 12 Reevesia species were compared, with the aim of better understanding the species circumscription and phylogenetic relationships within the genus and among genera in the family Malvaceae s.l. The 11 newly sequenced plastomes range between 161,532 and 161, 945 bp in length. The genomes contain 114 unique genes, 18 of which are duplicated in the inverted repeats (IRs). Gene content of these plastomes is nearly identical. All the protein-coding genes are under purifying selection in the Reevesia plastomes compared. The top ten hypervariable regions, SSRs, and the long repeats identified are potential molecular markers for future population genetic and phylogenetic studies. Phylogenetic analysis based on the whole plastomes confirmed the monophyly of Reevesia and a close relationship with Durio (traditional Bombacaceae) in subfamily Helicteroideae, but not with the morphologically similar genera Pterospermum and Sterculia (both of traditional Sterculiaceae). Phylogenetic relationships within Reevesia suggested that two species, R. pubescens and R. thyrsoidea, as newly defined, are not monophyletic. Six taxa, R. membranacea, R. xuefengensis, R. botingensis, R. lofouensis, R. longipetiolata and R. pycnantha, are suggested to be recognized.

Coming full circle: On the origin and evolution of the looping model for enhancer-promoter communication.

In mammalian organisms, enhancers can regulate transcription from great genomic distances. How enhancers affect distal gene expression has been a major question in the field of gene regulation. One model to explain how enhancers communicate with their target promoters, the chromatin looping model, posits that enhancers and promoters come in close spatial proximity to mediate communication. Chromatin looping has been broadly accepted as a means for enhancer-promoter communication, driven by accumulating in vitro and in vivo evidence. The genome is now known to be folded into a complex 3D arrangement, created and maintained in part by the interplay of the Cohesin complex and the DNA-binding protein CTCF. In the last few years, however, doubt over the relationship between looping and transcriptional activation has emerged, driven by studies finding that only a modest number of genes are perturbed with acute degradation of looping machinery components. In parallel, newer models describing distal enhancer action have also come to prominence. In this article, we explore the emergence and development of the looping model as a means for enhancer-promoter communication and review the contrasting evidence between historical gene-specific and current global data for the role of chromatin looping in transcriptional regulation. We also discuss evidence for alternative models to chromatin looping and their support in the literature. We suggest that, while there is abundant evidence for chromatin looping as a major mechanism for enhancer function, enhancer-promoter communication is likely mediated by more than one mechanism in an enhancer- and context-dependent manner.