Search

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.

BACKGROUND AND AIMS: Variation in micro-RNA (miRNA) levels in blood has been associated with alterations of physiological functions of the cardiovascular system. Circulating miRNA have the potential to become reliable biomarkers for risk stratification and early detection of cardiovascular events. Recurrent thrombotic events in patients with established coronary artery disease (CAD) demonstrate the need for personalized approaches to secondary prevention, especially in light of recent novel treatment approaches. METHODS: In a single center cohort study, whole blood samples were collected from 437 subjects undergoing cardiac catheterization, who were followed for recurrent cardiovascular events during a mean follow up of 1.5 years. We selected a case cohort (n = 22) with recurrent thrombotic events on standard medical therapy (stent thrombosis (n = 6) or spontaneous myocardial infarction (MI) (n = 16)) and a matched cohort with CAD, but uneventful clinical follow up (n = 26), as well as a control group with cardiovascular risk factors, but without angiographic CAD (n = 24). We performed complete miRNA next generation sequencing of RNA extracted from whole blood samples (including leukocytes and platelets). RESULTS: A differential pattern of miRNA expression was found among controls, CAD patients with no events, and CAD patients with recurrent events. MiRNA previously associated with MI, CAD, endothelial function, vascular smooth muscle cells, platelets, angiogenesis, heart failure, cardiac hypertrophy, arrhythmia, and stroke were found variably expressed in our case-control cohorts. Seventy miRNA (FDR < 0.05) were linked to the risk of recurrent myocardial infarction and future stent thrombosis, as compared to CAD patients with subsequently uneventful follow up. CONCLUSIONS: MiRNA next generation sequencing demonstrates altered fingerprint profile of whole blood miRNA expression among subjects with subsequent recurrent thrombotic events on standard medical therapy (‘non-responders’), as compared to subjects with no recurrent cardiovascular events. MiRNA profiling may be useful to identify high risk subjects and provide additional insights into disease mechanisms not currently attenuated with standard medical therapy used in CAD treatment.

siRNA stabilized for in vivo applications is filtered and reabsorbed in the renal proximal tubule (PT), reducing mRNA expression transiently. Prior siRNA efforts have successfully prevented upregulation of mRNA in response to injury. We proposed reducing constitutive gene and protein expression of LRP2 (megalin) in order to understand its molecular regulation in mice. Using siRNA targeting mouse LRP2 (siLRP2), reduction of LRP2 mRNA expression was compared to scrambled siRNA (siSCR) in mouse PT cells. Mice received siLRP2 administration optimized for dose, administration site, carrier solution, administration frequency, and administration duration. Kidney cortex was collected upon sacrifice. Renal gene and protein expression were compared by qRT-PCR, immunoblot, and immunohistochemistry (IHC). Compared to siSCR, siLRP2 reduced mRNA expression in PT cells to 16.6%±0.6%. In mouse kidney cortex, siLRP2 reduced mRNA expression to 74.8 ± 6.3% 3 h and 70.1 ± 6.3% 6 h after administration. mRNA expression rebounded at 12 h (160.6 ± 11.2%). No megalin renal protein expression reduction was observed by immunoblot or IHC, even after serial twice daily dosing for 3.5 days. Megalin is a constitutively expressed protein. Although LRP2 renal mRNA expression reduction was achieved, siRNA remains a costly and inefficient intervention to reduce in vivo megalin protein expression.

Next-generation sequencing (NGS) studies have identified large numbers of genetic variants that are predicted to alter miRNA–mRNA interactions. We developed a novel high-throughput bioassay, PASSPORT-seq, that can functionally test in parallel 100s of these variants in miRNA binding sites (mirSNPs). The results are highly reproducible across both technical and biological replicates. The utility of the bioassay was demonstrated by testing 100 mirSNPs in HEK293, HepG2, and HeLa cells. The results of several of the variants were validated in all three cell lines using traditional individual luciferase assays. Fifty-five mirSNPs were functional in at least one of three cell lines (FDR ≤ 0.05); 11, 36, and 27 of them were functional in HEK293, HepG2, and HeLa cells, respectively. Only four of the variants were functional in all three cell lines, which demonstrates the cell-type specific effects of mirSNPs and the importance of testing the mirSNPs in multiple cell lines. Using PASSPORT-seq, we functionally tested 111 variants in the 3′ UTR of 17 pharmacogenes that are predicted to alter miRNA regulation. Thirty-three of the variants tested were functional in at least one cell line.

CYP2B6*6 and CYP2B6*18 are the most clinically important variants causing reduced CYP2B6 protein expression and activity. However, these variants do not account for all variability in CYP2B6 activity. Emerging evidence has shown that genetic variants in the 3'UTR may explain variable drug response by altering microRNA regulation. Five 3'UTR variants were associated with significantly altered efavirenz AUC0-48 (8-OH-EFV/EFV) ratios in healthy human volunteers. The rs70950385 (AG>CA) variant, predicted to create a microRNA binding site for miR-1275, was associated with 33% decreased CYP2B6 activity among normal metabolizers (AG/AG vs. CA/CA [p