A porcine brain-wide RNA editing landscape

AbstractAdenosine-to-inosine (A-to-I) RNA editing, catalyzed by ADAR enzymes, is an essential post-transcriptional modification. Although hundreds of thousands of RNA editing sites have been reported in mammals, brain-wide analysis of the RNA editing in the mammalian brain remains rare. Here, a genome-wide RNA-editing investigation is performed in 119 samples, representing 30 anatomically defined subregions in the pig brain. We identify a total of 682,037 A-to-I RNA editing sites of which 97% are not identified before. Within the pig brain, cerebellum and olfactory bulb are regions with most edited transcripts. The editing level of sites residing in protein-coding regions are similar across brain regions, whereas region-distinct editing is observed in repetitive sequences. Highly edited conserved recoding events in pig and human brain are found in neurotransmitter receptors, demonstrating the evolutionary importance of RNA editing in neurotransmission functions. Although potential data biases caused by age, sex or health status are not considered, this study provides a rich resource to better understand the evolutionary importance of post-transcriptional RNA editing.

AbstractAutism spectrum disorder (ASD) is a genetically complex, clinically heterogeneous neurodevelopmental disease. Recently, our understanding of the molecular abnormalities in ASD has been expanded through transcriptomic analyses of postmortem brains. However, a crucial molecular pathway involved in synaptic development, RNA editing, has not yet been studied on a genome-wide scale. Here, we profiled the global patterns of adenosine-to-inosine (A-to-I) editing in a large cohort of post-mortem ASD brains. Strikingly, we observed a global bias of hypo-editing in ASD brains, common to different brain regions and involving many genes with known neurobiological functions. Through genome-wide protein-RNA binding analyses and detailed molecular assays, we show that the Fragile X proteins, FMRP and FXR1P, interact with ADAR proteins and modulate A-to-I editing. Furthermore, we observed convergent patterns of RNA editing alterations in ASD and Fragile X syndrome, thus establishing RNA editing as a molecular link underlying these two highly related diseases. Our findings support a role for RNA editing dysregulation in ASD and highlight novel mechanisms for RNA editing regulation.

Abstract Motivation Genome-wide association studies (GWAS) are a powerful method to detect even weak associations between variants and phenotypes; however, many of the identified associated variants are in non-coding regions, and presumably influence gene expression regulation. Identifying potential drug targets, i.e. causal protein-coding genes, therefore, requires crossing the genetics results with functional data. Results We present a novel data integration pipeline that analyses GWAS results in the light of experimental epigenetic and cis-regulatory datasets, such as ChIP-Seq, Promoter-Capture Hi-C or eQTL, and presents them in a single report, which can be used for inferring likely causal genes. This pipeline was then fed into an interactive data resource. Availability and implementation The analysis code is available at www.github.com/Ensembl/postgap and the interactive data browser at postgwas.opentargets.io.

AbstractCerebral amyloid angiopathy (CAA) contributes to accelerated cognitive decline in Alzheimer’s disease (AD) dementia and is a common finding at autopsy. The APOEε4 allele and male sex have previously been reported to associate with increased CAA in AD. To inform biomarker and therapeutic target discovery, we aimed to identify additional genetic risk factors and biological pathways involved in this vascular component of AD etiology. We present a genome-wide association study of CAA pathology in AD cases and report sex- and APOE-stratified assessment of this phenotype. Genome-wide genotypes were collected from 853 neuropathology-confirmed AD cases scored for CAA across five brain regions, and imputed to the Haplotype Reference Consortium panel. Key variables and genome-wide genotypes were tested for association with CAA in all individuals and in sex and APOEε4 stratified subsets. Pathway enrichment was run for each of the genetic analyses. Implicated loci were further investigated for functional consequences using brain transcriptome data from 1,186 samples representing seven brain regions profiled as part of the AMP-AD consortium. We confirmed association of male sex, AD neuropathology and APOEε4 with increased CAA, and identified a novel locus, LINC-PINT, associated with lower CAA amongst APOEε4-negative individuals (rs10234094-C, beta = −3.70 [95% CI −0.49—−0.24]; p = 1.63E-08). Transcriptome profiling revealed higher LINC-PINT expression levels in AD cases, and association of rs10234094-C with altered LINC-PINT splicing. Pathway analysis indicates variation in genes involved in neuronal health and function are linked to CAA in AD patients. Further studies in additional and diverse cohorts are needed to assess broader translation of our findings.

Genodermatosis such as hair disorders mostly follow a monogenic mode of inheritance. Congenital hypotrichosis (HY) belong to this group of disorders and is characterized by abnormally reduced hair since birth. The purpose of this study was to characterize the clinical phenotype of a breed-specific non-syndromic form of HY in Belted Galloway cattle and to identify the causative genetic variant for this recessive disorder. An affected calf born in Switzerland presented with multiple small to large areas of alopecia on the limbs and on the dorsal part of the head, neck, and back. A genome-wide association study using Swiss and US Belted Galloway cattle encompassing 12 cases and 61 controls revealed an association signal on chromosome 29. Homozygosity mapping in a subset of cases refined the HY locus to a 1.5 Mb critical interval and subsequent Sanger sequencing of protein-coding exons of positional candidate genes revealed a stop gain variant in the HEPHL1 gene that encodes a multi-copper ferroxidase protein so-called hephaestin like 1 (c.1684A>T; p.Lys562*). A perfect concordance between the homozygous presence of this most likely pathogenic loss-of-function variant and the HY phenotype was found. Genotyping of more than 700 purebred Swiss and US Belted Galloway cattle showed the global spread of the mutation. This study provides a molecular test that will permit the avoidance of risk matings by systematic genotyping of relevant breeding animals. This rare recessive HEPHL1-related form of hypotrichosis provides a novel large animal model for similar human conditions. The results have been incorporated in the Online Mendelian Inheritance in Animals (OMIA) database (OMIA 002230-9913).

Abstract Amborella trichopoda is placed close to the base of the angiosperm lineage (basal angiosperm). By genome-wide RNA sequencing, we identified 184C-to-U RNA editing sites in the plastid genome of Amborella. This number is much higher than that observed in other angiosperms including maize (44 sites), rice (39 sites) and grape (115 sites). Despite the high frequency of RNA editing, the biased distribution of RNA editing sites in the genome, target codon preference and nucleotide preference adjacent to the edited cytidine are similar to that in other angiosperms, suggesting a common editing machinery. Consistent with this idea, the Amborella nuclear genome encodes 2–3 times more of the E- and DYW-subclass members of pentatricopeptide repeat proteins responsible for RNA editing site recognition in plant organelles. Among 165 editing sites in plastid protein coding sequences in Amborella, 100 sites were conserved at least in one out of 38 species selected to represent key branching points of the angiosperm phylogenetic tree. We assume these 100 sites represent at least a subset of the sites in the plastid editotype of ancestral angiosperms. We then mapped the loss and gain of editing sites on the phylogenetic tree of angiosperms. Our results support the idea that the evolution of angiosperms has led to the loss of RNA editing sites in plastids.

Recently, the prevalence and importance of RNA editing have been illuminated in mammals. However, studies on RNA editing of pigs, a widely used biomedical model animal, are rare. Here we collected RNA sequencing data across 11 tissues and identified more than 490,000 RNA editing sites. We annotated their biological features, detected flank sequence characteristics of A-to-I editing sites and the impact of A-to-I editing on miRNA–mRNA interactions, and identified RNA editing quantitative trait loci (edQTL). Sus scrofa RNA editing sites showed high enrichment in repetitive regions with a median editing level as 15.38%. Expectedly, 96.3% of the editing sites located in non-coding regions including intron, 3′ UTRs, intergenic, and gene proximal regions. There were 2233 editing sites located in the coding regions and 980 of them caused missense mutation. Our results indicated that to an A-to-I editing site, the adjacent four nucleotides, two before it and two after it, have a high impact on the editing occurrences. A commonly observed editing motif is CCAGG. We found that 4552 A-to-I RNA editing sites could disturb the original binding efficiencies of miRNAs and 4176 A-to-I RNA editing sites created new potential miRNA target sites. In addition, we performed edQTL analysis and found that 1134 edQTLs that significantly affected the editing levels of 137 RNA editing sites. Finally, we constructed PRESDB, the first pig RNA editing sites database. The site provides necessary functions associated with Sus scrofa RNA editing study.

AbstractBioinformatic analyses and tools make extensive use of k-mers (fixed contiguous strings of k nucleotides) as an informational unit. K-mer analyses are both useful and fast, but are strongly affected by single nucleotide polymorphisms or sequencing errors, effectively hindering direct-analyses of whole regions and decreasing their usability between evolutionary distant samples. Q-grams or spaced seeds, subsequences generated with a pattern of used-and-skipped nucleotides, overcome many of these limitations but introduce larger complexity which hinders their wider adoption.We introduce a concept of skip-mers, a cyclic pattern of used-and-skipped positions of k nucleotides spanning a region of size S ≥ k, and show how analyses are improved by using this simple subset of q-grams as a replacement for k-mers. The entropy of skip-mers increases with the larger span, capturing information from more distant positions and increasing the specificity, and uniqueness, of larger span skip-mers within a genome. In addition, skip-mers constructed in cycles of 1 or 2 nucleotides in every 3 (or a multiple of 3) lead to increased sensitivity in the coding regions of genes, by grouping together the more conserved nucleotides of the protein-coding regions.We implemented a set of tools to count and intersect skip-mers between different datasets, a simple task given that the properties of skip-mers make them a direct substitute for k-mers. We used these tools to show how skip-mers have advantages over k-mers in terms of entropy and increased sensitivity to detect conserved coding sequence, allowing better identification of genic matches between evolutionarily distant species. We then show benefits for multi-genome analyses provided by increased and better correlated coverage of conserved skip-mers across multiple samples.Software availabilitythe skm-tools implementing the methods described in this manuscript are available under MIT license at http://github.com/bioinfologics/skm-tools/

AbstractIn order to identify host factors that impact Bovine Herpes Virus Type 1 (BHV-1) infection we previously applied a genome wide CRISPR knockout screen with a library covering all bovine protein coding genes. We compiled a list of both pro-viral and anti-viral proteins involved in BHV-1 replication; here we provide further analysis of those that are potentially involved in viral entry into the host cell. These entry related factors include the cell surface proteins PVR and PVRL2, a group of enzymes directly or indirectly associated with the biosynthesis of Heparan Sulfate Proteoglycans (HSPG), and proteins that reside in the Golgi apparatus engaging in intra-Golgi trafficking. For the first time, we provide evidence that PVRL2 serves a receptor for BHV-1, mediating more efficient entry than the previously identified PVR. By knocking out two enzymes that catalyze HSPG chain elongation, HST2ST1 and GLCE, we demonstrated the significance of HSPG in BHV-1 entry. Another intriguing cluster of genes, COG1, COG2 and COG4-7 encodes for six subunits of the conserved oligomeric Golgi (COG) complex. MDBK cells lacking COG6 were less infectable by BHV-1 but release newly produced virions more efficiently as evidenced by fewer but bigger plaques compared to control cells, suggesting impaired HSPG biosynthesis. To facilitate candidate validation, we devised a one-step multiplex CRISPR interference (CRISPRi) system named CRISPR3i that enables quick and simultaneous deployment of three CRISPRs for efficient gene inactivation. Using CRISPR3i, we verified an additional 23 candidates, with many implicated in cellular entry.