PASSPORT-seq: A Novel High-Throughput Bioassay to Functionally Test Polymorphisms in Micro-RNA Target Sites

Massively parallel screen uncovers many rare 3' UTR variants regulating mRNA abundance of cancer driver genes

Understanding the function of rare non-coding variants represents a significant challenge. Using MapUTR, a screening method, we studied the function of rare 3' UTR variants affecting mRNA abundance post-transcriptionally. Among 17,301 rare gnomAD variants, an average of 24.5% were functional, with 70% in cancer-related genes, many in critical cancer pathways. This observation motivated an interrogation of 11,929 somatic mutations, uncovering 3928 (33%) functional mutations in 155 cancer driver genes. Functional MapUTR variants were enriched in microRNA- or protein-binding sites and may underlie outlier gene expression in tumors. Further, we introduce untranslated tumor mutational burden (uTMB), a metric reflecting the amount of somatic functional MapUTR variants of a tumor and show its potential in predicting patient survival. Through prime editing, we characterized three variants in cancer-relevant genes (MFN2, FOSL2, and IRAK1), demonstrating their cancer-driving potential. Our study elucidates the function of tens of thousands of non-coding variants, nominates non-coding cancer driver mutations, and demonstrates their potential contributions to cancer.

The Search for and Functional Analysis of Genetic Variants in microRNA-Binding Sites using Massively Parallel Reporter Assay

A large-scale search for the genetic variants with a bias in the representation of alleles in transcriptome data (AE SNPs) and the binding sites in microRNA 3'-UTRs was performed and their functional significance was assessed using massively parallel reporter assay (MPRA). Of the 629,559 associated "SNP-gene" pairs (eQTLs) discovered in the human liver tissue according to the GTEx Analysis V8 data, 4394 polymorphic positions in the 3'-UTRs of the genes, which represent the eQTLs for these genes were selected. The TargetScanHuman 7.0 algorithm and PolymiRTS database were searched for the potential microRNA-binding sites. Of the predicted microRNA sites affected by eQTL-SNPs, we selected 51 sites with the best evidence of functionality according to Ago2-CLIP-seq, CLEAR-CLIP, and eCLIP-seq for RNA-binding proteins. For MPRA, a library of the plasmids carrying the main and alternative alleles for each AE SNP (in total, 102 constructs) was created. Allele-specific expression for 6 SNPs was detected by transfection of the HepG2 cell line with the constructed plasmid library and sequencing of target DNA and RNA sequences using the Illumina (MiSeq) platform.

Functional 3'-UTR Variants Identify Regulatory Mechanisms Impacting Alcohol Use Disorder and Related Traits

Although genome-wide association studies (GWAS) have identified loci associated with alcohol consumption and alcohol use disorder (AUD), they do not identify which variants are functional. To approach this, we evaluated the impact of variants in 3' untranslated regions (3'-UTRs) of genes in loci associated with substance use and neurological disorders using a massively parallel reporter assay (MPRA) in neuroblastoma and microglia cells. Functionally impactful variants explained a higher proportion of heritability of alcohol traits than non-functional variants. We identified genes whose 3'UTR activities are associated with AUD and alcohol consumption by combining variant effects from MPRA with GWAS results. We examined their effects by evaluating gene expression after CRISPR inhibition of neuronal cells and stratifying brain tissue samples by MPRA-derived 3'-UTR activity. A pathway analysis of differentially expressed genes identified inflammation response pathways. These analyses suggest that variation in response to inflammation contributes to the propensity to increase alcohol consumption.

5. Collaborative Study on the Genetics of Alcoholism: Functional genomics

Alcohol Use Disorder is a complex genetic disorder, involving genetic, neural, and environmental factors, and their interactions. The Collaborative Study on the Genetics of Alcoholism (COGA) has been investigating these factors and identified putative alcohol use disorder risk genes through genome-wide association studies. In this review, we describe advances made by COGA in elucidating the functional changes induced by alcohol use disorder risk genes using multimodal approaches with human cell lines and brain tissue. These studies involve investigating gene regulation in lymphoblastoid cells from COGA participants and in post-mortem brain tissues. High throughput reporter assays are being used to identify single nucleotide polymorphisms in which alternate alleles differ in driving gene expression. Specific single nucleotide polymorphisms (both coding or noncoding) have been modeled using induced pluripotent stem cells derived from COGA participants to evaluate the effects of genetic variants on transcriptomics, neuronal excitability, synaptic physiology, and the response to ethanol in human neurons from individuals with and without alcohol use disorder. We provide a perspective on future studies, such as using polygenic risk scores and populations of induced pluripotent stem cell-derived neurons to identify signaling pathways related with responses to alcohol. Starting with genes or loci associated with alcohol use disorder, COGA has demonstrated that integration of multimodal data within COGA participants and functional studies can reveal mechanisms linking genomic variants with alcohol use disorder, and potential targets for future treatments.

SNPs in 3'UTR miRNA Target Sequences Associated with Individual Drug Susceptibility

The complementary interaction of microRNAs (miRNAs) with their binding sites in the 3'untranslated regions (3'UTRs) of target gene mRNAs represses translation, playing a leading role in gene expression control. MiRNA recognition elements (MREs) in the 3'UTRs of genes often contain single nucleotide polymorphisms (SNPs), which can change the binding affinity for target miRNAs leading to dysregulated gene expression. Accumulated data suggest that these SNPs can be associated with various human pathologies (cancer, diabetes, neuropsychiatric disorders, and cardiovascular diseases) by disturbing the interaction of miRNAs with their MREs located in mRNA 3'UTRs. Numerous data show the role of SNPs in 3'UTR MREs in individual drug susceptibility and drug resistance mechanisms. In this review, we brief the data on such SNPs focusing on the most rigorously proven cases. Some SNPs belong to conventional genes from the drug-metabolizing system (in particular, the genes coding for cytochromes P450 (CYP 450), phase II enzymes (SULT1A1 and UGT1A), and ABCB3 transporter and their expression regulators (PXR and GATA4)). Other examples of SNPs are related to the genes involved in DNA repair, RNA editing, and specific drug metabolisms. We discuss the gene-by-gene studies and genome-wide approaches utilized or potentially utilizable to detect the MRE SNPs associated with individual response to drugs.

Mapping the miRNA-mRNA Interactome in Human Hepatocytes and Identification of Functional mirSNPs in Pharmacogenes

MiRNAs regulate the expression of hepatic genes involved in pharmacokinetics and pharmacodynamics. Genetic variants affecting miRNA binding (mirSNPs) have been associated with altered drug response, but previously used methods to identify miRNA binding sites and functional mirSNPs in pharmacogenes are indirect and limited by low throughput. We utilized the high-throughput chimeric-eCLIP assay to directly map thousands of miRNA-mRNA interactions and define the miRNA binding sites in primary hepatocytes. We then used the high-throughput PASSPORT-seq assay to functionally test 262 potential mirSNPs with coordinates overlapping the identified miRNA binding sites. Using chimeric-eCLIP, we identified a network of 448 miRNAs that collectively target 11,263 unique genes in primary hepatocytes pooled from 100 donors. Our data provide an extensive map of miRNA binding of each gene, including pharmacogenes, expressed in primary hepatocytes. For example, we identified the hsa-mir-27b-DPYD interaction at a previously validated binding site. A second example is our identification of 19 unique miRNAs that bind to CYP2B6 across 20 putative binding sites on the transcript. Using PASSPORT-seq, we then identified 24 mirSNPs that functionally impacted reporter mRNA levels. To our knowledge, this is the most comprehensive identification of miRNA binding sites in pharmacogenes. Combining chimeric-eCLIP with PASSPORT-seq successfully identified functional mirSNPs in pharmacogenes that may affect transcript levels through altered miRNA binding. These results provide additional insights into potential mechanisms contributing to interindividual variability in drug response.

Allele-specific expression and high-throughput reporter assay reveal functional genetic variants associated with alcohol use disorders

Genome-wide association studies (GWAS) of complex traits, such as alcohol use disorders (AUD), usually identify variants in non-coding regions and cannot by themselves distinguish whether the associated variants are functional or in linkage disequilibrium with the functional variants. Transcriptome studies can identify genes whose expression differs between alcoholics and controls. To test which variants associated with AUD may cause expression differences, we integrated data from deep RNA-seq and GWAS of four postmortem brain regions from 30 subjects with AUD and 30 controls to analyze allele-specific expression (ASE). We identified 88 genes with differential ASE in subjects with AUD compared to controls. Next, to test one potential mechanism contributing to the differential ASE, we analyzed single nucleotide polymorphisms (SNPs) in the 3' untranslated regions (3'UTR) of these genes. Of the 88 genes with differential ASE, 61 genes contained 437 SNPs in the 3'UTR with at least one heterozygote among the subjects studied. Using a modified PASSPORT-seq (parallel assessment of polymorphisms in miRNA target-sites by sequencing) assay, we identified 25 SNPs that affected RNA levels in a consistent manner in two neuroblastoma cell lines, SH-SY5Y and SK-N-BE(2). Many of these SNPs are in binding sites of miRNAs and RNA-binding proteins, indicating that these SNPs are likely causal variants of AUD-associated differential ASE. In sum, we demonstrate that a combination of computational and experimental approaches provides a powerful strategy to uncover functionally relevant variants associated with the risk for AUD.

Identification of Functional Genetic Variants Associated With Alcohol Dependence and Related Phenotypes Using a High-Throughput Assay

BACKGROUND

Genome-wide association studies (GWAS) of alcohol dependence (AD) and related phenotypes have identified multiple loci, but the functional variants underlying the loci have in most cases not been identified. Noncoding variants can influence phenotype by affecting gene expression; for example, variants in the 3' untranslated regions (3'UTR) can affect gene expression posttranscriptionally.

METHODS

We adapted a high-throughput assay known as PASSPORT-seq (parallel assessment of polymorphisms in miRNA target sites by sequencing) to identify among variants associated with AD and related phenotypes those that cause differential expression in neuronal cell lines. Based upon meta-analyses of alcohol-related traits in African American and European Americans in the Collaborative Study on the Genetics of Alcoholism, we tested 296 single nucleotide polymorphisms (SNPs with meta-analysis p values ≤ 0.001) that were located in 3'UTRs.

RESULTS

We identified 60 SNPs that affected gene expression (false discovery rate [FDR] < 0.05) in SH-SY5Y cells and 92 that affected expression in SK-N-BE(2) cells. Among these, 30 SNPs altered RNA levels in the same direction in both cell lines. Many of these SNPs reside in the binding sites of miRNAs and RNA-binding proteins and are expression quantitative trait loci of genes including KIF6,FRMD4A,CADM2,ADD2,PLK2, and GAS7.

CONCLUSION

The SNPs identified in the PASSPORT-seq assay are functional variants that might affect the risk for AD and related phenotypes. Our study provides insights into gene regulation in AD and demonstrates the value of PASSPORT-seq as a tool to screen genetic variants in GWAS loci for one potential mechanism of action.

The Impact of Population Variation in the Analysis of microRNA Target Sites

The impact of population variation in the analysis of regulatory interactions is an underdeveloped area. MicroRNA target recognition occurs via pairwise complementarity. Consequently, a number of computational prediction tools have been developed to identify potential target sites that can be further validated experimentally. However, as microRNA target predictions are done mostly considering a reference genome sequence, target sites showing variation among populations are neglected. Here, we studied the variation at microRNA target sites in human populations and quantified their impact in microRNA target prediction. We found that African populations carry a significant number of potential microRNA target sites that are not detectable in the current human reference genome sequence. Some of these targets are conserved in primates and only lost in Out-of-Africa populations. Indeed, we identified experimentally validated microRNA/transcript interactions that are not detected in standard microRNA target prediction programs, yet they have segregating target alleles abundant in non-European populations. In conclusion, we show that ignoring population diversity may leave out regulatory elements essential to understand disease and gene expression, particularly neglecting populations of African origin.

Prediction of drug response and adverse drug reactions: From twin studies to Next Generation Sequencing

Understanding and predicting inter-individual differences related to the success of drug therapy is of tremendous importance, both during drug development and for clinical applications. Importantly, while seminal twin studies indicate that the majority of inter-individual differences in drug disposition are driven by hereditary factors, common genetic polymorphisms explain only less than half of this genetically encoded variability. Recent progress in Next Generation Sequencing (NGS) technologies has for the first time allowed to comprehensively map the genetic landscape of human pharmacogenes. Importantly, these projects have unveiled vast numbers of rare genetic variants, which are estimated to contribute substantially to the missing heritability of drug metabolism phenotypes. However, functional interpretation of these rare variants remains challenging and constitutes one of the important frontiers of contemporary pharmacogenomics. Furthermore, NGS technologies face challenges in the interrogation of genes residing in complex genomic regions, such as CYP2D6 and HLA genes. We here provide an update of the implementation of pharmacogenomic variations in the clinical setting and present emerging strategies that facilitate the translation of NGS data into clinically useful information. Importantly, we anticipate that these developments will soon result in a paradigm shift of pre-emptive genotyping away from the interrogation to candidate variants and towards the comprehensive profiling of an individuals genotype, thus allowing for a true individualization of patient drug treatment regimens.