ZNF91 deletion in human embryonic stem cells leads to ectopic activation of SVA retrotransposons and up-regulation of KRAB zinc finger gene clusters

The arms race of ray-finned fish against the derepression of LTR retroelements

Transposable elements (TEs) are dynamic components of eukaryotic genomes that are important in speciation. Evolutionary pressure has led to the evolution of silencing mechanisms, such as the involvement of Krüppel box-associated zinc finger proteins (KRAB-ZFPs) that regulate TE transcription via TRIM28-mediated NuRD complex recruitment in sarcopterygians. The transcriptional activity of genes that encode proteins of the NuRD complex was also detected in actinopterygians, even though genes encoding crucial proteins such as TRIM28 and KRAB-ZNF are missing in this evolutionary lineage. This study hypothesized that TRIM33 could serve as a replacement for TRIM28 in this lineage and suggested an evolutionary relationship between the sarcopterygian KRAB-ZFPs, the actinopterygian KRAB-like ZNFs and the fish N-terminal zinc-associated domain (FiNZ) ZNFs that are specific to cyprinids. Through coimmunoprecipitation analyses conducted in the basal teleost Anguilla marmorata and the cyprinid Danio rerio, we evaluated the interaction of TRIM33 with the KRAB-like and FiNZ ZFPs. Moreover, we explored the expression profiles of TEs and genes that encode proteins that are potentially involved in their silencing during zebrafish development. Our results may provide new evidence for more thoroughly understanding the evolutionary mechanisms underlying the adaptation of ray-finned fish.

G-quadruplexes in an SVA retrotransposon cause aberrant TAF1 gene expression in X-linked dystonia parkinsonism

G-quadruplexes (G4s) are non-canonical nucleic acid structures that form in guanine (G)-rich genomic regions. X-linked dystonia parkinsonism (XDP) is an inherited neurodegenerative disease in which a SINE-VNTR-Alu (SVA) retrotransposon, characterised by amplification of a G-rich repeat, is inserted into the coding sequence of TAF1, a key partner of RNA polymerase II. XDP SVA alters TAF1 expression, but the cause of this outcome in XDP remains unknown. To assess whether G4s form in XDP SVA and affect TAF1 expression, we first characterised bioinformatically predicted XDP SVA G4s in vitro. We next showed that highly stable G4s can form and stop polymerase amplification at the SVA region from patient-derived fibroblasts and neural progenitor cells. Using chromatin immunoprecipitazion (ChIP) with an anti-G4 antibody coupled to sequencing or quantitative PCR, we showed that XDP SVA G4s are folded even when embedded in a chromatin context in patient-derived cells. Using the G4 ligands BRACO-19 and quarfloxin and total RNA-sequencing analysis, we showed that stabilisation of the XDP SVA G4s reduces TAF1 transcripts downstream and around the SVA, and increases upstream transcripts, while destabilisation using the G4 unfolder PhpC increases TAF1 transcripts. Our data indicate that G4 formation in the XDP SVA is a major cause of aberrant TAF1 expression, opening the way for the development of strategies to unfold G4s and potentially target the disease.

Landscape of evolutionary arms races between transposable elements and KRAB-ZFP family

Transposable elements (TEs) are mobile parasitic sequences that have expanded within the host genome. It has been hypothesized that host organisms have expanded the Krüppel-associated box-containing zinc finger proteins (KRAB-ZFPs), which epigenetically suppress TEs, to counteract disorderly TE transpositions. This process is referred to as the evolutionary arms race. However, the extent to which this evolutionary arms race occurred across various TE families remains unclear. In the present study, we systematically explored the evolutionary arms race between TE families and human KRAB-ZFPs using public ChIP-seq data. We discovered and characterized new instances of evolutionary arms races with KRAB-ZFPs in endogenous retroviruses. Furthermore, we found that the regulatory landscape shaped by this arms race contributed to the gene regulatory networks. In summary, our results provide insight into the impact of the evolutionary arms race on TE families, the KRAB-ZFP family, and host gene regulatory networks.

Mini-heterochromatin domains constrain the cis-regulatory impact of SVA transposons in human brain development and disease

SVA (SINE (short interspersed nuclear element)-VNTR (variable number of tandem repeats)-Alu) retrotransposons remain active in humans and contribute to individual genetic variation. Polymorphic SVA alleles harbor gene regulatory potential and can cause genetic disease. However, how SVA insertions are controlled and functionally impact human disease is unknown. Here we dissect the epigenetic regulation and influence of SVAs in cellular models of X-linked dystonia parkinsonism (XDP), a neurodegenerative disorder caused by an SVA insertion at the TAF1 locus. We demonstrate that the KRAB zinc finger protein ZNF91 establishes H3K9me3 and DNA methylation over SVAs, including polymorphic alleles, in human neural progenitor cells. The resulting mini-heterochromatin domains attenuate the cis-regulatory impact of SVAs. This is critical for XDP pathology; removal of local heterochromatin severely aggravates the XDP molecular phenotype, resulting in increased TAF1 intron retention and reduced expression. Our results provide unique mechanistic insights into how human polymorphic transposon insertions are recognized and how their regulatory impact is constrained by an innate epigenetic defense system.

ZNF91 is an endogenous repressor of the molecular phenotype associated with X-linked dystonia-parkinsonism (XDP)

X-linked dystonia-parkinsonism (XDP) is a severe neurodegenerative disorder resulting from an inherited intronic SINE-Alu-VNTR (SVA) retrotransposon in the TAF1 gene that causes dysregulation of TAF1 transcription. The specific mechanism underlying this dysregulation remains unclear, but it is hypothesized to involve the formation of G-quadruplexes (G4) structures within the XDP-SVA that impede transcription. In this study, we show that ZNF91, a critical repressor of SVA retrotransposons, specifically binds to G4-forming DNA sequences. Further, we found that genetic deletion of ZNF91 exacerbates the molecular phenotype associated with the XDP-SVA insertion in patient cells, while no difference was observed when ZNF91 was deleted from isogenic control cells. Additionally, we observed a significant age-related reduction in ZNF91 expression in whole blood and brain, indicating a progressive loss of repression of the XDP-SVA in XDP. These findings indicate that ZNF91 plays a crucial role in controlling the molecular phenotype associated with XDP. Since ZNF91 binds to G4-forming DNA sequences in SVAs, this suggests that interactions between ZNF91 and G4-forming sequences in the XDP-SVA minimize the severity of the molecular phenotype. Our results showing that ZNF91 expression levels significantly decrease with age provide a potential explanation for the age-related progressive neurodegenerative character of XDP. Collectively, our study provides important insights into the protective role of ZNF91 in XDP pathogenesis and suggests that restoring ZNF91 expression, destabilization of G4s, or targeted repression of the XDP-SVA could be future therapeutic strategies to prevent or treat XDP.

Prevalence of and gene regulatory constraints on transcriptional adaptation in single cells

BACKGROUND

Cells and tissues have a remarkable ability to adapt to genetic perturbations via a variety of molecular mechanisms. Nonsense-induced transcriptional compensation, a form of transcriptional adaptation, has recently emerged as one such mechanism, in which nonsense mutations in a gene trigger upregulation of related genes, possibly conferring robustness at cellular and organismal levels. However, beyond a handful of developmental contexts and curated sets of genes, no comprehensive genome-wide investigation of this behavior has been undertaken for mammalian cell types and conditions. How the regulatory-level effects of inherently stochastic compensatory gene networks contribute to phenotypic penetrance in single cells remains unclear.

RESULTS

We analyze existing bulk and single-cell transcriptomic datasets to uncover the prevalence of transcriptional adaptation in mammalian systems across diverse contexts and cell types. We perform regulon gene expression analyses of transcription factor target sets in both bulk and pooled single-cell genetic perturbation datasets. Our results reveal greater robustness in expression of regulons of transcription factors exhibiting transcriptional adaptation compared to those of transcription factors that do not. Stochastic mathematical modeling of minimal compensatory gene networks qualitatively recapitulates several aspects of transcriptional adaptation, including paralog upregulation and robustness to mutation. Combined with machine learning analysis of network features of interest, our framework offers potential explanations for which regulatory steps are most important for transcriptional adaptation.

CONCLUSIONS

Our integrative approach identifies several putative hits-genes demonstrating possible transcriptional adaptation-to follow-up on experimentally and provides a formal quantitative framework to test and refine models of transcriptional adaptation.

Causes and Consequences of Varying Transposable Element Activity: An Evolutionary Perspective

Transposable elements (TEs) are genomic parasites found in nearly all eukaryotes, including humans. This evolutionary success of TEs is due to their replicative activity, involving insertion into new genomic locations. TE activity varies at multiple levels, from between taxa to within individuals. The rapidly accumulating evidence of the influence of TE activity on human health, as well as the rapid growth of new tools to study it, motivated an evaluation of what we know about TE activity thus far. Here, we discuss why TE activity varies, and the consequences of this variation, from an evolutionary perspective. By studying TE activity in nonhuman organisms in the context of evolutionary theories, we can shed light on the factors that affect TE activity. While the consequences of TE activity are usually deleterious, some have lasting evolutionary impacts by conferring benefits on the host or affecting other evolutionary processes.

Genomic profiling of six human somatic histone H1 variants denotes that H1X accumulates at recently incorporated transposable elements

Histone H1, a vital component in chromatin structure, binds to linker DNA and regulates nuclear processes. We have investigated the distribution of histone H1 variants in a breast cancer cell line using ChIP-Seq. Two major groups of variants are identified: H1.2, H1.3, H1.5 and H1.0 are abundant in low GC regions (B compartment), while H1.4 and H1X preferentially localize in high GC regions (A compartment). Examining their abundance within transposable elements (TEs) reveals that H1X and H1.4 are enriched in recently-incorporated TEs (SVA and SINE-Alu), while H1.0/H1.2/H1.3/H1.5 are more abundant in older elements. Notably, H1X is particularly enriched in SVA families, while H1.4 shows the highest abundance in young AluY elements. Although low GC variants are generally enriched in LINE, LTR and DNA repeats, H1X and H1.4 are also abundant in a subset of recent LINE-L1 and LTR repeats. H1X enrichment at SVA and Alu is consistent across multiple cell lines. Further, H1X depletion leads to TE derepression, suggesting its role in maintaining TE repression. Overall, this study provides novel insights into the differential distribution of histone H1 variants among repetitive elements, highlighting the potential involvement of H1X in repressing TEs recently incorporated within the human genome.

Prevalence of and gene regulatory constraints on transcriptional adaptation in single cells

Cells and tissues have a remarkable ability to adapt to genetic perturbations via a variety of molecular mechanisms. Nonsense-induced transcriptional compensation, a form of transcriptional adaptation, has recently emerged as one such mechanism, in which nonsense mutations in a gene can trigger upregulation of related genes, possibly conferring robustness at cellular and organismal levels. However, beyond a handful of developmental contexts and curated sets of genes, to date, no comprehensive genome-wide investigation of this behavior has been undertaken for mammalian cell types and contexts. Moreover, how the regulatory-level effects of inherently stochastic compensatory gene networks contribute to phenotypic penetrance in single cells remains unclear. Here we combine computational analysis of existing datasets with stochastic mathematical modeling and machine learning to uncover the widespread prevalence of transcriptional adaptation in mammalian systems and the diverse single-cell manifestations of minimal compensatory gene networks. Regulon gene expression analysis of a pooled single-cell genetic perturbation dataset recapitulates important model predictions. Our integrative approach uncovers several putative hits-genes demonstrating possible transcriptional adaptation-to follow up on experimentally, and provides a formal quantitative framework to test and refine models of transcriptional adaptation.

Structural Evolution of Gene Promoters Driven by Primate-Specific KRAB Zinc Finger Proteins

Krüppel-associated box (KRAB) zinc finger proteins (KZNFs) recognize and repress transposable elements (TEs); TEs are DNA elements that are capable of replicating themselves throughout our genomes with potentially harmful consequences. However, genes from this family of transcription factors have a much wider potential for genomic regulation. KZNFs have become integrated into gene-regulatory networks through the control of TEs that function as enhancers and gene promoters; some KZNFs also bind directly to gene promoters, suggesting an additional, more direct layer of KZNF co-option into gene-regulatory networks. Binding site analysis of ZNF519, ZNF441, and ZNF468 suggests the structural evolution of KZNFs to recognize TEs can result in coincidental binding to gene promoters independent of TE sequences. We show a higher rate of sequence turnover in gene promoter KZNF binding sites than neighboring regions, implying a selective pressure is being applied by the binding of a KZNF. Through CRISPR/Cas9 mediated genetic deletion of ZNF519, ZNF441, and ZNF468, we provide further evidence for genome-wide co-option of the KZNF-mediated gene-regulatory functions; KZNF knockout leads to changes in expression of KZNF-bound genes in neuronal lineages. Finally, we show that the opposite can be established upon KZNF overexpression, further strengthening the support for the role of KZNFs as bona-fide gene regulators. With no eminent role for ZNF519 in controlling its TE target, our study may provide a snapshot into the early stages of the completed co-option of a KZNF, showing the lasting, multilayered impact that retrovirus invasions and host response mechanisms can have upon the evolution of our genomes.

Evolutionary dynamics between transposable elements and their host genomes: mechanisms of suppression and escape

Transposable elements (TEs) are ubiquitous among eukaryotic species. Their evolutionary persistence is likely due to a combination of tolerogenic, evasive/antagonistic, and cooperative interactions with their host genomes. Here, we focus on metazoan species and review recent advances related to the harmful effects of TE insertions, including how epigenetic effects and TE-derived RNAs can damage host cells. We discuss new findings related to host pathways that silence TEs, such as the piRNA pathway and the APOBEC3 and Kruppel-associated box zinc finger gene families. Finally, we summarize novel strategies used by TEs to evade host silencing, including the Y chromosome as a permissive niche for TE mobilization and TE counterdefense strategies to block host silencing factors.

CRISPR deletion of a SINE-VNTR-Alu (SVA_67) retrotransposon demonstrates its ability to differentially modulate gene expression at the MAPT locus

Background

SINE-VNTR-Alu (SVA) retrotransposons are hominid-specific elements which have been shown to play important roles in processes such as chromatin structure remodelling and regulation of gene expression demonstrating that these repetitive elements exert regulatory functions. We have previously shown that the presence or absence of a specific SVA element, termed SVA_67, was associated with differential expression of several genes at the MAPT locus, a locus associated with Parkinson's Disease (PD) and frontotemporal dementia. However, we were not able to demonstrate that causation of differential gene expression was directed by the SVA due to lack of functional validation.

Methods

We performed CRISPR to delete SVA_67 in the HEK293 cell line. Quantification of target gene expression was performed using qPCR to assess the effects on expression in response to the deletion of SVA_67. Differences between CRISPR edit and control cell lines were analysed using two-tailed t-test with a minimum 95% confidence interval to determine statistical significance.

Results

In this study, we provide data highlighting the SVA-specific effect on differential gene expression. We demonstrate that the hemizygous deletion of the endogenous SVA_67 in CRISPR edited cell lines was associated with differential expression of several genes at the MAPT locus associated with neurodegenerative diseases including KANSL1, MAPT and LRRC37A.

Discussion

This data is consistent with our previous bioinformatic work of differential gene expression analysis using transcriptomic data from the Parkinson's Progression Markers Initiative (PPMI) cohort. As SVAs have regulatory influences on gene expression, and insertion polymorphisms contribute to interpersonal differences in expression patterns, these results highlight the potential contribution of these elements to complex diseases with potentially many genetic components, such as PD.

Transposable elements drive the evolution of metazoan zinc finger genes

Cys2-His2 zinc finger genes (ZNFs) form the largest family of transcription factors in metazoans. ZNF evolution is highly dynamic and characterized by the rapid expansion and contraction of numerous subfamilies across the animal phylogeny. The forces and mechanisms underlying rapid ZNF evolution remain poorly understood, but there is growing evidence that, in tetrapods, the targeting and repression of lineage-specific transposable elements (TEs) plays a critical role in the evolution of the Krüppel-associated box ZNF (KZNF) subfamily. Currently, it is unknown whether this function and coevolutionary relationship is unique to KZNFs or is a broader feature of metazoan ZNFs. Here, we present evidence that genomic conflict with TEs has been a central driver of the diversification of ZNFs in animals. Sampling from 3221 genome assemblies, we show that the copy number of retroelements correlates with that of ZNFs across at least 750 million years of metazoan evolution. Using computational predictions, we show that ZNFs preferentially bind TEs in diverse animal species. We further investigate the largest ZNF subfamily found in cyprinid fish, which is characterized by a conserved sequence we dubbed the fish N-terminal zinc finger-associated (FiNZ) domain. Zebrafish possess approximately 700 FiNZ-ZNFs, many of which are evolving adaptively under positive selection. Like mammalian KZNFs, most zebrafish FiNZ-ZNFs are expressed at the onset of zygotic genome activation, and blocking their translation using morpholinos during early embryogenesis results in derepression of transcriptionally active TEs. Together, these data suggest that ZNF diversification has been intimately connected to TE expansion throughout animal evolution.

Genetic features and genomic targets of human KRAB-zinc finger proteins

Krüppel-associated box (KRAB) domain-containing zinc finger proteins (KZFPs) are one of the largest groups of transcription factors encoded by tetrapods, with 378 members in human alone. KZFP genes are often grouped in clusters reflecting amplification by gene and segment duplication since the gene family first emerged more than 400 million years ago. Previous work has revealed that many KZFPs recognize transposable element (TE)-embedded sequences as genomic targets, and that KZFPs facilitate the co-option of the regulatory potential of TEs for the benefit of the host. Here, we present a comprehensive survey of the genetic features and genomic targets of human KZFPs, notably completing past analyses by adding data on close to a hundred family members. General principles emerge from our study of the TE-KZFP regulatory system, which point to multipronged evolutionary mechanisms underlaid by highly complex and combinatorial modes of action with strong influences on human speciation.

Hominoid SVA-lncRNA AK057321 targets human-specific SVA retrotransposons in SCN8A and CDK5RAP2 to initiate neuronal maturation

SINE-VNTR-Alu (SVA) retrotransposons arose and expanded in the genome of hominoid primates concurrent with the slowing of brain maturation. We report genes with intronic SVA transposons are enriched for neurodevelopmental disease and transcribed into long non-coding SVA-lncRNAs. Human-specific SVAs in microcephaly CDK5RAP2 and epilepsy SCN8A gene introns repress their expression via transcription factor ZNF91 to delay neuronal maturation. Deleting the SVA in CDK5RAP2 initiates multi-dimensional and in SCN8A selective sodium current neuronal maturation by upregulating these genes. SVA-lncRNA AK057321 forms RNA:DNA heteroduplexes with the genomic SVAs and upregulates these genes to initiate neuronal maturation. SVA-lncRNA AK057321 also promotes species-specific cortex and cerebellum-enriched expression upregulating human genes with intronic SVAs (e.g., HTT, CHAF1B and KCNJ6) but not mouse orthologs. The diversity of neuronal genes with intronic SVAs suggest this hominoid-specific SVA transposon-based gene regulatory mechanism may act at multiple steps to specialize and achieve neoteny of the human brain.

Primate-specific transposable elements shape transcriptional networks during human development

The human genome contains more than 4.5 million inserts derived from transposable elements (TEs), the result of recurrent waves of invasion and internal propagation throughout evolution. For new TE copies to be inherited, they must become integrated in the genome of the germline or pre-implantation embryo, which requires that their source TE be expressed at these stages. Accordingly, many TEs harbor DNA binding sites for the pluripotency factors OCT4, NANOG, SOX2, and KLFs and are transiently expressed during embryonic genome activation. Here, we describe how many primate-restricted TEs have additional binding sites for lineage-specific transcription factors driving their expression during human gastrulation and later steps of fetal development. These TE integrants serve as lineage-specific enhancers fostering the transcription, amongst other targets, of KRAB-zinc finger proteins (KZFPs) of comparable evolutionary age, which in turn corral the activity of TE-embedded regulatory sequences in a similarly lineage-restricted fashion. Thus, TEs and their KZFP controllers play broad roles in shaping transcriptional networks during early human development.

Genomic features underlie the co-option of SVA transposons as cis-regulatory elements in human pluripotent stem cells

Domestication of transposable elements (TEs) into functional cis-regulatory elements is a widespread phenomenon. However, the mechanisms behind why some TEs are co-opted as functional enhancers while others are not are underappreciated. SINE-VNTR-Alus (SVAs) are the youngest group of transposons in the human genome, where ~3,700 copies are annotated, nearly half of which are human-specific. Many studies indicate that SVAs are among the most frequently co-opted TEs in human gene regulation, but the mechanisms underlying such processes have not yet been thoroughly investigated. Here, we leveraged CRISPR-interference (CRISPRi), computational and functional genomics to elucidate the genomic features that underlie SVA domestication into human stem-cell gene regulation. We found that ~750 SVAs are co-opted as functional cis-regulatory elements in human induced pluripotent stem cells. These SVAs are significantly closer to genes and harbor more transcription factor binding sites than non-co-opted SVAs. We show that a long DNA motif composed of flanking YY1/2 and OCT4 binding sites is enriched in the co-opted SVAs and that these two transcription factors bind consecutively on the TE sequence. We used CRISPRi to epigenetically repress active SVAs in stem cell-like NCCIT cells. Epigenetic perturbation of active SVAs strongly attenuated YY1/OCT4 binding and influenced neighboring gene expression. Ultimately, SVA repression resulted in ~3,000 differentially expressed genes, 131 of which were the nearest gene to an annotated SVA. In summary, we demonstrated that SVAs modulate human gene expression, and uncovered that location and sequence composition contribute to SVA domestication into gene regulatory networks.

Active endogenous retroviral elements in human pluripotent stem cells play a role in regulating host gene expression

Human endogenous retroviruses, also called LTR elements, can be bound by transcription factors and marked by different histone modifications in different biological contexts. Recently, individual LTR or certain subclasses of LTRs such as LTR7/HERVH and LTR5_Hs/HERVK families have been identified as cis-regulatory elements. However, there are still many LTR elements with unknown functions. Here, we dissected the landscape of histone modifications and regulatory map of LTRs by integrating 98 ChIP-seq data in human embryonic stem cells (ESCs), and annotated the active LTRs enriching enhancer/promoter-related histone marks. Notably, we found that MER57E3 functionally acted as proximal regulatory element to activate respective ZNF gene. Additionally, HERVK transcript could mainly function in nucleus to activate the adjacent genes. Since LTR5_Hs/LTR5 was bound by many early embryo-specific transcription factors, we further investigated the expression dynamics in different pluripotent states. LTR5_Hs/LTR5/HERVK exhibited higher expression level in naïve ESCs and extended pluripotent stem cells (EPSCs). Functionally, the LTR5_Hs/LTR5 with high activity could serve as a distal enhancer to regulate the host genes. Ultimately, our study not only provides a comprehensive regulatory map of LTRs in human ESCs, but also explores the regulatory models of MER57E3 and LTR5_Hs/LTR5 in host genome.

DNA sequence-dependent formation of heterochromatin nanodomains

The mammalian epigenome contains thousands of heterochromatin nanodomains (HNDs) marked by di- and trimethylation of histone H3 at lysine 9 (H3K9me2/3), which have a typical size of 3–10 nucleosomes. However, what governs HND location and extension is only partly understood. Here, we address this issue by introducing the chromatin hierarchical lattice framework (ChromHL) that predicts chromatin state patterns with single-nucleotide resolution. ChromHL is applied to analyse four HND types in mouse embryonic stem cells that are defined by histone methylases SUV39H1/2 or GLP, transcription factor ADNP or chromatin remodeller ATRX. We find that HND patterns can be computed from PAX3/9, ADNP and LINE1 sequence motifs as nucleation sites and boundaries that are determined by DNA sequence (e.g. CTCF binding sites), cooperative interactions between nucleosomes as well as nucleosome-HP1 interactions. Thus, ChromHL rationalizes how patterns of H3K9me2/3 are established and changed via the activity of protein factors in processes like cell differentiation.

The ability to predict epigenetic regulation is an important challenge in biology. Here the authors describe heterochromatin nanodomains (HNDs) and compare four different types of H3K9me2/3-marked HNDs in mouse embryonic stem cells. They further develop a computational framework to predict genome-wide HND maps from DNA sequence and protein concentrations, at single-nucleotide resolution.

A hidden layer of structural variation in transposable elements reveals potential genetic modifiers in human disease-risk loci

Genome-wide association studies (GWAS) have been highly informative in discovering disease-associated loci but are not designed to capture all structural variations in the human genome. Using long-read sequencing data, we discovered widespread structural variation within SINE-VNTR-Alu (SVA) elements, a class of great ape-specific transposable elements with gene-regulatory roles, which represents a major source of structural variability in the human population. We highlight the presence of structurally variable SVAs (SV-SVAs) in neurological disease-associated loci, and we further associate SV-SVAs to disease-associated SNPs and differential gene expression using luciferase assays and expression quantitative trait loci data. Finally, we genetically deleted SV-SVAs in the BIN1 and CD2AP Alzheimer's disease-associated risk loci and in the BCKDK Parkinson's disease-associated risk locus and assessed multiple aspects of their gene-regulatory influence in a human neuronal context. Together, this study reveals a novel layer of genetic variation in transposable elements that may contribute to identification of the structural variants that are the actual drivers of disease associations of GWAS loci.

Transcriptional states of retroelement-inserted regions and specific KRAB zinc finger protein association are correlated with DNA methylation of retroelements in human male germ cells

DNA methylation, repressive histone modifications, and PIWI-interacting RNAs are essential for controlling retroelement silencing in mammalian germ lines. Dysregulation of retroelement silencing is associated with male sterility. Although retroelement silencing mechanisms have been extensively studied in mouse germ cells, little progress has been made in humans. Here, we show that the Krüppel-associated box domain zinc finger proteins are associated with DNA methylation of retroelements in human primordial germ cells. Further, we show that the hominoid-specific retroelement SINE-VNTR-Alus (SVA) is subjected to transcription-directed de novo DNA methylation during human spermatogenesis. The degree of de novo DNA methylation in SVAs varies among human individuals, which confers significant inter-individual epigenetic variation in sperm. Collectively, our results highlight potential molecular mechanisms for the regulation of retroelements in human male germ cells.