Disease-Associated Short Tandem Repeats Co-localize with Chromatin Domain Boundaries

The role of GATA4 in mesenchymal stem cell senescence: A new frontier in regenerative medicine

The Mesenchymal Stem Cell (MSC) is a multipotent progenitor cell with known differentiation potential towards various cell lineage, making it an appealing candidate for regenerative medicine. One major contributing factor to age-related MSC dysfunction is cellular senescence, which is the hallmark of relatively irreversible growth arrest and changes in functional properties. GATA4, a zinc-finger transcription factor, emerges as a critical regulator in MSC biology. Originally identified as a key regulator of heart development and specification, GATA4 has since been connected to several aspects of cellular processes, including stem cell proliferation and differentiation. Accumulating evidence suggests that the involvement of GATA4-nuclear signalizing in the process of MSC senescence-related traits may contribute to age-induced alterations in MSC behavior. GATA4 emerged as the central player in MSC senescence, interacting with several signaling pathways. Studies have shown that GATA4 expression is reduced with age in MSCs, which is associated with increased expression levels of senescence markers and impaired regenerative potential. At the mechanistic level, GATA4 regulates the expression of genes involved in cell cycle regulation, DNA repair, and oxidative stress response, thereby influencing the senescence phenotype in MSCs. The findings underscore the critical function of GATA4 in MSC homeostasis and suggest a promising new target to restore stem cell function during aging and disease. A better understanding of the molecular mechanisms that underlie GATA4 mediated modulation of MSC senescence would provide an opportunity to develop new therapies to revitalize old MSCs to increase their regenerative function for therapeutic purposes in regenerative medicine.

Gestational trophoblastic disease: STR genotyping for precision diagnosis

INTRODUCTION

Gestational trophoblastic disease (GTD) encompasses a constellation of rare to common gynecologic conditions stemming from aberrant gestations with distinct genetic backgrounds and variable degrees of trophoblast proliferation of either neoplastic or non-neoplastic nature. GTD is categorized into hydatidiform moles and gestational trophoblastic neoplasms, and their clinical outcomes vary widely across different subtypes. Prompt and accurate diagnosis plays a pivotal role in the effective management and prognostication of patients. Short tandem repeats (STRs) are repetitive DNA sequences dispersed throughout the human genome and inherit a tremendous genetic polymorphism among individuals. Widely recognized for its applications in forensic identity and paternity testing, the relevance of STR genotyping in the diagnosis of GTD has emerged as an essential ancillary test in the classification and management of GTD of both non-neoplastic hydatidiform moles and gestational trophoblastic tumors.

AREA COVERED

This review discusses fundamental principles, laboratory operation, and diagnostic interpretations of STR genotyping in the context of diagnosis and differential diagnosis of GTD. PubMed was searched for all references up to 2024.

EXPERT OPINION

STR genotyping is the gold standard in the diagnosis and subclassification of hydatidiform moles and has an important application in diagnostic workup and risk stratifications of gestational trophoblastic tumors as well.

Machine and Deep Learning Methods for Predicting 3D Genome Organization

Three-dimensional (3D) chromatin interactions, such as enhancer-promoter interactions (EPIs), loops, topologically associating domains (TADs), and A/B compartments, play critical roles in a wide range of cellular processes by regulating gene expression. Recent development of chromatin conformation capture technologies has enabled genome-wide profiling of various 3D structures, even with single cells. However, current catalogs of 3D structures remain incomplete and unreliable due to differences in technology, tools, and low data resolution. Machine learning methods have emerged as an alternative to obtain missing 3D interactions and/or improve resolution. Such methods frequently use genome annotation data (ChIP-seq, DNAse-seq, etc.), DNA sequencing information (k-mers and transcription factor binding site (TFBS) motifs), and other genomic properties to learn the associations between genomic features and chromatin interactions. In this review, we discuss computational tools for predicting three types of 3D interactions (EPIs, chromatin interactions, and TAD boundaries) and analyze their pros and cons. We also point out obstacles to the computational prediction of 3D interactions and suggest future research directions.

Nanopore sequencing as a novel method of characterising anorexia nervosa risk loci

BACKGROUND

Anorexia nervosa (AN) is a polygenic, severe metabopsychiatric disorder with poorly understood aetiology. Eight significant loci have been identified by genome-wide association studies (GWAS) and single nucleotide polymorphism (SNP)-based heritability was estimated to be ~ 11-17, yet causal variants remain elusive. It is therefore important to define the full spectrum of genetic variants in the wider regions surrounding these significantly associated loci. The hypothesis we evaluate here is that unrecognised or relatively unexplored variants in these regions exist and are promising targets for future functional analyses. To test this hypothesis, we implemented a novel approach with targeted nanopore sequencing (Oxford Nanopore Technologies) for 200 kb regions centred on each of the eight AN-associated loci in 10 AN case samples. Our bioinformatics pipeline entailed base-calling and alignment with Dorado and minimap2 software, followed by variant calling with four separate tools, Sniffles2, Clair3, Straglr, and NanoVar. We then leveraged publicly available databases to characterise these loci in putative functional context and prioritise a subset of potentially relevant variants.

RESULTS

Targeted nanopore sequencing effectively enriched the target regions (average coverage 14.64x). To test our hypothesis, we curated a list of 20 prioritised variants in non-coding regions, poorly represented in the current human reference genome but that may have functional consequences in AN pathology. Notably, we identified a polymorphic SINE-VNTR-Alu like sub-family D element (SVA-D), intergenic with IP6K2 and PRKAR2A, and a poly-T short tandem repeat (STR) in the 3'UTR of FOXP1.

CONCLUSIONS

Our results highlight the potential of targeted nanopore sequencing for characterising poorly resolved or complex variation, which may be initially obscured in risk-associated regions detected by GWAS. Some of the variants identified in this way, such as the polymorphic SVA-D and poly-T STR, could contribute to mechanisms of phenotypic risk, through regulation of several neighbouring genes implicated in AN biology, and affect post-transcriptional processing of FOXP1, respectively. This exploratory investigation was not powered to detect functional effects, however, the variants we observed using this method are poorly represented in the current human reference genome and accompanying databases, and further examination of these may provide new opportunities for improved understanding of genetic risk mechanisms of AN.

Beyond the Synapse: FMR1 and FMRP Molecular Mechanisms in the Nucleus

FMR1 (Fragile X messenger ribonucleoprotein 1), located on the X-chromosome, encodes the multi-functional FMR1 protein (FMRP), critical to brain development and function. Trinucleotide CGG repeat expansions at this locus cause a range of neurological disorders, collectively referred to as Fragile X-related conditions. The most well-known of these is Fragile X syndrome, a neurodevelopmental disorder associated with syndromic facial features, autism, intellectual disabilities, and seizures. However, CGG expansions of different sizes also confer a risk of neuropsychiatric and neurodegenerative disorders throughout the lifespan, through distinct molecular mechanisms. Although Fragile X syndrome is associated with downstream synaptic deficits and neuronal hyperexcitability, work in the past decade has demonstrated that both the causative FMR1 trinucleotide repeat expansion and FMRP itself play important roles in nuclear function and regulation, including non-canonical nucleic acid structure formation and chromatin dynamics. These effects are critical to cellular pathophysiology, although the full extent of their contribution to clinical phenotypes is only just emerging. Here, we present a focused review on some of the nuclear consequences of FMR1/FMRP dysregulation, including parallels in other repeat expansion disorders, ranging from studies in model systems to human cells and tissues.

DRED: A Comprehensive Database of Genes Related to Repeat Expansion Diseases

Expansion of tandem repeats in genes often causes severe diseases, such as fragile X syndrome, Huntington's disease, and spinocerebellar ataxia. However, information on genes associated with repeat expansion diseases is scattered throughout the literature, systematic prediction of potential genes that may cause diseases via repeat expansion is also lacking. Here, we develop DRED, a Database of genes related to Repeat Expansion Diseases, as a manually-curated database that covers all known 61 genes related to repeat expansion diseases reported in PubMed and OMIM, along with the detailed repeat information for each gene. DRED also includes 516 genes with the potential to cause diseases via repeat expansion, which were predicted based on their repeat composition, genetic variations, genomic features, and disease associations. Various types of information on repeat expansion diseases and their corresponding genes/repeats are presented in DRED, together with links to external resources, such as NCBI and ClinVar. DRED provides user-friendly interfaces with comprehensive functions, and can serve as a central data resource for basic research and repeat expansion disease-related medical diagnosis. DRED is freely accessible at http://omicslab.genetics.ac.cn/dred, and will be frequently updated to include newly reported genes related to repeat expansion diseases.

Identification of responsible sequences which mutations cause maternal H19-ICR hypermethylation with Beckwith-Wiedemann syndrome-like overgrowth

Beckwith-Wiedemann syndrome (BWS) is caused by a gain of methylation (GOM) at the imprinting control region within the Igf2-H19 domain on the maternal allele (H19-ICR GOM). Mutations in the binding sites of several transcription factors are involved in H19-ICR GOM and BWS. However, the responsible sequence(s) for H19-ICR GOM with BWS-like overgrowth has not been identified in mice. Here, we report that a mutation in the SOX-OCT binding site (SOBS) causes partial H19-ICR GOM, which does not extend beyond CTCF binding site 3 (CTS3). Moreover, simultaneously mutating both SOBS and CTS3 causes complete GOM of the entire H19-ICR, leading to the misexpression of the imprinted genes, and frequent BWS-like overgrowth. In addition, CTS3 is critical for CTCF/cohesin-mediated chromatin conformation. These results indicate that SOBS and CTS3 are the sequences in which mutations cause H19-ICR GOM leading to BWS-like overgrowth and are essential for maintaining the unmethylated state of maternal H19-ICR.

De novo structural variants in autism spectrum disorder disrupt distal regulatory interactions of neuronal genes

Three-dimensional genome organization plays a critical role in gene regulation, and disruptions can lead to developmental disorders by altering the contact between genes and their distal regulatory elements. Structural variants (SVs) can disturb local genome organization, such as the merging of topologically associating domains upon boundary deletion. Testing large numbers of SVs experimentally for their effects on chromatin structure and gene expression is time and cost prohibitive. To address this, we propose a computational approach to predict SV impacts on genome folding, which can help prioritize causal hypotheses for functional testing. We developed a weighted scoring method that measures chromatin contact changes specifically affecting regions of interest, such as regulatory elements or promoters, and implemented it in the SuPreMo-Akita software (Gjoni and Pollard 2024). With this tool, we ranked hundreds of de novo SVs (dnSVs) from autism spectrum disorder (ASD) individuals and their unaffected siblings based on predicted disruptions to nearby neuronal regulatory interactions. This revealed that putative cis-regulatory element interactions (CREints) are more disrupted by dnSVs from ASD probands versus unaffected siblings. We prioritized candidate variants that disrupt ASD CREints and validated our top-ranked locus using isogenic excitatory neurons with and without the dnSV, confirming accurate predictions of disrupted chromatin contacts. This study establishes disrupted genome folding as a potential genetic mechanism in ASD and provides a general strategy for prioritizing variants predicted to disrupt regulatory interactions across tissues.

C9orf72 expansion creates the unstable folate-sensitive fragile site FRA9A

The hyper-unstable Chr9p21 locus, harbouring the interferon gene cluster, oncogenes and C9orf72, is linked to multiple diseases. C9orf72 (GGGGCC)n expansions ( C9orf72 Exp) are associated with incompletely penetrant amyotrophic lateral sclerosis, frontotemporal dementia and autoimmune disorders. C9orf72 Exp patients display hyperactive cGAS-STING-linked interferon immune and DNA damage responses, but the source of immuno-stimulatory or damaged DNA is unknown. Here, we show C9orf72 Exp in pre-symptomatic and ALS-FTD patient cells and brains cause the folate-sensitive chromosomal fragile site, FRA9A. FRA9A centers on >33kb of C9orf72 as highly-compacted chromatin embedded in an 8.2Mb fragility zone spanning 9p21, encompassing 46 genes, making FRA9A one of the largest fragile sites. C9orf72 Exp cells show chromosomal instability, heightened global- and Chr9p-enriched sister-chromatid exchanges, truncated-Chr9s, acentric-Chr9s and Chr9-containing micronuclei, providing endogenous sources of damaged and immunostimulatory DNA. Cells from one C9orf72 Exp patient contained highly-rearranged FRA9A-expressing Chr9 with Chr9-wide dysregulated gene expression. Somatic C9orf72 Exp repeat instability and chromosomal fragility are sensitive to folate-deficiency. Age-dependent repeat instability, chromosomal fragility, and chromosomal instability can be transferred to CNS and peripheral tissues of transgenic C9orf72 Exp mice, implicating C9orf72 Exp as the source. Our results highlight unappreciated effects of C9orf72 expansions that trigger vitamin-sensitive chromosome fragility, adding structural variations to the disease-enriched 9p21 locus, and likely elsewhere.

C9orf72 repeat expansion creates the unstable folate-sensitive fragile site FRA9A

The hyper-unstable Chr9p21 locus, harbouring the interferon gene cluster, oncogenes and C9orf72, is linked to multiple diseases. C9orf72 (GGGGCC)n expansions (C9orf72Exp) are associated with incompletely penetrant amyotrophic lateral sclerosis, frontotemporal dementia and autoimmune disorders. C9orf72Exp patients display hyperactive cGAS-STING-linked interferon immune and DNA damage responses, but the source of immunostimulatory or damaged DNA is unknown. Here, we show C9orf72Exp in pre-symptomatic and amyotrophic lateral sclerosis-frontotemporal dementia patient cells and brains cause the folate-sensitive chromosomal fragile site, FRA9A. FRA9A centers on >33 kb of C9orf72 as highly compacted chromatin embedded in an 8.2 Mb fragility zone spanning 9p21, encompassing 46 genes, making FRA9A one of the largest fragile sites. C9orf72Exp cells show chromosomal instability, heightened global- and Chr9p-enriched sister-chromatid exchanges, truncated-Chr9s, acentric-Chr9s and Chr9-containing micronuclei, providing endogenous sources of damaged and immunostimulatory DNA. Cells from one C9orf72Exp patient contained a highly rearranged FRA9A-expressing Chr9 with Chr9-wide dysregulated gene expression. Somatic C9orf72Exp repeat instability and chromosomal fragility are sensitive to folate deficiency. Age-dependent repeat instability, chromosomal fragility and chromosomal instability can be transferred to CNS and peripheral tissues of transgenic C9orf72Exp mice, implicating C9orf72Exp as the source. Our results highlight unappreciated effects of C9orf72 expansions that trigger vitamin-sensitive chromosome fragility, adding structural variations to the disease-enriched 9p21 locus, and likely elsewhere.

Cell-type-specific loops linked to RNA polymerase II elongation in human neural differentiation

DNA is folded into higher-order structures that shape and are shaped by genome function. The role of long-range loops in the establishment of new gene expression patterns during cell fate transitions remains poorly understood. Here, we investigate the link between cell-specific loops and RNA polymerase II (RNA Pol II) during neural lineage commitment. We find thousands of loops decommissioned or gained de novo upon differentiation of human induced pluripotent stem cells (hiPSCs) to neural progenitor cells (NPCs) and post-mitotic neurons. During hiPSC-to-NPC and NPC-to-neuron transitions, genes changing from RNA Pol II initiation to elongation are >4-fold more likely to anchor cell-specific loops than repressed genes. Elongated genes exhibit significant mRNA upregulation when connected in cell-specific promoter-enhancer loops but not invariant promoter-enhancer loops or promoter-promoter loops or when unlooped. Genes transitioning from repression to RNA Pol II initiation exhibit a slight mRNA increase independent of loop status. Our data link cell-specific loops and robust RNA Pol II-mediated elongation during neural cell fate transitions.

Loops are geometric catalysts for DNA integration

The insertion of DNA elements within genomes underpins both genetic diversity and disease when unregulated. Most of DNA insertions are not random and the physical mechanisms underlying the integration site selection are poorly understood. Here, we perform Molecular Dynamics simulations to study the insertion of DNA elements, such as viral DNA or transposons, into naked DNA or chromatin substrates. More specifically, we explore the role of loops within the polymeric substrate and discover that they act as 'geometric catalysts' for DNA integration by reducing the energy barrier for substrate deformation. Additionally, we discover that the 1D pattern and 3D conformation of loops have a marked effect on the distribution of integration sites. Finally, we show that loops may compete with nucleosomes to attract DNA integrations. These results may be tested in vitro and they may help to understand patterns of DNA insertions with implications in genome evolution and engineering.

Tau beyond Tangles: DNA Damage Response and Cytoskeletal Protein Crosstalk on Neurodegeneration

Neurons in the brain are continuously exposed to various sources of DNA damage. Although the mechanisms of DNA damage repair in mitotic cells have been extensively characterized, the repair pathways in post-mitotic neurons are still largely elusive. Moreover, inaccurate repair can result in deleterious mutations, including deletions, insertions, and chromosomal translocations, ultimately compromising genomic stability. Since neurons are terminally differentiated cells, they cannot employ homologous recombination (HR) for double-strand break (DSB) repair, suggesting the existence of neuron-specific repair mechanisms. Our research has centered on the microtubule-associated protein tau (MAPT), a crucial pathological protein implicated in neurodegenerative diseases, and its interplay with neurons' DNA damage response (DDR). This review aims to provide an updated synthesis of the current understanding of the complex interplay between DDR and cytoskeletal proteins in neurons, with a particular focus on the role of tau in neurodegenerative disorders.

Intertwining roles of R-loops and G-quadruplexes in DNA repair, transcription and genome organization

R-loops are three-stranded nucleic acid structures that are abundant and widespread across the genome and that have important physiological roles in many nuclear processes. Their accumulation is observed in cancers and neurodegenerative disorders. Recent studies have implicated a function for R-loops and G-quadruplex (G4) structures, which can form on the displaced single strand of R-loops, in three-dimensional genome organization in both physiological and pathological contexts. Here we discuss the interconnected functions of DNA:RNA hybrids and G4s within R-loops, their impact on DNA repair and gene regulatory networks, and their emerging roles in genome organization during development and disease.

FISHnet: Detecting chromatin domains in single-cell sequential Oligopaints imaging data

Sequential Oligopaints DNA FISH is an imaging technique that measures higher-order genome folding at single-allele resolution via multiplexed, probe-based tracing. Currently there is a paucity of algorithms to identify 3D genome features in sequential Oligopaints data. Here, we present FISHnet, a graph theory method based on optimization of network modularity to detect chromatin domains and boundaries in pairwise distance matrices. FISHnet uncovers cell type-specific domain-like folding patterns on single alleles, thus enabling future studies aiming to elucidate the role for single-cell folding variation on genome function.

STRAS:a snakemake pipeline for genome-wide short tandem repeats annotation and score

High-throughput whole genome sequencing (WGS) is clinically used in finding single nucleotide variants and small indels. Several bioinformatics tools are developed to call short tandem repeats (STRs) copy numbers from WGS data, such as ExpansionHunter denovo, GangSTR and HipSTR. However, expansion disorders are rare and it is hard to find candidate expansions in single patient sequencing data with ~ 800,000 STRs calls. In this paper I describe a snakemake pipeline for genome-wide STRs Annotation and Score (STRAS) using a Random Forest (RF) model to predict pathogenicity. The predictor was validated by benchmark data from Clinvar and PUBMED. True positive rate was 93.8%. True negative rate was 98.0%.Precision was 98.6% and recall rate was 93.8%. F1-score was 0.961. Sensitivity was 93.8% and specificity was 99.6%. These results showed STRAS could be a useful tool for clinical researchers to find STR loci of interest and filter out neutral STRs. STRAS is freely available at https://github.com/fancheyu5/STRAS .

Multiscale 3D genome organization underlies duck fatty liver with no adipose inflammation or serious injury

Non-alcoholic fatty liver disease (NAFLD) is the most common form of chronic liver disease. Little is known about how gene expression and chromatin structure are regulated in NAFLD due to lack of suitable model. Ducks naturally develop fatty liver similar to serious human non-alcoholic fatty liver (NAFL) without adipose inflammation and liver fibrosis, thus serves as a good model for investigating molecular mechanisms of adipose metabolism and anti-inflammation. Here, we constructed a NAFLD model without adipose inflammation and liver fibrosis in ducks. By performing dynamic pathological and transcriptomic analyses, we identified critical genes involving in regulation of the NF-κB and MHCII signaling, which usually lead to adipose inflammation and liver fibrosis. We further generated dynamic three-dimensional chromatin maps during liver fatty formation and recovery. This showed that ducks enlarged hepatocyte cell nuclei to reduce inter-chromosomal interaction, decompress chromatin structure, and alter strength of intra-TAD and loop interactions during fatty liver formation. These changes partially contributed to the tight control the NF-κB and the MHCII signaling. Our analysis uncovers duck chromatin reorganization might be advantageous to maintain liver regenerative capacity and reduce adipose inflammation. These findings shed light on new strategies for NAFLD control.

Stratifying TAD boundaries pinpoints focal genomic regions of regulation, damage, and repair

Advances in chromatin mapping have exposed the complex chromatin hierarchical organization in mammals, including topologically associating domains (TADs) and their substructures, yet the functional implications of this hierarchy in gene regulation and disease progression are not fully elucidated. Our study delves into the phenomenon of shared TAD boundaries, which are pivotal in maintaining the hierarchical chromatin structure and regulating gene activity. By integrating high-resolution Hi-C data, chromatin accessibility, and DNA double-strand breaks (DSBs) data from various cell lines, we systematically explore the complex regulatory landscape at high-level TAD boundaries. Our findings indicate that these boundaries are not only key architectural elements but also vibrant hubs, enriched with functionally crucial genes and complex transcription factor binding site-clustered regions. Moreover, they exhibit a pronounced enrichment of DSBs, suggesting a nuanced interplay between transcriptional regulation and genomic stability. Our research provides novel insights into the intricate relationship between the 3D genome structure, gene regulation, and DNA repair mechanisms, highlighting the role of shared TAD boundaries in maintaining genomic integrity and resilience against perturbations. The implications of our findings extend to understanding the complexities of genomic diseases and open new avenues for therapeutic interventions targeting the structural and functional integrity of TAD boundaries.

A comprehensive benchmarking with interpretation and operational guidance for the hierarchy of topologically associating domains

Topologically associating domains (TADs), megabase-scale features of chromatin spatial architecture, are organized in a domain-within-domain TAD hierarchy. Within TADs, the inner and smaller subTADs not only manifest cell-to-cell variability, but also precisely regulate transcription and differentiation. Although over 20 TAD callers are able to detect TAD, their usability in biomedicine is confined by a disagreement of outputs and a limit in understanding TAD hierarchy. We compare 13 computational tools across various conditions and develop a metric to evaluate the similarity of TAD hierarchy. Although outputs of TAD hierarchy at each level vary among callers, data resolutions, sequencing depths, and matrices normalization, they are more consistent when they have a higher similarity of larger TADs. We present comprehensive benchmarking of TAD hierarchy callers and operational guidance to researchers of life science researchers. Moreover, by simulating the mixing of different types of cells, we confirm that TAD hierarchy is generated not simply from stacking Hi-C heatmaps of heterogeneous cells. Finally, we propose an air conditioner model to decipher the role of TAD hierarchy in transcription.

SAFB restricts contact domain boundaries associated with L1 chimeric transcription

Long interspersed element-1 (LINE-1 or L1) comprises 17% of the human genome, continuously generates genetic variations, and causes disease in certain cases. However, the regulation and function of L1 remain poorly understood. Here, we uncover that L1 can enrich RNA polymerase IIs (RNA Pol IIs), express L1 chimeric transcripts, and create contact domain boundaries in human cells. This impact of L1 is restricted by a nuclear matrix protein scaffold attachment factor B (SAFB) that recognizes transcriptionally active L1s by binding L1 transcripts to inhibit RNA Pol II enrichment. Acute inhibition of RNA Pol II transcription abolishes the domain boundaries associated with L1 chimeric transcripts, indicating a transcription-dependent mechanism. Deleting L1 impairs domain boundary formation, and L1 insertions during evolution have introduced species-specific domain boundaries. Our data show that L1 can create RNA Pol II-enriched regions that alter genome organization and that SAFB regulates L1 and RNA Pol II activity to preserve gene regulation.

3D epigenomics and 3D epigenopathies

Mammalian genomes are intricately compacted to form sophisticated 3-dimensional structures within the tiny nucleus, so called 3D genome folding. Despite their shapes reminiscent of an entangled yarn, the rapid development of molecular and next-generation sequencing technologies (NGS) has revealed that mammalian genomes are highly organized in a hierarchical order that delicately affects transcription activities. An increasing amount of evidence suggests that 3D genome folding is implicated in diseases, giving us a clue on how to identify novel therapeutic approaches. In this review, we will study what 3D genome folding means in epigenetics, what types of 3D genome structures there are, how they are formed, and how the technologies have developed to explore them. We will also discuss the pathological implications of 3D genome folding. Finally, we will discuss how to leverage 3D genome folding and engineering for future studies. [BMB Reports 2024; 57(5): 216-231].

Enhancer-promoter specificity in gene transcription: molecular mechanisms and disease associations

Although often located at a distance from their target gene promoters, enhancers are the primary genomic determinants of temporal and spatial transcriptional specificity in metazoans. Since the discovery of the first enhancer element in simian virus 40, there has been substantial interest in unraveling the mechanism(s) by which enhancers communicate with their partner promoters to ensure proper gene expression. These research efforts have benefited considerably from the application of increasingly sophisticated sequencing- and imaging-based approaches in conjunction with innovative (epi)genome-editing technologies; however, despite various proposed models, the principles of enhancer-promoter interaction have still not been fully elucidated. In this review, we provide an overview of recent progress in the eukaryotic gene transcription field pertaining to enhancer-promoter specificity. A better understanding of the mechanistic basis of lineage- and context-dependent enhancer-promoter engagement, along with the continued identification of functional enhancers, will provide key insights into the spatiotemporal control of gene expression that can reveal therapeutic opportunities for a range of enhancer-related diseases.

Single-cell multiplex chromatin and RNA interactions in ageing human brain

Dynamically organized chromatin complexes often involve multiplex chromatin interactions and sometimes chromatin-associated RNA1-3. Chromatin complex compositions change during cellular differentiation and ageing, and are expected to be highly heterogeneous among terminally differentiated single cells4-7. Here we introduce the multinucleic acid interaction mapping in single cells (MUSIC) technique for concurrent profiling of multiplex chromatin interactions, gene expression and RNA-chromatin associations within individual nuclei. When applied to 14 human frontal cortex samples from older donors, MUSIC delineated diverse cortical cell types and states. We observed that nuclei exhibiting fewer short-range chromatin interactions were correlated with both an 'older' transcriptomic signature and Alzheimer's disease pathology. Furthermore, the cell type exhibiting chromatin contacts between cis expression quantitative trait loci and a promoter tends to be that in which these cis expression quantitative trait loci specifically affect the expression of their target gene. In addition, female cortical cells exhibit highly heterogeneous interactions between XIST non-coding RNA and chromosome X, along with diverse spatial organizations of the X chromosomes. MUSIC presents a potent tool for exploration of chromatin architecture and transcription at cellular resolution in complex tissues.

Machine and deep learning methods for predicting 3D genome organization

Three-Dimensional (3D) chromatin interactions, such as enhancer-promoter interactions (EPIs), loops, Topologically Associating Domains (TADs), and A/B compartments play critical roles in a wide range of cellular processes by regulating gene expression. Recent development of chromatin conformation capture technologies has enabled genome-wide profiling of various 3D structures, even with single cells. However, current catalogs of 3D structures remain incomplete and unreliable due to differences in technology, tools, and low data resolution. Machine learning methods have emerged as an alternative to obtain missing 3D interactions and/or improve resolution. Such methods frequently use genome annotation data (ChIP-seq, DNAse-seq, etc.), DNA sequencing information (k-mers, Transcription Factor Binding Site (TFBS) motifs), and other genomic properties to learn the associations between genomic features and chromatin interactions. In this review, we discuss computational tools for predicting three types of 3D interactions (EPIs, chromatin interactions, TAD boundaries) and analyze their pros and cons. We also point out obstacles of computational prediction of 3D interactions and suggest future research directions.

Single-cell multiplex chromatin and RNA interactions in aging human brain

The dynamically organized chromatin complexes often involve multiplex chromatin interactions and sometimes chromatin-associated RNA (caRNA) 1-3. Chromatin complex compositions change during cellular differentiation and aging, and are expected to be highly heterogeneous among terminally differentiated single cells 4-7. Here we introduce the Multi-Nucleic Acid Interaction Mapping in Single Cell (MUSIC) technique for concurrent profiling of multiplex chromatin interactions, gene expression, and RNA-chromatin associations within individual nuclei. Applied to 14 human frontal cortex samples from elderly donors, MUSIC delineates diverse cortical cell types and states. We observed the nuclei exhibiting fewer short-range chromatin interactions are correlated with an "older" transcriptomic signature and with Alzheimer's pathology. Furthermore, the cell type exhibiting chromatin contacts between cis expression quantitative trait loci (cis eQTLs) and a promoter tends to be the cell type where these cis eQTLs specifically affect their target gene's expression. Additionally, the female cortical cells exhibit highly heterogeneous interactions between the XIST non-coding RNA and Chromosome X, along with diverse spatial organizations of the X chromosomes. MUSIC presents a potent tool for exploring chromatin architecture and transcription at cellular resolution in complex tissues.

Computational methods for analysing multiscale 3D genome organization

Recent progress in whole-genome mapping and imaging technologies has enabled the characterization of the spatial organization and folding of the genome in the nucleus. In parallel, advanced computational methods have been developed to leverage these mapping data to reveal multiscale three-dimensional (3D) genome features and to provide a more complete view of genome structure and its connections to genome functions such as transcription. Here, we discuss how recently developed computational tools, including machine-learning-based methods and integrative structure-modelling frameworks, have led to a systematic, multiscale delineation of the connections among different scales of 3D genome organization, genomic and epigenomic features, functional nuclear components and genome function. However, approaches that more comprehensively integrate a wide variety of genomic and imaging datasets are still needed to uncover the functional role of 3D genome structure in defining cellular phenotypes in health and disease.

Epigenetic regulatory layers in the 3D nucleus

Nearly 7 decades have elapsed since Francis Crick introduced the central dogma of molecular biology, as part of his ideas on protein synthesis, setting the fundamental rules of sequence information transfer from DNA to RNAs and proteins. We have since learned that gene expression is finely tuned in time and space, due to the activities of RNAs and proteins on regulatory DNA elements, and through cell-type-specific three-dimensional conformations of the genome. Here, we review major advances in genome biology and discuss a set of ideas on gene regulation and highlight how various biomolecular assemblies lead to the formation of structural and regulatory features within the nucleus, with roles in transcriptional control. We conclude by suggesting further developments that will help capture the complex, dynamic, and often spatially restricted events that govern gene expression in mammalian cells.

RExPRT: a machine learning tool to predict pathogenicity of tandem repeat loci

Expansions of tandem repeats (TRs) cause approximately 60 monogenic diseases. We expect that the discovery of additional pathogenic repeat expansions will narrow the diagnostic gap in many diseases. A growing number of TR expansions are being identified, and interpreting them is a challenge. We present RExPRT (Repeat EXpansion Pathogenicity pRediction Tool), a machine learning tool for distinguishing pathogenic from benign TR expansions. Our results demonstrate that an ensemble approach classifies TRs with an average precision of 93% and recall of 83%. RExPRT's high precision will be valuable in large-scale discovery studies, which require prioritization of candidate loci for follow-up studies.

LUSTR: a new customizable tool for calling genome-wide germline and somatic short tandem repeat variants

BACKGROUND

Short tandem repeats (STRs) are widely distributed across the human genome and are associated with numerous neurological disorders. However, the extent that STRs contribute to disease is likely under-estimated because of the challenges calling these variants in short read next generation sequencing data. Several computational tools have been developed for STR variant calling, but none fully address all of the complexities associated with this variant class.

RESULTS

Here we introduce LUSTR which is designed to address some of the challenges associated with STR variant calling by enabling more flexibility in defining STR loci, allowing for customizable modules to tailor analyses, and expanding the capability to call somatic and multiallelic STR variants. LUSTR is a user-friendly and easily customizable tool for targeted or unbiased genome-wide STR variant screening that can use either predefined or novel genome builds. Using both simulated and real data sets, we demonstrated that LUSTR accurately infers germline and somatic STR expansions in individuals with and without diseases.

CONCLUSIONS

LUSTR offers a powerful and user-friendly approach that allows for the identification of STR variants and can facilitate more comprehensive studies evaluating the role of pathogenic STR variants across human diseases.

Spatially coordinated heterochromatinization of long synaptic genes in fragile X syndrome

Short tandem repeat (STR) instability causes transcriptional silencing in several repeat expansion disorders. In fragile X syndrome (FXS), mutation-length expansion of a CGG STR represses FMR1 via local DNA methylation. Here, we find megabase-scale H3K9me3 domains on autosomes and encompassing FMR1 on the X chromosome in FXS patient-derived iPSCs, iPSC-derived neural progenitors, EBV-transformed lymphoblasts, and brain tissue with mutation-length CGG expansion. H3K9me3 domains connect via inter-chromosomal interactions and demarcate severe misfolding of TADs and loops. They harbor long synaptic genes replicating at the end of S phase, replication-stress-induced double-strand breaks, and STRs prone to stepwise somatic instability. CRISPR engineering of the mutation-length CGG to premutation length reverses H3K9me3 on the X chromosome and multiple autosomes, refolds TADs, and restores gene expression. H3K9me3 domains can also arise in normal-length iPSCs created with perturbations linked to genome instability, suggesting their relevance beyond FXS. Our results reveal Mb-scale heterochromatinization and trans interactions among loci susceptible to instability.

TAD evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find 14% of all human TAD boundaries to be shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons compared to species-specific boundaries. CRISPR-Cas9 knockouts of an ultraconserved boundary in a mouse model lead to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in the upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations and showcases the functional importance of TAD evolution.

Genome-wide analysis of the interplay between chromatin-associated RNA and 3D genome organization in human cells

The interphase genome is dynamically organized in the nucleus and decorated with chromatin-associated RNA (caRNA). It remains unclear whether the genome architecture modulates the spatial distribution of caRNA and vice versa. Here, we generate a resource of genome-wide RNA-DNA and DNA-DNA contact maps in human cells. These maps reveal the chromosomal domains demarcated by locally transcribed RNA, hereafter termed RNA-defined chromosomal domains. Further, the spreading of caRNA is constrained by the boundaries of topologically associating domains (TADs), demonstrating the role of the 3D genome structure in modulating the spatial distribution of RNA. Conversely, stopping transcription or acute depletion of RNA induces thousands of chromatin loops genome-wide. Activation or suppression of the transcription of specific genes suppresses or creates chromatin loops straddling these genes. Deletion of a specific caRNA-producing genomic sequence promotes chromatin loops that straddle the interchromosomal target sequences of this caRNA. These data suggest a feedback loop where the 3D genome modulates the spatial distribution of RNA, which in turn affects the dynamic 3D genome organization.

Neuronal DNA double-strand breaks lead to genome structural variations and 3D genome disruption in neurodegeneration

Persistent DNA double-strand breaks (DSBs) in neurons are an early pathological hallmark of neurodegenerative diseases including Alzheimer's disease (AD), with the potential to disrupt genome integrity. We used single-nucleus RNA-seq in human postmortem prefrontal cortex samples and found that excitatory neurons in AD were enriched for somatic mosaic gene fusions. Gene fusions were particularly enriched in excitatory neurons with DNA damage repair and senescence gene signatures. In addition, somatic genome structural variations and gene fusions were enriched in neurons burdened with DSBs in the CK-p25 mouse model of neurodegeneration. Neurons enriched for DSBs also had elevated levels of cohesin along with progressive multiscale disruption of the 3D genome organization aligned with transcriptional changes in synaptic, neuronal development, and histone genes. Overall, this study demonstrates the disruption of genome stability and the 3D genome organization by DSBs in neurons as pathological steps in the progression of neurodegenerative diseases.

Repetitive DNA sequence detection and its role in the human genome

Repetitive DNA sequences playing critical roles in driving evolution, inducing variation, and regulating gene expression. In this review, we summarized the definition, arrangement, and structural characteristics of repeats. Besides, we introduced diverse biological functions of repeats and reviewed existing methods for automatic repeat detection, classification, and masking. Finally, we analyzed the type, structure, and regulation of repeats in the human genome and their role in the induction of complex diseases. We believe that this review will facilitate a comprehensive understanding of repeats and provide guidance for repeat annotation and in-depth exploration of its association with human diseases.

A landscape of complex tandem repeats within individual human genomes

Markedly expanded tandem repeats (TRs) have been correlated with ~60 diseases. TR diversity has been considered a clue toward understanding missing heritability. However, haplotype-resolved long TRs remain mostly hidden or blacked out because their complex structures (TRs composed of various units and minisatellites containing >10-bp units) make them difficult to determine accurately with existing methods. Here, using a high-precision algorithm to determine complex TR structures from long, accurate reads of PacBio HiFi, an investigation of 270 Japanese control samples yields several genome-wide findings. Approximately 322,000 TRs are difficult to impute from the surrounding single-nucleotide variants. Greater genetic divergence of TR loci is significantly correlated with more events of younger replication slippage. Complex TRs are more abundant than single-unit TRs, and a tendency for complex TRs to consist of <10-bp units and single-unit TRs to be minisatellites is statistically significant at loci with ≥500-bp TRs. Of note, 8909 loci with extended TRs (>100b longer than the mode) contain several known disease-associated TRs and are considered candidates for association with disorders. Overall, complex TRs and minisatellites are found to be abundant and diverse, even in genetically small Japanese populations, yielding insights into the landscape of long TRs.

Spatial and temporal organization of the genome: Current state and future aims of the 4D nucleome project

The four-dimensional nucleome (4DN) consortium studies the architecture of the genome and the nucleus in space and time. We summarize progress by the consortium and highlight the development of technologies for (1) mapping genome folding and identifying roles of nuclear components and bodies, proteins, and RNA, (2) characterizing nuclear organization with time or single-cell resolution, and (3) imaging of nuclear organization. With these tools, the consortium has provided over 2,000 public datasets. Integrative computational models based on these data are starting to reveal connections between genome structure and function. We then present a forward-looking perspective and outline current aims to (1) delineate dynamics of nuclear architecture at different timescales, from minutes to weeks as cells differentiate, in populations and in single cells, (2) characterize cis-determinants and trans-modulators of genome organization, (3) test functional consequences of changes in cis- and trans-regulators, and (4) develop predictive models of genome structure and function.

Three-dimensional genome structure and function

Linear DNA undergoes a series of compression and folding events, forming various three-dimensional (3D) structural units in mammalian cells, including chromosomal territory, compartment, topologically associating domain, and chromatin loop. These structures play crucial roles in regulating gene expression, cell differentiation, and disease progression. Deciphering the principles underlying 3D genome folding and the molecular mechanisms governing cell fate determination remains a challenge. With advancements in high-throughput sequencing and imaging techniques, the hierarchical organization and functional roles of higher-order chromatin structures have been gradually illuminated. This review systematically discussed the structural hierarchy of the 3D genome, the effects and mechanisms of cis-regulatory elements interaction in the 3D genome for regulating spatiotemporally specific gene expression, the roles and mechanisms of dynamic changes in 3D chromatin conformation during embryonic development, and the pathological mechanisms of diseases such as congenital developmental abnormalities and cancer, which are attributed to alterations in 3D genome organization and aberrations in key structural proteins. Finally, prospects were made for the research about 3D genome structure, function, and genetic intervention, and the roles in disease development, prevention, and treatment, which may offer some clues for precise diagnosis and treatment of related diseases.

Short Tandem Repeat DNA Profiling Using Perylene-Oligonucleotide Fluorescence Assay

We report an amplification-free genotyping method to determine the number of human short tandem repeats (STRs). DNA-based STR profiling is a robust method for genetic identification purposes such as forensics and biobanking and for identifying specific molecular subtypes of cancer. STR detection requires polymerase amplification, which introduces errors that obscure the correct genotype. We developed a new method that requires no polymerase. First, we synthesized perylene-nucleoside reagents and incorporated them into oligonucleotide probes that recognize five common human STRs. Using these probes and a bead-based hybridization approach, accurate STR detection was achieved in only 1.5 h, including DNA preparation steps, with up to a 1000-fold target DNA enrichment. This method was comparable to PCR-based assays. Using standard fluorometry, the limit of detection was 2.00 ± 0.07 pM for a given target. We used this assay to accurately identify STRs from 50 human subjects, achieving >98% consensus with sequencing data for STR genotyping.

DNA-initiated epigenetic cascades driven by C9orf72 hexanucleotide repeat

The C9orf72 hexanucleotide repeat expansion (HRE) is the most frequent genetic cause of the neurodegenerative diseases amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Here, we describe the pathogenic cascades that are initiated by the C9orf72 HRE DNA. The HRE DNA binds to its protein partner DAXX and promotes its liquid-liquid phase separation, which is capable of reorganizing genomic structures. An HRE-dependent nuclear accumulation of DAXX drives chromatin remodeling and epigenetic changes such as histone hypermethylation and hypoacetylation in patient cells. While regulating global gene expression, DAXX plays a key role in the suppression of basal and stress-inducible expression of C9orf72 via chromatin remodeling and epigenetic modifications of the promoter of the major C9orf72 transcript. Downregulation of DAXX or rebalancing the epigenetic modifications mitigates the stress-induced sensitivity of C9orf72-patient-derived motor neurons. These studies reveal a C9orf72 HRE DNA-dependent regulatory mechanism for both local and genomic architectural changes in the relevant diseases.

3D genome alterations and editing in pathology

The human genome is folded into a multi-level 3D structure that controls many nuclear functions including gene expression. Recently, alterations in 3D genome organization were associated with several genetic diseases and cancer. As a consequence, experimental approaches are now being developed to modify the global 3D genome organization and that of specific loci. Here, we discuss emerging experimental approaches of 3D genome editing that may prove useful in biomedicine.

The motif composition of variable number tandem repeats impacts gene expression

Understanding the impact of DNA variation on human traits is a fundamental question in human genetics. Variable number tandem repeats (VNTRs) make up ∼3% of the human genome but are often excluded from association analysis owing to poor read mappability or divergent repeat content. Although methods exist to estimate VNTR length from short-read data, it is known that VNTRs vary in both length and repeat (motif) composition. Here, we use a repeat-pangenome graph (RPGG) constructed on 35 haplotype-resolved assemblies to detect variation in both VNTR length and repeat composition. We align population-scale data from the Genotype-Tissue Expression (GTEx) Consortium to examine how variations in sequence composition may be linked to expression, including cases independent of overall VNTR length. We find that 9422 out of 39,125 VNTRs are associated with nearby gene expression through motif variations, of which only 23.4% are accessible from length. Fine-mapping identifies 174 genes to be likely driven by variation in certain VNTR motifs and not overall length. We highlight two genes, CACNA1C and RNF213, that have expression associated with motif variation, showing the utility of RPGG analysis as a new approach for trait association in multiallelic and highly variable loci.

The 3D genome and its impacts on human health and disease

Eukaryotic genomes are highly compacted in the cell nucleus. Two loci separated by a long linear distance can be brought into proximity in space through DNA-binding proteins and RNAs, which contributes profoundly to the regulation of gene expression. Recent technology advances have enabled the development and application of the chromosome conformation capture (3C) technique and a host of 3C-based methods that enable genome-scale investigations into changes in chromatin high-order structures during diverse physiological processes and diseases. In this review, we introduce 3C-based technologies and discuss how they can be utilized to glean insights into the impacts of three-dimensional (3D) genome organization in normal physiological and disease processes.

Insights into Enhancer RNAs: Biogenesis and Emerging Role in Brain Diseases

Enhancers are cis-acting elements that control the transcription of target genes and are transcribed into a class of noncoding RNAs (ncRNAs) termed enhancer RNAs (eRNAs). eRNAs have shorter half-lives than mRNAs and long noncoding RNAs; however, the frequency of transcription of eRNAs is close to that of mRNAs. eRNA expression is associated with a high level of histone mark H3K27ac and a low level of H3K27me3. Although eRNAs only account for a small proportion of ncRNAs, their functions are important. eRNAs can not only increase enhancer activity by promoting the formation of enhancer-promoter loops but also regulate transcriptional activation. Increasing numbers of studies have found that eRNAs play an important role in the occurrence and development of brain diseases; however, further research into eRNAs is required. This review discusses the concept, characteristics, classification, function, and potential roles of eRNAs in brain diseases.

Native functions of short tandem repeats

Over a third of the human genome is comprised of repetitive sequences, including more than a million short tandem repeats (STRs). While studies of the pathologic consequences of repeat expansions that cause syndromic human diseases are extensive, the potential native functions of STRs are often ignored. Here, we summarize a growing body of research into the normal biological functions for repetitive elements across the genome, with a particular focus on the roles of STRs in regulating gene expression. We propose reconceptualizing the pathogenic consequences of repeat expansions as aberrancies in normal gene regulation. From this altered viewpoint, we predict that future work will reveal broader roles for STRs in neuronal function and as risk alleles for more common human neurological diseases.

TAD Evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species, and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find that only 14% of all human TAD boundaries are shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons, compared to species-specific boundaries. CRISPR-Cas9 knockouts of two ultraconserved boundaries in mouse models leads to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations, and showcase the functional importance of TAD evolution.

Emerging regulatory mechanisms of noncoding RNAs in topologically associating domains

Topologically associating domains (TADs) are integral to spatial genome organization, instructing gene expression, and cell fate. Recently, several advances have uncovered roles for noncoding RNAs (ncRNAs) in the regulation of the form and function of mammalian TADs. Phase separation has also emerged as a potential arbiter of ncRNAs in the regulation of TADs. In this review we discuss the implications of these novel findings in relation to how ncRNAs might structurally and functionally regulate TADs from two perspectives: moderating loop extrusion through interactions with architectural proteins, and facilitating TAD phase separation. Additionally, we propose future studies and directions to investigate these phenomena.

CasSABER for Programmable In Situ Visualization of Low and Nonrepetitive Gene Loci

Site-specific imaging of target genes using CRISPR probes is essential for understanding the molecular mechanisms of gene function and engineering tools to modulate its downstream pathways. Herein, we develop CRISPR/Cas9-mediated signal amplification by exchange reaction (CasSABER) for programmable in situ imaging of low and nonrepetitive regions of the target gene in the cell nucleus. The presynthesized primer-exchange reaction (PER) probe is able to hybridize multiple fluorophore-bearing imager strands to specifically light up dCas9/sgRNA target-bound gene loci, enabling in situ imaging of fixed cellular gene loci with high specificity and signal-to-noise ratio. In combination with a multiround branching strategy, we successfully detected nonrepetitive gene regions using a single sgRNA. As an intensity-codable and orthogonal probe system, CasSABER enables the adjustable amplification of local signals in fixed cells, resulting in the simultaneous visualization of multicopy and single-copy gene loci with similar fluorescence intensity. Owing to avoiding the complexity of controlling in situ mutistep enzymatic reactions, CasSABER shows good reliability, sensitivity, and ease of implementation, providing a rapid and cost-effective molecular toolkit for studying multigene interaction in fundamental research and gene diagnosis.

Mutation and selection processes regulating short tandem repeats give rise to genetic and phenotypic diversity across species

Short tandem repeats (STRs) are units of 1-6 bp that repeat in a tandem fashion in DNA. Along with single nucleotide polymorphisms and large structural variations, they are among the major genomic variants underlying genetic, and likely phenotypic, divergence. STRs experience mutation rates that are orders of magnitude higher than other well-studied genotypic variants. Frequent copy number changes result in a wide range of alleles, and provide unique opportunities for modulating complex phenotypes through variation in repeat length. While classical studies have identified key roles of individual STR loci, the advent of improved sequencing technology, high-quality genome assemblies for diverse species, and bioinformatics methods for genome-wide STR analysis now enable more systematic study of STR variation across wide evolutionary ranges. In this review, we explore mutation and selection processes that affect STR copy number evolution, and how these processes give rise to varying STR patterns both within and across species. Finally, we review recent examples of functional and adaptive changes linked to STRs.