Variable patterns of retrotransposition in different HeLa strains provide mechanistic insights into SINE RNA mobilization processes.

Alu elements are non-autonomous Short INterspersed Elements (SINEs) derived from the 7SL RNA gene that are present at over one million copies in human genomic DNA. Alu mobilizes by a mechanism known as retrotransposition, which requires the Long INterspersed Element-1 (LINE-1 or L1) ORF2 -encoded protein (ORF2p). Here, we demonstrate that HeLa strains differ in their capacity to support Alu retrotransposition. Human Alu elements retrotranspose efficiently in HeLa-HA and HeLa-CCL2 ( Alu -permissive) strains, but not in HeLa-JVM or HeLa-H1 ( Alu -nonpermissive) strains. A similar pattern of retrotransposition was observed for other 7SL RNA -derived SINEs and tRNA -derived SINEs. In contrast, mammalian LINE-1s, a zebrafish LINE, a human SINE-VNTR - Alu ( SVA ) element, and an L1 ORF1 -containing messenger RNA can retrotranspose in all four HeLa strains. Using an in vitro reverse transcriptase-based assay, we show that Alu RNAs associate with ORF2p and are converted into cDNAs in both Alu -permissive and Alu -nonpermissive HeLa strains, suggesting that 7SL - and tRNA -derived SINE RNAs use strategies to 'hijack' L1 ORF2p that are distinct from those used by SVA elements and ORF1 -containing mRNAs. These data further suggest ORF2p associates with the Alu RNA poly(A) tract in both Alu -permissive and Alu -nonpermissive HeLa strains, but that Alu retrotransposition is blocked after this critical step in Alu -nonpermissive HeLa strains.

Biallelic loss-of-function variants of SLC12A9 cause lysosome dysfunction and a syndromic neurodevelopmental disorder.

PURPOSE: Pathogenic variants of FIG4 generate enlarged lysosomes and neurological and developmental disorders. To identify additional genes regulating lysosomal volume, we carried out a genome-wide activation screen to detect suppression of enlarged lysosomes in FIG4-/- cells. METHODS: The CRISPR-a gene activation screen utilized sgRNAs from the promoters of protein-coding genes. Fluorescence-activated cell sorting separated cells with correction of the enlarged lysosomes from uncorrected cells. Patient variants of SLC12A9 were identified by exome or genome sequencing and studied by segregation analysis and clinical characterization. RESULTS: Overexpression of SLC12A9, a solute co-transporter, corrected lysosomal swelling in FIG4-/- cells. SLC12A9 (NP_064631.2) colocalized with LAMP2 at the lysosome membrane. Biallelic variants of SLC12A9 were identified in 3 unrelated probands with neurodevelopmental disorders. Common features included intellectual disability, skeletal and brain structural abnormalities, congenital heart defects, and hypopigmented hair. Patient 1 was homozygous for nonsense variant p.(Arg615∗), patient 2 was compound heterozygous for p.(Ser109Lysfs∗20) and a large deletion, and proband 3 was compound heterozygous for p.(Glu290Glyfs∗36) and p.(Asn552Lys). Fibroblasts from proband 1 contained enlarged lysosomes that were corrected by wild-type SLC12A9 cDNA. Patient variant p.(Asn552Lys) failed to correct the lysosomal defect. CONCLUSION: Impaired function of SLC12A9 results in enlarged lysosomes and a recessive disorder with a recognizable neurodevelopmental phenotype.

Benchmarking splice variant prediction algorithms using massively parallel splicing assays.

BACKGROUND: Variants that disrupt mRNA splicing account for a sizable fraction of the pathogenic burden in many genetic disorders, but identifying splice-disruptive variants (SDVs) beyond the essential splice site dinucleotides remains difficult. Computational predictors are often discordant, compounding the challenge of variant interpretation. Because they are primarily validated using clinical variant sets heavily biased to known canonical splice site mutations, it remains unclear how well their performance generalizes. RESULTS: We benchmark eight widely used splicing effect prediction algorithms, leveraging massively parallel splicing assays (MPSAs) as a source of experimentally determined ground-truth. MPSAs simultaneously assay many variants to nominate candidate SDVs. We compare experimentally measured splicing outcomes with bioinformatic predictions for 3,616 variants in five genes. Algorithms' concordance with MPSA measurements, and with each other, is lower for exonic than intronic variants, underscoring the difficulty of identifying missense or synonymous SDVs. Deep learning-based predictors trained on gene model annotations achieve the best overall performance at distinguishing disruptive and neutral variants, and controlling for overall call rate genome-wide, SpliceAI and Pangolin have superior sensitivity. Finally, our results highlight two practical considerations when scoring variants genome-wide: finding an optimal score cutoff, and the substantial variability introduced by differences in gene model annotation, and we suggest strategies for optimal splice effect prediction in the face of these issues. CONCLUSION: SpliceAI and Pangolin show the best overall performance among predictors tested, however, improvements in splice effect prediction are still needed especially within exons.

High-Throughput Splicing Assays Identify Known and Novel WT1 Exon 9 Variants in Nephrotic Syndrome.

Introduction: Frasier syndrome (FS) is a rare Mendelian form of nephrotic syndrome (NS) caused by variants which disrupt the proper splicing of WT1. This key transcription factor gene is alternatively spliced at exon 9 to produce 2 isoforms ("KTS+" and "KTS-"), which are normally expressed in the kidney at a ∼2:1 (KTS+:KTS-) ratio. FS results from variants that reduce this ratio by disrupting the splice donor of the KTS+ isoform. FS is extremely rare, and it is unclear whether any variants beyond the 8 already known could cause FS. Methods: To prospectively identify other splicing-disruptive variants, we leveraged a massively parallel splicing assay. We tested every possible single nucleotide variant (n = 519) in and around WT1 exon 9 for effects upon exon inclusion and KTS+/- ratio. Results: Splice disruptive variants (SDVs) made up 11% of the tested point variants overall and were tightly concentrated near the canonical acceptor and the KTS+/- alternate donors. Our map successfully identified all 8 known FS or focal segmental glomerulosclerosis (FSGS) variants and 16 additional novel variants which were comparably disruptive to these known pathogenic variants. We also identified 19 variants that, conversely, increased the KTS+/KTS- ratio, of which 2 are observed in unrelated individuals with 46,XX ovotesticular disorder of sex development (46,XX OTDSD). Conclusion: This splicing effect map can serve as functional evidence to guide the clinical interpretation of newly observed variants in and around WT1 exon 9.

Using a modular massively parallel reporter assay to discover context-specific regulatory grammars in type 2 diabetes.

Recent genome-wide association studies have established that most complex disease-associated loci are found in noncoding regions where defining their function is nontrivial. In this study, we leverage a modular massively parallel reporter assay (MPRA) to uncover sequence features linked to context-specific regulatory activity. We screened enhancer activity across a panel of 198-bp fragments spanning over 10k type 2 diabetes- and metabolic trait-associated variants in the 832/13 rat insulinoma cell line, a relevant model of pancreatic beta cells. We explored these fragments' context sensitivity by comparing their activities when placed up-or downstream of a reporter gene, and in combination with either a synthetic housekeeping promoter (SCP1) or a more biologically relevant promoter corresponding to the human insulin gene ( INS ). We identified clear effects of MPRA construct design on measured fragment enhancer activity. Specifically, a subset of fragments (n = 702/11,656) displayed positional bias, evenly distributed across up- and downstream preference. A separate set of fragments exhibited promoter bias (n = 698/11,656), mostly towards the cell-specific INS promoter (73.4%). To identify sequence features associated with promoter preference, we used Lasso regression with 562 genomic annotations and discovered that fragments with INS promoter-biased activity are enriched for HNF1 motifs. HNF1 family transcription factors are key regulators of glucose metabolism disrupted in maturity onset diabetes of the young (MODY), suggesting genetic convergence between rare coding variants that cause MODY and common T2D-associated regulatory variants. We designed a follow-up MPRA containing HNF1 motif-enriched fragments and observed several instances where deletion or mutation of HNF1 motifs disrupted the INS promoter-biased enhancer activity, specifically in the beta cell model but not in a skeletal muscle cell line, another diabetes-relevant cell type. Together, our study suggests that cell-specific regulatory activity is partially influenced by enhancer-promoter compatibility and indicates that careful attention should be paid when designing MPRA libraries to capture context-specific regulatory processes at disease-associated genetic signals.

Benchmarking splice variant prediction algorithms using massively parallel splicing assays.

Background: Variants that disrupt mRNA splicing account for a sizable fraction of the pathogenic burden in many genetic disorders, but identifying splice-disruptive variants (SDVs) beyond the essential splice site dinucleotides remains difficult. Computational predictors are often discordant, compounding the challenge of variant interpretation. Because they are primarily validated using clinical variant sets heavily biased to known canonical splice site mutations, it remains unclear how well their performance generalizes. Results: We benchmarked eight widely used splicing effect prediction algorithms, leveraging massively parallel splicing assays (MPSAs) as a source of experimentally determined ground-truth. MPSAs simultaneously assay many variants to nominate candidate SDVs. We compared experimentally measured splicing outcomes with bioinformatic predictions for 3,616 variants in five genes. Algorithms' concordance with MPSA measurements, and with each other, was lower for exonic than intronic variants, underscoring the difficulty of identifying missense or synonymous SDVs. Deep learning-based predictors trained on gene model annotations achieved the best overall performance at distinguishing disruptive and neutral variants. Controlling for overall call rate genome-wide, SpliceAI and Pangolin also showed superior overall sensitivity for identifying SDVs. Finally, our results highlight two practical considerations when scoring variants genome-wide: finding an optimal score cutoff, and the substantial variability introduced by differences in gene model annotation, and we suggest strategies for optimal splice effect prediction in the face of these issues. Conclusion: SpliceAI and Pangolin showed the best overall performance among predictors tested, however, improvements in splice effect prediction are still needed especially within exons.

A maternal germline mutator phenotype in a family affected by heritable colorectal cancer.

Variation in DNA repair genes can increase cancer risk by elevating the rate of oncogenic mutation. Defects in one such gene, MUTYH, are known to elevate the incidence of colorectal cancer in a recessive Mendelian manner, and some evidence has also linked MUTYH to elevated incidence of other cancers as well as elevated mutation rates in normal somatic and germline cells. Here, we use whole genome sequencing to measure germline de novo mutation rates in a large extended family affected by pathogenic MUTYH variation and a history of colorectal cancer. Although this family's genotype, p.Y179C/V234M (c.536A>G/700G>A on transcript NM_001128425), contains a variant with conflicting functional interpretations, we use an in vitro cell line assay to determine that it partially attenuates MUTYH's function. In the children of mothers affected by the Y179C/V234M genotype, we identify an elevation of the C>A mutation rate that is weaker than mutator effects previously reported to be caused by other pathogenic MUTYH genotypes, suggesting that mutation rates in normal tissues may be useful for classifying cancer-associated variation along a continuum of severity. Surprisingly, we detect no significant elevation of the C>A mutation rate in children born to a father with the same biallelic MUTYH genotype, despite calculating that we should have adequate power to detect such a mutator effect. This suggests that the oxidative stress repaired by MUTYH may contribute more to female reproductive aging than male reproductive aging in the general population.

Population-scale skeletal muscle single-nucleus multi-omic profiling reveals extensive context specific genetic regulation.

Skeletal muscle, the largest human organ by weight, is relevant to several polygenic metabolic traits and diseases including type 2 diabetes (T2D). Identifying genetic mechanisms underlying these traits requires pinpointing the relevant cell types, regulatory elements, target genes, and causal variants. Here, we used genetic multiplexing to generate population-scale single nucleus (sn) chromatin accessibility (snATAC-seq) and transcriptome (snRNA-seq) maps across 287 frozen human skeletal muscle biopsies representing 456,880 nuclei. We identified 13 cell types that collectively represented 983,155 ATAC summits. We integrated genetic variation to discover 6,866 expression quantitative trait loci (eQTL) and 100,928 chromatin accessibility QTL (caQTL) (5% FDR) across the five most abundant cell types, cataloging caQTL peaks that atlas-level snATAC maps often miss. We identified 1,973 eGenes colocalized with caQTL and used mediation analyses to construct causal directional maps for chromatin accessibility and gene expression. 3,378 genome-wide association study (GWAS) signals across 43 relevant traits colocalized with sn-e/caQTL, 52% in a cell-specific manner. 77% of GWAS signals colocalized with caQTL and not eQTL, highlighting the critical importance of population-scale chromatin profiling for GWAS functional studies. GWAS-caQTL colocalization showed distinct cell-specific regulatory paradigms. For example, a C2CD4A/B T2D GWAS signal colocalized with caQTL in muscle fibers and multiple chromatin loop models nominated VPS13C, a glucose uptake gene. Sequence of the caQTL peak overlapping caSNP rs7163757 showed allelic regulatory activity differences in a human myocyte cell line massively parallel reporter assay. These results illuminate the genetic regulatory architecture of human skeletal muscle at high-resolution epigenomic, transcriptomic, and cell state scales and serve as a template for population-scale multi-omic mapping in complex tissues and traits.

Saturation-scale functional evidence supports clinical variant interpretation in Lynch syndrome.

BACKGROUND: Lynch syndrome (LS) is a cancer predisposition syndrome affecting more than 1 in every 300 individuals worldwide. Clinical genetic testing for LS can be life-saving but is complicated by the heavy burden of variants of uncertain significance (VUS), especially missense changes. RESULT: To address this challenge, we leverage a multiplexed analysis of variant effect (MAVE) map covering >94% of the 17,746 possible missense variants in the key LS gene MSH2. To establish this map's utility in large-scale variant reclassification, we overlay it on clinical databases of >15,000 individuals with LS gene variants uncovered during clinical genetic testing. We validate these functional measurements in a cohort of individuals with paired tumor-normal test results and find that MAVE-based function scores agree with the clinical interpretation for every one of the MSH2 missense variants with an available classification. We use these scores to attempt reclassification for 682 unique missense VUS, among which 34 scored as deleterious by our function map, in line with previously published rates for other cancer predisposition genes. Combining functional data and other evidence, ten missense VUS are reclassified as pathogenic/likely pathogenic, and another 497 could be moved to benign/likely benign. Finally, we apply these functional scores to paired tumor-normal genetic tests and identify a subset of patients with biallelic somatic loss of function, reflecting a sporadic Lynch-like Syndrome with distinct implications for treatment and relatives' risk. CONCLUSION: This study demonstrates how high-throughput functional assays can empower scalable VUS resolution and prospectively generate strong evidence for variant classification.

Gene amplification acts as a molecular foothold to facilitate cross-species adaptation and evasion of multiple antiviral pathways.

Cross-species spillover events are responsible for many of the pandemics in human history including COVID-19; however, the evolutionary mechanisms that enable these events are poorly understood. We have previously modeled this process using a chimeric vaccinia virus expressing the rhesus cytomegalovirus-derived protein kinase R (PKR) antagonist RhTRS1 in place of its native PKR antagonists: E3L and K3L (VACVΔEΔK + RhTRS1). Using this virus, we demonstrated that gene amplification of rhtrs1 occurred early during experimental evolution and was sufficient to fully rescue virus replication in partially resistant African green monkey (AGM) fibroblasts. Notably, this rapid gene amplification also allowed limited virus replication in otherwise completely non-permissive human fibroblasts, suggesting that gene amplification may act as a 'molecular foothold' to facilitate viral adaptation to multiple species. In this study, we demonstrate that there are multiple barriers to VACVΔEΔK + RhTRS1 replication in human cells, mediated by both PKR and ribonuclease L (RNase L). We experimentally evolved three AGM-adapted virus populations in human fibroblasts. Each population adapted to human cells bimodally, via an initial 10-fold increase in replication after only two passages followed by a second 10-fold increase in replication by passage 9. Using our Illumina-based pipeline, we found that some single nucleotide polymorphisms (SNPs) which had evolved during the prior AGM adaptation were rapidly lost, while thirteen single-base substitutions and short indels increased over time, including two SNPs unique to human foreskin fibroblast (HFF)-adapted populations. Many of these changes were associated with components of the viral RNA polymerase, although no variant was shared between all three populations. Taken together, our results demonstrate that rhtrs1 amplification was sufficient to increase viral tropism after passage in an 'intermediate species' and subsequently enabled the virus to adopt different, species-specific adaptive mechanisms to overcome distinct barriers to viral replication in AGM and human cells.

GCAF(TMEM251) regulates lysosome biogenesis by activating the mannose-6-phosphate pathway.

The mannose-6-phosphate (M6P) biosynthetic pathway for lysosome biogenesis has been studied for decades and is considered a well-understood topic. However, whether this pathway is regulated remains an open question. In a genome-wide CRISPR/Cas9 knockout screen, we discover TMEM251 as the first regulator of the M6P modification. Deleting TMEM251 causes mistargeting of most lysosomal enzymes due to their loss of M6P modification and accumulation of numerous undigested materials. We further demonstrate that TMEM251 localizes to the Golgi and is required for the cleavage and activity of GNPT, the enzyme that catalyzes M6P modification. In zebrafish, TMEM251 deletion leads to severe developmental defects including heart edema and skeletal dysplasia, which phenocopies Mucolipidosis Type II. Our discovery provides a mechanism for the newly discovered human disease caused by TMEM251 mutations. We name TMEM251 as GNPTAB cleavage and activity factor (GCAF) and its related disease as Mucolipidosis Type V.

Gene amplification acts as a molecular foothold to facilitate cross-species adaptation and evasion of multiple antiviral pathways.

Cross-species spillover events are responsible for many of the pandemics in human history including COVID-19; however, the evolutionary mechanisms that enable these events are poorly understood. We have previously modeled this process using a chimeric vaccinia virus expressing the rhesus cytomegalovirus-derived PKR antagonist RhTRS1 in place of its native PKR antagonists; E3L and K3L (VACVΔEΔK+RhTRS1). Using this virus, we demonstrated that gene amplification of rhtrs1 occurred early during experimental evolution and was sufficient to fully rescue virus replication in partially resistant African green monkey (AGM) fibroblasts. Notably, this rapid gene amplification also allowed limited virus replication in otherwise completely non-permissive human fibroblasts, suggesting that gene amplification may act as a "molecular foothold" to facilitate viral adaptation to multiple species. In this study, we demonstrate that there are multiple barriers to VACVΔEΔK+RhTRS1 replication in human cells, mediated by both PKR and RNase L. We experimentally evolved three AGM-adapted virus populations in human fibroblasts. Each population adapted to human cells bimodally, via an initial 10-fold increase in replication after only two passages followed by a second 10-fold increase in replication by passage nine. Using our Illumina-based pipeline, we found that some SNPs which had evolved during the prior AGM adaptation were rapidly lost, while 13 single-base substitutions and short indels increased over time, including two SNPs unique to HFF adapted populations. Many of these changes were associated with components of the viral RNA polymerase, although no variant was shared between all three populations. Taken together, our results demonstrate that rhtrs1 amplification was sufficient to increase viral tropism after passage in an "intermediate species" and subsequently enabled the virus to adopt different, species-specific adaptive mechanisms to overcome distinct barriers to viral replication in AGM and human cells.

Longitudinal profiling of clonal hematopoiesis provides insight into clonal dynamics.

BACKGROUND: Clonal hematopoiesis of indeterminate potential (CHIP), the age-related expansion of mutant hematopoietic stem cells, confers risk for multiple diseases of aging including hematologic cancer and cardiovascular disease. Whole-exome or genome sequencing can detect CHIP, but due to those assays' high cost, most population studies have been cross-sectional, sequencing only a single timepoint per individual. RESULTS: We developed and validated a cost-effective single molecule molecular inversion probe sequencing (smMIPS) assay for detecting CHIP, targeting the 11 most frequently mutated genes in CHIP along with 4 recurrent mutational hotspots. We sequenced 548 multi-timepoint samples collected from 182 participants in the Women's Health Initiative cohort, across a median span of 16 years. We detected 178 driver mutations reaching variant allele frequency ≥ 2% in at least one timepoint, many of which were detectable well below this threshold at earlier timepoints. The majority of clonal mutations (52.1%) expanded over time (with a median doubling period of 7.43 years), with the others remaining static or decreasing in size in the absence of any cytotoxic therapy. CONCLUSIONS: Targeted smMIPS sequencing can sensitively measure clonal dynamics in CHIP. Mutations that reached the conventional threshold for CHIP (2% frequency) tended to continue growing, indicating that after CHIP is acquired, it is generally not lost. The ability to cost-effectively profile CHIP longitudinally will enable future studies to investigate why some CHIP clones expand, and how their dynamics relate to health outcomes at a biobank scale.

High-throughput splicing assays identify missense and silent splice-disruptive POU1F1 variants underlying pituitary hormone deficiency.

Pituitary hormone deficiency occurs in ∼1:4,000 live births. Approximately 3% of the cases are due to mutations in the alpha isoform of POU1F1, a pituitary-specific transcriptional activator. We found four separate heterozygous missense variants in unrelated individuals with hypopituitarism that were predicted to affect a minor isoform, POU1F1 beta, which can act as a transcriptional repressor. These variants retain repressor activity, but they shift splicing to favor the expression of the beta isoform, resulting in dominant-negative loss of function. Using a high-throughput splicing reporter assay, we tested 1,070 single-nucleotide variants in POU1F1. We identified 96 splice-disruptive variants, including 14 synonymous variants. In separate cohorts, we found two additional synonymous variants nominated by this screen that co-segregate with hypopituitarism. This study underlines the importance of evaluating the impact of variants on splicing and provides a catalog for interpretation of variants of unknown significance in POU1F1.

A Transcription Start Site Map in Human Pancreatic Islets Reveals Functional Regulatory Signatures.

Identifying the tissue-specific molecular signatures of active regulatory elements is critical to understand gene regulatory mechanisms. Here, we identify transcription start sites (TSS) using cap analysis of gene expression (CAGE) across 57 human pancreatic islet samples. We identify 9,954 reproducible CAGE tag clusters (TCs), ∼20% of which are islet specific and occur mostly distal to known gene TSS. We integrated islet CAGE data with histone modification and chromatin accessibility profiles to identify epigenomic signatures of transcription initiation. Using a massively parallel reporter assay, we validated the transcriptional enhancer activity for 2,279 of 3,378 (∼68%) tested islet CAGE elements (5% false discovery rate). TCs within accessible enhancers show higher enrichment to overlap type 2 diabetes genome-wide association study (GWAS) signals than existing islet annotations, which emphasizes the utility of mapping CAGE profiles in disease-relevant tissue. This work provides a high-resolution map of transcriptional initiation in human pancreatic islets with utility for dissecting active enhancers at GWAS loci.

Chromatin information content landscapes inform transcription factor and DNA interactions.

Interactions between transcription factors and chromatin are fundamental to genome organization and regulation and, ultimately, cell state. Here, we use information theory to measure signatures of organized chromatin resulting from transcription factor-chromatin interactions encoded in the patterns of the accessible genome, which we term chromatin information enrichment (CIE). We calculate CIE for hundreds of transcription factor motifs across human samples and identify two classes: low and high CIE. The 10-20% of common and tissue-specific high CIE transcription factor motifs, associate with higher protein-DNA residence time, including different binding site subclasses of the same transcription factor, increased nucleosome phasing, specific protein domains, and the genetic control of both chromatin accessibility and gene expression. These results show that variations in the information encoded in chromatin architecture reflect functional biological variation, with implications for cell state dynamics and memory.

Massively parallel functional testing of MSH2 missense variants conferring Lynch syndrome risk.

The lack of functional evidence for the majority of missense variants limits their clinical interpretability and poses a key barrier to the broad utility of carrier screening. In Lynch syndrome (LS), one of the most highly prevalent cancer syndromes, nearly 90% of clinically observed missense variants are deemed "variants of uncertain significance" (VUS). To systematically resolve their functional status, we performed a massively parallel screen in human cells to identify loss-of-function missense variants in the key DNA mismatch repair factor MSH2. The resulting functional effect map is substantially complete, covering 94% of the 17,746 possible variants, and is highly concordant (96%) with existing functional data and expert clinicians' interpretations. The large majority (89%) of missense variants were functionally neutral, perhaps unexpectedly in light of its evolutionary conservation. These data provide ready-to-use functional evidence to resolve the ∼1,300 extant missense VUSs in MSH2 and may facilitate the prospective classification of newly discovered variants in the clinic.

SOX10-regulated promoter use defines isoform-specific gene expression in Schwann cells.

BACKGROUND: Multicellular organisms adopt various strategies to tailor gene expression to cellular contexts including the employment of multiple promoters (and the associated transcription start sites (TSSs)) at a single locus that encodes distinct gene isoforms. Schwann cells-the myelinating cells of the peripheral nervous system (PNS)-exhibit a specialized gene expression profile directed by the transcription factor SOX10, which is essential for PNS myelination. SOX10 regulates promoter elements associated with unique TSSs and gene isoforms at several target loci, implicating SOX10-mediated, isoform-specific gene expression in Schwann cell function. Here, we report on genome-wide efforts to identify SOX10-regulated promoters and TSSs in Schwann cells to prioritize genes and isoforms for further study. RESULTS: We performed global TSS analyses and mined previously reported ChIP-seq datasets to assess the activity of SOX10-bound promoters in three models: (i) an adult mammalian nerve; (ii) differentiating primary Schwann cells, and (iii) cultured Schwann cells with ablated SOX10 function. We explored specific characteristics of SOX10-dependent TSSs, which provides confidence in defining them as SOX10 targets. Finally, we performed functional studies to validate our findings at four previously unreported SOX10 target loci: ARPC1A, CHN2, DDR1, and GAS7. These findings suggest roles for the associated SOX10-regulated gene products in PNS myelination. CONCLUSIONS: In sum, we provide comprehensive computational and functional assessments of SOX10-regulated TSS use in Schwann cells. The data presented in this study will stimulate functional studies on the specific mRNA and protein isoforms that SOX10 regulates, which will improve our understanding of myelination in the peripheral nerve.