Analysis of REST binding sites with canonical and non-canonical motifs in human cell lines.

BACKGROUND: Repressor element 1 (RE1) silencing transcription factor (REST) is a transcriptional repressor abundantly expressed in aging human brains. It is known to regulate genes associated with oxidative stress, inflammation, and neurological disorders by binding to a canonical form of sequence motif and its non-canonical variations. Although analysis of genomic sequence motifs is crucial to understand transcriptional regulation by transcription factors (TFs), a comprehensive characterization of various forms of RE1 motifs in human cell lines has not been performed. RESULTS: Here, we analyzed 23 ENCODE REST ChIP-seq datasets from diverse human cell lines and identified a non-redundant set of 68,975 loci with ChIP-seq peaks. Our systematic characterization of these binding sites revealed that the canonical form of REST binding motif was found primarily in ChIP-seq peaks shared across multiple cell lines, while non-canonical forms of motifs were identified in both cell-line-specific binding sites and those shared across cell lines. Remarkably, we observed a notable prevalence of non-canonical motifs that corresponded to half segments of the canonical motif. Furthermore, our analysis unveiled the presence of cell-line-specific REST binding patterns, as evidenced by the clustering of ChIP-seq experiments according to their respective cell lines. This observation underscores the cell-line specificity of REST binding at certain genomic loci, implying intricate cell-line-specific regulatory mechanisms. CONCLUSIONS: Overall, our study provides a comprehensive characterization of REST binding motifs in human cell lines and genome-wide RE1 motif profiles. These findings contribute to a deeper understanding of REST-mediated transcriptional regulation and highlight the importance of considering cell-line-specific effects in future investigations.

Somatic cancer driver mutations are enriched and associated with inflammatory states in Alzheimer's disease microglia.

Alzheimer's disease (AD) is an age-associated neurodegenerative disorder characterized by progressive neuronal loss and pathological accumulation of the misfolded proteins amyloid-ß and tau1,2. Neuroinflammation mediated by microglia and brain-resident macrophages plays a crucial role in AD pathogenesis1-5, though the mechanisms by which age, genes, and other risk factors interact remain largely unknown. Somatic mutations accumulate with age and lead to clonal expansion of many cell types, contributing to cancer and many non-cancer diseases6,7. Here we studied somatic mutation in normal aged and AD brains by three orthogonal methods and in three independent AD cohorts. Analysis of bulk RNA sequencing data from 866 samples from different brain regions revealed significantly higher (~two-fold) overall burdens of somatic single-nucleotide variants (sSNVs) in AD brains compared to age-matched controls. Molecular-barcoded deep (>1000X) gene panel sequencing of 311 prefrontal cortex samples showed enrichment of sSNVs and somatic insertions and deletions (sIndels) in cancer driver genes in AD brain compared to control, with recurrent, and often multiple, mutations in genes implicated in clonal hematopoiesis (CH)8,9. Pathogenic sSNVs were enriched in CSF1R+ microglia of AD brains, and the high proportion of microglia (up to 40%) carrying some sSNVs in cancer driver genes suggests mutation-driven microglial clonal expansion (MiCE). Analysis of single-nucleus RNA sequencing (snRNAseq) from temporal neocortex of 62 additional AD cases and controls exhibited nominally increased mosaic chromosomal alterations (mCAs) associated with CH10,11. Microglia carrying mCA showed upregulated pro-inflammatory genes, resembling the transcriptomic features of disease-associated microglia (DAM) in AD. Our results suggest that somatic driver mutations in microglia are common with normal aging but further enriched in AD brain, driving MiCE with inflammatory and DAM signatures. Our findings provide the first insights into microglial clonal dynamics in AD and identify potential new approaches to AD diagnosis and therapy.

ATM-deficiency-induced microglial activation promotes neurodegeneration in ataxia-telangiectasia.

While ATM loss of function has long been identified as the genetic cause of ataxia-telangiectasia (A-T), how it leads to selective and progressive degeneration of cerebellar Purkinje and granule neurons remains unclear. ATM expression is enriched in microglia throughout cerebellar development and adulthood. Here, we find evidence of microglial inflammation in the cerebellum of patients with A-T using single-nucleus RNA sequencing. Pseudotime analysis revealed that activation of A-T microglia preceded upregulation of apoptosis-related genes in granule and Purkinje neurons and that microglia exhibited increased neurotoxic cytokine signaling to granule and Purkinje neurons in A-T. To confirm these findings experimentally, we performed transcriptomic profiling of A-T induced pluripotent stem cell (iPSC)-derived microglia, which revealed cell-intrinsic microglial activation of cytokine production and innate immune response pathways compared to controls. Furthermore, A-T microglia co-culture with either control or A-T iPSC-derived neurons was sufficient to induce cytotoxicity. Taken together, these studies reveal that cell-intrinsic microglial activation may promote neurodegeneration in A-T.

Somatic Mosaicism in Amyotrophic Lateral Sclerosis and Frontotemporal Dementia Reveals Widespread Degeneration from Focal Mutations.

Although mutations in dozens of genes have been implicated in familial forms of amyotrophic lateral sclerosis (fALS) and frontotemporal degeneration (fFTD), most cases of these conditions are sporadic (sALS and sFTD), with no family history, and their etiology remains obscure. We tested the hypothesis that somatic mosaic mutations, present in some but not all cells, might contribute in these cases, by performing ultra-deep, targeted sequencing of 88 genes associated with neurodegenerative diseases in postmortem brain and spinal cord samples from 404 individuals with sALS or sFTD and 144 controls. Known pathogenic germline mutations were found in 20.6% of ALS, and 26.5% of FTD cases. Predicted pathogenic somatic mutations in ALS/FTD genes were observed in 2.7% of sALS and sFTD cases that did not carry known pathogenic or novel germline mutations. Somatic mutations showed low variant allele fraction (typically <2%) and were often restricted to the region of initial discovery, preventing detection through genetic screening in peripheral tissues. Damaging somatic mutations were preferentially enriched in primary motor cortex of sALS and prefrontal cortex of sFTD, mirroring regions most severely affected in each disease. Somatic mutation analysis of bulk RNA-seq data from brain and spinal cord from an additional 143 sALS cases and 23 controls confirmed an overall enrichment of somatic mutations in sALS. Two adult sALS cases were identified bearing pathogenic somatic mutations in DYNC1H1 and LMNA, two genes associated with pediatric motor neuron degeneration. Our study suggests that somatic mutations in fALS/fFTD genes, and in genes associated with more severe diseases in the germline state, contribute to sALS and sFTD, and that mosaic mutations in a small fraction of cells in focal regions of the nervous system can ultimately result in widespread degeneration.

Impaired neural stress resistance and loss of REST in bipolar disorder.

Neurodevelopmental changes and impaired stress resistance have been implicated in the pathogenesis of bipolar disorder (BD), but the underlying regulatory mechanisms are unresolved. Here we describe a human cerebral organoid model of BD that exhibits altered neural development, elevated neural network activity, and a major shift in the transcriptome. These phenotypic changes were reproduced in cerebral organoids generated from iPS cell lines derived in different laboratories. The BD cerebral organoid transcriptome showed highly significant enrichment for gene targets of the transcriptional repressor REST. This was associated with reduced nuclear REST and REST binding to target gene recognition sites. Reducing the oxygen concentration in organoid cultures to a physiological range ameliorated the developmental phenotype and restored REST expression. These effects were mimicked by treatment with lithium. Reduced nuclear REST and derepression of REST targets genes were also observed in the prefrontal cortex of BD patients. Thus, an impaired cellular stress response in BD cerebral organoids leads to altered neural development and transcriptional dysregulation associated with downregulation of REST. These findings provide a new model and conceptual framework for exploring the molecular basis of BD.

Pan-cancer analysis reveals multifaceted roles of retrotransposon-fusion RNAs.

Transposon-derived transcripts are abundant in RNA sequences, yet their landscape and function, especially for fusion transcripts derived from unannotated or somatically acquired transposons, remains underexplored. Here, we developed a new bioinformatic tool to detect transposon-fusion transcripts in RNA-sequencing data and performed a pan-cancer analysis of 10,257 cancer samples across 34 cancer types as well as 3,088 normal tissue samples. We identified 52,277 cancer-specific fusions with ~30 events per cancer and hotspot loci within transposons vulnerable to fusion formation. Exonization of intronic transposons was the most prevalent genic fusions, while somatic L1 insertions constituted a small fraction of cancer-specific fusions. Source L1s and HERVs, but not Alus showed decreased DNA methylation in cancer upon fusion formation. Overall cancer-specific L1 fusions were enriched in tumor suppressors while Alu fusions were enriched in oncogenes, including recurrent Alu fusions in EZH2 predictive of patient survival. We also demonstrated that transposon-derived peptides triggered CD8+ T-cell activation to the extent comparable to EBV viruses. Our findings reveal distinct epigenetic and tumorigenic mechanisms underlying transposon fusions across different families and highlight transposons as novel therapeutic targets and the source of potent neoantigens.

Contribution and therapeutic implications of retroelement insertions in ataxia telangiectasia.

Certain classes of genetic variation still escape detection in clinical sequencing analysis. One such class is retroelement insertion, which has been reported as a cause of Mendelian diseases and may offer unique therapeutic implications. Here, we conducted retroelement profiling on whole-genome sequencing data from a cohort of 237 individuals with ataxia telangiectasia (A-T). We found 15 individuals carrying retroelement insertions in ATM, all but one of which integrated in noncoding regions. Systematic functional characterization via RNA sequencing, RT-PCR, and/or minigene splicing assays showed that 12 out of 14 intronic insertions led or contributed to ATM loss of function by exon skipping or activating cryptic splice sites. We also present proof-of-concept antisense oligonucleotides that suppress cryptic exonization caused by a deep intronic retroelement insertion. These results provide an initial systematic estimate of the contribution of retroelements to the genetic architecture of recessive Mendelian disorders as ∼2.1%-5.5%. Our study highlights the importance of retroelement insertions as causal variants and therapeutic targets in genetic diseases.

A framework for individualized splice-switching oligonucleotide therapy.

Splice-switching antisense oligonucleotides (ASOs) could be used to treat a subset of individuals with genetic diseases1, but the systematic identification of such individuals remains a challenge. Here we performed whole-genome sequencing analyses to characterize genetic variation in 235 individuals (from 209 families) with ataxia-telangiectasia, a severely debilitating and life-threatening recessive genetic disorder2,3, yielding a complete molecular diagnosis in almost all individuals. We developed a predictive taxonomy to assess the amenability of each individual to splice-switching ASO intervention; 9% and 6% of the individuals had variants that were 'probably' or 'possibly' amenable to ASO splice modulation, respectively. Most amenable variants were in deep intronic regions that are inaccessible to exon-targeted sequencing. We developed ASOs that successfully rescued mis-splicing and ATM cellular signalling in patient fibroblasts for two recurrent variants. In a pilot clinical study, one of these ASOs was used to treat a child who had been diagnosed with ataxia-telangiectasia soon after birth, and showed good tolerability without serious adverse events for three years. Our study provides a framework for the prospective identification of individuals with genetic diseases who might benefit from a therapeutic approach involving splice-switching ASOs.

Contribution of Somatic Ras/Raf/Mitogen-Activated Protein Kinase Variants in the Hippocampus in Drug-Resistant Mesial Temporal Lobe Epilepsy.

Importance: Mesial temporal lobe epilepsy (MTLE) is the most common focal epilepsy subtype and is often refractory to antiseizure medications. While most patients with MTLE do not have pathogenic germline genetic variants, the contribution of postzygotic (ie, somatic) variants in the brain is unknown. Objective: To test the association between pathogenic somatic variants in the hippocampus and MTLE. Design, Setting, and Participants: This case-control genetic association study analyzed the DNA derived from hippocampal tissue of neurosurgically treated patients with MTLE and age-matched and sex-matched neurotypical controls. Participants treated at level 4 epilepsy centers were enrolled from 1988 through 2019, and clinical data were collected retrospectively. Whole-exome and gene-panel sequencing (each genomic region sequenced more than 500 times on average) were used to identify candidate pathogenic somatic variants. A subset of novel variants was functionally evaluated using cellular and molecular assays. Patients with nonlesional and lesional (mesial temporal sclerosis, focal cortical dysplasia, and low-grade epilepsy-associated tumors) drug-resistant MTLE who underwent anterior medial temporal lobectomy were eligible. All patients with available frozen tissue and appropriate consents were included. Control brain tissue was obtained from neurotypical donors at brain banks. Data were analyzed from June 2020 to August 2022. Exposures: Drug-resistant MTLE. Main Outcomes and Measures: Presence and abundance of pathogenic somatic variants in the hippocampus vs the unaffected temporal neocortex. Results: Of 105 included patients with MTLE, 53 (50.5%) were female, and the median (IQR) age was 32 (26-44) years; of 30 neurotypical controls, 11 (36.7%) were female, and the median (IQR) age was 37 (18-53) years. Eleven pathogenic somatic variants enriched in the hippocampus relative to the unaffected temporal neocortex (median [IQR] variant allele frequency, 1.92 [1.5-2.7] vs 0.3 [0-0.9]; P = .01) were detected in patients with MTLE but not in controls. Ten of these variants were in PTPN11, SOS1, KRAS, BRAF, and NF1, all predicted to constitutively activate Ras/Raf/mitogen-activated protein kinase (MAPK) signaling. Immunohistochemical studies of variant-positive hippocampal tissue demonstrated increased Erk1/2 phosphorylation, indicative of Ras/Raf/MAPK activation, predominantly in glial cells. Molecular assays showed abnormal liquid-liquid phase separation for the PTPN11 variants as a possible dominant gain-of-function mechanism. Conclusions and Relevance: Hippocampal somatic variants, particularly those activating Ras/Raf/MAPK signaling, may contribute to the pathogenesis of sporadic, drug-resistant MTLE. These findings may provide a novel genetic mechanism and highlight new therapeutic targets for this common indication for epilepsy surgery.

A cell type-aware framework for nominating non-coding variants in Mendelian regulatory disorders.

Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN) cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generated single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. Seventy-five percent of elements (44 of 59) validated in an in vivo transgenic reporter assay, demonstrating that single cell accessibility is a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieved significant reduction in our variant search space and nominated candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1, and EBF3 - as well as new candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work provides novel non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

Prevalence and mechanisms of somatic deletions in single human neurons during normal aging and in DNA repair disorders.

Replication errors and various genotoxins cause DNA double-strand breaks (DSBs) where error-prone repair creates genomic mutations, most frequently focal deletions, and defective repair may lead to neurodegeneration. Despite its pathophysiological importance, the extent to which faulty DSB repair alters the genome, and the mechanisms by which mutations arise, have not been systematically examined reflecting ineffective methods. Here, we develop PhaseDel, a computational method to detect focal deletions and characterize underlying mechanisms in single-cell whole genome sequences (scWGS). We analyzed high-coverage scWGS of 107 single neurons from 18 neurotypical individuals of various ages, and found that somatic deletions increased with age and in highly expressed genes in human brain. Our analysis of 50 single neurons from DNA repair-deficient diseases with progressive neurodegeneration (Cockayne syndrome, Xeroderma pigmentosum, and Ataxia telangiectasia) reveals elevated somatic deletions compared to age-matched controls. Distinctive mechanistic signatures and transcriptional associations suggest roles for somatic deletions in neurodegeneration.

Somatic mutations in single human cardiomyocytes reveal age-associated DNA damage and widespread oxidative genotoxicity.

The accumulation of somatic DNA mutations over time is a hallmark of aging in many dividing and nondividing cells but has not been studied in postmitotic human cardiomyocytes. Using single-cell whole-genome sequencing, we identified and characterized the landscape of somatic single-nucleotide variants (sSNVs) in 56 single cardiomyocytes from 12 individuals (aged from 0.4 to 82 years). Cardiomyocyte sSNVs accumulate with age at rates that are faster than in many dividing cell types and nondividing neurons. Cardiomyocyte sSNVs show distinctive mutational signatures that implicate failed nucleotide excision repair and base excision repair of oxidative DNA damage, and defective mismatch repair. Since age-accumulated sSNVs create many damaging mutations that disrupt gene functions, polyploidization in cardiomyocytes may provide a mechanism of genetic compensation to minimize the complete knockout of essential genes during aging. Age-related accumulation of cardiac mutations provides a paradigm to understand the influence of aging on cardiac dysfunction.

A neurodevelopmental disorder caused by a novel de novo SVA insertion in exon 13 of the SRCAP gene.

Pathogenic variants in the SRCAP (SNF2-related CREBBP activator protein) gene, which encodes a chromatin-remodeling ATPase, cause neurodevelopmental disorders including Floating Harbor syndrome (FLHS). Here, we report the discovery of a de novo transposon insertion in SRCAP exon 13 from trio genome sequencing in a 28-year-old female with failure to thrive, developmental delay, mood disorder and seizure disorder. The insertion was a full-length (~2.8 kb), antisense-oriented SVA insertion relative to the SRCAP transcript, bearing a 5' transduction and hallmarks of target-primed reverse transcription. The 20-bp 5' transduction allowed us to trace the source SVA element to an intron of a long non-coding RNA on chromosome 12, which is highly expressed in testis. RNA sequencing and qRT-PCR confirmed significant depletion of SRCAP expression and low-level exon skipping in the proband. This case highlights a novel disease-causing structural variant and the importance of transposon analysis in a clinical diagnostic setting.

Somatic genomic changes in single Alzheimer's disease neurons.

Dementia in Alzheimer's disease progresses alongside neurodegeneration1-4, but the specific events that cause neuronal dysfunction and death remain poorly understood. During normal ageing, neurons progressively accumulate somatic mutations5 at rates similar to those of dividing cells6,7 which suggests that genetic factors, environmental exposures or disease states might influence this accumulation5. Here we analysed single-cell whole-genome sequencing data from 319 neurons from the prefrontal cortex and hippocampus of individuals with Alzheimer's disease and neurotypical control individuals. We found that somatic DNA alterations increase in individuals with Alzheimer's disease, with distinct molecular patterns. Normal neurons accumulate mutations primarily in an age-related pattern (signature A), which closely resembles 'clock-like' mutational signatures that have been previously described in healthy and cancerous cells6-10. In neurons affected by Alzheimer's disease, additional DNA alterations are driven by distinct processes (signature C) that highlight C>A and other specific nucleotide changes. These changes potentially implicate nucleotide oxidation4,11, which we show is increased in Alzheimer's-disease-affected neurons in situ. Expressed genes exhibit signature-specific damage, and mutations show a transcriptional strand bias, which suggests that transcription-coupled nucleotide excision repair has a role in the generation of mutations. The alterations in Alzheimer's disease affect coding exons and are predicted to create dysfunctional genetic knockout cells and proteostatic stress. Our results suggest that known pathogenic mechanisms in Alzheimer's disease may lead to genomic damage to neurons that can progressively impair function. The aberrant accumulation of DNA alterations in neurodegeneration provides insight into the cascade of molecular and cellular events that occurs in the development of Alzheimer's disease.

Persistent DNA damage associated with ATM kinase deficiency promotes microglial dysfunction.

The autosomal recessive genome instability disorder Ataxia-telangiectasia, caused by mutations in ATM kinase, is characterized by the progressive loss of cerebellar neurons. We find that DNA damage associated with ATM loss results in dysfunctional behaviour of human microglia, immune cells of the central nervous system. Microglial dysfunction is mediated by the pro-inflammatory RELB/p52 non-canonical NF-κB transcriptional pathway and leads to excessive phagocytic clearance of neuronal material. Activation of the RELB/p52 pathway in ATM-deficient microglia is driven by persistent DNA damage and is dependent on the NIK kinase. Activation of non-canonical NF-κB signalling is also observed in cerebellar microglia of individuals with Ataxia-telangiectasia. These results provide insights into the underlying mechanisms of aberrant microglial behaviour in ATM deficiency, potentially contributing to neurodegeneration in Ataxia-telangiectasia.

Rates and Patterns of Clonal Oncogenic Mutations in the Normal Human Brain.

Although oncogenic mutations have been found in nondiseased, proliferative nonneural tissues, their prevalence in the human brain is unknown. Targeted sequencing of genes implicated in brain tumors in 418 samples derived from 110 individuals of varying ages, without tumor diagnoses, detected oncogenic somatic single-nucleotide variants (sSNV) in 5.4% of the brains, including IDH1 R132H. These mutations were largely present in subcortical white matter and enriched in glial cells and, surprisingly, were less common in older individuals. A depletion of high-allele frequency sSNVs representing macroscopic clones with age was replicated by analysis of bulk RNA sequencing data from 1,816 nondiseased brain samples ranging from fetal to old age. We also describe large clonal copy number variants and that sSNVs show mutational signatures resembling those found in gliomas, suggesting that mutational processes of the normal brain drive early glial oncogenesis. This study helps understand the origin and early evolution of brain tumors. SIGNIFICANCE: In the nondiseased brain, clonal oncogenic mutations are enriched in white matter and are less common in older individuals. We revealed early steps in acquiring oncogenic variants, which are essential to understanding brain tumor origins and building new mutational baselines for diagnostics.This article is highlighted in the In This Issue feature, p. 1.

Whole-genome analysis reveals the contribution of non-coding de novo transposon insertions to autism spectrum disorder.

BACKGROUND: Retrotransposons have been implicated as causes of Mendelian disease, but their role in autism spectrum disorder (ASD) has not been systematically defined, because they are only called with adequate sensitivity from whole genome sequencing (WGS) data and a large enough cohort for this analysis has only recently become available. RESULTS: We analyzed WGS data from a cohort of 2288 ASD families from the Simons Simplex Collection by establishing a scalable computational pipeline for retrotransposon insertion detection. We report 86,154 polymorphic retrotransposon insertions-including > 60% not previously reported-and 158 de novo retrotransposition events. The overall burden of de novo events was similar between ASD individuals and unaffected siblings, with 1 de novo insertion per 29, 117, and 206 births for Alu, L1, and SVA respectively, and 1 de novo insertion per 21 births total. However, ASD cases showed more de novo L1 insertions than expected in ASD genes. Additionally, we observed exonic insertions in loss-of-function intolerant genes, including a likely pathogenic exonic insertion in CSDE1, only in ASD individuals. CONCLUSIONS: These findings suggest a modest, but important, impact of intronic and exonic retrotransposon insertions in ASD, show the importance of WGS for their analysis, and highlight the utility of specific bioinformatic tools for high-throughput detection of retrotransposon insertions.

Molecular Signature Expands the Landscape of Driver Negative Thyroid Cancers.

Thyroid cancer is the most common endocrine malignancy. However, the cytological diagnosis of follicular thyroid carcinoma (FTC), Hürthle cell carcinoma (HCC), and follicular variant of papillary thyroid carcinoma (FVPTC) and their benign counterparts is a challenge for preoperative diagnosis. Nearly 20-30% of biopsied thyroid nodules are classified as having indeterminate risk of malignancy and incur costs to the health care system. Based on that, 120 patients were screened for the main driver mutations previously described in thyroid cancer. Subsequently, 14 mutation-negative cases that are the main source of diagnostic errors (FTC, HCC, or FVPTC) underwent RNA-Sequencing analysis. Somatic variants in candidate driver genes (ECD, NUP98,LRP1B, NCOR1, ATM, SOS1, and SPOP) and fusions were described. NCOR1 and SPOP variants underwent validation. Moreover, expression profiling of driver-negative samples was compared to 16 BRAF V600E, RAS, or PAX8-PPARg positive samples. Negative samples were separated in two clusters, following the expression pattern of the RAS/PAX8-PPARg or BRAF V600E positive samples. Both negative groups showed distinct BRS, ERK, and TDS scores, tumor mutation burden, signaling pathways and immune cell profile. Altogether, here we report novel gene variants and describe cancer-related pathways that might impact preoperative diagnosis and provide insights into thyroid tumor biology.

Genomic approaches to trace the history of human brain evolution with an emerging opportunity for transposon profiling of ancient humans.

Transposable elements (TEs) significantly contribute to shaping the diversity of the human genome, and lines of evidence suggest TEs as one of driving forces of human brain evolution. Existing computational approaches, including cross-species comparative genomics and population genetic modeling, can be adapted for the study of the role of TEs in evolution. In particular, diverse ancient and archaic human genome sequences are increasingly available, allowing reconstruction of past human migration events and holding the promise of identifying and tracking TEs among other evolutionarily important genetic variants at an unprecedented spatiotemporal resolution. However, highly degraded short DNA templates and other unique challenges presented by ancient human DNA call for major changes in current experimental and computational procedures to enable the identification of evolutionarily important TEs. Ancient human genomes are valuable resources for investigating TEs in the evolutionary context, and efforts to explore ancient human genomes will potentially provide a novel perspective on the genetic mechanism of human brain evolution and inspire a variety of technological and methodological advances. In this review, we summarize computational and experimental approaches that can be adapted to identify and validate evolutionarily important TEs, especially for human brain evolution. We also highlight strategies that leverage ancient genomic data and discuss unique challenges in ancient transposon genomics.