Ekaterina E. Khrameeva

Active chromatin and transcription play a key role in chromosome partitioning into topologically associating domains

Recent advances enabled by the Hi-C technique have unraveled many principles of chromosomal folding that were subsequently linked to disease and gene regulation. In particular, Hi-C revealed that chromosomes of animals are organized into topologically associating domains (TADs), evolutionary conserved compact chromatin domains that influence gene expression. Mechanisms that underlie partitioning of the genome into TADs remain poorly understood. To explore principles of TAD folding in Drosophila melanogaster, we performed Hi-C and poly(A)(+) RNA-seq in four cell lines of various origins (S2, Kc167, DmBG3-c2, and OSC). Contrary to previous studies, we find that regions between TADs (i.e., the inter-TADs and TAD boundaries) in Drosophila are only weakly enriched with the insulator protein dCTCF, while another insulator protein Su(Hw) is preferentially present within TADs. However, Drosophila inter-TADs harbor active chromatin and constitutively transcribed (housekeeping) genes. Accordingly, we find that binding of insulator proteins dCTCF and Su(Hw) predicts TAD boundaries much worse than active chromatin marks do. Interestingly, inter-TADs correspond to decompacted inter-bands of polytene chromosomes, whereas TADs mostly correspond to densely packed bands. Collectively, our results suggest that TADs are condensed chromatin domains depleted in active chromatin marks, separated by regions of active chromatin. We propose the mechanism of TAD self-assembly based on the ability of nucleosomes from inactive chromatin to aggregate, and lack of this ability in acetylated nucleosomal arrays. Finally, we test this hypothesis by polymer simulations and find that TAD partitioning may be explained by different modes of inter-nucleosomal interactions for active and inactive chromatin.

RBM24 is a major regulator of muscle-specific alternative splicing.

Cell-type-specific splicing generates numerous alternatively spliced transcripts playing important roles for organ development and homeostasis, but only a few tissue-specific splicing factors have been identified. We found that RBM24 governs a large number of muscle-specific splicing events that are critically involved in cardiac and skeletal muscle development and disease. Targeted inactivation of RBM24 in mice disrupted cardiac development and impaired sarcomerogenesis in striated muscles. In vitro splicing assays revealed that recombinant RBM24 is sufficient to promote muscle-specific exon inclusion in nuclear extracts of nonmuscle cells. Furthermore, we demonstrate that binding of RBM24 to an intronic splicing enhancer (ISE) is essential and sufficient to overcome repression of exon inclusion by an exonic splicing silencer (ESS) containing PTB and hnRNP A1/A2 binding sites. Introduction of ESS and ISE converted a constitutive exon into an RMB24-dependent alternative exon. We reason that RBM24 is a major regulator of alternative splicing in striated muscles.

Crosslinking-immunoprecipitation (iCLIP) analysis reveals global regulatory roles of hnRNP L

Heterogeneous nuclear ribonucleoprotein L (hnRNP L) is a multifunctional RNA-binding protein that is involved in many different processes, such as regulation of transcription, translation, and RNA stability. We have previously characterized hnRNP L as a global regulator of alternative splicing, binding to CA-repeat, and CA-rich RNA elements. Interestingly, hnRNP L can both activate and repress splicing of alternative exons, but the precise mechanism of hnRNP L-mediated splicing regulation remained unclear. To analyze activities of hnRNP L on a genome-wide level, we performed individual-nucleotide resolution crosslinking-immunoprecipitation in combination with deep-sequencing (iCLIP-Seq). Sequence analysis of the iCLIP crosslink sites showed significant enrichment of C/A motifs, which perfectly agrees with the in vitro binding consensus obtained earlier by a SELEX approach, indicating that in vivo hnRNP L binding targets are mainly determined by the RNA-binding activity of the protein. Genome-wide mapping of hnRNP L binding revealed that the protein preferably binds to introns and 3′ UTR. Additionally, position-dependent splicing regulation by hnRNP L was demonstrated: The protein represses splicing when bound to intronic regions upstream of alternative exons, and in contrast, activates splicing when bound to the downstream intron. These findings shed light on the longstanding question of differential hnRNP L-mediated splicing regulation. Finally, regarding 3′ UTR binding, hnRNP L binding preferentially overlaps with predicted microRNA target sites, indicating global competition between hnRNP L and microRNA binding. Translational regulation by hnRNP L was validated for a subset of predicted target 3′UTRs.

Neanderthal ancestry drives evolution of lipid catabolism in contemporary Europeans

Although Neanderthals are extinct, fragments of their genomes persist in contemporary humans. Here we show that while the genome-wide frequency of Neanderthal-like sites is approximately constant across all contemporary out-of-Africa populations, genes involved in lipid catabolism contain more than threefold excess of such sites in contemporary humans of European descent. Evolutionally, these genes show significant association with signatures of recent positive selection in the contemporary European, but not Asian or African populations. Functionally, the excess of Neanderthal-like sites in lipid catabolism genes can be linked with a greater divergence of lipid concentrations and enzyme expression levels within this pathway, seen in contemporary Europeans, but not in the other populations. We conclude that sequence variants that evolved in Neanderthals may have given a selective advantage to anatomically modern humans that settled in the same geographical areas.

A Novel Intra-U1 snRNP Cross-Regulation Mechanism: Alternative Splicing Switch Links U1C and U1-70K Expression

The U1 small nuclear ribonucleoprotein (snRNP)-specific U1C protein participates in 5′ splice site recognition and regulation of pre-mRNA splicing. Based on an RNA-Seq analysis in HeLa cells after U1C knockdown, we found a conserved, intra-U1 snRNP cross-regulation that links U1C and U1-70K expression through alternative splicing and U1 snRNP assembly. To investigate the underlying regulatory mechanism, we combined mutational minigene analysis, in vivo splice-site blocking by antisense morpholinos, and in vitro binding experiments. Alternative splicing of U1-70K pre-mRNA creates the normal (exons 7–8) and a non-productive mRNA isoform, whose balance is determined by U1C protein levels. The non-productive isoform is generated through a U1C-dependent alternative 3′ splice site, which requires an adjacent cluster of regulatory 5′ splice sites and binding of intact U1 snRNPs. As a result of nonsense-mediated decay (NMD) of the non-productive isoform, U1-70K mRNA and protein levels are down-regulated, and U1C incorporation into the U1 snRNP is impaired. U1-70K/U1C-deficient particles are assembled, shifting the alternative splicing balance back towards productive U1-70K splicing, and restoring assembly of intact U1 snRNPs. Taken together, we established a novel feedback regulation that controls U1-70K/U1C homeostasis and ensures correct U1 snRNP assembly and function.

Genomic study of the Ket: a Paleo-Eskimo-related ethnic group with significant ancient North Eurasian ancestry

The Kets, an ethnic group in the Yenisei River basin, Russia, are considered the last nomadic hunter-gatherers of Siberia, and Ket language has no transparent affiliation with any language family. We investigated connections between the Kets and Siberian and North American populations, with emphasis on the Mal’ta and Paleo-Eskimo ancient genomes, using original data from 46 unrelated samples of Kets and 42 samples of their neighboring ethnic groups (Uralic-speaking Nganasans, Enets, and Selkups). We genotyped over 130,000 autosomal SNPs, identified mitochondrial and Y-chromosomal haplogroups, and performed high-coverage genome sequencing of two Ket individuals. We established that Nganasans, Kets, Selkups, and Yukaghirs form a cluster of populations most closely related to Paleo-Eskimos in Siberia (not considering indigenous populations of Chukotka and Kamchatka). Kets are closely related to modern Selkups and to some Bronze and Iron Age populations of the Altai region, with all these groups sharing a high degree of Mal’ta ancestry. Implications of these findings for the linguistic hypothesis uniting Ket and Na-Dene languages into a language macrofamily are discussed.

Spatial Proximity and Similarity of the Epigenetic State of Genome Domains

Recent studies demonstrate that the organization of the chromatin within the nuclear space might play a crucial role in the regulation of gene expression. The ongoing progress in determination of the 3D structure of the nuclear chromatin allows one to study correlations between spatial proximity of genome domains and their epigenetic state. We combined the data on three-dimensional architecture of the whole human genome with results of high-throughput studies of the chromatin functional state and observed that fragments of different chromosomes that are spatially close tend to have similar patterns of histone modifications, methylation state, DNAse sensitivity, expression level, and chromatin states in general. Moreover, clustering of genome regions by spatial proximity produced compact clusters characterized by the high level of histone modifications and DNAse sensitivity and low methylation level, and loose clusters with the opposite characteristics. We also associated the spatial proximity data with previously detected chimeric transcripts and the results of RNA-seq experiments and observed that the frequency of formation of chimeric transcripts from fragments of two different chromosomes is higher among spatially proximal genome domains. A fair fraction of these chimeric transcripts seems to arise post-transcriptionally via trans-splicing.

Lipidome determinants of maximal lifespan in mammals

Maximal lifespan of mammalian species, even if closely related, may differ more than 10-fold, however the nature of the mechanisms that determine this variability is unresolved. Here, we assess the relationship between maximal lifespan duration and concentrations of more than 20,000 lipid compounds, measured in 669 tissue samples from 6 tissues of 35 species representing three mammalian clades: primates, rodents and bats. We identify lipids associated with species’ longevity across the three clades, uncoupled from other parameters, such as basal metabolic rate, body size, or body temperature. These lipids clustered in specific lipid classes and pathways, and enzymes linked to them display signatures of greater stabilizing selection in long-living species, and cluster in functional groups related to signaling and protein-modification processes. These findings point towards the existence of defined molecular mechanisms underlying variation in maximal lifespan among mammals.

Inter-Chromosomal Contact Networks Provide Insights into Mammalian Chromatin Organization

The recent advent of conformation capture techniques has provided unprecedented insights into the spatial organization of chromatin. We present a large-scale investigation of the inter-chromosomal segment and gene contact networks in embryonic stem cells of two mammalian organisms: humans and mice. Both interaction networks are characterized by a high degree of clustering of genome regions and the existence of hubs. Both genomes exhibit similar structural characteristics such as increased flexibility of certain Y chromosome regions and co-localization of centromere-proximal regions. Spatial proximity is correlated with the functional similarity of genes in both species. We also found a significant association between spatial proximity and the co-expression of genes in the human genome. The structural properties of chromatin are also species specific, including the presence of two highly interactive regions in mouse chromatin and an increased contact density on short, gene-rich human chromosomes, thereby indicating their central nuclear position. Trans-interacting segments are enriched in active marks in human and had no distinct feature profile in mouse. Thus, in contrast to interactions within individual chromosomes, the inter-chromosomal interactions in human and mouse embryonic stem cells do not appear to be conserved.

Mutants of monomeric red fluorescent protein mRFP1 at residue 66: Structure modeling by molecular dynamics and search for correlations with spectral properties

To study the interrelation between the spectral and structural properties of fluorescent proteins, structures of mutants of monomeric red fluorescent protein mRFP1 with all possible point mutations of Glu66 (except replacement by Pro) were simulated by molecular dynamics. A global search for correlations between geometrical structure parameters and some spectral characteristics (absorption maximum wavelength, integral extinction coefficient at the absorption maximum, excitation maximum wavelength, emission maximum wavelength, and quantum yield) was performed for the chromophore and its 6 environment in mRFP1, Q66A, Q66L, Q66S, Q66C, Q66H, and Q66N. The correlation coefficients (0.81–0.87) were maximal for torsion angles in phenolic and imidazolidine rings as well as for torsion angles in the regions of connection between these rings and chromophore attachment to β-barrel. The data can be used to predict the spectral properties of fluorescent proteins based on their structures and to reveal promising positions for directed mutagenesis.