Evidence of influence of genomic DNA sequence on human X chromosome inactivation

X chromosome inactivation and embryonic stem cells

X chromosome inactivation (XCI) is a process required to equalize the dosage of X-encoded genes between female and male cells. XCI is initiated very early during female embryonic development or upon differentiation of female embryonic stem (ES) cells and results in inactivation of one X chromosome in every female somatic cell. The regulation of XCI involves factors that also play a crucial role in ES cell maintenance and differentiation and the XCI process therefore provides a beautiful paradigm to study ES cell biology. In this chapter we describe the important cis and trans acting regulators of XCI and introduce the models that have been postulated to explain initiation of XCI in female cells only. We also discuss the proteins involved in the establishment of the inactive X chromosome and describe the different chromatin modifications associated with the inactivation process. Finally, we describe the potential of mouse and human ES and induced pluripotent stem (iPS) cells as model systems to study the XCI process.

Characterisation of inactivation domains and evolutionary strata in human X chromosome through Markov segmentation

Markov segmentation is a method of identifying compositionally different subsequences in a given symbolic sequence. We have applied this technique to the DNA sequence of the human X chromosome to analyze its compositional structure. The human X chromosome is known to have acquired DNA through distinct evolutionary events and is believed to be composed of five evolutionary strata. In addition, in female mammals all copies of X chromosome in excess of one are transcriptionally inactivated. The location of a gene is correlated with its ability to undergo inactivation, but correlations between evolutionary strata and inactivation domains are less clear. Our analysis provides an accurate estimate of the location of stratum boundaries and gives a high-resolution map of compositionally different regions on the X chromosome. This leads to the identification of a novel stratum, as well as segments wherein a group of genes either undergo inactivation or escape inactivation in toto. We identify oligomers that appear to be unique to inactivation domains alone.

Genomic analyses of sex chromosome evolution

Sex chromosomes have evolved multiple times in many taxa. The recent explosion in the availability of whole genome sequences from a variety of organisms makes it possible to investigate sex chromosome evolution within and across genomes. Comparative genomic studies have shown that quite distant species may share fundamental properties of sex chromosome evolution, while very similar species can evolve unique sex chromosome systems. Furthermore, within-species genomic analyses can illuminate chromosome-wide sequence and expression polymorphisms. Here, we explore recent advances in the study of vertebrate sex chromosomes achieved using genomic analyses.

Genomic features that predict allelic imbalance in humans suggest patterns of constraint on gene expression variation

Variation in gene expression is an important contributor to phenotypic diversity within and between species. Although this variation often has a genetic component, identification of the genetic variants driving this relationship remains challenging. In particular, measurements of gene expression usually do not reveal whether the genetic basis for any observed variation lies in cis or in trans to the gene, a distinction that has direct relevance to the physical location of the underlying genetic variant, and which may also impact its evolutionary trajectory. Allelic imbalance measurements identify cis-acting genetic effects by assaying the relative contribution of the two alleles of a cis-regulatory region to gene expression within individuals. Identification of patterns that predict commonly imbalanced genes could therefore serve as a useful tool and also shed light on the evolution of cis-regulatory variation itself. Here, we show that sequence motifs, polymorphism levels, and divergence levels around a gene can be used to predict commonly imbalanced genes in a human data set. Reduction of this feature set to four factors revealed that only one factor significantly differentiated between commonly imbalanced and nonimbalanced genes. We demonstrate that these results are consistent between the original data set and a second published data set in humans obtained using different technical and statistical methods. Finally, we show that variation in the single allelic imbalance-associated factor is partially explained by the density of genes in the region of a target gene (allelic imbalance is less probable for genes in gene-dense regions), and, to a lesser extent, the evenness of expression of the gene across tissues and the magnitude of negative selection on putative regulatory regions of the gene. These results suggest that the genomic distribution of functional cis-regulatory variants in the human genome is nonrandom, perhaps due to local differences in evolutionary constraint.

The role of LINEs and CpG islands in dosage compensation on the chicken Z chromosome

Most avian Z genes are expressed more highly in ZZ males than ZW females, suggesting that chromosome-wide mechanisms of dosage compensation have not evolved. Nevertheless, a small percentage of Z genes are expressed at similar levels in males and females, an indication that a yet unidentified mechanism compensates for the sex difference in copy number. Primary DNA sequences are thought to have a role in determining chromosome gene inactivation status on the mammalian X chromosome. However, it is currently unknown whether primary DNA sequences also mediate chicken Z gene compensation status. Using a combination of chicken DNA sequences and Z gene compensation profiles of 310 genes, we explored the relationship between Z gene compensation status and primary DNA sequence features. Statistical analysis of different Z chromosomal features revealed that long interspersed nuclear elements (LINEs) and CpG islands are enriched on the Z chromosome compared with 329 other DNA features. Linear support vector machine (SVM) classifiers, using primary DNA sequences, correctly predict the Z compensation status for >60% of all Z-linked genes. CpG islands appear to be the most accurate classifier and alone can correctly predict compensation of 63% of Z genes. We also show that LINE CR1 elements are enriched 2.7-fold on the chicken Z chromosome compared with autosomes and that chicken chromosomal length is highly correlated with percentage LINE content. However, the position of LINE elements is not significantly associated with dosage compensation status of Z genes. We also find a trend for a higher proportion of CpG islands in the region of the Z chromosome with the fewest dosage-compensated genes compared with the region containing the greatest concentration of compensated genes. Comparison between chicken and platypus genomes shows that LINE elements are not enriched on sex chromosomes in platypus, indicating that LINE accumulation is not a feature of all sex chromosomes. Our results suggest that CpG islands are not randomly distributed on the Z chromosome and may influence Z gene dosage compensation status.

X chromosome inactivation and Xist evolution in a rodent lacking LINE-1 activity

Dosage compensation in eutherian mammals occurs by inactivation of one X chromosome in females. Silencing of that X chromosome is initiated by Xist, a large non-coding RNA, whose coating of the chromosome extends in cis from the X inactivation center. LINE-1 (L1) retrotransposons have been implicated as possible players for propagation of the Xist signal, but it has remained unclear whether they are essential components. We previously identified a group of South American rodents in which L1 retrotransposition ceased over 8 million years ago and have now determined that at least one species of these rodents, Oryzomys palustris, still retains X inactivation. We have also isolated and analyzed the majority of the Xist RNA from O. palustris and a sister species retaining L1 activity, Sigmodon hispidus, to determine if evolution in these sequences has left signatures that might suggest a critical role for L1 elements in Xist function. Comparison of rates of Xist evolution in the two species fails to support L1 involvement, although other explanations are possible. Similarly, comparison of known repeats and potential RNA secondary structures reveals no major differences with the exception of a new repeat in O. palustris that has potential to form new secondary structures.

X-chromosome upregulation and inactivation: two sides of the dosage compensation mechanism in mammals

Mammals have a very complex, tightly controlled, and developmentally regulated process of dosage compensation. One form of the process equalizes expression of the X-linked genes, present as a single copy in males (XY) and as two copies in females (XX), by inactivation of one of the two X-chromosomes in females. The second form of the process leads to balanced expression between the X-linked and autosomal genes by transcriptional upregulation of the active X in males and females. However, not all X-linked genes are absolutely balanced. This review is focused on the recent advances in studying the dosage compensation phenomenon in mammals.

The impact of local genome sequence on defining heterochromatin domains

Characterizing how genomic sequence interacts with trans-acting regulatory factors to implement a program of gene expression in eukaryotic organisms is critical to understanding genome function. One means by which patterns of gene expression are achieved is through the differential packaging of DNA into distinct types of chromatin. While chromatin state exerts a major influence on gene expression, the extent to which cis-acting DNA sequences contribute to the specification of chromatin state remains incompletely understood. To address this, we have used a fission yeast sequence element (L5), known to be sufficient to nucleate heterochromatin, to establish de novo heterochromatin domains in the Schizosaccharomyces pombe genome. The resulting heterochromatin domains were queried for the presence of H3K9 di-methylation and Swi6p, both hallmarks of heterochromatin, and for levels of gene expression. We describe a major effect of genomic sequences in determining the size and extent of such de novo heterochromatin domains. Heterochromatin spreading is antagonized by the presence of genes, in a manner that can occur independent of strength of transcription. Increasing the dosage of Swi6p results in increased heterochromatin proximal to the L5 element, but does not result in an expansion of the heterochromatin domain, suggesting that in this context genomic effects are dominant over trans effects. Finally, we show that the ratio of Swi6p to H3K9 di-methylation is sequence-dependent and correlates with the extent of gene repression. Taken together, these data demonstrate that the sequence content of a genomic region plays a significant role in shaping its response to encroaching heterochromatin and suggest a role of DNA sequence in specifying chromatin state.

Epigenetic packaging of DNA sequence into chromatin is a major force in shaping the function of complex genomes. Different types of chromatin have distinct effects on gene expression, and thus chromatin state imparts distinct features on the associated genomic DNA. Our study focuses on the transition between two opposing chromatin states: euchromatin, which generally correlates with gene expression, and heterochromatin, which is typically refractive to gene expression. While heterochromatin is capable of spreading into euchromatic domains, the parameters that influence such spreading are unknown. We established heterochromatin at ectopic sites in the genome and evaluated whether specific DNA sequences affected the extent of heterochromatin spreading and the transition between heterochromatin and euchromatin. We found that the nature of the genomic DNA neighboring the heterochromatic sequence dramatically affected the extent of heterochromatin spreading. In particular, the presence of genes antagonized the spread of heterochromatin, whereas neutral sequence elements were incorporated into the domain. This study demonstrates that genome sequence and chromatin identity are inextricably linked; features of both interact to determine the structural and functional fate of underlying DNA sequences.

EpiGRAPH: user-friendly software for statistical analysis and prediction of (epi)genomic data

The EpiGRAPH web service http://epigraph.mpi-inf.mpg.de/ enables biologists to uncover hidden associations in vertebrate genome and epigenome datasets. Users can upload sets of genomic regions and EpiGRAPH will test multiple attributes (including DNA sequence, chromatin structure, epigenetic modifications and evolutionary conservation) for enrichment or depletion among these regions. Furthermore, EpiGRAPH learns to predictively identify similar genomic regions. This paper demonstrates EpiGRAPH's practical utility in a case study on monoallelic gene expression and describes its novel approach to reproducible bioinformatic analysis.

[Different effects on reporter gene expression by distinct L1-ORF2 segements]

An intact L1 define element is 6 kb in length in human genome. The majority of the L1s is truncated and has direction difference, implying that it is interesting to study the effects of different length and directions of L1s on gene. In this work, 7 different segments were obtained from L1-open reading frame 2 (ORF2), each of which was 280 bp in length. Each segment was connected into 8 repeats in head and tail tandem manner and was inserted to downstream of GFP gene in different directions in pEGFP-C1. The inserted ORF2 segments in the sense orientation caused much stronger inhibition on gene transcription and protein expression than antisense sequences did. Among all segments, the first and ninth 280 bp segments of ORF2 in both orientations induced weaker inhibition on gene transcription than other segments in the same orientations and did not induce transcriptional elongation. The distribution of Alu in most regions of genome was inverse ratio with L1. The inserted Alus in both orientations inhibited GFP gene expression, but the inhibition in antisense orientation was stronger than that in sense orientation and the sense Alu was the sequence of inducing transcription elongation. A-rich of ORF2 was probably the molecular basis of its sense orientation with stronger inhibition on gene expression.

Dosage compensation and gene expression on the mammalian X chromosome: one plus one does not always equal two

Counting chromosomes is not just simple math. Although normal males and females differ in sex chromosome content (XY vs. XX), X chromosome imbalance is tolerated because dosage compensation mechanisms have evolved to ensure functional equivalence. In mammals this is accomplished by two processes--X chromosome inactivation that silences most genes on one X chromosome in females, leading to functional X monosomy for most genes in both sexes, and X chromosome upregulation that results in increased gene expression on the single active X in males and females, equalizing dosage relative to autosomes. This review focuses on genes on the X chromosome, and how gene content, organization and expression levels can be influenced by these two processes. Special attention is given to genes that are not X inactivated, and are not necessarily fully dosage compensated. These genes that "escape" X inactivation are of medical importance as they explain phenotypes in individuals with sex chromosome aneuploidies and may impact normal traits and disorders that differ between men and women. Moreover, escape genes give insight into how X chromosome inactivation is spread and maintained on the X.

X chromosome inactivation: heterogeneity of heterochromatin

The silent X chromosome in mammalian females is a classic example of facultative heterochromatin, the term highlighting the compacted and inactive nature of the chromosome. However, it is now clear that the heterochromatin of the inactive X is not homogeneous--as indeed, not all genes on the inactive X are silenced. We summarize known features and events of X inactivation in different mouse and human model systems, and highlight the heterogeneity of chromatin along the inactive X. Characterizing this heterogeneity is likely to provide insight into the cis-acting sequences involved in X chromosome inactivation.

Escape from X chromosome inactivation is an intrinsic property of the Jarid1c locus

Although most genes on one X chromosome in mammalian females are silenced by X inactivation, some "escape" X inactivation and are expressed from both active and inactive Xs. How these escape genes are transcribed from a largely inactivated chromosome is not fully understood, but underlying genomic sequences are likely involved. We developed a transgene approach to ask whether an escape locus is autonomous or is instead influenced by X chromosome location. Two BACs carrying the mouse Jarid1c gene and adjacent X-inactivated transcripts were randomly integrated into mouse XX embryonic stem cells. Four lines with single-copy, X-linked transgenes were identified, and each was inserted into regions that are normally X-inactivated. As expected for genes that are normally subject to X inactivation, transgene transcripts Tspyl2 and Iqsec2 were X-inactivated. However, allelic expression and RNA/DNA FISH indicate that transgenic Jarid1c escapes X inactivation. Therefore, transgenes at 4 different X locations recapitulate endogenous inactive X expression patterns. We conclude that escape from X inactivation is an intrinsic feature of the Jarid1c locus and functionally delimit this escape domain to the 112-kb maximum overlap of the BACs tested. Additionally, although extensive chromatin differences normally distinguish active and inactive loci, unmodified BACs direct proper inactive X expression patterns, establishing that primary DNA sequence alone, in a chromosome position-independent manner, is sufficient to determine X chromosome inactivation status. This transgene approach will enable further dissection of key elements of escape domains and allow rigorous testing of specific genomic sequences on inactive X expression.

Analysis of transposon interruptions suggests selection for L1 elements on the X chromosome

It has been hypothesised that the massive accumulation of L1 transposable elements on the X chromosome is due to their function in X inactivation, and that the accumulation of Alu elements near genes is adaptive. We tested the possible selective advantage of these two transposable element (TE) families with a novel method, interruption analysis. In mammalian genomes, a large number of TEs interrupt other TEs due to the high overall abundance and age of repeats, and these interruptions can be used to test whether TEs are selectively neutral. Interruptions of TEs, which are beneficial for the host, are expected to be deleterious and underrepresented compared with neutral ones. We found that L1 elements in the regions of the X chromosome that contain the majority of the inactivated genes are significantly less frequently interrupted than on the autosomes, while L1s near genes that escape inactivation are interrupted with higher frequency, supporting the hypothesis that L1s on the X chromosome play a role in its inactivation. In addition, we show that TEs are less frequently interrupted in introns than in intergenic regions, probably due to selection against the expansion of introns, but the insertion pattern of Alus is comparable to other repeats.

Recent experimental findings (for example the ENCODE project) show that many functional non-coding regions of genomes are not conserved across species, making the in-silico discovery of such regions challenging. Transposable elements (TEs), which represent 45 percent of the human genome and typically show no sequence conservation, are particularly intriguing from this point of view, because the highly nonrandom genomic distribution of many TE families in genomes has led to hypotheses that their presence is adaptive and have an epigenetic (regulatory) function. We use a novel approach based on the analysis of interrupted TEs to investigate if repeats are under selection that does not rely on sequence conservation. L1 elements, the most active transposable elements of the human genome, are highly overrepresented on the X-chromosome and were suggested to enhance its inactivation in mammals. We find that the interruption pattern of L1 repeats indicates a function for L1 elements in the inactivation of the mammalian X chromosome. Additionally, we show that a considerable fraction of TEs in introns are under selection for integrity, possibly due to selection on intron size or on TEs themselves.

Sequence information encoded in DNA that may influence long-range chromatin structure correlates with human chromosome functions

Little is known about the possible function of the bulk of the human genome. We have recently shown that long-range regular oscillation in the motif non-T, A/T, G (VWG) existing at ten-nucleotide multiples influences large-scale nucleosome array formation. In this work, we have determined the locations of all 100 kb regions that are predicted to form distinctive chromatin structures throughout each human chromosome (except Y). Using these data, we found that a significantly greater fraction of 300 kb sequences lacked annotated transcripts in genomic DNA regions > or = 300 kb that contained nearly continuous chromatin organizing signals than in control regions. We also found a relationship between the meiotic recombination frequency and the presence of strong VWG chromatin organizing signals. Large (> or = 300 kb) genomic DNA regions having low average recombination frequency are enriched in chromatin organizing signals. As additional controls, we show using chromosome 1 that the VWG motif signals are not enriched in randomly selected DNA regions having the mean size of the recombination coldspots, and that non-VWG motif sets do not generate signals that are enriched in recombination coldspots. We also show that tandemly repeated alpha satellite DNA contains strong VWG signals for the formation of distinctive nucleosome arrays, consistent with the low recombination activity of centromeres. Our correlations cannot be explained simply by variations in the GC content. Our findings suggest that a specific set of periodic DNA motifs encoded in genomic DNA, which provide signals for chromatin organization, influence human chromosome function.

Computational epigenetics

Epigenetic research aims to understand heritable gene regulation that is not directly encoded in the DNA sequence. Epigenetic mechanisms such as DNA methylation and histone modifications modulate the packaging of the DNA in the nucleus and thereby influence gene expression. Patterns of epigenetic information are faithfully propagated over multiple cell divisions, which makes epigenetic regulation a key mechanism for cellular differentiation and cell fate decisions. In addition, incomplete erasure of epigenetic information can lead to complex patterns of non-Mendelian inheritance. Stochastic and environment-induced epigenetic defects are known to play a major role in cancer and ageing, and they may also contribute to mental disorders and autoimmune diseases. Recent technical advances such as ChIP-on-chip and ChIP-seq have started to convert epigenetic research into a high-throughput endeavor, to which bioinformatics is expected to make significant contributions. Here, we review pioneering computational studies that have contributed to epigenetic research. In addition, we give a brief introduction into epigenetics-targeted at bioinformaticians who are new to the field-and we outline future challenges in computational epigenetics.

Global patterns of sequence evolution in Drosophila

Background

Sequencing of the genomes of several Drosophila allows for the first precise analyses of how global sequence patterns change among multiple, closely related animal species. A basic question is whether there are characteristic features that differentiate chromosomes within a species or between different species.

Results

We explored the euchromatin of the chromosomes of seven Drosophila species to establish their global patterns of DNA sequence diversity. Between species, differences in the types and amounts of simple sequence repeats were found. Within each species, the autosomes have almost identical oligonucleotide profiles. However, X chromosomes and autosomes have, in all species, a qualitatively different composition. The X chromosomes are less complex than the autosomes, containing both a higher amount of simple DNA sequences and, in several cases, chromosome-specific repetitive sequences. Moreover, we show that the right arm of the X chromosome of Drosophila pseudoobscura, which evolved from an autosome 10 – 18 millions of years ago, has a composition which is identical to that of the original, left arm of the X chromosome.

Conclusion

The consistent differences among species, differences among X chromosomes and autosomes and the convergent evolution of X and neo-X chromosomes demonstrate that strong forces are acting on drosophilid genomes to generate peculiar chromosomal landscapes. We discuss the relationships of the patterns observed with differential recombination and mutation rates and with the process of dosage compensation.

Faced with inequality: chicken do not have a general dosage compensation of sex-linked genes

Background

The contrasting dose of sex chromosomes in males and females potentially introduces a large-scale imbalance in levels of gene expression between sexes, and between sex chromosomes and autosomes. In many organisms, dosage compensation has thus evolved to equalize sex-linked gene expression in males and females. In mammals this is achieved by X chromosome inactivation and in flies and worms by up- or down-regulation of X-linked expression, respectively. While otherwise widespread in systems with heteromorphic sex chromosomes, the case of dosage compensation in birds (males ZZ, females ZW) remains an unsolved enigma.

Results

Here, we use a microarray approach to show that male chicken embryos generally express higher levels of Z-linked genes than female birds, both in soma and in gonads. The distribution of male-to-female fold-change values for Z chromosome genes is wide and has a mean of 1.4–1.6, which is consistent with absence of dosage compensation and sex-specific feedback regulation of gene expression at individual loci. Intriguingly, without global dosage compensation, the female chicken has significantly lower expression levels of Z-linked compared to autosomal genes, which is not the case in male birds.

Conclusion

The pronounced sex difference in gene expression is likely to contribute to sexual dimorphism among birds, and potentially has implication to avian sex determination. Importantly, this report, together with a recent study of sex-biased expression in somatic tissue of chicken, demonstrates the first example of an organism with a lack of global dosage compensation, providing an unexpected case of a viable system with large-scale imbalance in gene expression between sexes.

The disappearing Barr body in breast and ovarian cancers

Interest has recently reawakened in whether loss of the heterochromatic X chromosome (Barr body) is prevalent in certain breast and ovarian cancers, and new insights into the mechanisms involved have emerged. Mitotic segregation errors commonly explain the loss of the inactive X chromosome (Xi), but compromise of Xi heterochromatin in some cancers may signal broader deficits of nuclear heterochromatin. The debated link between BRCA1 and Xi might reflect a general relationship between BRCA1 and heterochromatin, which could connect BRCA1 to both epigenetic and genetic instability. We suggest that heterochromatic instability is a common but largely unexplored mechanism, leading to widespread genomic misregulation and the evolution of some cancers.

Sex chromosomes of basal placental mammals

Placental (eutherian) mammals are currently classified into four superordinal clades (Afrotheria, Xenarthra, Laurasiatheria and Supraprimates) of which one, the Afrotheria (a unique lineage of African origin), is generally considered to be basal. Therefore, Afrotheria provide a pivotal evolutionary link for studying fundamental differences between the sex chromosomes of human/mouse (both representatives of Supraprimates and the index species for studies of sex chromosomes) and those of the distantly related marsupials. In this study, we use female fibroblasts to investigate classical features of X chromosome inactivation including replication timing of the X chromosomes and Barr body formation. We also examine LINE-1 accumulation on the X chromosomes of representative afrotherians and look for evidence of a pseudoautosomal region (PAR). Our results demonstrate that asynchronous replication of the X chromosomes is common to Afrotheria, as with other mammals, and Barr body formation is observed across all Placentalia, suggesting that mechanisms controlling this evolved before their radiation. Finally, we provide evidence of a PAR (which marsupials lack) and demonstrate that LINE1 is accumulated on the afrotherian and xenarthran X, although this is probably not due to transposition events in a common ancestor, but rather ongoing selection to retain recently inserted LINE1 on the X.

Genomic environment predicts expression patterns on the human inactive X chromosome

What genomic landmarks render most genes silent while leaving others expressed on the inactive X chromosome in mammalian females? To date, signals determining expression status of genes on the inactive X remain enigmatic despite the availability of complete genomic sequences. Long interspersed repeats (L1s), particularly abundant on the X, are hypothesized to spread the inactivation signal and are enriched in the vicinity of inactive genes. However, both L1s and inactive genes are also more prevalent in ancient evolutionary strata. Did L1s accumulate there because of their role in inactivation or simply because they spent more time on the rarely recombining X? Here we utilize an experimentally derived inactivation profile of the entire human X chromosome to uncover sequences important for its inactivation, and to predict expression status of individual genes. Focusing on Xp22, where both inactive and active genes reside within evolutionarily young strata, we compare neighborhoods of genes with different inactivation states to identify enriched oligomers. Occurrences of such oligomers are then used as features to train a linear discriminant analysis classifier. Remarkably, expression status is correctly predicted for 84% and 91% of active and inactive genes, respectively, on the entire X, suggesting that oligomers enriched in Xp22 capture most of the genomic signal determining inactivation. To our surprise, the majority of oligomers associated with inactivated genes fall within L1 elements, even though L1 frequency in Xp22 is low. Moreover, these oligomers are enriched in parts of L1 sequences that are usually underrepresented in the genome. Thus, our results strongly support the role of L1s in X inactivation, yet indicate that a chromatin microenvironment composed of multiple genomic sequence elements determines expression status of X chromosome genes.