Giorgio Bernardi

Surfperches (Family Embiotocidae) are viviparous temperate reef fishes that brood their young. This life history trait translates into limited dispersal, strong population structure, and an unusually strong potential for local adaptation in a marine fish. As part of the California Conservation Genomics Project (CCGP), we sequenced the genome of the Black Surfperch, Embiotoca jacksoni, to establish a genomic model for understanding phylogeographic patterns of marine organisms in California. These patterns, in turn, may inform the design of marine protected areas using dispersal models based on genomic data. The genome of E. jacksoni is typical of marine fishes at less than 1Gb (genome size = 635 Mb), and our assembly is near-chromosome level (contig N50 = 6.5Mb, scaffold N50 = 15.5 Mb, BUSCO = 98.1%). Within the context of the CCGP, the genome will be used as a reference for future whole genome resequencing projects aimed at enhancing our knowledge of the population structure of the species, and efficacy of Marine Protected Areas across the state.

The chromatin of the human genome was analyzed at three DNA size levels. At the first, compartment level, two “gene spaces” were found many years ago: A GC-rich, gene-rich “genome core” and a GC-poor, gene-poor “genome desert”, the former corresponding to open chromatin centrally located in the interphase nucleus, the latter to closed chromatin located peripherally. This bimodality was later confirmed and extended by the discoveries (1) of LADs, the Lamina-Associated Domains, and InterLADs; (2) of two “spatial compartments”, A and B, identified on the basis of chromatin interactions; and (3) of “forests and prairies” characterized by high and low CpG islands densities. Chromatin compartments were shown to be associated with the compositionally different, flat and single- or multi-peak DNA structures of the two, GC-poor and GC-rich, “super-families” of isochores. At the second, sub-compartment, level, chromatin corresponds to flat isochores and to isochore loops (due to compositional DNA gradients) that are susceptible to extrusion. Finally, at the short-sequence level, two sets of sequences, GC-poor and GC-rich, define two different nucleosome spacings, a short one and a long one. In conclusion, chromatin structures are moulded according to a “genomic code” by DNA sequences that pervade the genome and leave no room for “junk”.

Recent investigations have revealed 1) that the isochores of the human genome group into two super‐families characterized by two different long‐range 3D structures, and 2) that these structures, essentially based on the distribution and topology of short sequences, mold primary chromatin domains (and define nucleosome binding). More specifically, GC‐poor, gene‐poor isochores are low‐heterogeneity sequences with oligo‐A spikes that mold the lamina‐associated domains (LADs), whereas GC‐rich, gene‐rich isochores are characterized by single or multiple GC peaks that mold the topologically associating domains (TADs). The formation of these “primary TADs” may be followed by extrusion under the action of cohesin and CTCF. Finally, the genomic code, which is responsible for the pervasive encoding and molding of primary chromatin domains (LADs and primary TADs, namely the “gene spaces”/“spatial compartments”) resolves the longstanding problems of “non‐coding DNA,” “junk DNA,” and “selfish DNA” leading to a new vision of the genome as shaped by DNA sequences.

The general model presented here for the formation of chromatin domains, LADs and TADs, is primarily based on the 3-D structures of the corresponding DNA sequences, the GC-poor and GC-rich isochores. Indeed, the low-heterogeneity GC-poor isochores locally are intrinsically stiff and curved because of the presence of interspersed oligo-Adenines. In contrast, the high-heterogeneity GC-rich isochores are in the shape of peaks characterized by increasing levels of GC and of interspersed oligo-Guanines. In LADs, oligo-Adenines induce local nucleosome depletions leading to structures that are well suited for the attachment to (and embedding in) the lamina. In TADs, the gradients of GC and of oligo-Guanines are responsible for a decreasing nucleosome density, decreasing supercoiling and increasing accessibility. This “moulding step” shapes the “primary TADs” into loops that lack self-interactions, being CTCF/cohesin-free structures. The cohesin complex then binds to the tips of “primary TADs” and slides down the loops, thanks to Nipbl, an essential factor for loading cohesin and for stimulating its ATPase activity and its translocation. This “extruding step” leads to closer contacts and to self-interactions in the loops and stops at the CTCF binding sites located at the base of the loops that are thus closed and insulated.

Recent investigations have shown that isochores are characterized by a 3-D structure which is primarily responsible for the topology of chromatin domains. More precisely, an analysis of human chromosome 21 demonstrated that low-heterogeneity, GC-poor isochores are characterized by the presence of oligo-Adenines that are intrinsically stiff, curved and unfavorable for nucleosome binding. This leads to a structure of the corresponding chromatin domains, the Lamina Associated Domains, or LADs, which is well suited for interaction with the lamina. In contrast, the high-heterogeneity GC-rich isochores are in the form of compositional peaks and valleys characterized by increasing gradients of oligo-Guanines in the peaks and oligo-Adenines in the valleys that lead to increasing nucleosome depletions in the corresponding chromatin domains, the Topological Associating Domains, or TADs. These results encouraged us to investigate in detail the di- and tri-nucleotide profiles of 100 Kb segments of chromosome 21, as well as those of the di- to octa-Adenines and di- to octa-Guanines in some representative regions of the chromosome. The results obtained show that the 3-D structures of isochores and chromatin domains depend not only upon oligo-Adenines and oligo-Guanines but also, to a lower but definite extent, upon the majority of di- and tri-nucleotides. This conclusion has strong implications for the biological role of non-coding sequences.

Recent investigations have shown that isochores are characterized by a 3-D structure which is primarily responsible for the topology of chromatin domains. More precisely, an analysis of human chromosome 21 demonstrated that GC-poor isochores are low-heterogeneity sequences characterized by the presence of oligo-Adenines that are intrinsically stiff, curved and unfavorable for nucleosome binding. This leads to a structure of the corresponding chromatin domains, the Lamina Associated Domains, or LADs, which is well suited for interaction with lamina. In contrast, the high-heteorogeneity GC-rich isochores are in the form of compositional peaks characterized by gradients of oligo-Guanines that lead to increasing nucleosome depletions in the corresponding chromatin domains, the Topological Associating Domains, or TADs. These results encouraged us to investigate in detail the di- and tri-nucleotide profiles of 100Kb segments of chromosome 21, as well as those of the di- to octa-Adenines and di- to octa-Guanines in several regions of the chromosome. The results obtained show that the 3-D structures of isochores and chromatin domains depend not only upon oligo-Adenines and oligo-Guanines but also, to a lower but definite extent, upon the majority of di- and tri-nucleotides. This conclusion, which applies to all human chromosome, has strong implications for the biological role of non-coding sequences.

In spite of the recent advances in the field of chromatin architecture ^1,2 , the formation mechanism of chromatin domains, TADs, the topologically associating domains, and LADs, the lamina associated domains, is still an open problem. While previous models only dealt with TADs and essentially relied on the architectural proteins CTCF and cohesin, the model presented here concerns both TADs and LADs and is primarily based on the corresponding DNA sequences, the GC-rich and GC-poor isochores, more specifically on their newly discovered 3-D structures. Indeed, the compositionally homogeneous GC-poor isochores were shown to be locally stiff because of the presence of interspersed oligo- Adenines ^4,5 , whereas the compositionally heterogeneous GC-rich isochores were found to be peak-shaped and characterized by increasing gradients of GC and of interspersed oligo- Guanines. In LADs, oligo-Adenines induce local nucleosome depletions ^4,5 that are responsible for a wavy structure well adapted for interaction with the lamina. In TADs, the increasing GC levels and increasing oligo-Guanines of the isochore peaks are responsible for a decreasing nucleosome density ^5,6 , a decreasing supercoiling ⁷ and an increasing accessibility ⁸ . These factors mould the loops of “primary TADs”, that lack self-interactions since they are CTCF/cohesin-free, yet transcriptionally functional structures ^9-11 . This “moulding step” is followed by a second step, in which the cohesin rings bind to the tips of the “primary TADs” and slide down the loops. This process is very likely due to Scc2/Nipbl, an essential factor not only for loading cohesin, but also for stimulating its translocation ¹² and its ATPase activity ¹³ . This “sliding step” creates self-interactions in the loops and stops at the CTCF binding sites located at the base of the loops that are thus closed and insulated.

A recent investigation showed the existence of correlations between the architectural features of mammalian interphase chromosomes and the compositional properties of isochores. This result prompted us to compare maps of the Topologically Associating Domains (TADs) and of the Lamina Associated Domains (LADs) with the corresponding isochore maps of mouse and human chromosomes. This approach revealed that: 1) TADs and LADs correspond to isochores, i.e., isochores are the genomic units that underlie chromatin domains; 2) the conservation of TADs and LADs in mammalian genomes is explained by the evolutionary conservation of isochores; 3) chromatin domains corresponding to GC-poor isochores (e.g., LADs) show not only self-interactions but also intrachromosomal interactions with other domains also corresponding to GC-poor isochores even if located far away; in contrast, chromatin domains corresponding to GC-rich isochores (e.g., TADs) show more localized chromosomal interactions, many of which are inter-chromosomal. In conclusion, this investigation establishes a link between DNA sequences and chromatin architecture, explains the evolutionary conservation of TADs and LADs and provides new information on the spatial distribution of GC-poor/gene-poor and GC-rich/gene-rich chromosomal regions in the interphase nucleus.

The formation of mammalian chromatin domains was investigated by analyzing the domain/isochore connection. This showed that LADs correspond to GC-poor isochores and are compositionally flat, flexible chromatin structures because of the local nucleosome depletions associated with the presence of oligo-A’s. In contrast, TADs correspond to GC-rich isochores that consist of single or (much more frequently) multiple, GC peaks that shape the single or multiple, loops of TADs. Indeed, the increasing nucleosome depletions linked to the GC gradients of isochore peaks lead to an increasing chromatin flexibility (accompanied by an increasing accessibility and decreasing supercoiling). In conclusion, isochores not only encode but also mold chromatin architecture; while architectural proteins play a role in closing and insulating TAD loops. An extension of this model concerns the encoding of open and closed chromosome compartments by alternating GC-rich and GC-poor isochores, the interactions among compartments defining the 3-D chromosome folding.

A recent investigation showed the existence of correlations between the architectural features of mammalian interphase chromosomes and the compositional properties of isochores. This result prompted us to compare maps of the Topologically Associating Domains (TADs) and of the Lamina Associated Domains (LADs) with the corresponding isochore maps of mouse and human chromosomes. This approach revealed that: 1) TADs and LADs correspond to isochores, i.e., isochores are the genomic units that underlie chromatin domains; 2) the conservation of TADs and LADs in mammalian genomes is explained by the evolutionary conservation of isochores; 3) chromatin domains corresponding to GC-poor isochores ( e.g., LADs) interact with other domains also corresponding to GC-poor isochores even if located far away on the chromosomes; in contrast, chromatin domains corresponding to GC-rich isochores ( e.g., TADs) show more localized chromosomal interactions, many of which are inter-chromosomal. In conclusion, this investigation establishes a link between DNA sequences and chromatin architecture, explains the evolutionary conservation of TADs and LADs and provides new information on the spatial distribution of GC-poor/gene-poor and GC-rich/gene-rich chromosomal regions in the interphase nucleus.