Tandem repeat polymorphisms: modulators of disease susceptibility and candidates for ‘missing heritability’

Identification of an (AC)n microsatellite in the Six1 gene promoter and its effect on production traits in Pietrain × Duroc × Landrace × Yorkshire pigs

The Sine oculis homeobox 1 (Six1) gene is important for skeletal muscle growth and fiber specification; therefore, it is considered as a promising candidate gene that may influence porcine growth and meat quality traits. Nevertheless, the association of Six1 with these processes and the mechanisms regulating its expression remain unclear. The objectives of this study were to identify variant sites of Six1 in different pig breeds, conduct association analysis to evaluate the relationship between polymorphisms of these variants and porcine production traits in Pietrain × Duroc × Landrace × Yorkshire commercial pigs, and explore the potential regulatory mechanisms of Six1 affecting production traits. A total of 12 variants were identified, including 10 single- nucleotide variations (SNVs), 1 insertion- deletion (Indel), and 1 (AC)n microsatellite. Association analysis demonstrated that the SNV, g.1595A>G, was significantly associated with meat color (redness, a*); individuals with the G allele had greater a* values (P < 0.05). Moreover, our results demonstrated that the (AC)n polymorphism in the Six1 promoter was significantly associated with weaning weight (P < 0.05), carcass weight (P < 0.05), and thoracic and lumbar back fat (P < 0.01).In addition, we found that the (AC)n variant was closely related with Six1 expression levels and demonstrated this polymorphism on promoter activity by in vitro experiments. Overall, this study provides novel evidence for elucidating the effects of Six1 on porcine production traits as promising candidate and describes two variants with these traits, which are potential reference markers for pig molecular breeding. In addition, our data on the relationship between porcine Six1 expression and the polymorphic (AC)n microsatellite in its promoter may facilitate similar studies in other species.

Massive variation of short tandem repeats with functional consequences across strains of Arabidopsis thaliana

Short tandem repeat (STR) mutations may comprise more than half of the mutations in eukaryotic coding DNA, yet STR variation is rarely examined as a contributor to complex traits. We assessed this contribution across a collection of 96 strains of Arabidopsis thaliana, genotyping 2046 STR loci each, using highly parallel STR sequencing with molecular inversion probes. We found that 95% of examined STRs are polymorphic, with a median of six alleles per STR across these strains. STR expansions (large copy number increases) are found in most strains, several of which have evident functional effects. These include three of six intronic STR expansions we found to be associated with intron retention. Coding STRs were depleted of variation relative to noncoding STRs, and we detected a total of 56 coding STRs (11%) showing low variation consistent with the action of purifying selection. In contrast, some STRs show hypervariable patterns consistent with diversifying selection. Finally, we detected 133 novel STR-phenotype associations under stringent criteria, most of which could not be detected with SNPs alone, and validated some with follow-up experiments. Our results support the conclusion that STRs constitute a large, unascertained reservoir of functionally relevant genomic variation.

Investigation of short tandem repeats in major depression using whole-genome sequencing data

BACKGROUND

Major depressive disorder (MDD) is a leading contributor to global disease burden. Recent studies have shown that genetic factors play significant roles in the susceptibility to this condition; however, the underlying genetic basis currently remains largely unknown. Short tandem repeat (STR) has been proposed as an explanatory factor in the "missing heritability" of complex diseases or traits.

METHODS

We investigated STR variations from 15 MDD patients and 10 ethnically matched healthy controls based on their deep whole-genome sequencing (WGS) data. The lobSTR software was used to computationally determine STRs.

RESULTS

The results of the Mexican-American sample showed that STRs are significantly richer in healthy controls than in MDD cases on each of the 23 chromosomes (all false discovery rates, FDR P-values < 0.0062); while for the Australian of European-ancestry sample, there was no statistically significant STRs difference between MDD cases and controls.

LIMITATIONS

High quality WGS costs limited obtaining larger datasets.

CONCLUSIONS

This preliminary work is the first study that STR variations are applied to investigate MDD based on WGS data. The results on Mexican-American population may imply that within the same ancestry, targeted sequencing on a specific chromosome or region of genome would be sufficient for examining the relationship between STR and MDD. Further studies should examine larger sequencing datasets on other ethnic groups.

Effect of Trinucleotide Repeats in the Huntington's Gene on Intelligence

Background

Huntington's Disease (HD) is caused by an abnormality in the HTT gene. This gene includes trinucleotide repeats ranging from 10 to 35, and when expanded beyond 39, causes HD. We previously reported that CAG repeats in the normal range had a direct and beneficial effect on brain development with higher repeats being associated with higher cognitive function. The current study now expands this line of inquiry to evaluate the effects of CAG repeat throughout the entire spectrum of repeats from 15 to 58.

Methods

We evaluated brain function in children ages 6–18 years old. DNA samples were processed to quantify the number of CAG repeats within HTT. Linear regression was used to determine if number of CAG repeats predicted measures of brain function.

Findings

The number of repeats in HTT, had a non-linear effect on a measure of general intelligence with an inverted U shape pattern. Increasing repeat length was associated with higher GAI scores up until roughly 40–41 repeats. After this peak, increasing repeat length was associated with declining GAI scores.

Interpretation

HTT may confer an advantage or a disadvantage depending upon the repeat length, playing a key role in the determination of intelligence, or causing a uniquely human brain disease.

•

The HTT gene includes trinucleotide repeats ranging from 10 to 35 in the normal population.

•

Huntington’s disease is caused when the trinucleotide repeats in the HTT gene expand beyond the normal range.

•

The gene is vital for brain growth in which each repeat contributes to the development of brain function, measured as intelligence.

•

Increasing repeats are beneficial to a certain point, then becomes detrimental to intelligence, forming an inverted U relationship.

A Faustian Bargain: Could the key to the evolution of the human brain be found in a dreadful disease? Huntington's Disease is a fatal progressive brain disease caused by a single gene. Results from the current study show that the same gene is important for brain development and intelligence. Therefore, this gene may have a simultaneous advantage and disadvantage: a role in the evolution of a superior human brain, yet the disadvantage of a uniquely human brain disease.

The effects of transcription and recombination on mutational dynamics of short tandem repeats

Short tandem repeats (STR) are ubiquitous components of the genomic architecture of most living organisms. Recent work has highlighted the widespread functional significance of such repeats, particularly around gene regulation, but the mutational processes underlying the evolution of these highly abundant and highly variable sequences are not fully understood. Traditional models assume that strand misalignment during replication is the predominant mechanism, but empirical data suggest the involvement of other processes including recombination and transcription. Despite this evidence, the relative influences of these processes have not previously been tested experimentally on a genome-wide scale. Using deep sequencing, we identify mutations at >200 microsatellites, across 700 generations in replicated populations of two otherwise identical sexual and asexual Saccharomyces cerevisiae strains. Using generalized linear models, we investigate correlates of STR mutability including the nature of the mutation, STR composition and contextual factors including recombination, transcription and replication origins. Sexual capability was not a significant predictor of microsatellite mutability, but, intriguingly, we identify transcription as a significant positive predictor. We also find that STR density is substantially increased in regions neighboring, but not within, recombination hotspots.

Tandem repeats mediating genetic plasticity in health and disease

Accumulating evidence suggests that many classes of DNA repeats exhibit attributes that distinguish them from other genetic variants, including the fact that they are more liable to mutation; this enables them to mediate genetic plasticity. The expansion of tandem repeats, particularly of short tandem repeats, can cause a range of disorders (including Huntington disease, various ataxias, motor neuron disease, frontotemporal dementia, fragile X syndrome and other neurological disorders), and emerging data suggest that tandem repeat polymorphisms (TRPs) can also regulate gene expression in healthy individuals. TRPs in human genomes may also contribute to the missing heritability of polygenic disorders. A better understanding of tandem repeats and their associated repeatome, as well as their capacity for genetic plasticity via both germline and somatic mutations, is needed to transform our understanding of the role of TRPs in health and disease.

Huntingtin gene repeat size variations affect risk of lifetime depression

Huntington disease (HD) is a severe neuropsychiatric disorder caused by a cytosine-adenine-guanine (CAG) repeat expansion in the HTT gene. Although HD is frequently complicated by depression, it is still unknown to what extent common HTT CAG repeat size variations in the normal range could affect depression risk in the general population. Using binary logistic regression, we assessed the association between HTT CAG repeat size and depression risk in two well-characterized Dutch cohorts─the Netherlands Study of Depression and Anxiety and the Netherlands Study of Depression in Older Persons─including 2165 depressed and 1058 non-depressed persons. In both cohorts, separately as well as combined, there was a significant non-linear association between the risk of lifetime depression and HTT CAG repeat size in which both relatively short and relatively large alleles were associated with an increased risk of depression (β = -0.292 and β = 0.006 for the linear and the quadratic term, respectively; both P < 0.01 after adjustment for the effects of sex, age, and education level). The odds of lifetime depression were lowest in persons with a HTT CAG repeat size of 21 (odds ratio: 0.71, 95% confidence interval: 0.52 to 0.98) compared to the average odds in the total cohort. In conclusion, lifetime depression risk was higher with both relatively short and relatively large HTT CAG repeat sizes in the normal range. Our study provides important proof-of-principle that repeat polymorphisms can act as hitherto unappreciated but complex genetic modifiers of depression.

Regulation of antilipopolysaccharide factors, ALFPm3 and ALFPm6, in Penaeus monodon

ALFPm6, a member of antimicrobial peptide in the antilipopolysaccharide factor (ALF) family from Penaeus monodon, plays important roles in shrimp immunity against pathogens. However, its antimicrobial activity and underlying mechanism have not been reported. The synthetic cyclic ALFPm6#29–52 peptide (cALFPm6#29–52) corresponding to the ALFPm6 LPS-binding domain can agglutinate and exhibited bacterial killing activity toward a Gram-negative bacterium, Escherichia coli 363 and Gram-positive bacteria, Bacillus megaterium, Aerococcus viridans, and Micrococcus luteus, with MIC values of 25–50 μM. Specifically, ALFPm6 and ALFPm3, the most abundant ALF isoforms, are different in terms of gene expression patterns upon pathogen infections. Herein, the regulation of ALFPm3 and ALFPm6 gene expression was studied. The 5′-upstream and promoter sequences were identified and the putative transcription factor (TF)-binding sites were predicted. The narrow down assay indicated that the ALFPm3 promoter and partial promoter of the ALFPm6 active regions were located at nucleotide positions (−814/+302) and (−282/+85), respectively. Mutagenesis of selected TF-binding sites revealed that Rel/NF-κB (−280/−270) of ALFPm3 and C/EBPβ (−88/−78) and Sp1 (−249/−238) sites of ALFPm6 were the activator-binding sites. Knockdown of the PmMyD88 and PmRelish genes in V. harveyi-infected shrimp suggested that the ALFPm3 gene was regulated by Toll and IMD pathways, while the ALFPm6 gene was regulated by the Toll pathway.

Large normal-range TBP and ATXN7 CAG repeat lengths are associated with increased lifetime risk of depression

Depression is one of the most prevalent and debilitating psychiatric disorders worldwide. Recently, we showed that both relatively short and relatively long cytosine-adenine-guanine (CAG) repeats in the huntingtin gene (HTT) are associated with an increased risk of lifetime depression. However, to what extent the variations in CAG repeat length in the other eight polyglutamine disease-associated genes (PDAGs) are associated with depression is still unknown. We determined the CAG repeat sizes of ATXN1, ATXN2, ATXN3, CACNA1A, ATXN7, TBP, ATN1 and AR in two well-characterized Dutch cohorts-the Netherlands Study of Depression and Anxiety and the Netherlands Study of Depression in Older Persons-including 2165 depressed and 1058 non-depressed individuals-aged 18-93 years. The association between PDAG CAG repeat size and the risk for depression was assessed via binary logistic regression. We found that the odds ratio (OR) for lifetime depression was significantly higher for individuals with >10, compared with subjects with ≤10, CAG repeats in both ATXN7 alleles (OR=1.90, confidence interval (CI) 1.26-2.85). For TBP we found a similar association: A CAG repeat length exceeding the median in both alleles was associated with an increased risk for lifetime depression (OR=1.33, CI 1.00-1.76). In conclusion, we observed that carriers of either ATXN7 or TBP alleles with relatively large CAG repeat sizes in both alleles had a substantially increased risk of lifetime depression. Our findings provide critical evidence for the notion that repeat polymorphisms can act as complex genetic modifiers of depression.

Genetic variations at the human growth hormone receptor (GHR) gene locus are associated with idiopathic short stature

GH plays an essential role in the growing child by binding to the growth hormone receptor (GHR) on target cells and regulating multiple growth promoting and metabolic effects. Mutations in the GHR gene coding regions result in GH insensitivity (dwarfism) due to a dysfunctional receptor protein. However, children with idiopathic short stature (ISS) show growth impairment without GH or GHR defects. We hypothesized that decreased expression of the GHR gene may be involved. To test this, we investigated whether common genetic variants (microsatellites, SNPs) in regulatory regions of the GHR gene region were associated with the ISS phenotype. Genotyping of a GT‐repeat microsatellite in the GHR 5′UTR in a Montreal ISS cohort (n = 37 ISS, n = 105 controls) revealed that the incidence of the long/short (L/S) genotype was 3.3× higher in ISS children than controls (P = 0.04, OR = 3.85). In an Italian replication cohort (n = 143 ISS, n = 282 controls), the medium/short (M/S) genotype was 1.9× more frequent in the male ISS than controls (P = 0.017, OR = 2.26). In both ISS cohorts, logistic regression analysis of 27 SNPs showed an association of ISS with rs4292454, while haplotype analysis revealed specific risk haplotypes in the 3′ haploblocks. In contrast, there were no differences in GT genotype frequencies in a cohort of short stature (SS) adults versus controls (CARTaGENE: n = 168 SS, n = 207 controls) and the risk haplotype in the SS cohort was located in the most 5′ haploblock. These data suggest that the variants identified are potentially genetic markers specifically associated with the ISS phenotype.

CRISPR-Cas9-targeted fragmentation and selective sequencing enable massively parallel microsatellite analysis

Microsatellites are multi-allelic and composed of short tandem repeats (STRs) with individual motifs composed of mononucleotides, dinucleotides or higher including hexamers. Next-generation sequencing approaches and other STR assays rely on a limited number of PCR amplicons, typically in the tens. Here, we demonstrate STR-Seq, a next-generation sequencing technology that analyses over 2,000 STRs in parallel, and provides the accurate genotyping of microsatellites. STR-Seq employs in vitro CRISPR-Cas9-targeted fragmentation to produce specific DNA molecules covering the complete microsatellite sequence. Amplification-free library preparation provides single molecule sequences without unique molecular barcodes. STR-selective primers enable massively parallel, targeted sequencing of large STR sets. Overall, STR-Seq has higher throughput, improved accuracy and provides a greater number of informative haplotypes compared with other microsatellite analysis approaches. With these new features, STR-Seq can identify a 0.1% minor genome fraction in a DNA mixture composed of different, unrelated samples.

The mechanism of transactivation regulation due to polymorphic short tandem repeats (STRs) using IGF1 promoter as a model

Functional short tandem repeats (STR) are polymorphic in the population, and the number of repeats regulates the expression of nearby genes (known as expression STR, eSTR). STR in IGF1 promoter has been extensively studied for its association with IGF1 concentration in blood and various clinical traits and represents an important eSTR. We previously used an in-vitro luciferase reporter model to examine the interaction between STRs and SNPs in IGF1 promoter. Here, we further explored the mechanism how the number of repeats of the STR regulates gene transcription. An inverse correlation between the number of repeats and the extent of transactivation was found in a haplotype consisting of three promoter SNPs (C-STR-T-T). We showed that these adjacent SNPs located outside the STR were required for the STR to function as eSTR. The C allele of rs35767 provides a binding site for CCAAT/enhancer-binding-protein δ (C/EBPD), which is essential for the gradational transactivation property of eSTR and FOXA3 may also be involved. Therefore, we propose a mechanism in which the gradational transactivation by the eSTR is caused by the interaction of one or more transcriptional complexes located outside the STR, rather than by direct binding to a repeat motif of the STR.

Single Amino Acid Repeats in the Proteome World: Structural, Functional, and Evolutionary Insights

Microsatellites or simple sequence repeats (SSR) are abundant, highly diverse stretches of short DNA repeats present in all genomes. Tandem mono/tri/hexanucleotide repeats in the coding regions contribute to single amino acids repeats (SAARs) in the proteome. While SSRs in the coding region always result in amino acid repeats, a majority of SAARs arise due to a combination of various codons representing the same amino acid and not as a consequence of SSR events. Certain amino acids are abundant in repeat regions indicating a positive selection pressure behind the accumulation of SAARs. By analysing 22 proteomes including the human proteome, we explored the functional and structural relationship of amino acid repeats in an evolutionary context. Only ~15% of repeats are present in any known functional domain, while ~74% of repeats are present in the disordered regions, suggesting that SAARs add to the functionality of proteins by providing flexibility, stability and act as linker elements between domains. Comparison of SAAR containing proteins across species reveals that while shorter repeats are conserved among orthologs, proteins with longer repeats, >15 amino acids, are unique to the respective organism. Lysine repeats are well conserved among orthologs with respect to their length and number of occurrences in a protein. Other amino acids such as glutamic acid, proline, serine and alanine repeats are generally conserved among the orthologs with varying repeat lengths. These findings suggest that SAARs have accumulated in the proteome under positive selection pressure and that they provide flexibility for optimal folding of functional/structural domains of proteins. The insights gained from our observations can help in effective designing and engineering of proteins with novel features.

Variability in a Short Tandem Repeat Mediates Complex Epistatic Interactions in Arabidopsis thaliana

Short tandem repeats (STRs) are hypervariable genetic elements that occur frequently in coding regions. Their high mutation rate readily generates genetic variation, contributing to adaptive evolution and human diseases. We previously reported that natural ELF3 polyglutamine variants cause reciprocal genetic incompatibilities in two divergent Arabidopsis thaliana backgrounds. Here, we dissect the genetic architecture of this incompatibility, revealing as many as four loci putatively interacting with ELF3 We were able to specifically identify one such ELF3-interacting gene, LSH9 We further used a yeast two-hybrid strategy to identify proteins whose physical interactions with ELF3 were affected by polyglutamine tract length. We found two proteins for which this was the case, ELF4 and AtGLDP1. Using these two approaches, we identify specific genetic interactions and physical mechanisms by which the ELF3 polyglutamine tract may mediate the observed genetic incompatibilities. Our work elucidates how STR variation, which is generally underascertained in population-scale sequencing, can contribute to phenotypic variation. Furthermore, our results support our proposal that highly variable STR loci can contribute to the epistatic component of heritability.

Genotypic Context and Epistasis in Individuals and Populations

Genes encode components of coevolved and interconnected networks. The effect of genotype on phenotype therefore depends on genotypic context through gene interactions known as epistasis. Epistasis is important in predicting phenotype from genotype for an individual. It is also examined in population studies to identify genetic risk factors in complex traits and to predict evolution under selection. Paradoxically, the effects of genotypic context in individuals and populations are distinct and sometimes contradictory. We argue that predicting genotype from phenotype for individuals based on population studies is difficult and, especially in human genetics, likely to result in underestimating the effects of genotypic context.

A novel polymorphic repeat in the upstream regulatory region of the estrogen-induced gene EIG121 is not associated with the risk of developing breast or endometrial cancer

BACKGROUND

The estrogen-induced gene 121 (EIG121) has been associated with breast and endometrial cancers, but its mechanism of action remains unknown. In a genome-wide search for tandem repeats, we found that EIG121 contains a short tandem repeat (STR) in its upstream regulatory region which has the potential to alter gene expression. The presence of this STR has not previously been analysed in relation to breast or endometrial cancer risk.

RESULTS

In this study, the lengths of this STR were determined by PCR, fragment analysis and sequencing using DNA from 223 breast cancer patients, 204 endometrial cancer patients and 220 healthy controls to determine if they were associated with the risk of developing breast or endometrial cancer. We found this repeat to be highly variable with the number of copies of the AG motif ranging from 27 to 72 and having a bimodal distribution. No statistically significant association was identified between the length of this STR and the risk of developing breast or endometrial cancer or age at diagnosis.

CONCLUSIONS

The STR in the upstream regulatory region of EIG121 is highly polymorphic, but is not associated with the risk of developing breast or endometrial cancer in the cohorts analysed here. While this polymorphic STR in the regulatory region of EIG121 appears to have no impact on the risk of developing breast or endometrial cancer, its association with disease recurrence or overall survival remains to be determined.

A polymorphic repeat in the IGF1 promoter influences the risk of endometrial cancer

Due to the lack of high-throughput genetic assays for tandem repeats, there is a paucity of knowledge about the role they may play in disease. A polymorphic CA repeat in the promoter region of the insulin-like growth factor 1 gene (IGF1 has been studied extensively over the past 10 years for association with the risk of developing breast cancer, among other cancers, with variable results. The aim of this study was to determine if this CA repeat is associated with the risk of developing breast cancer and endometrial cancer. Using a case-control design, we analysed the length of this CA repeat in a series of breast cancer and endometrial cancer cases and compared this with a control population. Our results showed an association when both alleles were considered in breast and endometrial cancers (P=0.029 and 0.011, respectively), but this did not pass our corrected threshold for significance due to multiple testing. When the allele lengths were analysed categorically against the most common allele length of 19 CA repeats, an association was observed with the risk of endometrial cancer due to a reduction in the number of long alleles (P=0.013). This was confirmed in an analysis of the long alleles separately for endometrial cancer risk (P=0.0012). Our study found no association between the length of this polymorphic CA repeat and breast cancer risk. The significant association observed between the CA repeat length and the risk of developing endometrial cancer has not been previously reported.

Survival of the Curviest: Noise-Driven Selection for Synergistic Epistasis

A major goal of human genetics is to elucidate the genetic architecture of human disease, with the goal of fueling improvements in diagnosis and the understanding of disease pathogenesis. The degree to which epistasis, or non-additive effects of risk alleles at different loci, accounts for common disease traits is hotly debated, in part because the conditions under which epistasis evolves are not well understood. Using both theory and evolutionary simulation, we show that the occurrence of common diseases (i.e. unfit phenotypes with frequencies on the order of 1%) can, under the right circumstances, be expected to be driven primarily by synergistic epistatic interactions. Conditions that are necessary, collectively, for this outcome include a strongly non-linear phenotypic landscape, strong (but not too strong) selection against the disease phenotype, and "noise" in the genotype-phenotype map that is both environmental (extrinsic, time-correlated) and developmental (intrinsic, uncorrelated) and, in both cases, neither too little nor too great. These results suggest ways in which geneticists might identify, a priori, those disease traits for which an "epistatic explanation" should be sought, and in the process better focus ongoing searches for risk alleles.

Cognitive endophenotypes, gene-environment interactions and experience-dependent plasticity in animal models of schizophrenia

Schizophrenia is a devastating brain disorder caused by a complex and heterogeneous combination of genetic and environmental factors. In order to develop effective new strategies to prevent and treat schizophrenia, valid animal models are required which accurately model the disorder, and ideally provide construct, face and predictive validity. The cognitive deficits in schizophrenia represent some of the most debilitating symptoms and are also currently the most poorly treated. Therefore it is crucial that animal models are able to capture the cognitive dysfunction that characterizes schizophrenia, as well as the negative and psychotic symptoms. The genomes of mice have, prior to the recent gene-editing revolution, proven the most easily manipulable of mammalian laboratory species, and hence most genetic targeting has been performed using mouse models. Importantly, when key environmental factors of relevance to schizophrenia are experimentally manipulated, dramatic changes in the phenotypes of these animal models are often observed. We will review recent studies in rodent models which provide insight into gene-environment interactions in schizophrenia. We will focus specifically on environmental factors which modulate levels of experience-dependent plasticity, including environmental enrichment, cognitive stimulation, physical activity and stress. The insights provided by this research will not only help refine the establishment of optimally valid animal models which facilitate development of novel therapeutics, but will also provide insight into the pathogenesis of schizophrenia, thus identifying molecular and cellular targets for future preclinical and clinical investigations.

Understanding the genetic architectonics of complex CNS traits: Lost by the association, but found in the interaction?

Recent evidence supports the value of endophenotypes and genome-wide association studies in psychiatric genetics, and their importance for dissecting the neural pathways and molecular mechanisms of complex neuropsychiatric disorders. Continuing this important discussion, here we outline three new mechanisms by which novel classes of genes may facilitate CNS pathogenesis without directly worsening its individual 'established' endophenotypes. These putative genetic mechanisms can apply to other human disorders in general, and may also be used for designing novel effective CNS drug treatments.

Association of the FAM46A gene VNTRs and BAG6 rs3117582 SNP with non small cell lung cancer (NSCLC) in Croatian and Norwegian populations

We analyzed for associations between a variable number of tandem repeat (VNTR) polymorphism in the Family with sequence similarity 46, member A (FAM46A) gene and a single nucleotide polymorphism (rs3117582) in the BCL2-Associated Athanogene 6 (BAG6) with non small cell lung cancer in Croatian and Norwegian subjects. A total of 503 (262 Croatian and 241Norwegian) non small cell lung cancer patients and 897 controls (568 Croatian and 329 Norwegian) were analyzed. We found that the frequency of allele b (three VNTR repeats) of FAM46A gene was significantly increased in the patients compared to the healthy controls in the Croatian and the combined Croatian and Norwegian subjects. Genotype frequencies of cd (four and five VNTR repeats) and cc (four VNTR repeats homozygote) of the FAM46A gene were significantly decreased in the patients compared to the healthy controls in the Croatian and Norwegian subjects, respectively. Logistic regression analyses revealed FAM46A genotype cc to be an independent predictive factor for non small cell lung cancer risk in the Norwegian subjects after adjustment for age, gender and smoking status. This is the first study to suggest an association between the FAM46A gene VNTR polymorphisms and non small cell lung cancer. We found also that BAG6 rs3117582 SNP was associated with non small cell lung cancer in the Norwegian subjects and the combined Croatian-Norwegian subjects corroborating the earlier finding that BAG6 rs3117582 SNP was associated with lung cancer in Europeans. Logistic regression analyses revealed that genotypes and alleles of BAG6 were independent predictive factor for non small cell lung cancer risk in the Norwegian and combined Croatian-Norwegian subjects, after adjustment for age and gender.

Statistical approaches to detecting and analyzing tandem repeats in genomic sequences

Tandem repeats (TRs) are frequently observed in genomes across all domains of life. Evidence suggests that some TRs are crucial for proteins with fundamental biological functions and can be associated with virulence, resistance, and infectious/neurodegenerative diseases. Genome-scale systematic studies of TRs have the potential to unveil core mechanisms governing TR evolution and TR roles in shaping genomes. However, TR-related studies are often non-trivial due to heterogeneous and sometimes fast evolving TR regions. In this review, we discuss these intricacies and their consequences. We present our recent contributions to computational and statistical approaches for TR significance testing, sequence profile-based TR annotation, TR-aware sequence alignment, phylogenetic analyses of TR unit number and order, and TR benchmarks. Importantly, all these methods explicitly rely on the evolutionary definition of a tandem repeat as a sequence of adjacent repeat units stemming from a common ancestor. The discussed work has a focus on protein TRs, yet is generally applicable to nucleic acid TRs, sharing similar features.

Susceptibility to large-joint osteoarthritis (hip and knee) is associated with BAG6 rs3117582 SNP and the VNTR polymorphism in the second exon of the FAM46A gene on chromosome 6

Family with sequence similarity 46, member A (FAM46A) gene VNTR and BCL2-Associated Athanogene 6 (BAG6) gene rs3117582 polymorphisms were genotyped in a case-control study with 474 large-joint (hip and knee) osteoarthritis (OA) patients and 568 controls in Croatian population by candidate-gene approach for association with OA. We found that BAG6 rs3117582 SNP genotypes were associated with protection (major allele homozygote) and susceptibility (major-minor allele heterozygote) to OA. BAG6 rs3117582 major allele (A) was associated with reduced risk to OA while the minor allele (C) was associated with increased risk to OA. We identified 6 alleles harboring 2 to 7 repeats making 20 genotypes for FAM46A. A rare FAM46A VNTR genotype comprising VNTR alleles with four and seven repeats (c/f) was associated with increased OA risk in both genders. The genotype with four and six repeats (c/e) was also associated with increased risk to OA in males. A polymorphic FAM46A allele with six repeats (e) was associated with reduced risk to OA in females. Our results suggest association between the FAM46A gene, BAG6 gene and OA in Croatian population, respectively. This is the first study to show associations between these genetic loci and OA.

The accuracy, feasibility and challenges of sequencing short tandem repeats using next-generation sequencing platforms

To date we have little knowledge of how accurate next-generation sequencing (NGS) technologies are in sequencing repetitive sequences beyond known limitations to accurately sequence homopolymers. Only a handful of previous reports have evaluated the potential of NGS for sequencing short tandem repeats (microsatellites) and no empirical study has compared and evaluated the performance of more than one NGS platform with the same dataset. Here we examined yeast microsatellite variants from both long-read (454-sequencing) and short-read (Illumina) NGS platforms and compared these to data derived through Sanger sequencing. In addition, we investigated any locus-specific biases and differences that might have resulted from variability in microsatellite repeat number, repeat motif or type of mutation. Out of 112 insertion/deletion variants identified among 45 microsatellite amplicons in our study, we found 87.5% agreement between the 454-platform and Sanger sequencing in frequency of variant detection after Benjamini-Hochberg correction for multiple tests. For a subset of 21 microsatellite amplicons derived from Illumina sequencing, the results of short-read platform were highly consistent with the other two platforms, with 100% agreement with 454-sequencing and 93.6% agreement with the Sanger method after Benjamini-Hochberg correction. We found that the microsatellite attributes copy number, repeat motif and type of mutation did not have a significant effect on differences seen between the sequencing platforms. We show that both long-read and short-read NGS platforms can be used to sequence short tandem repeats accurately, which makes it feasible to consider the use of these platforms in high-throughput genotyping. It appears the major requirement for achieving both high accuracy and rare variant detection in microsatellite genotyping is sufficient read depth coverage. This might be a challenge because each platform generates a consistent pattern of non-uniform sequence coverage, which, as our study suggests, may affect some types of tandem repeats more than others.

Research review: Polygenic methods and their application to psychiatric traits

BACKGROUND

Despite evidence from twin and family studies for an important contribution of genetic factors to both childhood and adult onset psychiatric disorders, identifying robustly associated specific DNA variants has proved challenging. In the pregenomics era the genetic architecture (number, frequency and effect size of risk variants) of complex genetic disorders was unknown. Empirical evidence for the genetic architecture of psychiatric disorders is emerging from the genetic studies of the last 5 years.

METHODS AND SCOPE

We review the methods investigating the polygenic nature of complex disorders. We provide mini-guides to genomic profile (or polygenic) risk scoring and to estimation of variance (or heritability) from common SNPs; a glossary of key terms is also provided. We review results of applications of the methods to psychiatric disorders and related traits and consider how these methods inform on missing heritability, hidden heritability and still-missing heritability.

FINDINGS

Genome-wide genotyping and sequencing studies are providing evidence that psychiatric disorders are truly polygenic, that is they have a genetic architecture of many genetic variants, including risk variants that are both common and rare in the population. Sample sizes published to date are mostly underpowered to detect effect sizes of the magnitude presented by nature, and these effect sizes may be constrained by the biological validity of the diagnostic constructs.

CONCLUSIONS

Increasing the sample size for genome wide association studies of psychiatric disorders will lead to the identification of more associated genetic variants, as already found for schizophrenia. These loci provide the starting point of functional analyses that might eventually lead to new prevention and treatment options and to improved biological validity of diagnostic constructs. Polygenic analyses will contribute further to our understanding of complex genetic traits as sample sizes increase and as sample resources become richer in phenotypic descriptors, both in terms of clinical symptoms and of nongenetic risk factors.

The overdue promise of short tandem repeat variation for heritability

Short tandem repeat (STR) variation has been proposed as a major explanatory factor in the heritability of complex traits in humans and model organisms. However, we still struggle to incorporate STR variation into genotype-phenotype maps. We review here the promise of STRs in contributing to complex trait heritability and highlight the challenges that STRs pose due to their repetitive nature. We argue that STR variants are more likely than single-nucleotide variants to have epistatic interactions, reiterate the need for targeted assays to genotype STRs accurately, and call for more appropriate statistical methods in detecting STR-phenotype associations. Lastly, we suggest that somatic STR variation within individuals may serve as a read-out of disease susceptibility, and is thus potentially a valuable covariate for future association studies.

The conserved PFT1 tandem repeat is crucial for proper flowering in Arabidopsis thaliana

It is widely appreciated that short tandem repeat (STR) variation underlies substantial phenotypic variation in organisms. Some propose that the high mutation rates of STRs in functional genomic regions facilitate evolutionary adaptation. Despite their high mutation rate, some STRs show little to no variation in populations. One such STR occurs in the Arabidopsis thaliana gene PFT1 (MED25), where it encodes an interrupted polyglutamine tract. Although the PFT1 STR is large (∼270 bp), and thus expected to be extremely variable, it shows only minuscule variation across A. thaliana strains. We hypothesized that the PFT1 STR is under selective constraint, due to previously undescribed roles in PFT1 function. We investigated this hypothesis using plants expressing transgenic PFT1 constructs with either an endogenous STR or synthetic STRs of varying length. Transgenic plants carrying the endogenous PFT1 STR generally performed best in complementing a pft1 null mutant across adult PFT1-dependent traits. In stark contrast, transgenic plants carrying a PFT1 transgene lacking the STR phenocopied a pft1 loss-of-function mutant for flowering time phenotypes and were generally hypomorphic for other traits, establishing the functional importance of this domain. Transgenic plants carrying various synthetic constructs occupied the phenotypic space between wild-type and pft1 loss-of-function mutants. By varying PFT1 STR length, we discovered that PFT1 can act as either an activator or repressor of flowering in a photoperiod-dependent manner. We conclude that the PFT1 STR is constrained to its approximate wild-type length by its various functional requirements. Our study implies that there is strong selection on STRs not only to generate allelic diversity, but also to maintain certain lengths pursuant to optimal molecular function.

A systematic evaluation of short tandem repeats in lipid candidate genes: riding on the SNP-wave

Structural genetic variants as short tandem repeats (STRs) are not targeted in SNP-based association studies and thus, their possible association signals are missed. We systematically searched for STRs in gene regions known to contribute to total cholesterol, HDL cholesterol, LDL cholesterol and triglyceride levels in two independent studies (KORA F4, n = 2553 and SAPHIR, n = 1648), resulting in 16 STRs that were finally evaluated. In a combined dataset of both studies, the sum of STR alleles was regressed on each phenotype, adjusted for age and sex. The association analyses were repeated for SNPs in a 200 kb region surrounding the respective STRs in the KORA F4 Study. Three STRs were significantly associated with total cholesterol (within LDLR, the APOA1/C3/A4/A5/BUD13 gene region and ABCG5/8), five with HDL cholesterol (3 within CETP, one in LPL and one inAPOA1/C3/A4/A5/BUD13), three with LDL cholesterol (LDLR, ABCG5/8 and CETP) and two with triglycerides (APOA1/C3/A4/A5/BUD13 and LPL). None of the investigated STRs, however, showed a significant association after adjusting for the lead or adjacent SNPs within that gene region. The evaluated STRs were found to be well tagged by the lead SNP within the respective gene regions. Therefore, the STRs reflect the association signals based on surrounding SNPs. In conclusion, none of the STRs contributed additionally to the SNP-based association signals identified in GWAS on lipid traits.

Microsatellite interruptions stabilize primate genomes and exist as population-specific single nucleotide polymorphisms within individual human genomes

Interruptions of microsatellite sequences impact genome evolution and can alter disease manifestation. However, human polymorphism levels at interrupted microsatellites (iMSs) are not known at a genome-wide scale, and the pathways for gaining interruptions are poorly understood. Using the 1000 Genomes Phase-1 variant call set, we interrogated mono-, di-, tri-, and tetranucleotide repeats up to 10 units in length. We detected ∼26,000–40,000 iMSs within each of four human population groups (African, European, East Asian, and American). We identified population-specific iMSs within exonic regions, and discovered that known disease-associated iMSs contain alleles present at differing frequencies among the populations. By analyzing longer microsatellites in primate genomes, we demonstrate that single interruptions result in a genome-wide average two- to six-fold reduction in microsatellite mutability, as compared with perfect microsatellites. Centrally located interruptions lowered mutability dramatically, by two to three orders of magnitude. Using a biochemical approach, we tested directly whether the mutability of a specific iMS is lower because of decreased DNA polymerase strand slippage errors. Modeling the adenomatous polyposis coli tumor suppressor gene sequence, we observed that a single base substitution interruption reduced strand slippage error rates five- to 50-fold, relative to a perfect repeat, during synthesis by DNA polymerases α, β, or η. Computationally, we demonstrate that iMSs arise primarily by base substitution mutations within individual human genomes. Our biochemical survey of human DNA polymerase α, β, δ, κ, and η error rates within certain microsatellites suggests that interruptions are created most frequently by low fidelity polymerases. Our combined computational and biochemical results demonstrate that iMSs are abundant in human genomes and are sources of population-specific genetic variation that may affect genome stability. The genome-wide identification of iMSs in human populations presented here has important implications for current models describing the impact of microsatellite polymorphisms on gene expression.

Microsatellites are short tandem repeat DNA sequences located throughout the human genome that display a high degree of inter-individual variation. This characteristic makes microsatellites an attractive tool for population genetics and forensics research. Some microsatellites affect gene expression, and mutations within such microsatellites can cause disease. Interruption mutations disrupt the perfect repeated array and are frequently associated with altered disease risk, but they have not been thoroughly studied in human genomes. We identified interrupted mono-, di-, tri- and tetranucleotide MSs (iMS) within individual genomes from African, European, Asian and American population groups. We show that many iMSs, including some within disease-associated genes, are unique to a single population group. By measuring the conservation of microsatellites between human and chimpanzee genomes, we demonstrate that interruptions decrease the probability of microsatellite mutations throughout the genome. We demonstrate that iMSs arise in the human genome by single base changes within the DNA, and provide biochemical data suggesting that these stabilizing changes may be created by error-prone DNA polymerases. Our genome-wide study supports the model in which iMSs act to stabilize individual genomes, and suggests that population-specific differences in microsatellite architecture may be an avenue by which genetic ancestry impacts individual disease risk.

Differential distributions of mononucleotide repeat sequences in 256 viral genomes and its potential implications

Mononucleotide repeats (MNRs) have been systematically investigated in the genomes of eukaryotic and prokaryotic organisms. However, detailed information on the distribution of MNRs in viral genomes is limited. In this study, we examined the distributions of MNRs in 256 fully sequenced virus genomes which showed extensive variations across viral genomes, and is significantly influenced by both genome size and CG content. Furthermore, the ratio of the observed to the expected number of MNRs (O/E ratio) appears to be influenced by both the host range and genome type of a particular virus. Additionally, the densities and frequencies of MNRs in genic regions are lower than in non-coding regions, suggesting that selective pressure acts on viral genomes. We also discuss the potential functional roles that these MNR loci could play in virus genomes. To our knowledge, this is the first analysis focusing on MNRs in viruses, and our study could have potential implications for a deeper understanding of virus genome stability and the co-evolution that occurs between a virus and its host.

Identifying novel interventional strategies for psychiatric disorders: integrating genomics, 'enviromics' and gene-environment interactions in valid preclinical models

Psychiatric disorders affect a substantial proportion of the population worldwide. This high prevalence, combined with the chronicity of the disorders and the major social and economic impacts, creates a significant burden. As a result, an important priority is the development of novel and effective interventional strategies for reducing incidence rates and improving outcomes. This review explores the progress that has been made to date in establishing valid animal models of psychiatric disorders, while beginning to unravel the complex factors that may be contributing to the limitations of current methodological approaches. We propose some approaches for optimizing the validity of animal models and developing effective interventions. We use schizophrenia and autism spectrum disorders as examples of disorders for which development of valid preclinical models, and fully effective therapeutics, have proven particularly challenging. However, the conclusions have relevance to various other psychiatric conditions, including depression, anxiety and bipolar disorders. We address the key aspects of construct, face and predictive validity in animal models, incorporating genetic and environmental factors. Our understanding of psychiatric disorders is accelerating exponentially, revealing extraordinary levels of genetic complexity, heterogeneity and pleiotropy. The environmental factors contributing to individual, and multiple, disorders also exhibit breathtaking complexity, requiring systematic analysis to experimentally explore the environmental mediators and modulators which constitute the 'envirome' of each psychiatric disorder. Ultimately, genetic and environmental factors need to be integrated via animal models incorporating the spatiotemporal complexity of gene-environment interactions and experience-dependent plasticity, thus better recapitulating the dynamic nature of brain development, function and dysfunction.

Association of variable number of tandem repeats in the coding region of the FAM46A gene, FAM46A rs11040 SNP and BAG6 rs3117582 SNP with susceptibility to tuberculosis

We analyzed for association between the Family with sequence similarity 46, member A (FAM46A) gene (located on chromosome 6q14.1), BCL2-Associated Athanogene 6 (BAG6) gene (located on chromosome 6p21.3) and tuberculosis in Croatian Caucasian. We genotyped the FAM46A rs11040 SNP, FAM46A VNTR and BAG6 rs3117582 polymorphisms in a case-control study with 257 tuberculosis patients and 493 healthy individuals in a Croatian Caucasian population. We found that genotype FAM46A 3/3 (three VNTR repeats homozygote) was associated with susceptibility to tuberculosis (p<0.0015, P_corr.<0.029, Odds ratio = 2.42, 95% Confidence Interval = 1.34–4.3). This association suggests that the protein domain encoded by the VNTR might be important for the function of the FAM46A protein, which, in turn, could be relevant in developing tuberculosis. In addition, we found that FAM46A rs11040 SNP:FAM46A VNTR:BAG6 haplotype 132 (G-3-C) is associated with susceptibility to tuberculosis (p<0.012, p_corr.<0.024, Odds ratio 3.45, 95% Confidence Interval = 1.26–9.74). This may suggests that the interaction between the FAM46A and BAG6 proteins may be involved in tuberculosis etiology. We found also that infection of human macrophages with heat-killed M. tuberculosis (H37Rv) led to over-expression of FAM46A (VNTR 3/4) transcript. This is the first study to show associations between the FAM46A gene VNTR polymorphisms, FAM46A rs11040 SNP:FAM46A VNTR:BAG6 haplotypes and any disease.

STaRRRT: a table of short tandem repeats in regulatory regions of the human genome

Background

Tandem repeats (TRs) are unstable regions commonly found within genomes that have consequences for evolution and disease. In humans, polymorphic TRs are known to cause neurodegenerative and neuromuscular disorders as well as being associated with complex diseases such as diabetes and cancer. If present in upstream regulatory regions, TRs can modify chromatin structure and affect transcription; resulting in altered gene expression and protein abundance. The most common TRs are short tandem repeats (STRs), or microsatellites. Promoter located STRs are considerably more polymorphic than coding region STRs. As such, they may be a common driver of phenotypic variation. To study STRs located in regulatory regions, we have performed genome-wide analysis to identify all STRs present in a region that is 2 kilobases upstream and 1 kilobase downstream of the transcription start sites of genes.

Results

The Short Tandem Repeats in Regulatory Regions Table, STaRRRT, contains the results of the genome-wide analysis, outlining the characteristics of 5,264 STRs present in the upstream regulatory region of 4,441 human genes. Gene set enrichment analysis has revealed significant enrichment for STRs in cellular, transcriptional and neurological system gene promoters and genes important in ion and calcium homeostasis. The set of enriched terms has broad similarity to that seen in coding regions, suggesting that regulatory region STRs are subject to similar evolutionary pressures as STRs in coding regions and may, like coding region STRs, have an important role in controlling gene expression.

Conclusions

STaRRRT is a readily-searchable resource for investigating potentially polymorphic STRs that could influence the expression of any gene of interest. The processes and genes enriched for regulatory region STRs provide potential novel targets for diagnosing and treating disease, and support a role for these STRs in the evolution of the human genome.

Evolution and homoplasy at the Bem6 microsatellite locus in three sweetpotato whitefly (Bemisia tabaci) cryptic species

BACKGROUND

The evolution of individual microsatellite loci is often complex and homoplasy is common but often goes undetected. Sequencing alleles at a microsatellite locus can provide a more complete picture of the common evolutionary mechanisms occurring at that locus and can reveal cases of homoplasy. Within species homoplasy can lead to an underestimate of differentiation among populations and among species homoplasy can produce a misleading interpretation regarding shared alleles and hybridization. This is especially problematic with cryptic species.

RESULTS

By sequencing alleles from three cryptic species of the sweetpotato whitefly (Bemisia tabaci), designated MEAM1, MED, and NW, the evolution of the putatively dinucleotide Bem6 (CA₈)imp microsatellite locus is inferred as one of primarily stepwise mutation occurring at four distinct heptaucleotide tandem repeats. In two of the species this pattern yields a compound tandem repeat. Homoplasy was detected both among species and within species.

CONCLUSIONS

In the absence of sequencing, size homoplasious alleles at the Bem6 locus lead to an overestimate of alleles shared and hybridization among cryptic species of Bemisia tabaci. Furthermore, the compound heptanucleotide motif structure of a putative dinucleotide microsatellite has implications for the nomenclature of heptanucleotide tandem repeats with step-wise evolution.

Microsatellite tandem repeats are abundant in human promoters and are associated with regulatory elements

Tandem repeats are genomic elements that are prone to changes in repeat number and are thus often polymorphic. These sequences are found at a high density at the start of human genes, in the gene’s promoter. Increasing empirical evidence suggests that length variation in these tandem repeats can affect gene regulation. One class of tandem repeats, known as microsatellites, rapidly alter in repeat number. Some of the genetic variation induced by microsatellites is known to result in phenotypic variation. Recently, our group developed a novel method for measuring the evolutionary conservation of microsatellites, and with it we discovered that human microsatellites near transcription start sites are often highly conserved. In this study, we examined the properties of microsatellites found in promoters. We found a high density of microsatellites at the start of genes. We showed that microsatellites are statistically associated with promoters using a wavelet analysis, which allowed us to test for associations on multiple scales and to control for other promoter related elements. Because promoter microsatellites tend to be G/C rich, we hypothesized that G/C rich regulatory elements may drive the association between microsatellites and promoters. Our results indicate that CpG islands, G-quadruplexes (G4) and untranslated regulatory regions have highly significant associations with microsatellites, but controlling for these elements in the analysis does not remove the association between microsatellites and promoters. Due to their intrinsic lability and their overlap with predicted functional elements, these results suggest that many promoter microsatellites have the potential to affect human phenotypes by generating mutations in regulatory elements, which may ultimately result in disease. We discuss the potential functions of human promoter microsatellites in this context.

Characterizing social behavior in genetically targeted mouse models of brain disorders

Fragile X syndrome, the leading inherited cause of mental retardation and autism spectrum disorders worldwide, is caused by a tandem repeat expansion in the FMR1 (fragile X mental retardation 1) gene. It presents with a distinct behavioral phenotype which overlaps significantly with that of autism. Emerging evidence suggests that tandem repeat polymorphisms (TRPs) might also play a key role in modulating disease susceptibility for a range of common polygenic disorders, including the broader autism spectrum of disorders (ASD) and other forms of psychiatric illness such as schizophrenia, depression, and bipolar disorder [1]. In order to understand how TRPs and associated gene mutations mediate pathogenesis, various mouse models have been generated. A crucial step in such functional genomics is high-quality behavioral and cognitive phenotyping. This chapter presents a basic behavioral battery for standardized tests for assaying social phenotypes in mouse models of brain disorders, with a focus on aggression.

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.

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.