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) 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 mRNA 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 SINEs 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.

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.

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.

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.