Population-scale tissue transcriptomics maps long non-coding RNAs to complex disease.

Long non-coding RNA (lncRNA) genes have well-established and important impacts on molecular and cellular functions. However, among the thousands of lncRNA genes, it is still a major challenge to identify the subset with disease or trait relevance. To systematically characterize these lncRNA genes, we used Genotype Tissue Expression (GTEx) project v8 genetic and multi-tissue transcriptomic data to profile the expression, genetic regulation, cellular contexts, and trait associations of 14,100 lncRNA genes across 49 tissues for 101 distinct complex genetic traits. Using these approaches, we identified 1,432 lncRNA gene-trait associations, 800 of which were not explained by stronger effects of neighboring protein-coding genes. This included associations between lncRNA quantitative trait loci and inflammatory bowel disease, type 1 and type 2 diabetes, and coronary artery disease, as well as rare variant associations to body mass index.

A catalog of GWAS fine-mapping efforts in autoimmune disease.

Genome-wide association studies (GWASs) have enabled unbiased identification of genetic loci contributing to common complex diseases. Because GWAS loci often harbor many variants and genes, it remains a major challenge to move from GWASs' statistical associations to the identification of causal variants and genes that underlie these association signals. Researchers have applied many statistical and functional fine-mapping strategies to prioritize genetic variants and genes as potential candidates. There is no gold standard in fine-mapping approaches, but consistent results across different approaches can improve confidence in the fine-mapping findings. Here, we combined text mining with a systematic review and formed a catalog of 85 studies with evidence of fine mapping for at least one autoimmune GWAS locus. Across all fine-mapping studies, we compiled 230 GWAS loci with allelic heterogeneity estimates and predictions of causal variants and trait-relevant genes. These 230 loci included 455 combinations of locus-by-disease association signals with 15 autoimmune diseases. Using these estimates, we assessed the probability of mediating disease risk associations across genes in GWAS loci and identified robust signals of causal disease biology. We predict that this comprehensive catalog of GWAS fine-mapping efforts in autoimmune disease will greatly help distill the plethora of information in the field and inform therapeutic strategies.

Determination of sensed and paced atrial-ventricular delay in cardiac resynchronization therapy.

BACKGROUND: Optimization of atrial-ventricular delay (AVD) during atrial sensing (SAVD) and pacing (PAVD) provides the most effective cardiac resynchronization therapy (CRT). We demonstrate a novel electrocardiographic methodology for quantifying electrical synchrony and optimizing SAVD/PAVD. METHODS: We studied 40 CRT patients with LV activation delay. Atrial-sensed to RV-sensed (As-RVs) and atrial-paced to RV-sensed (Ap-RVs) intervals were measured from intracardiac electrograms (IEGM). LV-only pacing was performed over a range of SAVD/PAVD settings. Electrical dyssynchrony (cardiac resynchronization index; CRI) was measured at each setting using a multilead ECG system placed over the anterior and posterior torso. Biventricular pacing, which included multiple interventricular delays, was also conducted in a subset of 10 patients. RESULTS: When paced LV-only, peak CRI was similar (93 ± 5% vs. 92 ± 5%) during atrial sensing or pacing but optimal PAVD was 61 ± 31 ms greater than optimal SAVD. The difference between As-RVs and Ap-RVs intervals on IEGMs (62 ± 31 ms) was nearly identical. The slope of the correlation line (0.98) and the correlation coefficient r (0.99) comparing the 2 methods of assessing SAVD-PAVD offset were nearly 1 and the y-intercept (0.63 ms) was near 0. During simultaneous biventricular (BiV) pacing at short AVD, SAVD and PAVD programming did not affect CRI, but CRI was significantly (p < .05) lower during atrial sensing at long AVD. CONCLUSIONS: A novel methodology for measuring electrical dyssynchrony was used to determine electrically optimal SAVD/PAVD during LV-only pacing. When BiV pacing, shorter AVDs produce better electrical synchrony.

Effectiveness of Platelet-Rich Plasma in Reducing Pain and Increasing Function After Acute Lateral Ankle Sprain: A Critically Appraised Topic.

CLINICAL SCENARIO: Ankle sprains are one of the most common injuries in athletics, and many lead to recurrent sprains, chronic ankle instability, and persistent symptoms. Treatment improvements are needed. Platelet-rich plasma (PRP) involves formulating autologous plasma with higher platelet concentration to be injected in the desired tissue. There is currently high-quality evidence supporting the use of PRP with lateral epicondylitis and knee osteoarthritis to accelerate the healing process and decrease pain. CLINICAL QUESTION: Does the injection of PRP relieve pain faster and improve function compared with no injection or placebo in patients with a lateral ankle sprain? SUMMARY OF KEY FINDINGS: A computerized search yielded 191 studies; of these, 3 studies fit the inclusion and exclusion criteria. PRP injection reduces pain and increases function after lateral ankle sprain 5 to 8 weeks after intervention. CLINICAL BOTTOM LINE: The use of PRP after lateral ankle sprain to decrease pain and increase function is supported with moderate evidence. STRENGTH OF RECOMMENDATION: Based on the Strength of Recommendation Taxonomy, evidence from the included studies is considered as level B, reflecting limited quality patient-oriented evidence.

A Nucleophilic Activity-Based Probe Enables Profiling of PLP-Dependent Enzymes.

PLP-dependent enzymes represent an important class of highly "druggable" enzymes that perform a wide array of critical reactions to support all organisms. Inhibition of individual members of this family of enzymes has been validated as a therapeutic target for pathologies ranging from infection with Mycobacterium tuberculosis to epilepsy. Given the broad nature of the activities within this family of enzymes, we envisioned a universally acting probe to characterize existing and putative members of the family that also includes the necessary chemical moieties to enable activity-based protein profiling experiments. Hence, we developed a probe that contains an N-hydroxyalanine warhead that acts as a covalent inhibitor of PLP-dependent enzymes, a linear diazirine for UV crosslinking, and an alkyne moiety to enable enrichment of crosslinked proteins. Our molecule was used to study PLP-dependent enzymes in vitro as well as look at whole-cell lysates of M. tuberculosis and assess inhibitory activity. The probe was able to enrich and identify LysA, a PLP-dependent enzyme crucial for lysine biosynthesis, through mass spectrometry. Overall, our study shows the utility of this trifunctional first-generation probe. We anticipate further optimization of probes for PLP-dependent enzymes will enable the characterization of rationally designed covalent inhibitors of PLP-dependent enzymes, which will expedite the preclinical characterization of these important therapeutic targets.

Uncovering Active Bacterial Symbionts in Three Species of Pollen-feeding Beetles (Nitidulidae: Meligethinae).

Microbial symbionts enable many phytophagous insects to specialize on plant-based diets through a range of metabolic services. Pollen comprises one-plant tissue consumed by such herbivores. While rich in lipids and proteins, its nutrient content is often imbalanced and difficult-to-access due to a digestibly recalcitrant cell wall. Pollen quality can be further degraded by harmful allelochemicals. To identify microbes that may aid in palynivory, we performed cDNA-based 16S rRNA metabarcoding on three related pollen beetles (Nitidulidae: Meligethinae) exhibiting different dietary breadths: Brassicogethes aeneus, B. matronalis, and Meligethes atratus. Nine bacterial symbionts (i.e., 97% OTUs) exhibited high metabolic activity during active feeding. Subsequent PCR surveys revealed varying prevalence of those from three Rickettsialles genera-Lariskella, Rickettsia, and Wolbachia-within beetle populations. Our findings lay the groundwork for future studies on the influence of phylogeny and diet on palynivorous insect microbiomes, and roles of symbionts in the use of challenging diets.

Genetic effects on gene expression across human tissues.

Characterization of the molecular function of the human genome and its variation across individuals is essential for identifying the cellular mechanisms that underlie human genetic traits and diseases. The Genotype-Tissue Expression (GTEx) project aims to characterize variation in gene expression levels across individuals and diverse tissues of the human body, many of which are not easily accessible. Here we describe genetic effects on gene expression levels across 44 human tissues. We find that local genetic variation affects gene expression levels for the majority of genes, and we further identify inter-chromosomal genetic effects for 93 genes and 112 loci. On the basis of the identified genetic effects, we characterize patterns of tissue specificity, compare local and distal effects, and evaluate the functional properties of the genetic effects. We also demonstrate that multi-tissue, multi-individual data can be used to identify genes and pathways affected by human disease-associated variation, enabling a mechanistic interpretation of gene regulation and the genetic basis of disease.

Genomic profiling of human vascular cells identifies TWIST1 as a causal gene for common vascular diseases.

Genome-wide association studies have identified multiple novel genomic loci associated with vascular diseases. Many of these loci are common non-coding variants that affect the expression of disease-relevant genes within coronary vascular cells. To identify such genes on a genome-wide level, we performed deep transcriptomic analysis of genotyped primary human coronary artery smooth muscle cells (HCASMCs) and coronary endothelial cells (HCAECs) from the same subjects, including splicing Quantitative Trait Loci (sQTL), allele-specific expression (ASE), and colocalization analyses. We identified sQTLs for TARS2, YAP1, CFDP1, and STAT6 in HCASMCs and HCAECs, and 233 ASE genes, a subset of which are also GTEx eGenes in arterial tissues. Colocalization of GWAS association signals for coronary artery disease (CAD), migraine, stroke and abdominal aortic aneurysm with GTEx eGenes in aorta, coronary artery and tibial artery discovered novel candidate risk genes for these diseases. At the CAD and stroke locus tagged by rs2107595 we demonstrate colocalization with expression of the proximal gene TWIST1. We show that disrupting the rs2107595 locus alters TWIST1 expression and that the risk allele has increased binding of the NOTCH signaling protein RBPJ. Finally, we provide data that TWIST1 expression influences vascular SMC phenotypes, including proliferation and calcification, as a potential mechanism supporting a role for TWIST1 in CAD.

Genetic and Epigenetic Fine Mapping of Complex Trait Associated Loci in the Human Liver.

Deciphering the impact of genetic variation on gene regulation is fundamental to understanding common, complex human diseases. Although histone modifications are important markers of gene regulatory elements of the genome, any specific histone modification has not been assayed in more than a few individuals in the human liver. As a result, the effects of genetic variation on histone modification states in the liver are poorly understood. Here, we generate the most comprehensive genome-wide dataset of two epigenetic marks, H3K4me3 and H3K27ac, and annotate thousands of putative regulatory elements in the human liver. We integrate these findings with genome-wide gene expression data collected from the same human liver tissues and high-resolution promoter-focused chromatin interaction maps collected from human liver-derived HepG2 cells. We demonstrate widespread functional consequences of natural genetic variation on putative regulatory element activity and gene expression levels. Leveraging these extensive datasets, we fine-map a total of 74 GWAS loci that have been associated with at least one complex phenotype. Our results reveal a repertoire of genes and regulatory mechanisms governing complex disease development and further the basic understanding of genetic and epigenetic regulation of gene expression in the human liver tissue.

Genomic architecture of Shh-dependent cochlear morphogenesis.

The mammalian cochlea develops from a ventral outgrowth of the otic vesicle in response to Shh signaling. Mouse embryos lacking Shh or its essential signal transduction components display cochlear agenesis; however, a detailed understanding of the transcriptional network mediating this process is unclear. Here, we describe an integrated genomic approach to identify Shh-dependent genes and associated regulatory sequences that promote cochlear duct morphogenesis. A comparative transcriptome analysis of otic vesicles from mouse mutants exhibiting loss (Smoecko ) and gain (Shh-P1) of Shh signaling reveal a set of Shh-responsive genes partitioned into four expression categories in the ventral half of the otic vesicle. This target gene classification scheme provides novel insight into several unanticipated roles for Shh, including priming the cochlear epithelium for subsequent sensory development. We also mapped regions of open chromatin in the inner ear by ATAC-seq that, in combination with Gli2 ChIP-seq, identified inner ear enhancers in the vicinity of Shh-responsive genes. These datasets are useful entry points for deciphering Shh-dependent regulatory mechanisms involved in cochlear duct morphogenesis and establishment of its constituent cell types.

Functional Characterization of Organoids Derived From Irreversibly Damaged Liver of Patients With NASH.

BACKGROUND AND AIMS: NASH will soon become the leading cause of liver transplantation in the United States and is also associated with increased COVID-19 mortality. Currently, there are no Food and Drug Administration-approved drugs available that slow NASH progression or address NASH liver involvement in COVID-19. Because animal models cannot fully recapitulate human NASH, we hypothesized that stem cells isolated directly from end-stage liver from patients with NASH may address current knowledge gaps in human NASH pathology. APPROACH AND RESULTS: We devised methods that allow the derivation, proliferation, hepatic differentiation, and extensive characterization of bipotent ductal organoids from irreversibly damaged liver from patients with NASH. The transcriptomes of organoids derived from NASH liver, but not healthy liver, show significant up-regulation of proinflammatory and cytochrome p450-related pathways, as well as of known liver fibrosis and tumor markers, with the degree of up-regulation being patient-specific. Functionally, NASH liver organoids exhibit reduced passaging/growth capacity and hallmarks of NASH liver, including decreased albumin production, increased free fatty acid-induced lipid accumulation, increased sensitivity to apoptotic stimuli, and increased cytochrome P450 metabolism. After hepatic differentiation, NASH liver organoids exhibit reduced ability to dedifferentiate back to the biliary state, consistent with the known reduced regenerative ability of NASH livers. Intriguingly, NASH liver organoids also show strongly increased permissiveness to severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2) vesicular stomatitis pseudovirus as well as up-regulation of ubiquitin D, a known inhibitor of the antiviral interferon host response. CONCLUSION: Expansion of primary liver stem cells/organoids derived directly from irreversibly damaged liver from patients with NASH opens up experimental avenues for personalized disease modeling and drug development that has the potential to slow human NASH progression and to counteract NASH-related SARS-CoV-2 effects.

Relationship between QRS duration and resynchronization window for CRT optimization: Implications for CRT in narrow QRS patients.

AIMS: Cardiac resynchronization therapy (CRT) response is proportional to QRS duration (QRSd). We hypothesize that this is, in part, due to slower conduction velocity and hence wider range of programmed device settings that produce adequate electrical wavefront fusion and resynchronization in wider QRSd patients. METHODS: CRT patients (n = 122) with left ventricular (LV) conduction delay, sinus rhythm and intact atrioventricular node conduction were studied. Patients were categorized by QRSd: narrow (<120 ms; n = 20); moderate (120-150 ms, n = 37); and prolonged (≥150 ms; n = 65). Electrocardiographic data was acquired during native rhythm and LV-only pacing at varying atrioventricular delays (AVDs). Electrical synchrony was quantified as cardiac resynchronization index (CRI) using multi­lead electrocardiographic systems and a proprietary algorithm that quantified wavefront fusion. A Gaussian distribution equation was fitted to CRI response. RESULTS: Peak CRI was high (87.6 ± 6.3%) and similar (p = 0.716) across QRSd groups. The standard deviation of the Gaussian distribution significantly correlated with QRSd (R = 0.614, p < 0.001), and progressively and significantly (p < 0.001) increased as QRSd increased from narrow (34.8 ± 10.0 ms), to moderate (50.6 ± 8.4 ms), to prolonged (67.6 ± 18.3 ms). At AVDs 20 and 40 ms from optimal, CRI differed significantly (p < 0.001) between groups, with progressively higher CRI values as native QRSd increased. CONCLUSION: Electrical resynchronization with optimally programmed LV-only pacing was similar between patients with varying QRSd, including patients with narrow QRSd. The resynchronization window that corresponded with optimal electrical resynchronization decreased as native QRSd decreased. This finding provides one potential explanation for the lack of significant benefit of CRT in narrow QRSd patients in previous studies.

Electrical dyssynchrony mapping and cardiac resynchronization therapy.

PURPOSE: There is no clinical methodology for quantification or display of electrical dyssynchrony over a wide range of atrial-ventricular delays (AVD) and ventricular-ventricular delays (VVD) in patients with cardiac resynchronization therapy (CRT). This study aimed to develop a new methodology, based on wavefront fusion, for mapping electrical synchrony. METHODS: A cardiac resynchronization index (CRI) was measured at multiple device settings in 90 patients. Electrical dyssynchrony maps (EDM) were constructed for each patient to display CRI at any combination of AVD and VVD. An optimal synchrony line (OSL) depicted the AVD/VVD combinations producing the highest CRIs. Fusion of right ventricular paced (RVp), left ventricular paced (LVp), and native wavefront offsets were calculated. RESULTS: CRI significantly increased (p < 0.0001) from 58.0 ± 28.1% at baseline to 98.3 ± 1.7% at optimized settings. EDMs in patients with high-grade heart block (n = 20) had an OSL parallel to the simultaneous biventricular pacing (BiVPVV-SIM) line with leftward shift across all AVDs (RVp-LVpOFFSET = 50.5 ± 29.8 ms). EDMs in patients with intact AV node conduction (n = 64) had an OSL parallel to the BiVPVV-SIM line with leftward shift at short AVDs (RVp-LVpOFFSET = 33.4 ± 23.3 ms), curvilinear at intermediate AVDs (triple fusion), and vertical at long AVDs (native-LVpOFFSET = 85.2 ± 22.8 ms) in all patients except those with poor LV lead position (n = 6). CONCLUSION: A new methodology is described for quantifying and graphing electrical dyssynchrony over a physiologic range of AVDs/VVDs. This methodology offers a noninvasive, practical, clinical approach for measuring electrical synchrony that could be applied to optimization of CRT devices.

Identifying differential regulatory control of APOE ɛ4 on African versus European haplotypes as potential therapeutic targets.

We previously demonstrated that in Alzheimer's disease (AD) patients, European apolipoprotein E (APOE) ε4 carriers express significantly more APOE ε4 in their brains than African AD carriers. We examined single nucleotide polymorphisms near APOE with significant frequency differences between African and European/Japanese APOE ε4 haplotypes that could contribute to this difference in expression through regulation. Two enhancer massively parallel reporter assay (MPRA) approaches were performed, supplemented with single fragment reporter assays. We used Capture C analyses to support interactions with the APOE promoter. Introns within TOMM40 showed increased enhancer activity in the European/Japanese versus African haplotypes in astrocytes and microglia. This region overlaps with APOE promoter interactions as assessed by Capture C analysis. Single variant analyses pinpoints rs2075650/rs157581, and rs59007384 as functionally different on these haplotypes. Identification of the mechanisms for differential regulatory function for APOE expression between African and European/Japanese haplotypes could lead to therapeutic targets for APOE ε4 carriers.

The flow of sickle blood in glass capillaries: Fundamentals and potential applications.

We have characterized the imbibed horizontal flow of sickle blood into 100-µm-diameter glass capillaries. We find that blood containing sickled cells typically traverses the capillaries between three and four times as slowly as oxygenated cells from the same patient for all genotypes tested, including SS, AS, SC and Sß+ thalassemia blood. Blood from SS patients treated with hydroxyurea has a viscosity intermediate between the SS and AA values. Blood containing cells that are not rigidified, such as normal red cells or oxygenated sickle cells, follows a simple Lucas-Washburn flow throughout the length of the 3-cm capillary. By fitting the flexible-cell data to the Lucas-Washburn model, a viscosity can be derived that is in good agreement with previous measurements over a range of volume fractions and is obtained using an apparatus that is far more complex. Deoxygenation sickles and thus rigidifies the cells, and their flow begins as Lucas-Washburn, albeit with higher viscosity than flexible cells. However, the flow further slows as a dense mass of cells forms behind the meniscus and increases in length as flow progresses. By assuming that the dense mass of cells exerts a frictional force proportional to its length, we derive an equation that is formally equivalent to vertical imbibition, even though the flow is horizontal, and this equation reproduces the observed behavior well. We present a simple theory using activity coefficients that accounts for this viscosity and its variation without adjustable parameters. In the course of control experiments, we have found that deoxygenation increases the flexibility of normal human red cells, an observation only recently published for mouse cells and previously unreported for human erythrocytes. Together, these studies form the foundation for an inexpensive and rapid point-of-care device to diagnose sickle cell disease or to determine blood viscosity in resource-challenged settings.

High-throughput characterization of genetic effects on DNA-protein binding and gene transcription.

Many variants associated with complex traits are in noncoding regions and contribute to phenotypes by disrupting regulatory sequences. To characterize these variants, we developed a streamlined protocol for a high-throughput reporter assay, Biallelic Targeted STARR-seq (BiT-STARR-seq), that identifies allele-specific expression (ASE) while accounting for PCR duplicates through unique molecular identifiers. We tested 75,501 oligos (43,500 SNPs) and identified 2720 SNPs with significant ASE (FDR < 10%). To validate disruption of binding as one of the mechanisms underlying ASE, we developed a new high-throughput allele-specific binding assay for NFKB1. We identified 2684 SNPs with allele-specific binding (ASB) (FDR < 10%); 256 of these SNPs also had ASE (OR = 1.97, P-value = 0.0006). Of variants associated with complex traits, 1531 resulted in ASE, and 1662 showed ASB. For example, we characterized that the Crohn's disease risk variant for rs3810936 increases NFKB1 binding and results in altered gene expression.

Genetic-Variation-Driven Gene-Expression Changes Highlight Genes with Important Functions for Kidney Disease.

Chronic kidney disease (CKD) is a complex gene-environmental disease affecting close to 10% of the US population. Genome-wide association studies (GWASs) have identified sequence variants, localized to non-coding genomic regions, associated with kidney function. Despite these robust observations, the mechanism by which variants lead to CKD remains a critical unanswered question. Expression quantitative trait loci (eQTL) analysis is a method to identify genetic variation associated with gene expression changes in specific tissue types. We hypothesized that an integrative analysis combining CKD GWAS and kidney eQTL results can identify candidate genes for CKD. We performed eQTL analysis by correlating genotype with RNA-seq-based gene expression levels in 96 human kidney samples. Applying stringent statistical criteria, we detected 1,886 genes whose expression differs with the sequence variants. Using direct overlap and Bayesian methods, we identified new potential target genes for CKD. With respect to one of the target genes, lysosomal beta A mannosidase (MANBA), we observed that genetic variants associated with MANBA expression in the kidney showed statistically significant colocalization with variants identified in CKD GWASs, indicating that MANBA is a potential target gene for CKD. The expression of MANBA was significantly lower in kidneys of subjects with risk alleles. Suppressing manba expression in zebrafish resulted in renal tubule defects and pericardial edema, phenotypes typically induced by kidney dysfunction. Our analysis shows that gene-expression changes driven by genetic variation in the kidney can highlight potential new target genes for CKD development.

A Dementia-Associated Risk Variant near TMEM106B Alters Chromatin Architecture and Gene Expression.

Neurodegenerative diseases pose an extraordinary threat to the world's aging population, yet no disease-modifying therapies are available. Although genome-wide association studies (GWASs) have identified hundreds of risk loci for neurodegeneration, the mechanisms by which these loci influence disease risk are largely unknown. Here, we investigated the association between common genetic variants at the 7p21 locus and risk of the neurodegenerative disease frontotemporal lobar degeneration. We showed that variants associated with disease risk correlate with increased expression of the 7p21 gene TMEM106B and no other genes; co-localization analyses implicated a common causal variant underlying both association with disease and association with TMEM106B expression in lymphoblastoid cell lines and human brain. Furthermore, increases in the amount of TMEM106B resulted in increases in abnormal lysosomal phenotypes and cell toxicity in both immortalized cell lines and neurons. We then combined fine-mapping, bioinformatics, and bench-based approaches to functionally characterize all candidate causal variants at this locus. This approach identified a noncoding variant, rs1990620, that differentially recruits CTCF in lymphoblastoid cell lines and human brain to influence CTCF-mediated long-range chromatin-looping interactions between multiple cis-regulatory elements, including the TMEM106B promoter. Our findings thus provide an in-depth analysis of the 7p21 locus linked by GWASs to frontotemporal lobar degeneration, nominating a causal variant and causal mechanism for allele-specific expression and disease association at this locus. Finally, we show that genetic variants associated with risk of neurodegenerative diseases beyond frontotemporal lobar degeneration are enriched in CTCF-binding sites found in brain-relevant tissues, implicating CTCF-mediated gene regulation in risk of neurodegeneration more generally.

Transposable elements are the primary source of novelty in primate gene regulation.

Gene regulation shapes the evolution of phenotypic diversity. We investigated the evolution of liver promoters and enhancers in six primate species using ChIP-seq (H3K27ac and H3K4me1) to profile cis-regulatory elements (CREs) and using RNA-seq to characterize gene expression in the same individuals. To quantify regulatory divergence, we compared CRE activity across species by testing differential ChIP-seq read depths directly measured for orthologous sequences. We show that the primate regulatory landscape is largely conserved across the lineage, with 63% of the tested human liver CREs showing similar activity across species. Conserved CRE function is associated with sequence conservation, proximity to coding genes, cell-type specificity, and transcription factor binding. Newly evolved CREs are enriched in immune response and neurodevelopmental functions. We further demonstrate that conserved CREs bind master regulators, suggesting that while CREs contribute to species adaptation to the environment, core functions remain intact. Newly evolved CREs are enriched in young transposable elements (TEs), including Long-Terminal-Repeats (LTRs) and SINE-VNTR-Alus (SVAs), that significantly affect gene expression. Conversely, only 16% of conserved CREs overlap TEs. We tested the cis-regulatory activity of 69 TE subfamilies by luciferase reporter assays, spanning all major TE classes, and showed that 95.6% of tested TEs can function as either transcriptional activators or repressors. In conclusion, we demonstrated the critical role of TEs in primate gene regulation and illustrated potential mechanisms underlying evolutionary divergence among the primate species through the noncoding genome.

INFERNO: inferring the molecular mechanisms of noncoding genetic variants.

The majority of variants identified by genome-wide association studies (GWAS) reside in the noncoding genome, affecting regulatory elements including transcriptional enhancers. However, characterizing their effects requires the integration of GWAS results with context-specific regulatory activity and linkage disequilibrium annotations to identify causal variants underlying noncoding association signals and the regulatory elements, tissue contexts, and target genes they affect. We propose INFERNO, a novel method which integrates hundreds of functional genomics datasets spanning enhancer activity, transcription factor binding sites, and expression quantitative trait loci with GWAS summary statistics. INFERNO includes novel statistical methods to quantify empirical enrichments of tissue-specific enhancer overlap and to identify co-regulatory networks of dysregulated long noncoding RNAs (lncRNAs). We applied INFERNO to two large GWAS studies. For schizophrenia (36,989 cases, 113,075 controls), INFERNO identified putatively causal variants affecting brain enhancers for known schizophrenia-related genes. For inflammatory bowel disease (IBD) (12,882 cases, 21,770 controls), INFERNO found enrichments of immune and digestive enhancers and lncRNAs involved in regulation of the adaptive immune response. In summary, INFERNO comprehensively infers the molecular mechanisms of causal noncoding variants, providing a sensitive hypothesis generation method for post-GWAS analysis. The software is available as an open source pipeline and a web server.