Search

Determining the phenotypic effects of single nucleotide variants is critical for understanding the genome and interpreting clinical sequencing results. Base editors, including diversifying base editors that create C>N mutations, are potent tools for installing point mutations in mammalian genomes and studying their effect on cellular function. Numerous base editor options are available for such studies, but little information exists on how the composition of the editor (deaminase, recruitment method, and fusion architecture) affects editing. To address this knowledge gap, the effect of various design features, such as deaminase recruitment and delivery method (electroporation or lentiviral transduction), on editing was assessed across ∼200 synthetic target sites. The direct fusion of a hyperactive variant of activation-induced cytidine deaminase to the N-terminus of dCas9 (DivA-BE) produced the highest editing efficiency, ∼4-fold better than the previous CRISPR-X method. Additionally, DivA-BE mutagenized the DNA strand that anneals to the targeting sgRNA to create G>N mutations, which were absent when the deaminase was fused to the C-terminus of dCas9. The DivA-BE editors efficiently diversified their target sites, an ideal characteristic for discovering functional variants. These and other findings provide a comprehensive analysis of how design features influence the activity of several popular base editors.

5'-end modifications play key roles in determining RNA fates. Phospho-methylation is a noncanonical cap occurring on either 5'-PPP or 5'-P ends. We used ChemRAP, in which affinity purification of cellular proteins with chemically synthesized modified RNAs is coupled to quantitative proteomics, to identify 5'-Pme "readers". We show that 5'-Pme is directly recognized by EPRS, the central subunit of the multisynthetase complex (MSC), through its linker domain, which has previously been involved in key noncanonical EPRS and MSC functions. We further determine that the 5'-Pme writer BCDIN3D regulates the binding of EPRS to specific mRNAs, either at coding regions rich in MSC codons, or around start codons. In the case of LRPPRC (leucine-rich pentatricopeptide repeat containing), a nuclear-encoded mitochondrial protein associated with the French Canadian Leigh syndrome, BCDIN3D deficiency abolishes binding of EPRS around its mRNA start codon, increases its translation but ultimately results in LRPPRC mislocalization. Overall, our results suggest that BCDIN3D may regulate the translation of specific mRNA via RNA-5'-Pme., (© 2024. The Author(s).)

Noncoding variants of presumed regulatory function contribute to the heritability of neuropsychiatric disease. A total of 2,221 noncoding variants connected to risk for ten neuropsychiatric disorders, including autism spectrum disorder, attention deficit hyperactivity disorder, bipolar disorder, borderline personality disorder, major depression, generalized anxiety disorder, panic disorder, post-traumatic stress disorder, obsessive-compulsive disorder and schizophrenia, were studied in developing human neural cells. Integrating epigenomic and transcriptomic data with massively parallel reporter assays identified differentially-active single-nucleotide variants (daSNVs) in specific neural cell types. Expression-gene mapping, network analyses and chromatin looping nominated candidate disease-relevant target genes modulated by these daSNVs. Follow-up integration of daSNV gene editing with clinical cohort analyses suggested that magnesium transport dysfunction may increase neuropsychiatric disease risk and indicated that common genetic pathomechanisms may mediate specific symptoms that are shared across multiple neuropsychiatric diseases., (© 2023. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Associations between genetic variation and traits are often in noncoding regions with strong linkage disequilibrium (LD), where a single causal variant is assumed to underlie the association. We applied a massively parallel reporter assay (MPRA) to functionally evaluate genetic variants in high, local LD for independent cis-expression quantitative trait loci (eQTL). We found that 17.7% of eQTLs exhibit more than one major allelic effect in tight LD. The detected regulatory variants were highly and specifically enriched for activating chromatin structures and allelic transcription factor binding. Integration of MPRA profiles with eQTL/complex trait colocalizations across 114 human traits and diseases identified causal variant sets demonstrating how genetic association signals can manifest through multiple, tightly linked causal variants.

Type II diabetes (T2D) and specific cancers share many risk factors, however, the molecular mechanisms underlying these connections are often not well-understood. BCDIN3D is an RNA modifying enzyme that methylates specific precursor microRNAs and tRNA His . In addition to breast cancer, BCDIN3D may also be linked to metabolism, as its gene locus is associated with obesity and T2D. In order to uncover metabolic pathways regulated by BCDIN3D in cancer, we performed an unbiased analysis of the metabolome, transcriptome, and proteome of breast cancer cells depleted for BCDIN3D. Intersection of these analyses showed that BCDIN3D-depleted cells have increased levels of Fructose 1,6 Bisphosphate (F1,6-BP), the last six-carbon glycolytic intermediate accompanied by reduced glycolytic capacity. We further show that elevated F1,6-BP is due to downregulation of Aldolase C (ALDOC), an enzyme that cleaves F1,6-BP mainly in the brain, but whose high expression/amplification is associated with poor prognosis in breast cancer. BCDIN3D regulates ALDOC through a non-canonical mechanism involving the crucial let-7 microRNA family and its target site on the 3'UTR of ALDOC. Overall, our results reveal an important connection between BCDIN3D, let-7 and glycolysis that may be relevant to breast cancer, obesity, and T2D.

Recombinant adeno-associated virus (rAAV) vectors have the unique ability to promote targeted integration of transgenes via homologous recombination at specified genomic sites, reaching frequencies of 0.1%-1%. We studied genomic parameters that influence targeting efficiencies on a large scale. To do this, we generated more than 1,000 engineered, doxycycline-inducible target sites in the human HAP1 cell line and infected this polyclonal population with a library of AAV-DJ targeting vectors, with each carrying a unique barcode. The heterogeneity of barcode integration at each target site provided an assessment of targeting efficiency at that locus. We compared targeting efficiency with and without target site transcription for identical chromosomal positions. Targeting efficiency was enhanced by target site transcription, while chromatin accessibility was associated with an increased likelihood of targeting. ChromHMM chromatin states characterizing transcription and enhancers in wild-type K562 cells were also associated with increased AAV-HR efficiency with and without target site transcription, respectively. Furthermore, the amenability of a site to targeting was influenced by the endogenous transcriptional level of intersecting genes. These results define important parameters that may not only assist in designing optimal targeting vectors for genome editing, but also provide new insights into the mechanism of AAV-mediated homologous recombination., (Copyright © 2020 The American Society of Gene and Cell Therapy. Published by Elsevier Inc. All rights reserved.)

Type I interferons (IFNs) induce hundreds of IFN-stimulated genes (ISGs) in response to viral infection. Induction of these ISGs must be regulated for an efficient and controlled antiviral response, but post-transcriptional controls of these genes have not been well defined. Here, we identify a role for the RNA base modification N6-methyladenosine (m 6 A) in the regulation of ISGs. Using ribosome profiling and quantitative mass spectrometry, coupled with m 6 A-immunoprecipitation and sequencing, we identify a subset of ISGs, including IFITM1, whose translation is enhanced by m 6 A and the m 6 A methyltransferase proteins METTL3 and METTL14. We further determine that the m 6 A reader YTHDF1 increases the expression of IFITM1 in an m 6 A-binding-dependent manner. Importantly, we find that the m 6 A methyltransferase complex promotes the antiviral activity of type I IFN. Thus, these studies identify m 6 A as having a role in post-transcriptional control of ISG translation during the type I IFN response for antiviral restriction., Competing Interests: Declaration of interests C.E.M. is a cofounder and board member for Biotia and Onegevity Health and an advisor or compensated speaker for Abbvie, Acuamark Diagnostics, ArcBio, Bio-Rad, DNA Genotek, Genialis, Genpro, Illumina, New England Biolabs, QIAGEN, Whole Biome, and Zymo Research., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

Rare genetic variants are abundant across the human genome, and identifying their function and phenotypic impact is a major challenge. Measuring aberrant gene expression has aided in identifying functional, large-effect rare variants (RVs). Here, we expanded detection of genetically driven transcriptome abnormalities by analyzing gene expression, allele-specific expression, and alternative splicing from multitissue RNA-sequencing data, and demonstrate that each signal informs unique classes of RVs. We developed Watershed, a probabilistic model that integrates multiple genomic and transcriptomic signals to predict variant function, validated these predictions in additional cohorts and through experimental assays, and used them to assess RVs in the UK Biobank, the Million Veterans Program, and the Jackson Heart Study. Our results link thousands of RVs to diverse molecular effects and provide evidence to associate RVs affecting the transcriptome with human traits., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Background: Population structure among study subjects may confound genetic association studies, and lack of proper correction can lead to spurious findings. The Genotype-Tissue Expression (GTEx) project largely contains individuals of European ancestry, but the v8 release also includes up to 15% of individuals of non-European ancestry. Assessing ancestry-based adjustments in GTEx improves portability of this research across populations and further characterizes the impact of population structure on GWAS colocalization., Results: Here, we identify a subset of 117 individuals in GTEx (v8) with a high degree of population admixture and estimate genome-wide local ancestry. We perform genome-wide cis-eQTL mapping using admixed samples in seven tissues, adjusted by either global or local ancestry. Consistent with previous work, we observe improved power with local ancestry adjustment. At loci where the two adjustments produce different lead variants, we observe 31 loci (0.02%) where a significant colocalization is called only with one eQTL ancestry adjustment method. Notably, both adjustments produce similar numbers of significant colocalizations within each of two different colocalization methods, COLOC and FINEMAP. Finally, we identify a small subset of eQTL-associated variants highly correlated with local ancestry, providing a resource to enhance functional follow-up., Conclusions: We provide a local ancestry map for admixed individuals in the GTEx v8 release and describe the impact of ancestry and admixture on gene expression, eQTLs, and GWAS colocalization. While the majority of the results are concordant between local and global ancestry-based adjustments, we identify distinct advantages and disadvantages to each approach.

5' ends are important for determining the fate of RNA molecules. BCDIN3D is an RNA phospho-methyltransferase that methylates the 5' monophosphate of specific RNAs. In order to gain new insights into the molecular function of BCDIN3D, we performed an unbiased analysis of its interacting RNAs by Thermostable Group II Intron Reverse Transcriptase coupled to next generation sequencing (TGIRT-seq). Our analyses showed that BCDIN3D interacts with full-length phospho-methylated tRNAHis and miR-4454. Interestingly, we found that miR-4454 is not synthesized from its annotated genomic locus, which is a primer-binding site for an endogenous retrovirus, but rather by Dicer cleavage of mature tRNAHis. Sequence analysis revealed that miR-4454 is identical to the 3' end of tRNAHis. Moreover, we were able to generate this 'miRNA' in vitro through incubation of mature tRNAHis with Dicer. As found previously for several pre-miRNAs, a 5'P-tRNAHis appears to be a better substrate for Dicer cleavage than a phospho-methylated tRNAHis. Moreover, tRNAHis 3'-fragment/'miR-4454' levels increase in cells depleted for BCDIN3D. Altogether, our results show that in addition to microRNAs, BCDIN3D regulates tRNAHis 3'-fragment processing without negatively affecting tRNAHis's canonical function of aminoacylation., Competing Interests: I have read the journal's policy and the authors of this manuscript have the following competing interests: Thermostable Group II Intron Reverse Transcriptase (TGIRT) enzymes and methods for their use are the subject of patents and patent applications that have been licensed by the University of Texas and East Tennessee State University to InGex, LLC. AML and the University of Texas are minority equity holders in InGex, LLC, and AML and the University of Texas receive royalty payments from the sale of TGIRT enzymes and the licensing of intellectual property by InGex to other companies. The other authors declare no competing financial interests.

The retinal pigment epithelium (RPE) serves vital roles in ocular development and retinal homeostasis but has limited representation in large-scale functional genomics datasets. Understanding how common human genetic variants affect RPE gene expression could elucidate the sources of phenotypic variability in selected monogenic ocular diseases and pinpoint causal genes at genome-wide association study (GWAS) loci. We interrogated the genetics of gene expression of cultured human fetal RPE (fRPE) cells under two metabolic conditions and discovered hundreds of shared or condition-specific expression or splice quantitative trait loci (e/sQTLs). Co-localizations of fRPE e/sQTLs with age-related macular degeneration (AMD) and myopia GWAS data suggest new candidate genes, and mechanisms by which a common RDH5 allele contributes to both increased AMD risk and decreased myopia risk. Our study highlights the unique transcriptomic characteristics of fRPE and provides a resource to connect e/sQTLs in a critical ocular cell type to monogenic and complex eye disorders., Competing Interests: Competing interestsStephen Montgomery is on the Scientific Advisory Board of Prime Genomics. The remaining authors declare no competing interests.

RNAP II switching from the paused to the productive transcription elongation state is a pivotal regulatory step that requires specific phosphorylations catalyzed by the P-TEFb kinase. Nucleosolic P-TEFb activity is inhibited by its interaction with the ribonuclear protein complex built around the 7SK small nuclear RNA (7SK snRNP). MePCE is the RNA methyltransferase that methylates and stabilizes 7SK in the nucleosol. Here, we report that MePCE also binds chromatin through the histone H4 tail to serve as a P-TEFb activator at specific genes important for cellular identity. Notably, this histone binding abolishes MePCE's RNA methyltransferase activity toward 7SK, which explains why MePCE-bound P-TEFb on chromatin may not be associated with the full 7SK snRNP and is competent for RNAP II activation. Overall, our results suggest that crosstalk between the histone-binding and RNA methylation activities of MePCE regulates P-TEFb activation on chromatin in a 7SK- and Brd4-independent manner., (Copyright © 2018 The Author(s). Published by Elsevier Inc. All rights reserved.)

Enoxacin is a small molecule that stimulates RNA interference (RNAi) and acts as a growth inhibitor selectively in cancer but not in untransformed cells. Here, we used alkenox, a clickable enoxacin surrogate, coupled with quantitative mass spectrometry, to identify PIWIL3 as a mechanistic target of enoxacin. PIWIL3 is an Argonaute protein of the PIWI subfamily that is mainly expressed in the germline and that mediates RNAi through piRNAs. Our results suggest that cancer cells re-express PIWIL3 to repress RNAi through miRNAs and thus open a new opportunity for cancer-specific targeting.

Exosomes are small extracellular vesicles that carry heterogeneous cargo, including RNA, between cells. Increasing evidence suggests that exosomes are important mediators of intercellular communication and biomarkers of disease. Despite this, the variability of exosomal RNA between individuals has not been well quantified. To assess this variability, we sequenced the small RNA of cells and exosomes from a 17-member family. Across individuals, we show that selective export of miRNAs occurs not only at the level of specific transcripts, but that a cluster of 74 mature miRNAs on chromosome 14q32 is massively exported in exosomes while mostly absent from cells. We also observe more interindividual variability between exosomal samples than between cellular ones and identify four miRNA expression quantitative trait loci shared between cells and exosomes. Our findings indicate that genomically colocated miRNAs can be exported together and highlight the variability in exosomal miRNA levels between individuals as relevant for exosome use as diagnostics., (Copyright © 2017 Tsang et al.)

miRNAs play a central role in numerous pathologies including multiple cancer types. miR-191 has predominantly been studied as an oncogene, but the role of miR-191 in the proliferation of primary cells is not well characterized, and the miR-191 targetome has not been experimentally profiled. Here we utilized RNA induced silencing complex immunoprecipitations as well as gene expression profiling to construct a genome wide miR-191 target profile. We show that miR-191 represses proliferation in primary human fibroblasts, identify multiple proto-oncogenes as novel miR-191 targets, including CDK9, NOTCH2, and RPS6KA3, and present evidence that miR-191 extensively mediates target expression through coding sequence (CDS) pairing. Our results provide a comprehensive genome wide miR-191 target profile, and demonstrate miR-191's regulation of primary human fibroblast proliferation.

Background: The pathways regulating the transition of mammalian cells from quiescence to proliferation are mediated by multiple miRNAs. Despite significant improvements in our understanding of miRNA targeting, the majority of miRNA regulatory networks are still largely unknown and require experimental validation., Results: Here we identified miR-503, miR-103, and miR-494 as negative regulators of proliferation in primary human cells. We experimentally determined their genome wide target profiles using RNA-induced silencing complex (RISC) immunoprecipitations and gene expression profiling. Analysis of the genome wide target profiles revealed evidence of extensive regulation of gene expression through non-canonical target pairing by miR-503. We identified the proto-oncogene DDHD2 as a target of miR-503 that requires pairing outside of the canonical 5' seed region of miR-503, representing a novel mode of miRNA-target pairing. Further bioinformatics analysis implicated miR-503 and DDHD2 in breast cancer tumorigenesis., Conclusions: Our results provide an extensive genome wide set of targets for miR-503, miR-103, and miR-494, and suggest that miR-503 may act as a tumor suppressor in breast cancer by its direct non-canonical targeting of DDHD2.

The transition of mammalian cells from quiescence to proliferation is accompanied by the differential expression of several microRNAs (miRNAs) and transcription factors. However, the interplay between transcription factors and miRNAs in modulating gene regulatory networks involved in human cell proliferation is largely unknown. Here we show that the miRNA miR-22 promotes proliferation in primary human cells, and through a combination of Argonaute-2 immunoprecipitation and reporter assays, we identified multiple novel targets of miR-22, including several cell-cycle arrest genes that mediate the effects of the tumor-suppressor p53. In addition, we found that miR-22 suppresses interferon gene expression by directly targeting high mobility group box-1 and interferon regulatory factor (IRF)-5, preventing activation of IRF3 and NF-κB, which are activators of interferon genes. The expression of interferon genes is elevated in quiescent cells and their expression is inhibitory for cell proliferation. In addition, we find that miR-22 is activated by the transcription factor Myc when quiescent cells enter proliferation and that miR-22 inhibits the Myc transcriptional repressor MXD4, mediating a feed-forward loop to elevate Myc expression levels. Our results implicate miR-22 in downregulating the anti-proliferative p53 and interferon pathways and reveal a new transcription factor-miRNA network that regulates the transition of primary human cells from quiescence to proliferation.