Long tandem repeats as a form of genomic copy number variation: structure and length polymorphism of a chromosome 5p repeat in control and schizophrenia populations

Gene-Environment Interactions in Repeat Expansion Diseases: Mechanisms of Environmentally Induced Repeat Instability

Short tandem repeats (STRs) are units of 1-6 base pairs that occur in tandem repetition to form a repeat tract. STRs exhibit repeat instability, which generates expansions or contractions of the repeat tract. Over 50 diseases, primarily affecting the central nervous system and muscles, are characterized by repeat instability. Longer repeat tracts are typically associated with earlier age of onset and increased disease severity. Environmental exposures are suspected to play a role in the pathogenesis of repeat expansion diseases. Here, we review the current knowledge of mechanisms of environmentally induced repeat instability in repeat expansion diseases. The current evidence demonstrates that environmental factors modulate repeat instability via DNA damage and induction of DNA repair pathways, with distinct mechanisms for repeat expansion and contraction. Of particular note, oxidative stress is a key mediator of environmentally induced repeat instability. The preliminary evidence suggests epigenetic modifications as potential mediators of environmentally induced repeat instability. Future research incorporating an array of environmental exposures, new human cohorts, and improved model systems, with a continued focus on cell-types, tissues, and critical windows, will aid in identifying mechanisms of environmentally induced repeat instability. Identifying environmental modulators of repeat instability and their mechanisms of action will inform preventions, therapies, and public health measures.

Epigenetic patterns in a complete human genome

The completion of a telomere-to-telomere human reference genome, T2T-CHM13, has resolved complex regions of the genome, including repetitive and homologous regions. Here, we present a high-resolution epigenetic study of previously unresolved sequences, representing entire acrocentric chromosome short arms, gene family expansions, and a diverse collection of repeat classes. This resource precisely maps CpG methylation (32.28 million CpGs), DNA accessibility, and short-read datasets (166,058 previously unresolved chromatin immunoprecipitation sequencing peaks) to provide evidence of activity across previously unidentified or corrected genes and reveals clinically relevant paralog-specific regulation. Probing CpG methylation across human centromeres from six diverse individuals generated an estimate of variability in kinetochore localization. This analysis provides a framework with which to investigate the most elusive regions of the human genome, granting insights into epigenetic regulation.

Emerging roles of macrosatellite repeats in genome organization and disease development

Abundant repetitive DNA sequences are an enigmatic part of the human genome. Despite increasing evidence on the functionality of DNA repeats, their biologic role is still elusive and under frequent debate. Macrosatellites are the largest of the tandem DNA repeats, located on one or multiple chromosomes. The contribution of macrosatellites to genome regulation and human health was demonstrated for the D4Z4 macrosatellite repeat array on chromosome 4q35. Reduced copy number of D4Z4 repeats is associated with local euchromatinization and the onset of facioscapulohumeral muscular dystrophy. Although the role other macrosatellite families may play remains rather obscure, their diverse functionalities within the genome are being gradually revealed. In this review, we will outline structural and functional features of coding and noncoding macrosatellite repeats, and highlight recent findings that bring these sequences into the spotlight of genome organization and disease development.

Polycomb repressive complex 1 provides a molecular explanation for repeat copy number dependency in FSHD muscular dystrophy

Repression of repetitive elements is crucial to preserve genome integrity and has been traditionally ascribed to constitutive heterochromatin pathways. FacioScapuloHumeral Muscular Dystrophy (FSHD), one of the most common myopathies, is characterized by a complex interplay of genetic and epigenetic events. The main FSHD form is linked to a reduced copy number of the D4Z4 macrosatellite repeat on 4q35, causing loss of silencing and aberrant expression of the D4Z4-embedded DUX4 gene leading to disease. By an unknown mechanism, D4Z4 copy-number correlates with FSHD phenotype. Here we show that the DUX4 proximal promoter (DUX4p) is sufficient to nucleate the enrichment of both constitutive and facultative heterochromatin components and to mediate a copy-number dependent gene silencing. We found that both the CpG/GC dense DNA content and the repetitive nature of DUX4p arrays are important for their repressive ability. We showed that DUX4p mediates a copy number-dependent Polycomb Repressive Complex 1 (PRC1) recruitment, which is responsible for the copy-number dependent gene repression. Overall, we directly link genetic and epigenetic defects in FSHD by proposing a novel molecular explanation for the copy number-dependency in FSHD pathogenesis, and offer insight into the molecular functions of repeats in chromatin regulation.

Influence of Repressive Histone and DNA Methylation upon D4Z4 Transcription in Non-Myogenic Cells

We looked at a disease-associated macrosatellite array D4Z4 and focused on epigenetic factors influencing its chromatin state outside of the disease-context. We used the HCT116 cell line that contains the non-canonical polyadenylation (poly-A) signal required to stabilize somatic transcripts of the human double homeobox gene DUX4, encoded from D4Z4. In HCT116, D4Z4 is packaged into constitutive heterochromatin, characterized by DNA methylation and histone H3 tri-methylation at lysine 9 (H3K9me3), resulting in low basal levels of D4Z4-derived transcripts. However, a double knockout (DKO) of DNA methyltransferase genes, DNMT1 and DNMT3B, but not either alone, results in significant loss of DNA and H3K9 methylation. This is coupled with upregulation of transcript levels from the array, including DUX4 isoforms (DUX4-fl) that are abnormally expressed in somatic muscle in the disease Facioscapulohumeral muscular dystrophy (FSHD) along with DUX4 protein, as indicated indirectly by upregulation of bondafide targets of DUX4 in DKO but not HCT116 cells. Results from treatment with a chemical inhibitor of histone methylation in HCT116 suggest that in the absence of DNA hypomethylation, H3K9me3 loss alone is sufficient to facilitate DUX4-fl transcription. Additionally, characterization of a cell line from a patient with Immunodeficiency, Centromeric instability and Facial anomalies syndrome 1 (ICF1) possessing a non-canonical poly-A signal and DNA hypomethylation at D4Z4 showed DUX4 target gene upregulation in the patient when compared to controls in spite of retention of H3K9me3. Taken together, these data suggest that both DNA methylation and H3K9me3 are determinants of D4Z4 silencing. Moreover, we show that in addition to testis, there is appreciable expression of spliced and polyadenylated D4Z4 derived transcripts that contain the complete DUX4 open reading frame (ORF) along with DUX4 target gene expression in the thymus, suggesting that DUX4 may provide normal function in this somatic tissue.

Estimation of the RNU2 macrosatellite mutation rate by BRCA1 mutation tracing

Large tandem repeat sequences have been poorly investigated as severe technical limitations and their frequent absence from the genome reference hinder their analysis. Extensive allelotyping of this class of variation has not been possible until now and their mutational dynamics are still poorly known. In order to estimate the mutation rate of a macrosatellite, we analysed in detail the RNU2 locus, which displays at least 50 different alleles containing 5-82 copies of a 6.1 kb repeat unit. Mining data from the 1000 Genomes Project allowed us to precisely estimate copy numbers of the RNU2 repeat unit using read depth of coverage. This further revealed significantly different mean values in various recent modern human populations, favoring a scenario of fast evolution of this locus. Its proximity to a disease gene with numerous founder mutations, BRCA1, within the same linkage disequilibrium block, offered the unique opportunity to trace RNU2 arrays over a large timescale. Analysis of the transmission of RNU2 arrays associated with one ‘private’ mutation in an extended kindred and four founder mutations in multiple kindreds gave an estimation by maximum likelihood of 5 × 10⁻³ mutations per generation, which is close to that of microsatellites.

A novel tRNA variable number tandem repeat at human chromosome 1q23.3 is implicated as a boundary element based on conservation of a CTCF motif in mouse

The human genome contains numerous large tandem repeats, many of which remain poorly characterized. Here we report a novel transfer RNA (tRNA) tandem repeat on human chromosome 1q23.3 that shows extensive copy number variation with 9-43 repeat units per allele and displays evidence of meiotic and mitotic instability. Each repeat unit consists of a 7.3 kb GC-rich sequence that binds the insulator protein CTCF and bears the chromatin hallmarks of a bivalent domain in human embryonic stem cells. A tRNA containing tandem repeat composed of at least three 7.6-kb GC-rich repeat units reside within a syntenic region of mouse chromosome 1. However, DNA sequence analysis reveals that, with the exception of the tRNA genes that account for less than 6% of a repeat unit, the remaining 7.2 kb is not conserved with the notable exception of a 24 base pair sequence corresponding to the CTCF binding site, suggesting an important role for this protein at the locus.

Sequencing studies in human genetics: design and interpretation

Next-generation sequencing is becoming the primary discovery tool in human genetics. There have been many clear successes in identifying genes that are responsible for Mendelian diseases, and sequencing approaches are now poised to identify the mutations that cause undiagnosed childhood genetic diseases and those that predispose individuals to more common complex diseases. There are, however, growing concerns that the complexity and magnitude of complete sequence data could lead to an explosion of weakly justified claims of association between genetic variants and disease. Here, we provide an overview of the basic workflow in next-generation sequencing studies and emphasize, where possible, measures and considerations that facilitate accurate inferences from human sequencing studies.

Genome-wide analysis of macrosatellite repeat copy number variation in worldwide populations: evidence for differences and commonalities in size distributions and size restrictions

Background

Macrosatellite repeats (MSRs), usually spanning hundreds of kilobases of genomic DNA, comprise a significant proportion of the human genome. Because of their highly polymorphic nature, MSRs represent an extreme example of copy number variation, but their structure and function is largely understudied. Here, we describe a detailed study of six autosomal and two X chromosomal MSRs among 270 HapMap individuals from Central Europe, Asia and Africa. Copy number variation, stability and genetic heterogeneity of the autosomal macrosatellite repeats RS447 (chromosome 4p), MSR5p (5p), FLJ40296 (13q), RNU2 (17q) and D4Z4 (4q and 10q) and X chromosomal DXZ4 and CT47 were investigated.

Results

Repeat array size distribution analysis shows that all of these MSRs are highly polymorphic with the most genetic variation among Africans and the least among Asians. A mitotic mutation rate of 0.4-2.2% was observed, exceeding meiotic mutation rates and possibly explaining the large size variability found for these MSRs. By means of a novel Bayesian approach, statistical support for a distinct multimodal rather than a uniform allele size distribution was detected in seven out of eight MSRs, with evidence for equidistant intervals between the modes.

Conclusions

The multimodal distributions with evidence for equidistant intervals, in combination with the observation of MSR-specific constraints on minimum array size, suggest that MSRs are limited in their configurations and that deviations thereof may cause disease, as is the case for facioscapulohumeral muscular dystrophy. However, at present we cannot exclude that there are mechanistic constraints for MSRs that are not directly disease-related. This study represents the first comprehensive study of MSRs in different human populations by applying novel statistical methods and identifies commonalities and differences in their organization and function in the human genome.

The macrosatellite DXZ4 mediates CTCF-dependent long-range intrachromosomal interactions on the human inactive X chromosome

The human X-linked macrosatellite DXZ4 is a large tandem repeat located at Xq23 that is packaged into heterochromatin on the male X chromosome and female active X chromosome and, in response to X chromosome, inactivation is organized into euchromatin bound by the insulator protein CCCTC-binding factor (CTCF) on the inactive X chromosome (Xi). The purpose served by this unusual epigenetic regulation is unclear, but suggests a Xi-specific gain of function for DXZ4. Other less extensive bands of euchromatin can be observed on the Xi, but the identity of the underlying DNA sequences is unknown. Here, we report the identification of two novel human X-linked tandem repeats, located 58 Mb proximal and 16 Mb distal to the macrosatellite DXZ4. Both tandem repeats are entirely contained within the transcriptional unit of novel spliced transcripts. Like DXZ4, the tandem repeats are packaged into Xi-specific CTCF-bound euchromatin. These sequences undergo frequent CTCF-dependent interactions with DXZ4 on the Xi, implicating DXZ4 as an epigenetically regulated Xi-specific structural element and providing the first putative functional attribute of a macrosatellite in the human genome.

The mouse DXZ4 homolog retains Ctcf binding and proximity to Pls3 despite substantial organizational differences compared to the primate macrosatellite

Background

The X-linked macrosatellite DXZ4 is a large homogenous tandem repeat that in females adopts an alternative chromatin organization on the primate X chromosome in response to X-chromosome inactivation. It is packaged into heterochromatin on the active X chromosome but into euchromatin and bound by the epigenetic organizer protein CTCF on the inactive X chromosome. Because its DNA sequence diverges rapidly beyond the New World monkeys, the existence of DXZ4 outside the primate lineage is unknown.

Results

Here we extend our comparative genome analysis and report the identification and characterization of the mouse homolog of the macrosatellite. Furthermore, we provide evidence of DXZ4 in a conserved location downstream of the PLS3 gene in a diverse group of mammals, and reveal that DNA sequence conservation is restricted to the CTCF binding motif, supporting a central role for this protein at this locus. However, many features that characterize primate DXZ4 differ in mouse, including the overall size of the array, the mode of transcription, the chromatin organization and conservation between adjacent repeat units of DNA sequence and length. Ctcf binds Dxz4 but is not exclusive to the inactive X chromosome, as evidenced by association in some males and equal binding to both X chromosomes in trophoblast stem cells.

Conclusions

Characterization of Dxz4 reveals substantial differences in the organization of DNA sequence, chromatin packaging, and the mode of transcription, so the potential roles performed by this sequence in mouse have probably diverged from those on the primate X chromosome.

A long ncRNA links copy number variation to a polycomb/trithorax epigenetic switch in FSHD muscular dystrophy

Repetitive sequences account for more than 50% of the human genome. Facioscapulohumeral muscular dystrophy (FSHD) is an autosomal-dominant disease associated with reduction in the copy number of the D4Z4 repeat mapping to 4q35. By an unknown mechanism, D4Z4 deletion causes an epigenetic switch leading to de-repression of 4q35 genes. Here we show that the Polycomb group of epigenetic repressors targets D4Z4 in healthy subjects and that D4Z4 deletion is associated with reduced Polycomb silencing in FSHD patients. We identify DBE-T, a chromatin-associated noncoding RNA produced selectively in FSHD patients that coordinates de-repression of 4q35 genes. DBE-T recruits the Trithorax group protein Ash1L to the FSHD locus, driving histone H3 lysine 36 dimethylation, chromatin remodeling, and 4q35 gene transcription. This study provides insights into the biological function of repetitive sequences in regulating gene expression and shows how mutations of such elements can influence the progression of a human genetic disease.

Genome wide association studies (GWAS) and copy number variation (CNV) studies of the major psychoses: what have we learnt?

Schizophrenia (SZ) and bipolar disorder (BPD) have high heritabilities and are clinically and genetically complex. Genome wide association studies (GWAS) and studies of copy number variations (CNV) in SZ and BPD have allowed probing of their underlying genetic risks. In this systematic review, we assess extant genetic signals from published GWAS and CNV studies of SZ and BPD up till March 2011. Risk genes associated with SZ at genome wide significance level (p value<7.2 × 10(-8)) include zinc finger binding protein 804A (ZNF804A), major histocompatibility (MHC) region on chromosome 6, neurogranin (NRGN) and transcription factor 4 (TCF4). Risk genes associated with BPD include ankyrin 3, node of Ranvier (ANK3), calcium channel, voltage dependent, L type, alpha 1C subunit (CACNA1C), diacylglycerol kinase eta (DGKH), gene locus on chromosome 16p12, and polybromo-1 (PBRM1) and very recently neurocan gene (NCAN). Possible common genes underlying psychosis include ZNF804A, CACNA1C, NRGN and PBRM1. The CNV studies suggest that whilst CNVs are found in both SZ and BPD, the large deletions and duplications are more likely found in SZ rather than BPD. The validation of any genetic signal is likely confounded by genetic and phenotypic heterogeneities which are influenced by epistatic, epigenetic and gene-environment interactions. There is a pressing need to better integrate the multiple research platforms including systems biology computational models, genomics, cross disorder phenotyping studies, transcriptomics, proteomics, metabolomics, neuroimaging and clinical correlations in order to get us closer to a more enlightened understanding of the genetic and biological basis underlying these potentially crippling conditions.

Epigenetic regulation of the X-chromosomal macrosatellite repeat encoding for the cancer/testis gene CT47

Macrosatellite repeats (MSRs) present an extreme example of copy number variation, yet their epigenetic regulation in normal and malignant cells is largely understudied. The CT47 cancer/testis antigen located on human Xq24 is organized as an array of 4.8 kb large units. CT47 is expressed in the testis and in certain types of cancer, but not in non-malignant somatic tissue. We used CT47 as a model to study a possible correlation between copy number variation, epigenetic regulation and transcription originating from MSRs in normal and malignant cells. In lymphoblastoid cell lines and primary fibroblasts, CT47 expression was absent, consistent with the observed heterochromatic structure and DNA hypermethylation of the CT47 promoter. Heterochromatinization of CT47 occurs early during development as human embryonic stem cells show high levels of DNA methylation and repressive chromatin modifications in the absence of CT47 expression. In small-cell lung carcinoma cell lines with low levels of CT47 transcripts, we observed reduced levels of histone 3 lysine 9 trimethylation (H3K9me3) and trimethylated lysine 27 of histone H3 (H3K27me3) without concomitant increase in euchromatic histone modifications. DNA methylation levels in the promoter region of CT47 are also significantly reduced in these cells. This supports a model in which during oncogenic transformation, there is a relative loss of repressive chromatin markers resulting in leaky expression of CT47. We conclude that some MSRs, like CT47 and the autosomal MSRs TAF11-Like, PRR20, ZAV and D4Z4, the latter being involved in facioscapulohumeral muscular dystrophy, seem to be governed by common regulatory mechanisms with their abundant expression mostly being restricted to the germ line.

Variation in array size, monomer composition and expression of the macrosatellite DXZ4

Macrosatellites are some of the most polymorphic regions of the human genome, yet many remain uncharacterized despite the association of some arrays with disease susceptibility. This study sought to explore the polymorphic nature of the X-linked macrosatellite DXZ4. Four aspects of DXZ4 were explored in detail, including tandem repeat copy number variation, array instability, monomer sequence polymorphism and array expression. DXZ4 arrays contained between 12 and 100 3.0 kb repeat units with an average array containing 57. Monomers were confirmed to be arranged in uninterrupted tandem arrays by restriction digest analysis and extended fiber FISH, and therefore DXZ4 encompasses 36–288 kb of Xq23. Transmission of DXZ4 through three generations in three families displayed a high degree of meiotic instability (8.3%), consistent with other macrosatellite arrays, further highlighting the unstable nature of these sequences in the human genome. Subcloning and sequencing of complete DXZ4 monomers identified numerous single nucleotide polymorphisms and alleles for the three microsatellite repeats located within each monomer. Pairwise comparisons of DXZ4 monomer sequences revealed that repeat units from an array are more similar to one another than those originating from different arrays. RNA fluorescence in situ hybridization revealed significant variation in DXZ4 expression both within and between cell lines. DXZ4 transcripts could be detected originiating from both the active and inactive X chromosome. Expression levels of DXZ4 varied significantly between males, but did not relate to the size of the array, nor did inheritance of the same array result in similar expression levels. Collectively, these studies provide considerable insight into the polymorphic nature of DXZ4, further highlighting the instability and variation potential of macrosatellites in the human genome.

Characterization of DXZ4 conservation in primates implies important functional roles for CTCF binding, array expression and tandem repeat organization on the X chromosome

BACKGROUND

Comparative sequence analysis is a powerful means with which to identify functionally relevant non-coding DNA elements through conserved nucleotide sequence. The macrosatellite DXZ4 is a polymorphic, uninterrupted, tandem array of 3-kb repeat units located exclusively on the human X chromosome. While not obviously protein coding, its chromatin organization suggests differing roles for the array on the active and inactive X chromosomes.

RESULTS

In order to identify important elements within DXZ4, we explored preservation of DNA sequence and chromatin conformation of the macrosatellite in primates. We found that DXZ4 DNA sequence conservation beyond New World monkeys is limited to the promoter and CTCF binding site, although DXZ4 remains a GC-rich tandem array. Investigation of chromatin organization in macaques revealed that DXZ4 in males and on the active X chromosome is packaged into heterochromatin, whereas on the inactive X, DXZ4 was euchromatic and bound by CTCF.

CONCLUSIONS

Collectively, these data suggest an important conserved role for DXZ4 on the X chromosome involving expression, CTCF binding and tandem organization.

Expression, tandem repeat copy number variation and stability of four macrosatellite arrays in the human genome

Background

Macrosatellites are some of the largest variable number tandem repeats in the human genome, but what role these unusual sequences perform is unknown. Their importance to human health is clearly demonstrated by the 4q35 macrosatellite D4Z4 that is associated with the onset of the muscle degenerative disease facioscapulohumeral muscular dystrophy. Nevertheless, many other macrosatellite arrays in the human genome remain poorly characterized.

Results

Here we describe the organization, tandem repeat copy number variation, transmission stability and expression of four macrosatellite arrays in the human genome: the TAF11-Like array located on chromosomes 5p15.1, the SST1 arrays on 4q28.3 and 19q13.12, the PRR20 array located on chromosome 13q21.1, and the ZAV array at 9q32. All are polymorphic macrosatellite arrays that at least for TAF11-Like and SST1 show evidence of meiotic instability. With the exception of the SST1 array that is ubiquitously expressed, all are expressed at high levels in the testis and to a lesser extent in the brain.

Conclusions

Our results extend the number of characterized macrosatellite arrays in the human genome and provide the foundation for formulation of hypotheses to begin assessing their functional role in the human genome.